Search NASA⌕ Search

SEARCH · Search NASA

Results for “FUN3D”

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 271 records · Page 15

Design of a Distributed Propulsion Concept Using an Adjoint-Based Approach and Blade Element Theory to Minimize Power

The adjoint-based design capability in FUN3D is extended to allow efficient gradient-based optimization and design of concepts with highly integrated and distributed aero-propulsive systems. Previous work by the authors demonstrated the derivation and application of design sensitivities of flow power and vehicle forces with respect to design variables associated with actuator disk size, loading, and installation on the aircraft. In this work, calculations of propulsive power, shaft power, and propeller efficiency using blade element-based disk loading, along with sensitivity derivatives needed to perform adjoint-based design, have been implemented. This includes the derivation of additional design sensitivities for actuator disk variables with respect to the radial distributions of blade pitch angle and chord length. The blade element-based loading implementation allows us to tailor the actuator disk loading to provide greater design authority and calculate the torque imparted on the flow for the modeling of swirl effects. The design capability is demonstrated by the optimization of NASA's parallel hybrid electric PEGASUS aircraft concept. The optimization objective is the minimization of shaft power at the aerodynamic interface planes for the wing-mounted and tail-cone boundary layer ingestion propulsors, subject to vehicle performance and propulsive constraints.

Irian Ordaz↗

CFD Validation Study of a Hypersonic Cone-Slice-Flap Configuration

Model validation is the process of determining the degree of accuracy between the real world and the model. The result of model validation can either be used to improve the model through calibration or quantify the model-form uncertainty. This work focuses on providing the model-form uncertainty through an area metric for a hypersonic cone-slice-flap geometry configuration given uncertainty in both the simulation and experimental data. A procedure will be presented that can give an accurate representation of the model-form uncertainty using a small number of high-fidelity runs. The work will also assess the impact between using different high-fidelity solvers and different turbulence models.The goal of performing this comparison is to provide a quantifiable measurement for the accuracy of each solver and turbulence model for this type of design. A low-fidelity analysis will also be performed to get a model-form uncertainty for the lower fidelity analysis. The two high-fidelty CFD solvers used are VULCAN-CFD and FUN3D and the low-fidelity results will come from Cart3D. The experimental data came from the 20-inch Mach 6Tunnel located at NASA Langley Research Center.

Laura M. White↗

Flow Characterization of the NASA Langley Unitary Plan Wind Tunnel, Test Section 2: Computational Results

This is an abstract for an invited paper at the AIAA Aviation Conference, June 2021. The work described here is part of an effort of coordinated experiments in the Unitary Plan Wind Tunnel (UPWT) facility at the NASA Langley Research Center (LaRC) and matching CFD simulations. The primary goal of the work is to assess the productivity and true predictive accuracy of CFD, absent any guidance from experiment, in the high supersonic speed range as compared to experiments performed at the NASA LaRC’s UPWT facility. This report concerns CFD simulation of the primary flow-path in the empty wind tunnel, including the settling chamber, nozzle, test section, and some of the tunnel downstream of the test section. The Mach number in the test section ranges from M~2.4 to M~4.6, and the required area ratio variation is achieved by translation of a nozzle block which constricts the area of a loosely S-shaped throat. Flow past protuberances in the settling chamber and into this S-bend throat are predicted by CFD to generate streamwise vorticity that may, or may not, persist through the throat and into the test section as coherent vortices. Some flow conditions are notably unsteady at frequencies well below those of turbulence, due to unsteady separated flow ahead of the nozzle block. The bulk flow moves at velocities ranging from 'walking speed' in the settling chamber to M~4.6 in the test section. Heat transfer to the settling chamber walls and buoyancy are significant at high Mach number. Subtle variations in surface curvature in the nozzle generate Mach waves that propagate into the test section. CFD of the empty tunnel serves two purposes. Firstly, the full-tunnel simulations are used to provide upstream boundary conditions for CFD of vehicle aerodynamics which are generally performed in a domain confined to the wind tunnel test section; these companion studies are addressed in other papers. Secondly, it is a challenging test for CFD to resolve all of the empty tunnel flow phenomena relevant to flow in the test section. It requires a more complete definition of geometry than was originally anticipated. In addition, it requires good spatial and temporal accuracy, and turbulence modeling that performs well on specific phenomena such as corner flows. The boundary conditions and solution algorithms must perform well from incompressible to almost hypersonic speeds. The final state of the pre-test CFD was a product of an iterative self-improvement process. The initial simulations of the empty tunnel were deficient in many respects, but hints to those deficiencies were recognized in the solutions, and remedies were implemented. Possible further improvements will be studied in the post-test phase when comparisons with experimental data are possible. The CFD was performed by five separate collaborative teams using four different flow solvers: FUN3D, Overflow, Star-CCM+ and USM3D. The level of effort of these teams varied significantly, but each made important contributions to the goals of the work. A concerted effort to use uncertainty quantification methods (UQ) in CFD is also a goal of this work. To this end, variations in CFD results due to grid refinement, turbulence modeling, and boundary conditions have been characterized. Code-to-code variation is another means of assessing CFD uncertainty. All CFD solvers predict similarity among the primary flow features; these include the variation of Mach number due to changes in Reynolds number, and the bulk flow angularity due to tunnel-wall curvature. All CFD solvers also predict similar trends in secondary flows, such as the downwash in the side-wall boundary layers. Three of the high-spatial resolution simulations give similar predictions of a complex secondary flow phenomena, streamwise vortices generated in the S-bend throat that persist into the side-wall boundary layers of the test section. Two of the highest-resolution simulations, run with the same turbulence model in different CFD solvers, gave encouragingly similar predictions of a complex tertiary flow phenomenon, small transient "sprites" of upwash flows, resulting from vortices that presumably originate in the separated flow near the leading edge of the nozzle block. While the CFD was done in a "blind pre-test" mode, requests from the experimental team for CFD results pertaining to unsteadiness and total temperature variations in the test section resulted in CFD runs that included a cooled wall in settling chamber. This then led to a change in the standard practice for running the Overflow results, which would not have occurred without this "release" of this experimental information. CFD was also used to guide some measurements. The paper will focus on establishing the consensus among CFD results and understanding differences among those results. Initial findings from the efforts to characterize CFD uncertainty have been done and will be included in the paper.

Robert Edward Childs↗

Implicit Thermochemical Nonequilibrium Compressible Flow Simulations on Unstructured Grids Using GPUs

As next-generation exascale-class systems arrive, existing software must be updated accordingly to effectively utilize these systems. For high concurrency and energy efficiency, many of these systems utilize GPU architectures. In this work, we present a CUDA C++ implementation of FUN3D's thermochemical nonequilibrium capability for turbulent flows. Efficiency is demonstrated at scale using the Summit system at the Oak Ridge Leadership Computing Facility which is representative of future exascale systems. This work enables faster, higher fidelity, and scale-resolving simulations of thermochemical nonequilibrium flows including reentry, hypersonics, and combustion.

CFD, GPU, HPC, Hypersonics, Chemistry↗

In-Tunnel Simulations of the NASA Juncture Flow Model

The effect of wind tunnel walls on the development of flow separation in the wing-fuselage juncture area is assessed using FUN3D’s overset and wind tunnel controller capabilities. The simulation data obtained from the in-tunnel runs are compared with the free-air simulations. The simulation results are also evaluated using measurements from the recently completed Phase 2 of the Juncture Flow Experiment. In general, the simulation data compared well with the measurements, however there were no significant differences in the prediction of the separation location, size, and shape, between the free-air and the in-tunnel runs. The results of this study suggest that the effects of tunnel walls on the flow separation in the juncture region are relatively small.

Juncture Flow↗

CFD Validation Study of a Hypersonic Cone-Slice-Flap Configuration

Model validation is the process of determining the degree of accuracy between the real world and the model. The result of model validation can be used to either improve the model through calibration or quantify the model-form uncertainty. This work focuses on estimating the model-form uncertainty through an area metric for a hypersonic cone-slice-flap geometry configuration given the uncertainty in both the simulation and experimental data. A procedure that can give an accurate representation of the model-form uncertainty using a small number of high-fidelity runs is outlined. The work also assesses the impact of using different high-fidelity solvers and different turbulence models. The goal of performing this comparison is to provide a quantifiable measurement of the accuracy of each solver and turbulence model for this type of design. A low-fidelity analysis is also performed to get a model-form uncertainty for this type of analysis. The two high-fidelity CFD solvers used here are VULCAN-CFD and FUN3D, and the low-fidelity results comes from Cart3D.The experimental data comes from the 20-inch Mach 6 Tunnel located at NASA Langley Research Center. The work here shows that the models used tend to under-predict the drag and moment aerodynamic coefficients for this type of design leading to a model-form uncertainty estimate which can account for up to an 86% of the total uncertainty in the model predictions.

Laura M White↗

Reduced-Order Aerodynamic Modeling Based on CFD Frequency Responses from Multisine Inputs

A system identification analysis was performed to determine a reduced-order model (ROM) of a computational fluid dynamics (CFD) solver in support of linear aeroservoelastic model development and feedback control design. The approach was applied to the FUN3D code for the half-span wind tunnel test article used in the NASA-Boeing collaboration called the Integrated Adaptive Wing Technology Maturation (IAWTM) project. In a transonic flow condition, multiple inputs (11 structural mode displacements and 3 control surface deflections) were simultaneously excited with orthogonal phase-optimized multisines while multiple outputs (the corresponding 14 generalized aerodynamic forces) were recorded. From these recorded times series, the matrix of frequency responses was computed and subsequently fit using rational function approximations (RFAs). It was found that the entire (14 x 14) matrix of frequency responses could be determined from a single CFD run and that results generally followed trends predicted using other methods. Differences were attributed to the modeling fidelity and nonlinearities from structural mode and control surface interactions at higher reduced frequencies. More accurate fits of the RFAs to the frequency response data were obtained by making two CFD runs, one with only structural mode excitations and one with only control surface excitations, which reduced the degree of nonlinearity in the modeling data.

Aeroservoelasticity↗

The Influence of Computer Architecture on Performance and Scaling for Hypersonic Flow Simulations

It is critical to understand how hypersonic simulation tools perform on a range of computational platforms. This information will aid in the acquisition of appropriate hardware and the potential refactoring of hypersonic codes to run on different systems. In this paper, we consider two representative high-speed reacting flow cases: a model Mach 8 hypersonic waverider glide vehicle and a model hydrocarbon-fueled hypersonic ramjet propulsion system. In both scenarios, the flow fields are in chemical non-equilibrium and are modeled by the multi-species reacting Navier-Stokes equations. For these simulations we use several hypersonic simulation tools, including US3D, Kestrel, FUN3D, and the JENRER© flow solver. We explore several high performance computing systems containing IntelR© XeonR© Platinum processors, AMD EPYCTM7702 processors, and NVIDIAR© Tesla V100 devices. We compare performance and strong scaling between the different systems.

CPU↗

In-Tunnel Simulations of the NASA Juncture Flow Model

The effect of wind tunnel walls on the development of flow separation in the wing-fuselage juncture area is assessed using FUN3D’s overset and wind tunnel controller capabilities. The simulation data obtained from the in-tunnel runs are compared with the free-air simulations. The simulation results are also evaluated using measurements from the recently completed Phase 2 of the Juncture Flow Experiment. In general, the simulation data compared well with the measurements, however there were no significant differences in the prediction of the separation location, size, and shape, between the free-air and the in-tunnel runs. The results of this study suggest that the effects of tunnel walls on the flow separation in the juncture region are relatively small.

Juncture Flow↗

Study of Mach 0.8 Transonic Truss-Braced Wing Aircraft Wing-Strut Interference Effects

This paper presents computational study of transonic wing-strut interference effects of the Mach 0.8 Transonic Truss-Braced Wing (TTBW) aircraft using the high-fidelity CFD solver FUN3D. The study is conducted for the wing-strut and the wing-alone configurations at design Mach number0.8 and Reynolds number14.0 ́106.The interference effects are calculated by comparing the wing aerodynamics along the spanwise direction between the wing-strut and the wing-alone configurations. The presence of the strut underneath the wing induces a suction peak on the lower surface of the wing, which causes changes in aerodynamic forces and moments. The interference effects of the Mach 0.8 TTBW aircraft are compared with the Mach 0.745 TTBW aircraft. A transonic wing-strut interference aerodynamic correction model is developed for use in a lower-fidelity tool, VSPAERO, for rapid aerodynamic analysis of the Mach 0.8 TTBW aircraft.

TTBW↗

Aerodynamic Optimization of Mach 0.745 Transonic Truss-Braced Wing Aircraft with Variable-Camber Continuous Trailing-Edge Flap

This paper presents an aerodynamic optimization study of the Mach 0.745 Transonic Truss-Braced Wing (TTBW) aircraft with Variable-Camber Continuous Trailing-Edge Flap (VC-CTEF). The VCCTEF is a novel wing shaping control concept to improve aircraft aerodynamic efficiency. Drag reduction studies are conducted for three different VCCTEF configurations with 6-, 10-, and 21-spanwise sections, respectively. A simple VCCTEF actuator weight model is proposed to account the weight penalty of the actuator in the design. A vortex-lattice model of the TTBW model is developed with transonic small disturbance, integral boundary-layer, and wing-strut interference corrections for rapid aerodynamic performance evaluations. The VSPAERO model has been validated against wind tunnel test data. The optimization results show that the 21-spanwise sections VCCTEF has largest absolute drag coefficient reduction, however, the 6-spanwise sections VCCTEF provides a relatively better solution for drag reduction when the actuator weight penalty is considered. A high-fidelity CFD solver FUN3D is used to validate the VCCTEF optimization.

Aerodynamic Optimization↗

Transonic Correction to Theodorsen's Theory for Oscillating Airfoil in Pitch and Plunge Toward Flutter

This paper presents a transonic correction method for an oscillating airfoil in pitch and plunge. The proposed method applies correction functions to the Theodorsen’s theory to capture the transonic nonlinear aero- dynamics. These correction functions apply necessary corrections to the amplitudes and the phase angles of the unsteady lift and pitching moment coefficients to account for transonic aerodynamics. The proposed method also postulates a correction for the motion of the aerodynamic center which could be induced by moving shocks. A series of unsteady RANS CFD simulations of the airfoil at the mean aerodynamic chord of the Transonic Truss-Braced Wing aircraft are conducted using FUN3D to provide data to construct these transonic correction functions. The computed responses of the unsteady lift and pitching moment coefficients using these transonic correction functions match the CFD simulation results very well even when the pitching moment coefficient is highly nonlinear. A flutter analysis of an airfoil in pitch and plunge illustrates the potential use of the proposed transonic correction method.

Transonic Flutter↗

Transonic Correction Method for Flight Dynamic Stability Analysis of Mach 0.745 Transonic Truss-Braced Wing

This paper presents a transonic correction method for obtaining dynamic stability derivatives for flight dynamic stability analysis. The method provides a transonic correction to the Theodorsen's theory of unsteady aerodynamics using FUN3D CFD solver of unsteady Reynolds-averaged Navier-Stokes equations (RANS) for a series of wing sections of the Mach 0.745 Transonic Truss-Braced Wing in pitch and plunge oscillations. Unsteady lift and pitching moment coefficients are obtained and used to develop the correction terms in the Theodorsen's theory to account for transonic aerodynamics. The unsteady lift and pitching moment derivatives with respect to the unsteady angle of attack are obtained as frequency response functions of the reduced frequency. These frequency response functions are used to compute the dynamic stability derivatives of lift and pitching moment due to the angle of attack and pitch rate and the dynamic stability derivatives for the rolling moment and yawing moment with respect to the roll rate and yaw rate. A transonic correction is applied to steady-state stability derivatives computed by VSPAERO solver using transonic small disturbance code TSFOIL coupled to an integral boundary method. A dynamic stability analysis is conducted for longitudinal and lateral-directional motions. Without transonic corrections and dynamic stability derivatives, the analysis indicates an unstable phugoid mode. The transonic correction applied to the steady-state stability derivatives computed by VSPAERO shows a stable phugoid mode. This is due to the increase of the drag stability derivatives as a result of the additional wave drag contribution in transonic flow. The effect of the transonic dynamic stability derivatives is observed to be a significant contributor to the increase in the damping values of all the flight dynamic modes of the Mach 0.745 Transonic-Truss Braced Wing.

Transonic↗

Aeroelastic Trim Drag Optimization of Mach 0.745 Transonic Truss-Braced Wing Aircraft with Variable-Camber Continuous Trailing-Edge Flap

This paper presents an aeroelastic trim drag optimization study of the Mach 0.745 Transonic Truss-Braced Wing (TTBW) aircraft with the Variable-Camber Continuous Trailing-Edge Flap (VCCTEF). An aero-structural analysis solver VSPAERO with transonic small disturbance, integral boundary-layer, and wing-strut interference corrections coupled to finite-element code BEAM3D is developed to provide a rapid aircraft aeroelastic performance evaluation. The VSPAERO model has been validated against wind tunnel test data. Aeroelastic trim drag optimization studies are conducted for a VCCTEF configuration with 6-spanwise sections. Three different flight conditions corresponding to Mach 0.745 are selected for the TTBW aircraft aeroelastic trim drag optimization at the design and off-design cruise lift coefficients. The preliminary optimization results show that the TTBW aircraft with the optimized VCCTEF deflection achieves a drag reduction of about 9.1 counts, 5.9 counts, and 11.8 counts corresponding to the lift coefficients 0.68, 0.73, and 0.78, respectively. A high-fidelity CFD solver FUN3D is used to validate the aeroelastic trim drag optimization.

Aeroelastic trim drag↗

Aeroelastic Modeling and CFD Simulation of Wind-Tunnel Scale Aspect Ratio 13.5 Common Research Model

This paper presents an aeroelastic modeling and simulation study of an aspect ratio 13.5 wind-tunnel scale Common Research Model (CRM) with distributed flaps. A vortex-lattice VSPAERO model of the CRM model is developed. A transonic small disturbance/integral boundary layer correction method is implemented in the VSPAERO model to account for the transonic and viscous flow effects. The structural deformation of the CRM model is calculated using a NASTRAN equivalent beam model. The VSPAERO model is coupled to the NASTRAN equivalent beam model to provide a rapid aero-structural analysis. A validation of the VSPAERO aeroelastic model is conducted by comparing the results to FUN3D CFD aeroelastic simulation results. An aerodynamic database is generated using the developed VSPAERO aeroelastic model for the real-time drag optimization and maneuver load alleviation study of the wind-tunnel scale CRM model.

Aeroelastic Modeling↗

Parametric Mixer Design Methods for Acoustic Explorations of Internally Mixed Nozzles

As part of a jet noise reduction effort under the Commercial Supersonics Technology project, a process for the design and optimization of a parametric dual-stream lobe mixer is presented. The lobe mixer is fully defined from five independent parameters, including mixer length, lobe count, lobe penetration, lobe bias, and trailing edge slope. The mixers are inserted into an academic, constant-area nozzle with an external plug. The full geometry is meshed in Pointwise v18.4, and the RANS equations are solved in FUN3D v13.6. The results of a DoE study, involving ~80 different mixer designs, are presented, along with an associated parametric model for various performance response quantities. With a design space fully explored, observations and recommendations are given, and future work is discussed to further improve the design.

Dual-stream Mixer↗

Cross-Validation of Computational and Experimental Distributed Surface Pressures on the Space Launch System

This paper presents a new workflow for comparing experimental pressure-sensitive paint (PSP) data to computational fluid dynamic (CFD) simulations by way of mapping data from corresponding grids utilizing interpolation methods. In addition to generating quantitative and qualitative point-to-point comparisons between PSP and CFD data, this workflow extracts sectional loading data from both grids and generates lineload comparison charts for corresponding PSP and CFD runs. Experimental PSP data presented in this paper were taken from a 2016 NASA Ames Research Center Unitary Plan Wind Tunnel 11- by 11-Foot Transonic WindTunnel Facility test of the NASA Space Launch System. CFD simulation data for comparison purposes were generated using the FUN3D code. Overall, interpolation onto PSP grids versus CFD grids yields comparable surface pressure fields. However, lineload comparisons are easier to make on the CFD grid-mapped data due to the grid topology and the current capabilities of the lineload analysis tools at NASA Langley Research Center. This workflow is written using contemporary software (Python, Tecplot, PyTecplot), is compatible with existing tools at NASA Langley, and is developed to be adaptable depending on the situation.

SLS↗

Progress on Transonic Flutter and Shock Buffet Computationsin Support of the Third Aeroelastic Prediction Workshop

This paper reports on the progress of the NASA Langley team contributions to the third Aeroelastic Pre-diction Workshop’s (AePW-3) High Angle Working Group (HAWG). The primary objectives of HAWG is to predict the fluter dynamic pressure of the NASA Benchmark Supercritical Wing (BSCW) configuration at Mach 0.8 and 5◦angle of attack. The secondary objective is to determine if a shock-buffet onset is present at or near that flow condition. The computational results are obtained using FUN3D, an unstructured grid Reynolds-averaged Navier-Stokes solver developed at the NASA Langley Research Center. The preliminary analysis results show a computationally-obtained flutter dynamic pressure of approximately 120 psf. Initial results describing unforced BSCW unsteady flow environment at flutter condition are also presented.

Pawel Chwalowski↗