Results for “Linear systems solvers” — Page 8

NASA NTRS1999

Sonic Boom Prediction and Minimization of the Douglas Reference OPT5 Configuration

Conventional CFD methods and grids do not yield adequate resolution of the complex shock flow pattern generated by a real aircraft geometry. As a result, a unique grid topology and supersonic flow solver was developed at Northrop Grumman based on the characteristic behavior of supersonic wave patterns emanating from the aircraft. Using this approach, it was possible to compute flow fields with adequate resolution several body lengths below the aircraft. In this region, three-dimensional effects are diminished and conventional two-dimensional modified linear theory (MLT) can be applied to estimate ground pressure signatures or sonic booms. To accommodate real aircraft geometries and alleviate the burdensome grid generation task, an implicit marching multi-block, multi-grid finite-volume Euler code was developed as the basis for the sonic boom prediction methodology. The Thomas two-dimensional extrapolation method is built into the Euler code so that ground signatures can be obtained quickly and efficiently with minimum computational effort suitable to the aircraft design environment. The loudness levels of these signatures can then be determined using a NASA generated noise code. Since the Euler code is a three-dimensional flow field solver, the complete circumferential region below the aircraft is computed. The extrapolation of all this field data from a cylinder of constant radius leads to the definition of the entire boom corridor occurring directly below and off to the side of the aircraft's flight path yielding an estimate for the entire noise "annoyance" corridor in miles as well as its magnitude. An automated multidisciplinary sonic boom design optimization software system was developed during the latter part of HSR Phase 1. Using this system, it was found that sonic boom signatures could be reduced through optimization of a variety of geometric aircraft parameters. This system uses a gradient based nonlinear optimizer as the driver in conjunction with a computationally efficient Euler CFD solver (NIIM3DSB) for computing the three-dimensional near-field characteristics of the aircraft. The intent of the design system is to identify and optimize geometric design variables that have a beneficial impact on the ground sonic boom. The system uses a simple wave drag data format to specify the aircraft geometry. The geometry is internally enhanced and analytic methods are used to generate marching grids suitable for the multi-block Euler solver. The Thomas extrapolation method is integrated into this system, and hence, the aircraft's centerline ground sonic boom signature is also automatically computed for a specified cruise altitude and yields the parameters necessary to evaluate the design function. The entire design system has been automated since the gradient based optimization software requires many flow analyses in order to obtain the required sensitivity derivatives for each design variable in order to converge on an optimal solution. Hence, once the problem is defined which includes defining the objective function and geometric and aerodynamic constraints, the system will automatically regenerate the perturbed geometry, the necessary grids, the Euler solution, and finally the ground sonic boom signature at the request of the optimizer.

Siclari, Michael J.↗

NASA NTRS2023

Measurements and Computations of Natural Transition on the NASA Juncture-Flow Model with a Symmetric Wing

Experiments were performed in the 14- by 22-Foot Subsonic Tunnel to assess natural transition on the symmetric-airfoil wings of the NASA Juncture-Flow Model. Infrared thermography was used to visualize the heating on the upper surface of both wings of the full-span model, and on the fuselage, for angles of incidence ranging from -10° to 10° at a fixed Reynolds number of 2.4E6 based on the chord length at the wing planform break. The fuselage boundary layer transitioned well upstream of the wing-root leading edge for all conditions. Transition fronts were identified by a steep rise in the surface temperature, and the transition coordinates were transformed from an image-based to a body-fixed system. Additionally, the state of the boundary layer was estimated at pressure ports distributed on the wings through observation of the pressure coefficient as a function of the angle of incidence. For increasing angles of incidence, the transition front was observed to advance upstream, in a mostly spanwise-uniform fashion, from near midchord at α = 0°; however, for increasingly negative angles of incidence, the transition front first receded and then advanced in a nonuniform jagged manner that is typically observed with stationary crossflow. The transition wedges first appeared inboard of the wing break and then spread outboard to near the tip by α = -6°. The upstream shift in transition at positive angles of incidence and the outboard progression of crossflow-dominated transition at increasingly negative angles of incidence are consistent with trends identified in a computational assessment of the boundary-layer transition based on both linear stability analysis and Reynolds-averaged-Navier-Stokes-based transition models. The stability results obtained from the Langley Stability and Transition Analysis Code were used to recalibrate a dual N-factor criterion, which allowed for the prediction of transition fronts that showed excellent agreement with the experiment. The Reynolds-averaged-Navier-Stokes-based models, from the NASA OVERFLOW 2.3 solver, that accounted for the crossflow instability showed mixed results in comparison with the experiment, with the helicity-based Langtry-Menter model performing the best. The experimental data, particularly the cases involving strong influence from both Tollmien-Schlichting and crossflow instabilities, will be valuable for the continued validation and improvement of transition models.

Boundary layer transition↗

NASA NTRS2023

Measurements and Computations of Natural Transition on the NASA Juncture-Flow Model with a Symmetric Wing

Experiments were performed in the 14- by 22-Foot Subsonic Tunnel to assess natural transition on the symmetric-airfoil wings of the NASA Juncture-Flow Model. Infrared thermography was used to visualize the heating on the upper surface of both wings of the full-span model, and on the fuselage, for angles of incidence ranging from -10° to 10° at a fixed Reynolds number of 2.4E6 based on the chord length at the wing planform break. The fuselage boundary layer transitioned well upstream of the wing-root leading edge for all conditions. Transition fronts were identified by a steep rise in the surface temperature, and the transition coordinates were transformed from an image-based to a body-fixed system. Additionally, the state of the boundary layer was estimated at pressure ports distributed on the wings through observation of the pressure coefficient as a function of the angle of incidence. For increasing angles of incidence, the transition front was observed to advance upstream, in a mostly spanwise-uniform fashion, from near midchord at α = 0°; however, for increasingly negative angles of incidence, the transition front first receded and then advanced in a nonuniform jagged manner that is typically observed with stationary crossflow. The transition wedges first appeared inboard of the wing break and then spread outboard to near the tip by α = -6°. The upstream shift in transition at positive angles of incidence and the outboard progression of crossflow-dominated transition at increasingly negative angles of incidence are consistent with trends identified in a computational assessment of the boundary-layer transition based on both linear stability analysis and Reynolds-averaged-Navier-Stokes-based transition models. The stability results obtained from the Langley Stability and Transition Analysis Code were used to recalibrate a dual N-factor criterion, which allowed for the prediction of transition fronts that showed excellent agreement with the experiment. The Reynolds-averaged-Navier-Stokes-based models, from the NASA OVERFLOW 2.3 solver, that accounted for the crossflow instability showed mixed results in comparison with the experiment, with the helicity-based Langtry-Menter model performing the best. The experimental data, particularly the cases involving strong influence from both Tollmien-Schlichting and crossflow instabilities, will be valuable for the continued validation and improvement of transition models.

boundary layer transition↗

NASA NTRS2023

Measurements and Computations of Natural Transition on the NASA Juncture-Flow Model with a Symmetric Wing

Experiments were performed in the 14- by 22-Foot Subsonic Tunnel to assess natural transition on the symmetric-airfoil wings of the NASA Juncture-Flow Model. Infrared thermography was used to visualize the heating on the upper surface of both wings of the full-span model, and on the fuselage, for angles of incidence ranging from -10° to 10° at a fixed Reynolds number of 2.4E6 based on the chord length at the wing planform break. The fuselage boundary layer transitioned well upstream of the wing-root leading edge for all conditions. Transition fronts were identified by a steep rise in the surface temperature, and the transition coordinates were transformed from an image-based to a body-fixed system. Additionally, the state of the boundary layer was estimated at pressure ports distributed on the wings through observation of the pressure coefficient as a function of the angle of incidence. For increasing angles of incidence, the transition front was observed to advance upstream, in a mostly spanwise-uniform fashion, from near midchord at α = 0°; however, for increasingly negative angles of incidence, the transition front first receded and then advanced in a nonuniform jagged manner that is typically observed with stationary crossflow. The transition wedges first appeared inboard of the wing break and then spread outboard to near the tip by α = -6°. The upstream shift in transition at positive angles of incidence and the outboard progression of crossflow-dominated transition at increasingly negative angles of incidence are consistent with trends identified in a computational assessment of the boundary-layer transition based on both linear stability analysis and Reynolds-averaged-Navier-Stokes-based transition models. The stability results obtained from the Langley Stability and Transition Analysis Code were used to recalibrate a dual N-factor criterion, which allowed for the prediction of transition fronts that showed excellent agreement with the experiment. The Reynolds-averaged-Navier-Stokes-based models, from the NASA OVERFLOW 2.3 solver, that accounted for the crossflow instability showed mixed results in comparison with the experiment, with the helicity-based Langtry-Menter model performing the best. The experimental data, particularly the cases involving strong influence from both Tollmien-Schlichting and crossflow instabilities, will be valuable for the continued validation and improvement of transition models.

boundary layer transition↗

NASA NTRS2003

Aeroelastic System Development Using Proper Orthogonal Decomposition and Volterra Theory

This research combines Volterra theory and proper orthogonal decomposition (POD) into a hybrid methodology for reduced-order modeling of aeroelastic systems. The out-come of the method is a set of linear ordinary differential equations (ODEs) describing the modal amplitudes associated with both the structural modes and the POD basis functions for the uid. For this research, the structural modes are sine waves of varying frequency, and the Volterra-POD approach is applied to the fluid dynamics equations. The structural modes are treated as forcing terms which are impulsed as part of the uid model realization. Using this approach, structural and uid operators are coupled into a single aeroelastic operator. This coupling converts a free boundary uid problem into an initial value problem, while preserving the parameter (or parameters) of interest for sensitivity analysis. The approach is applied to an elastic panel in supersonic cross ow. The hybrid Volterra-POD approach provides a low-order uid model in state-space form. The linear uid model is tightly coupled with a nonlinear panel model using an implicit integration scheme. The resulting aeroelastic model provides correct limit-cycle oscillation prediction over a wide range of panel dynamic pressure values. Time integration of the reduced-order aeroelastic model is four orders of magnitude faster than the high-order solution procedure developed for this research using traditional uid and structural solvers.

Lucia, David J.↗

NASA NTRS2014

High-Lift OVERFLOW Analysis of the DLR-F11 Wind Tunnel Model

In response to the 2nd AIAA CFD High Lift Prediction Workshop, the DLR-F11 wind tunnel model is analyzed using the Reynolds-averaged Navier-Stokes flow solver OVERFLOW. A series of overset grids for a bracket-off landing configuration is constructed and analyzed as part of a general grid refinement study. This high Reynolds number (15.1 million) analysis is done at multiple angles-of-attack to evaluate grid resolution effects at operational lift levels as well as near stall. A quadratic constitutive relation recently added to OVERFLOW for improved solution accuracy is utilized for side-of-body separation issues at low angles-of-attack and outboard wing separation at stall angles. The outboard wing separation occurs when the slat brackets are added to the landing configuration and is a source of discrepancy between the predictions and experimental data. A detailed flow field analysis is performed at low Reynolds number (1.35 million) after pressure tube bundles are added to the bracket-on medium grid system with the intent of better understanding bracket/bundle wake interaction with the wing's boundary layer. Localized grid refinement behind each slat bracket and pressure tube bundle coupled with a time accurate analysis are exercised in an attempt to improve stall prediction capability. The results are inconclusive and suggest the simulation is missing a key element such as boundary layer transition. The computed lift curve is under-predicted through the linear range and over-predicted near stall, and the solution from the most complete configuration analyzed shows outboard wing separation occurring behind slat bracket 6 where the experiment shows it behind bracket 5. These results are consistent with most other participants of this workshop.

High-lift↗

NASA NTRS2018

Computational Investigation of Nominally-Orthogonal Pneumatic Active Flow Control for High-Lift Systems

We explore the feasibility of using nominally-orthogonal jets as active aerodynamic load control for multi-element high-lift systems, and whether the nominally-orthogonal jets can offer a variety of performance improvements. These nominally-orthogonal jets inject momentum normal to the airfoil surface near the flap trailing edge, where they create a vortex that entrains flow from the opposing side and change the airfoil circulation. Lift-enhancement opportunities of trailing edge nominally-orthogonal jets have previously been studied by Malavard et al. and Blaylock et al. on single-element airfoils; however, their effect on drag was not thoroughly investigated. In this study, we investigate two- dimensional nominally-orthogonal jet effects on both lift and drag on a two-element airfoil, NLR7301. We utilize Chimera Grid Tools to generate structured curvilinear overset grids, and the Reynolds-averaged Navier-Stokes solver OVERFLOW-2 to solve for the flow field around the airfoil. We perform various computational sensitivity studies on the baseline airfoil without a jet to validate computational results against benchmark experimental data. Using a Chimera overset grid topology, we demonstrate a similar lift-enhancement effect between a nominally-orthogonal jet and a nominally-orthogonal physical tab employed at the same location on the studied airfoil. After we introduce the nominally-orthogonal jet concept, we investigate nominally-orthogonal jets with various momentum coefficient settings, C(sub µ) = 0.00 − 0.04 and present a lift-enhancement relationship ∆C(sub l) ≃ 3.59 (C(sub µ)) for this airfoil. We discuss that utilizing a nominally-orthogonal jet with C(sub µ) = 0.01 can shift the linear region of the lift curve by a ∆C(sub l) = 0.36 for the pressure side jet and by a ∆C(sub l) = −0.27 for the suction side jet. Employing a nominally-orthogonal jet is also shown effective in altering the drag. To study the impact, we carry out a drag decomposition study in the form of drag polars. We show for a given C(sub l) = 2.50, a nominally-orthogonal jet with C(sub µ) = 0.01 on the pressure and suction side of the airfoil results in 113 drag count decrements and 41 drag count increments, respectively, compared to the baseline airfoil with no jet. These results show that large and controllable changes in aerodynamic performance can be achieved by relatively small active flow control inputs using the nominally-orthogonal jets presented in this study.

Nominally-Orthogonal↗

NASA NTRSDate not supplied

Preliminary Results Cycling GEOS-JEDI with GSI-based Background Errors

The first phase of transitioning the NASA GMAO GEOS atmospheric data assimilation capabilities to JEDI involves the replacement of the Grid-point Statistical Interpolation (GSI) with a corresponding JEDI analysis. This includes taking JEDI's Unified Observation Operator (UFO), its underlying dependencies, and the JEDI solver that enables a hybrid 4DEnVar strategy similar to what is used in the current GEOS-GSI system. Variational analysis involves at least two main components associated with the observation and background cost function terms. The first is directly related to the UFO, which is being carefully validated in a joint collaboration between GMAO and NCEP to demonstrate consistency with corresponding observations usage in GSI. The second component is the background term, which in a hybrid system involves the ability to set up both a climatologically-based term and an ensemble-based term. JEDI provides the means to implement both terms through its BUMP component. Use of BUMP would require a complete re-tune of both climatological and ensemble, which is a non-trivial exercise we would prefer to avoid. As an alternative, the work here studies the results of interfacing the GSI-background error capability (GSIBEC) into JEDI through SABER. With this, the exact same background error covariance formulation used in GSI can be employed in JEDI without need for re-tuning. This brief summary covers the work done to interface GSIBEC into JEDI and shows preliminary results where the background error covariances of the control (GEOS-GSI) and experiment (GEOS-JEDI) are identical in corresponding cycling experiments. The cycling exercise is obviously preliminary and so much can be expected from GEOS-JEDI when compared to GEOS-GSI. There is still a number of features that need closer attention and although in some cases in principle ready to cycle have been intentionally either turned off or not fully exercised (e.g., VarBC is applied but not cycled). Other features are still pending implementation, one such example is the implementation of the Tangent Linear Normal Mode Constraint. Still, results are quite encouraging as hopefully the discussion here illustrates.

Ricardo Todling↗

NASA NTRSDate not supplied

Aerial Vehicle Routing and Scheduling for UAS Traffic Management: A Monte Carlo Tree Search Approach

Numerous unmanned aircraft systems operating at low altitudes to deliver goods and services may one day become ubiquitous in our cities. In the Unmanned Aircraft Systems (UAS) Traffic Management (UTM) framework, such a concept is envisioned, where aerial vehicles operate beyond visual line of sight (BVLOS) within specifically reserved and time stamped “corridors” in the airspace. For example, these corridors or operational intent volumes can connect an aerial vehicle’s origin site to its destination site for package delivery operations. There may also be more than one corridor available for an aerial vehicle to choose from and often different corridors may intersect with one another. Thus, it is imperative to ensure flight trajectories belonging to different aerial vehicles are not in conflict. Per the UTM CONOPs, we assume that a vehicle almost always stays inside its corridor or operational volume. This work provides a framework for strategic deconfliction of UTM or package delivery drones, where we schedule the departure time of all vehicles subject to various temporal constraints (including the corridor deconfliction at the intersections). We present the “multi-route weighted package delivery problem” which serves as an exemplifying model for strategic deconfliction in UTM. In the multi-route weighted package delivery problem, a graph network is given which consists of a set of depots (source) and drop-off (destination) nodes, with multiple routes (defined as a sequence of waypoints) connecting the depots to drop-off nodes. In addition, routes are weighted by the associated ground risk and total travel distance for package delivery. The goal is for a known set of aerial vehicles to depart from the depots, choose a route and take off time, while avoiding conflicts with other aerial vehicles, and minimizing both risk and distance traveled. We provide a mixed integer linear programming (MILP) formulation of the problem, as well as a heuristic solution based on Monte Carlo Tree Search (MCTS) – a method used in game theory and artificial intelligence – to overcome limitations inherent to optimal solvers. Computational results show the advantages of using MCTS over the MILP formulation; the former can provide a sub-optimal solution quickly, and may sometimes even reach an optimal solution, whereas the latter may not even produce a solution in reasonable time. Furthermore, results from both the MILP formulation and MCTS methods were validated using a preliminary agent-based simulator implementing the UTM concept of operations. Thus, the MCTS method can be seen as a scalable solution to the complex multi-route weighted package delivery problem and may possibly be extended to similar complex optimization problems.

Kenny Chour↗

NASA NTRSDate not supplied

Scheduling Position, Navigation and Time Service Requests from Non-dedicated Lunar Constellations

This paper presents a centralized scheduler that satisfies user requests for Position, Navigation, and Time (PNT) services from an ad-hoc, non-dedicated orbital constellation around the Moon. Traditional, dedicated GNSS networks provide service 24/7, which allows users to acquire localization services at-will. For ad-hoc networks, a coordinated schedule is needed to ensure Quality of Service (QoS) guarantees for user localization, while satisfying non-dedicated assets’ usage constraints. This scheduler bridges this coordination gap by leveraging Mixed Integer-Linear Programming (MILP) to schedule this “as-needed” localization service while respecting the constraints on each asset. In upcoming decades there is expected to be a substantial increase in Lunar missions. Many of these missions will feature low-cost surface assets near the moon’s polar regions and small-sat science missions in orbit. Most missions need PNT capabilities to ensure safe operations and meet their science objectives, but low-cost missions may not be able to support the large power, mass, and weight that a weak GNSS or DSN based navigation solution would entail. Asset localization has been demonstrated using a decentralized extended Kalman Filter (DEKF) in the previously presented Lunar Autonomous PNT System (LAPS). Within the LAPS simulation environment, a module has been developed to generate the coordinated user-asset schedules described above; this Service Scheduler Module (SSM) allows for complete end-to-end testing of the entire system. Within SSM, a user service request consists of a location on the Lunar surface, a cumulative service duration, and a window in which service must occur. SSM takes as input these requests and the LAPS-predicted positional degree of precision as the QoS for each available set of orbital assets. A simple, baseline MILP model is formulated to provide the highest-precision service balanced across all requests. To reflect the non-dedicated nature of the constellation, this baseline model is augmented with additional asset-specific load capacity constraints or availability constraints. The load capacity constraints limit total time spent providing service, and the availability constraints reflect blockout times or availability windows when the assets are not otherwise occupied. SSM outputs two schedules: the user schedule to indicate their service times and expected QoS, and a satellite schedule to be transmitted to the orbiting constellation, describing when each non-dedicated asset provides PNT service. SSM is predominantly implemented in MATLAB and allows the use of any MILP solver to generate the resulting schedules. This paper describes the SSM - LAPS interface, how the output of LAPS is used to construct the MILP, and how SSM provides user localization service while satisfying constraints. It will also demonstrate the tool’s flexibility for formulating schedules for the end user and the constellation, focusing on scenarios that match real-world proposed missions. It will detail how SSM can be used to compare the addition of load capacity constraints, satellite availability constraints, and QoS guarantees for the users. Finally, we describe how SSM can be used to support the design of the ad-hoc constellation itself. The resulting integrated capability will support the design of future ad-hoc Lunar PNT networks, enabling high-quality, low-cost Lunar exploration

Swarm↗

NASA NTRS2023

Evaluation of Transport-Equation-Based Transition Models for High-Speed Boundary Layers Using OVERFLOW

Accurate modeling of laminar-turbulent transition is crucial for the design of hypersonic flight systems. However, the current transition models used in production CFD codes are insufficient for high-speed flows. Many extensions to low-speed models have been suggested; however, a thorough verification and validation effort is needed before these models can be used in design settings. Challenges include potentially missing details of the model implementation requirements and/or a complete specification of the input parameters needed to replicate the test findings. A meaningful assessment of the generalization capability of these models is also hindered by a lack of information regarding the specific flow configurations and associated grids employed for model calibration. As a key first step toward model verification, we present an independent assessment of two recently proposed models for high-speed transition, namely, a model within the SST-$\gamma$ framework and a model based on the SST-$\gamma-\nu_L$ equations. These models are implemented in the NASA OVERFLOW 2.3e solver and their performance in predicting first mode, second mode, and crossflow transition has been evaluated for several test cases in the supersonic and hypersonic regimes. Besides the test cases employed by the model developers, which could have also been used for model calibration, the present assessment includes supplementary configurations that contribute to an unbiased assessment of the models. The outcomes presented in this study indicate the potential for the models to be applied to high-speed flight configurations. Key steps toward future improvements to these models are also outlined.

High-speed flow↗