Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “weak grids”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS2017

Genetic Algorithm-Guided, Adaptive Model Order Reduction of Flexible Aircrafts

This paper presents a methodology for automated model order reduction (MOR) of flexible aircrafts to construct linear parameter-varying (LPV) reduced order models (ROM) for aeroservoelasticity (ASE) analysis and control synthesis in broad flight parameter space. The novelty includes utilization of genetic algorithms (GAs) to automatically determine the states for reduction while minimizing the trial-and-error process and heuristics requirement to perform MOR; balanced truncation for unstable systems to achieve locally optimal realization of the full model; congruence transformation for "weak" fulfillment of state consistency across the entire flight parameter space; and ROM interpolation based on adaptive grid refinement to generate a globally functional LPV ASE ROM. The methodology is applied to the X-56A MUTT model currently being tested at NASA/AFRC for flutter suppression and gust load alleviation. Our studies indicate that X-56A ROM with less than one-seventh the number of states relative to the original model is able to accurately predict system response among all input-output channels for pitch, roll, and ASE control at various flight conditions. The GA-guided approach exceeds manual and empirical state selection in terms of efficiency and accuracy. The adaptive refinement allows selective addition of the grid points in the parameter space where flight dynamics varies dramatically to enhance interpolation accuracy without over-burdening controller synthesis and onboard memory efforts downstream. The present MOR framework can be used by control engineers for robust ASE controller synthesis and novel vehicle design.

Numerical Analysi↗
NASA NTRS2020

The Evolution of a Circumplanetary Disc with a Dead Zone

We investigate whether the regular Galilean satellites could have formed in the dead zone of a circumplanetary disc. A deadzone is a region of weak turbulence in which the magnetorotational instability is suppressed, potentially an ideal environmentfor satellite formation. With the grid-based hydrodynamic codeFAR GO3D, we examine the evolution of a circumplanetary discmodel with a dead zone. Material accumulates in the dead zone of the disc leading to a higher total mass and but a similartemperature profile compared to a fully turbulent disc model. The tidal torque increases the rate of mass transport through thedead zone leading to a steady-state disc with a dead zone that does not undergo accretion outbursts. We explore a range of disc,dead zone, and mass inflow parameters and find that the maximum mass of the disc is around 0.001MJ. Since the total solidmass of such a disc is much lower, we find that there is not sufficient material in the disc forin situformation of the Galileansatellites and that external supplement is required.

Cheng Chen↗
NASA NTRS1992

Numerical investigation of tail buffet on F-18 aircraft

Numerical investigation of vortex induced tail buffet is conducted on the F-18 aircraft at high angles of attack. The Reynolds-averaged Navier-Stokes equations are integrated using a time-accurate, implicit procedure. A generalized overset zonal grid scheme is used to decompose the computational space around the complete aircraft with faired-over inlet. A weak coupling between the aerodynamics and structures is assumed to compute the structural oscillation of the flexible vertical tail. Time-accurate computations of the turbulent flow around the F-18 aircraft at 30 degrees angle of attack show the surface and off-surface flowfield details, including the unsteadiness created by the vortex burst and its interaction with the vertical twin tail which causes the tail buffet. The effect of installing a LEX fence on modifying the vortex structure upstream of the tail is also examined.

Rizk, Yehia M.↗
NASA NTRS2022

Updraft Dynamics and Microphysics: On the Added Value of The Cumulus Thermal Reference Frame in Simulations of Aerosol–Deep Convection Interactions

One fundamental question about atmospheric moist convection processes that remains debated is whether, or under which conditions, a relevant variability in background aerosol concentrations may have a significant dynamical impact on convective clouds and their associated precipitation. Furthermore, current climate models must parameterize both the microphysical and the cumulus convection processes, but this is usually implemented separately, whereas in nature there is a strong coupling between them. As a first step to improve our understanding of these two problems, we investigate how aerosol concentrations modify key properties of updrafts in eight large-eddy-permitting regional simulations of a case study of scattered convection over Houston, Texas, in which convection is explicitly simulated and microphysical processes are parameterized. Dynamical and liquid-phase microphysical responses are investigated using the following two different reference frames: static cloudy updraft grid cells versus tracked cumulus thermals. In both frameworks, we observe the expected microphysical responses to higher aerosol concentrations, such as higher cloud number concentrations and lower rain number concentrations. In terms of the dynamical responses, both frameworks indicate weak impacts of varying aerosol concentrations relative to the noise between simulations over the observationally derived range of aerosol variability for this case study. On the other hand, results suggest that thermals are more selective than cloudy updraft grid cells in terms of sampling the most active convective air masses. For instance, vertical velocity from thermals is significantly higher at upper levels than when sampled from cloudy updraft grid points, and several microphysical variables have higher average values in the cumulus thermal framework than in the cloudy updraft framework. In addition, the thermal analysis is seen to add rich quantitative information about the rates and covariability of microphysical processes spatially and throughout tracked thermal lifecycles, which can serve as a stronger foundation for improving subgrid-scale parameterizations.

aerosols↗
NASA NTRS2022

Updraft Dynamics and Microphysics: On the Added Value of the Cumulus Thermal Reference Frame in Simulations of Aerosol–Deep Convection Interactions

One fundamental question about atmospheric moist convection processes that remains debated is whether or under what conditions a relevant variability in background aerosol concentrations may have a significant dynamical impact on convective clouds and their associated precipitation. Furthermore, current climate models must parameterize both the microphysical and the cumulus convection processes, but this is usually implemented separately, whereas in nature there is a strong coupling between them. As a first step to improve our understanding of these two problems, we investigate how aerosol concentrations modify key properties of updrafts in eight large-eddy permitting regional simulations of a case study of scattered convection over Houston, Texas, in which convection is explicitly simulated and microphysical processes are parameterized. Dynamical and liquid-phase microphysical responses are investigated using two different reference frames: static cloudy-updraft grid cells versus tracked cumulus thermals. In both frameworks we observe the expected microphysical responses to higher aerosol concentrations, such as higher cloud number concentrations and lower rain number concentrations. In terms of the dynamical responses, both frameworks indicate weak impacts of varying aerosol concentrations relative to the noise between simulations over the observationally derived range of aerosol variability for this case study. On the other hand, results suggest that thermals are more selective than cloudy-updraft grid cells in terms of sampling the most active convective air masses. For instance, vertical velocity from thermals is significantly higher at upper levels than when sampled from cloudy-updraft grid points, and several microphysical variables have higher average values in the cumulus thermal framework than in the cloudy-updraft framework. In addition, the thermal analysis is seen to add rich quantitative information about the rates and covariability of microphysical processes spatially and throughout tracked thermal lifecycles, which can serve as a stronger foundation for improving subgrid-scale parameterizations.

Daniel Hernandez-Deckers↗
NASA NTRS2005

The Closure of the Ocean Mixed Layer Temperature Budget using Level-Coordinate Model Fields

Entrainment is an important element of the mixed layer mass, heat, and temperature budgets. Conventional procedures to estimate entrainment heat advection often do not permit the closure of heat and temperature budgets because of inaccuracies in its formulation. In this study a rigorous approach to evaluate the effect of entrainment using the output of a general circulation model (GCM) that does not have an explicit prognostic mixed layer model is described. The integral elements of the evaluation are 1) the rigorous estimates of the temperature difference between mixed layer water and entrained water at each horizontal grid point, 2) the formulation of the temperature difference such that the budget closes over a volume greater than one horizontal grid point, and 3) the apparent warming of the mixed layer during the mixed layer shoaling to account for the weak vertical temperature gradient within the mixed layer. This evaluation of entrainment heat advection is compared with the estimates by other commonly used ad hoc formulations by applying them in three regions: the north-central Pacific, the Kuroshio Extension, and the Nino-3 areas in the tropical Pacific. In all three areas the imbalance in the mixed layer temperature budget by the ad hoc estimates is significant, reaching a maximum of about 4 K yr(exp -1).

temperature gradient↗
NASA NTRS2015

Global Lightning Climatology from the Tropical Rainfall Measuring Mission (TRMM), Lightning Imaging Sensor (LIS) and the Optical Transient Detector (OTD)

The Tropical Rainfall Measuring Mission (TRMM) Lightning Imaging Sensor (LIS) has been collecting observations of total lightning in the global tropics and subtropics (roughly 38 deg S - 38 deg N) since December 1997. A similar instrument, the Optical Transient Detector, operated from 1995-2000 on another low earth orbit satellite that also saw high latitudes. Lightning data from these instruments have been used to create gridded climatologies and time series of lightning flash rate. These include a 0.5 deg resolution global annual climatology, and lower resolution products describing the annual cycle and the diurnal cycle. These products are updated annually. Results from the update through 2013 will be shown at the conference. The gridded products are publicly available for download. Descriptions of how each product can be used will be discussed, including strengths, weaknesses, and caveats about the smoothing and sampling used in various products.

Cecil, Daniel J.↗
NASA NTRS2008

Turbulence Model Behavior in Low Reynolds Number Regions of Aerodynamic Flowfields

The behaviors of the widely-used Spalart-Allmaras (SA) and Menter shear-stress transport (SST) turbulence models at low Reynolds numbers and under conditions conducive to relaminarization are documented. The flows used in the investigation include 2-D zero pressure gradient flow over a flat plate from subsonic to hypersonic Mach numbers, 2-D airfoil flow from subsonic to supersonic Mach numbers, 2-D subsonic sink-flow, and 3-D subsonic flow over an infinite swept wing (particularly its leading-edge region). Both models exhibit a range over which they behave transitionally in the sense that the flow is neither laminar nor fully turbulent, but these behaviors are different: the SST model typically has a well-defined transition location, whereas the SA model does not. Both models are predisposed to delayed activation of turbulence with increasing freestream Mach number. Also, both models can be made to achieve earlier activation of turbulence by increasing their freestream levels, but too high a level can disturb the turbulent solution behavior. The technique of maintaining freestream levels of turbulence without decay in the SST model, introduced elsewhere, is shown here to be useful in reducing grid-dependence of the model's transitional behavior. Both models are demonstrated to be incapable of predicting relaminarization; eddy viscosities remain weakly turbulent in accelerating or laterally-strained boundary layers for which experiment and direct simulations indicate turbulence suppression. The main conclusion is that these models are intended for fully turbulent high Reynolds number computations, and using them for transitional (e.g., low Reynolds number) or relaminarizing flows is not appropriate.

Rumsey, Christopher L.↗
NASA NTRS2008

Turbulence Model Behavior in Low Reynolds Number Regions of Aerodynamic Flowfields

The behaviors of the widely-used Spalart-Allmaras (SA) and Menter shear-stress transport (SST) turbulence models at low Reynolds numbers and under conditions conducive to relaminarization are documented. The flows used in the investigation include 2-D zero pressure gradient flow over a flat plate from subsonic to hypersonic Mach numbers, 2-D airfoil flow from subsonic to supersonic Mach numbers, 2-D subsonic sink-flow, and 3-D subsonic flow over an infinite swept wing (particularly its leading-edge region). Both models exhibit a range over which they behave 'transitionally' in the sense that the flow is neither laminar nor fully turbulent, but these behaviors are different: the SST model typically has a well-defined transition location, whereas the SA model does not. Both models are predisposed to delayed activation of turbulence with increasing freestream Mach number. Also, both models can be made to achieve earlier activation of turbulence by increasing their freestream levels, but too high a level can disturb the turbulent solution behavior. The technique of maintaining freestream levels of turbulence without decay in the SST model, introduced elsewhere, is shown here to be useful in reducing grid-dependence of the model's transitional behavior. Both models are demonstrated to be incapable of predicting relaminarization; eddy viscosities remain weakly turbulent in accelerating or laterally-strained boundary layers for which experiment and direct simulations indicate turbulence suppression. The main conclusion is that these models are intended for fully turbulent high Reynolds number computations, and using them for transitional (e.g., low Reynolds number) or relaminarizing flows is not appropriate.

Rumsey, Christopher L.↗
NASA NTRSDate not supplied

High-Lift Common Research Model: RANS, HRLES and WMLES Perspectives for CLmax Prediction Using LAVA

A unified assessment of three turbulence treatments: Reynolds Averaged Navier-Stokes(RANS), Hybrid RANS/LES (HRLES) and Equilibrium Wall-Modelled Large Eddy Simulation(WMLES) is presented for the High-Lift Common Research Model (CRM-HL). For the free-air configuration, steady-state RANS simulations show very accurate drag polar predictions in the low-𝛼linear regime. However, strong grid sensitivity is reported near the maximum lift-state(𝐶𝐿max), with finer-grids showing larger errors and predicting erroneous flow topologies on the wing. Our RANS simulations show that several corrections for the Spalart-Allmaras (SA)turbulence model widely used in the community lead to more erroneous results compared to the baseline closure, without exception. Both scale-resolving methods (HRLES and WMLES)address these drawbacks and predict an outboard separation pattern on the main element that is in good agreement with the oil flow photographs taken from the QinetiQ wind tunnel experiments, when LES-appropriate grids and numerical discretizations are used. While RANS simulations with the baseline SA closure do not show any wing-root separation post𝐶𝐿max, both HRLES and WMLES show onset of corner flow separation with varying degrees of progression, along with a weak pitch break in the wing-contribution of the overall pitching moment. This post-𝐶𝐿max pitch break seen in the free-air simulations is weaker than the break observed in experiments, with a weaker break reported in WMLES for each iteration of grid-refinement. In-tunnel simulations using both SA-baseline RANS and WMLES show a much stronger post-𝐶𝐿max break with the WMLES predictions showing excellent agreement with the experiment in terms of both the flow-topology observed and the pressure-coefficients at various spanwise stations. Sensitivity to the tunnel wall boundary layer is characterized via comparisons between viscous and inviscid treatments for the tunnel walls. WMLES predictions show moderate sensitivity at the predicted inboard flow-state at 𝐶𝐿max along with the progression towards a post-𝐶𝐿max stall; however, this stalled state at 𝛼≈20◦(inside the tunnel) obtained with both tunnel wall treatments appears to be largely identical.

TTT↗
NASA NTRS2015

Entropy Stable Wall Boundary Conditions for the Three-Dimensional Compressible Navier-Stokes Equations

Non-linear entropy stability and a summation-by-parts framework are used to derive entropy stable wall boundary conditions for the three-dimensional compressible Navier-Stokes equations. A semi-discrete entropy estimate for the entire domain is achieved when the new boundary conditions are coupled with an entropy stable discrete interior operator. The data at the boundary are weakly imposed using a penalty flux approach and a simultaneous-approximation-term penalty technique. Although discontinuous spectral collocation operators on unstructured grids are used herein for the purpose of demonstrating their robustness and efficacy, the new boundary conditions are compatible with any diagonal norm summation-by-parts spatial operator, including finite element, finite difference, finite volume, discontinuous Galerkin, and flux reconstruction/correction procedure via reconstruction schemes. The proposed boundary treatment is tested for three-dimensional subsonic and supersonic flows. The numerical computations corroborate the non-linear stability (entropy stability) and accuracy of the boundary conditions.

Parsani, Matteo↗
NASA NTRS2024

Spatiotemporal Evolution of Marine Heatwaves Globally

The spatiotemporal evolution of marine heatwaves (MHWs) is explored using a tracking algorithm termed Ocetrac that provides objective characterization of MHW spatiotemporal evolution. Candidate MHW grid points are defined in detrended gridded sea temperature data using a seasonally varying temperature threshold. Identified MHW points are collected into spatially distinct objects using edge detection with weak sensitivity to edge detection and size threshold criteria. These MHW objects are followed in space and time while allowing objects to split and merge. Ocetrac is applied to monthly satellite sea surface temperature data from September 1981 through January 2021. The resulting MHWs are characterized by their intensity, duration, and total area covered. The global analysis shows that MHWs in the Gulf of Maine and Mediterranean Sea evolve within a relatively small region, while major MHWs in the Pacific and Indian Oceans are linked in space and time. The largest and most long lasting MHW using this method lasts for 60 months from November 2013 to October 2018, encompassing previously identified MHW events including those in the Northeast Pacific (2014-2015), the Tasman Sea (2015-2016, 2017-2018), and the Great Barrier Reef (2016).

spatiotemporal↗
NASA NTRS1992

The Osher scheme for non-equilibrium reacting flows

An extension of the Osher upwind scheme to nonequilibrium reacting flows is presented. Owing to the presence of source terms, the Riemann problem is no longer self-similar and therefore its approximate solution becomes tedious. With simplicity in mind, a linearized approach which avoids an iterative solution is used to define the intermediate states and sonic points. The source terms are treated explicitly. Numerical computations are presented to demonstrate the feasibility, efficiency and accuracy of the proposed method. The test problems include a ZND (Zeldovich-Neumann-Doring) detonation problem for which spurious numerical solutions which propagate at mesh speed have been observed on coarse grids. With the present method, a change of limiter causes the solution to change from the physically correct CJ detonation solution to the spurious weak detonation solution.

Suresh, Ambady↗
NASA NTRS2009

Water Vapor Measurements by Howard University Raman Lidar during the WAVES 2006 Campaign

Retrieval of water vapor mixing ratio using the Howard University Raman Lidar is presented with emphasis on three aspects: i) performance of the lidar against collocated radiosondes and Raman lidar, ii) investigation of the atmospheric state variables when poor agreement between lidar and radiosondes values occurred and iii) a comparison with satellite-based measurements. The measurements were acquired during the Water Vapor Validation Experiment Sondes/Satellites 2006 field campaign. Ensemble averaging of water vapor mixing ratio data from ten night-time comparisons with Vaisala RS92 radiosondes shows on average an agreement within 10 % up to approx. 8 km. A similar analysis of lidar-to-lidar data of over 700 profiles revealed an agreement to within 20 % over the first 7 km (10 % below 4 km). A grid analysis, defined in the temperature - relative humidity space, was developed to characterize the lidar - radiosonde agreement and quantitatively localizes regions of strong and weak correlations as a function of altitude, temperature or relative humidity. Three main regions of weak correlation emerge: i) regions of low relative humidity and low temperature, ii) moderate relative humidity at low temperatures and iii) low relative humidity at moderate temperatures. Comparison of Atmospheric InfraRed Sounder and Tropospheric Emission Sounder satellites retrievals of moisture with that of Howard University Raman Lidar showed a general agreement in the trend but the formers miss a lot of the details in atmospheric structure due to their low resolution. A relative difference of about 20 % is usually found between lidar and satellites measurements.

Adam, M.↗
NASA NTRS1989

Application of a Two-dimensional Unsteady Viscous Analysis Code to a Supersonic Throughflow Fan Stage

The Rai ROTOR1 code for two-dimensional, unsteady viscous flow analysis was applied to a supersonic throughflow fan stage design. The axial Mach number for this fan design increases from 2.0 at the inlet to 2.9 at the outlet. The Rai code uses overlapped O- and H-grids that are appropriately packed. The Rai code was run on a Cray XMP computer; then data postprocessing and graphics were performed to obtain detailed insight into the stage flow. The large rotor wakes uniformly traversed the rotor-stator interface and dispersed as they passed through the stator passage. Only weak blade shock losses were computerd, which supports the design goals. High viscous effects caused large blade wakes and a low fan efficiency. Rai code flow predictions were essentially steady for the rotor, and they compared well with Chima rotor viscous code predictions based on a C-grid of similar density.

Steinke, Ronald J.↗
NASA NTRS1988

On the interpretation of integrated water vapor patterns in midlatitude cyclones derived from the Nimbus 7 scanning multichannel microwave radiometer

The present exploration of methods for the interpretation of the Nimbus 7 satellite's Scanning Multichannel Microwave Radiometer (SMMR) vapor patterns and the ways in which they relate to the dynamical structure of individual midlatitude storms employs gridded meteorological data from the First GARP Global Experiment special observing period in order to calculate diagnostic quantities. SMMR patterns for a storm at a weak stage determined from the diagnostic quantities are compared with SMMR patterns for the storm at a stronger stage. A more complete interpretation of the SMMR patterns emerges from these considerations.

Mcmurdie, Lynn A.↗
NASA NTRS1994

Network synthesis localization of two soft gamma repeaters

We introduce the method of 'network synthesis,' which allows the detection of very weak gamma-ray transient signals in the data of the Ulysses gamma-ray burst (GRB) experiment from repeating sources. It consists of defining a grid of alpha, delta values, and for each BATSE detection of a burst from a soft gamma repeater, predicting the arrival time of the burst at Ulysses and co-adding the Ulysses data rephased so that the burst signals are aligned in time and produce a detectable pulse. We demonstrate that this method identifies the position of the soft repeater SGR 1806-20, and apply it to the repeater B1900+14. We show that the counterpart to this burst source is probably in or in the vicinity of the Galactic supernova remnant G42.8+0.6.

Hurley, Kevin↗
NASA NTRS1996

An Investigation of Cavity Vortex Generators in Supersonic Flow

The purpose of this report is to document the results of experiments performed at the University of Kansas and at the NASA Langley Research Center (LaRC) into the use of shaped cavities to generate vortices in supersonic flow, as well as the progress made in simulating the observed flow using the PAB3D flow solver. The investigation was performed on 18 different cavity configurations installed in a convergent-divergent nozzle at the Jet Exit Facility at the LaRC. Pressure sensitive paint, static-pressure ports, focusing Schliern, and water tunnel flow visualization techniques were used to study the nature of the flow created by these cavities. The results of these investigations revealed that a shaped cavity can generate a pair of counter-rotating streamwise vortices in supersonic flow by creating weak, compression Mach waves and weak shocks. The PAB3D computer program, developed at the LaRC, was used to attempt to reproduce the experimental results. Unfortunately, due to problems with matching the grid blocks, no converged results were obtained. However, intermediate results, as well as a complete definition of the grid matching problems and suggested courses of actions are presented.

Hazlewood, Richard↗