Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “high order”

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 595 records · Page 33

Model reduction by weighted Component Cost Analysis

Component Cost Analysis considers any given system driven by a white noise process as an interconnection of different components, and assigns a metric called 'component cost' to each component. These component costs measure the contribution of each component to a predefined quadratic cost function. A reduced-order model of the given system may be obtained by deleting those components that have the smallest component costs. The theory of Component Cost Analysis is extended to include finite-bandwidth colored noises. The results also apply when actuators have dynamics of their own. Closed-form analytical expressions of component costs are also derived for a mechanical system described by its modal data. This is very useful to compute the modal costs of very high order systems. A numerical example for MINIMAST system is presented.

Kim, Jae H.↗

A normal stress subgrid-scale eddy viscosity model in large eddy simulation

The Smagorinsky subgrid-scale eddy viscosity model (SGS-EVM) is commonly used in large eddy simulations (LES) to represent the effects of the unresolved scales on the resolved scales. This model is known to be limited because its constant must be optimized in different flows, and it must be modified with a damping function to account for near-wall effects. The recent dynamic model is designed to overcome these limitations but is compositionally intensive as compared to the traditional SGS-EVM. In a recent study using direct numerical simulation data, Horiuti has shown that these drawbacks are due mainly to the use of an improper velocity scale in the SGS-EVM. He also proposed the use of the subgrid-scale normal stress as a new velocity scale that was inspired by a high-order anisotropic representation model. The testing of Horiuti, however, was conducted using DNS data from a low Reynolds number channel flow simulation. It was felt that further testing at higher Reynolds numbers and also using different flows (other than wall-bounded shear flows) were necessary steps needed to establish the validity of the new model. This is the primary motivation of the present study. The objective is to test the new model using DNS databases of high Reynolds number channel and fully developed turbulent mixing layer flows. The use of both channel (wall-bounded) and mixing layer flows is important for the development of accurate LES models because these two flows encompass many characteristic features of complex turbulent flows.

Horiuti, K.↗

High-Fidelity Simulations of a Quadrotor Vehicle for Urban Air Mobility

High-fidelity computational fluid dynamics simulations have been carried out to analyze NASA’s quadrotor air taxi concept for urban air mobility. High-order accurate schemes, dual-time stepping, and the delayed detached-eddy simulation model have been employed. The flow solver has been loosely coupled with the rotorcraft comprehensive analysis code. The vehicle simulated is a six-passenger quadrotor for air taxi operations. The quadrotor has been simulated in hover and cruise, with flight conditions representative of an air taxi mission. Initially, the quadcopter vehicle is analyzed using low-fidelity tools. Then, a complete analysis using high-fidelity simulations of the quadcopter is performed in conjunction with collective pitch control. Airloads, wake geometry, and performance are studied and compared. A comparison between isolated rotors and complete vehicle simulations is presented. NASA’s quadrotor is one of NASA’s concept vehicles designed in support of aircraft development for vertical take-off and landing air taxi operations.

ARMD↗

Application of CART3D to Complex Propulsion-Airframe Integration with Vehicle Sketch Pad

Vehicle Sketch Pad (VSP) is an easy-to-use modeler used to generate aircraft geometries for use in conceptual design and analysis. It has been used in the past to generate metageometries for aerodynamic analyses ranging from handbook methods to Navier-Stokes computational fluid dynamics (CFD). As desirable as it is to bring high order analyses, such as CFD, into the conceptual design process, this has been difficult and time consuming in practice due to the manual nature of both surface and volume grid generation. Over the last couple of years, VSP has had a major upgrade of its surface triangulation and export capability. This has enhanced its ability to work with Cart3D, an inviscid, three dimensional fluid flow toolset. The combination of VSP and Cart3D allows performing inviscid CFD on complex geometries with relatively high productivity. This paper will illustrate the use of VSP with Cart3D through an example case of a complex propulsion-airframe integration (PAI) of an over-wing nacelle (OWN) airliner configuration.

Hahn, Andrew S.↗

Head-on parallel blade-vortex interaction

An experimental and computational study was carried out to investigate the parallel head-on blade-vortex interaction (BVI) and its noise generation mechanism. A shock tube, with an enlarged test section, was used to generate a compressible starting vortex which interacted with a target airfoil. The dual-pulsed holographic interferometry (DPHI) technique and airfoil surface pressure measurements were employed to obtain quantitative flow data during the BVI. A thin-layer Navier-Stokes code (BV12D), with a high-order upwind-biased scheme and a multizonal grid, was also used to simulate numerically the phenomena occurring in the head-on BVI. The detailed structure of a convecting vortex was studied through independent measurements of density and pressure distributions across the vortex center. Results indicate that, in a strong head-on BVI, the opposite pressure peaks are generated on both sides of the leading edge as the vortex approaches. Then, as soon as the vortex passes by the leading edge, the high-pressure peak suddenly moves toward the low-peak-reducing in magnitude as it moves--simultaneously giving rise to the initial sound wave. In both experiment and computation, it is shown that the viscous effect plays a significant role in head-on BVIs.

Lee, Soogab↗

Computer extension and analytic continuation of Stokes' expansion for gravity waves

Stokes' infinitesimal-wave expansion for steady progressive free-surface waves has been extended to high order using a computer to perform the coefficient arithmetic. Stokes' expansion has been found to be incapable of yielding the highest wave for any value of the water depth since convergence is limited by a square-root branch-point some distance short of the maximum. By reformulating the problem using a different independent parameter, the highest waves are obtained correctly. Series summation and analytic continuation are facilitated by the use of Pade approximants. The method is valid in principle for any finite value of the wavelength and solutions of high accuracy can be obtained for most values of the wave height and water depth.

Schwartz, L. W.↗

Optical communications and a comparison of optical technologies for a high data rate return link from Mars

The important principles of direct- and heterodyne-detection optical free-space communications are reviewed. Signal-to-noise-ratio (SNR) and bit-error-rate (BER) expressions are derived for both the direct-detection and heterodyne-detection optical receivers. For the heterodyne system, performance degradation resulting from received-signal and local oscillator-beam misalignment and laser phase noise is analyzed. Determination of interfering background power from local and extended background sources is discussed. The BER performance of direct- and heterodyne-detection optical links in the presence of Rayleigh-distributed random pointing and tracking errors is described. Finally, several optical systems employing Nd:YAG, GaAs, and CO2 laser sources are evaluated and compared to assess their feasibility in providing high-data-rate (10- to 1000-Mbps) Mars-to-Earth communications. It is shown that the root mean square (rms) pointing and tracking accuracy is a critical factor in defining the system transmitting laser-power requirements and telescope size and that, for a given rms error, there is an optimum telescope aperture size that minimizes the required power. The results of the analysis conducted indicate that, barring the achievement of extremely small rms pointing and tracking errors (less than 0.2 microrad), the two most promising types of optical systems are those that use an Nd:YAG laser (lambda = 1.064 microns) and high-order pulse position modulator (PPM) and direct detection, and those that use a CO2 laser (lambda = 10.6 microns) and phase shifting keying homodyne modulation and coherent detection. For example, for a PPM order of M = 64 and an rms pointing accuracy of 0.4 microrad, an Nd:YAG system can be used to implement a 100-Mbps Mars link with a 40-cm transmitting telescope, a 20-W laser, and a 10-m receiving photon bucket. Under the same conditions, a CO2 system would require 3-m transmitting and receiving telescopes and a 32-W laser to implement such a link. Other types of optical systems, such as a semiconductor laser systems, are impractical in the presence of large rms pointing errors because of the high power requirements of the 100-Mbps Mars link, even when optimal-size telescopes are used.

Spence, Rodney L.↗

Method and technique for installing light-weight, fragile, high-temperature fiber insulation

A method of installing fragile, light weight, high temperature fiber insulation, particularly where the insulation is to be used as a seal strip providing a high order of thermal barrier insulation is discussed. The process is based on provision of a strip of the mineral batting cut oversize by a predetermined amount, saturated in a fugitive polymer solution, compressed in a mold, dried and cured to form a rigidized batting material which is machined to required shape. The machine dimensions would normally be at least nominally less than the dimensions of the cavity to be sealed. After insertion in the cavity, which may be a wire-mesh seal enclosure, the apparatus is subjected to baking at a temperature sufficiently high to cause the resin to burn off cleanly, leaving the batting substantially in its original condition and expanded into the cavity or seal enclosure.

Patel, B. C.↗

A method and technique for installing light-weight fragile, high-temperature fiber insulation

A method of installing fragile, light-weight, high-temperature fiber insulation, particularly where the insulation is to be used as a seal strip providing a high order of thermal barrier insulation is described. The process is based on provision of a strip of the mineral batting cut oversize by a predetermined amount, saturated in a fugitive polymer solution, compressed in a mold, dried and cured to form a rigidized batting material which may be machined to required shape. The machined dimensions would normally be at least nominally less than the dimensions of the cavity to be sealed. After insertion in the cavity, which may be a wire-mesh seal enclosure, the apparatus is subjected to baking at a temperature sufficiently high to cause the resin to burn off cleanly, leaving the batting substantially in its original condition and expanded into the cavity or seal enclosure.

Ballantine, T. J.↗

A time accurate finite volume high resolution scheme for three dimensional Navier-Stokes equations

A time accurate, three-dimensional, finite volume, high resolution scheme for solving the compressible full Navier-Stokes equations is presented. The present derivation is based on the upwind split formulas, specifically with the application of Roe's (1981) flux difference splitting. A high-order accurate (up to the third order) upwind interpolation formula for the inviscid terms is derived to account for nonuniform meshes. For the viscous terms, discretizations consistent with the finite volume concept are described. A variant of second-order time accurate method is proposed that utilizes identical procedures in both the predictor and corrector steps. Avoiding the definition of midpoint gives a consistent and easy procedure, in the framework of finite volume discretization, for treating viscous transport terms in the curvilinear coordinates. For the boundary cells, a new treatment is introduced that not only avoids the use of 'ghost cells' and the associated problems, but also satisfies the tangency conditions exactly and allows easy definition of viscous transport terms at the first interface next to the boundary cells. Numerical tests of steady and unsteady high speed flows show that the present scheme gives accurate solutions.

Liou, Meng-Sing↗

Analysis and design of numerical schemes for gas dynamics 1: Artificial diffusion, upwind biasing, limiters and their effect on accuracy and multigrid convergence

The theory of non-oscillatory scalar schemes is developed in this paper in terms of the local extremum diminishing (LED) principle that maxima should not increase and minima should not decrease. This principle can be used for multi-dimensional problems on both structured and unstructured meshes, while it is equivalent to the total variation diminishing (TVD) principle for one-dimensional problems. A new formulation of symmetric limited positive (SLIP) schemes is presented, which can be generalized to produce schemes with arbitrary high order of accuracy in regions where the solution contains no extrema, and which can also be implemented on multi-dimensional unstructured meshes. Systems of equations lead to waves traveling with distinct speeds and possibly in opposite directions. Alternative treatments using characteristic splitting and scalar diffusive fluxes are examined, together with modification of the scalar diffusion through the addition of pressure differences to the momentum equations to produce full upwinding in supersonic flow. This convective upwind and split pressure (CUSP) scheme exhibits very rapid convergence in multigrid calculations of transonic flow, and provides excellent shock resolution at very high Mach numbers.

Jameson, Antony↗

Further Development of a New, Flux-Conserving Newton Scheme for the Navier-Stokes Equations

This paper is one of a series of papers describing the development of a new numerical approach for solving the steady Navier-Stokes equations. The key features in the current development are (1) the discrete representation of the dependent variables by way of high order polynomial expansions, (2) the retention of all derivatives in the expansions as unknowns to be explicitly solved for, (3) the automatic balancing of fluxes at cell interfaces, and (4) the discrete simulation of both the integral and differential forms of the governing equations. The main purpose of this paper is, first, to provide a systematic and rigorous derivation of the conditions that are used to simulate the differential form of the Navier-Stokes equations, and second, to extend our previously-presented internal flow scheme to external flows and nonuniform grids. Numerical results are presented for high Reynolds number flow (Re = 100,000) around a finite flat plate, and detailed comparisons are made with the Blasius flat plate solution and Goldstein wake solution. It is shown that the error in the streamwise velocity decreases like r(sup alpha)(Delta)y(exp 2), where alpha approx. 0.25 and r = delta(y)/delta(x) is the grid aspect ratio.

Scott, James R.↗

3-D Structure of Sunspots using Imaging Spectroscopy

We use the Interferometric BIdimensional Spectrometer (IBIS) of the INAF/Arcetri Astrophysical Observatory and installed at the National Solar Observatory (NSO) Dunn Solar Telescope, to understand the structure of sunspots. Using the spectral lines FeI 6301.5 A, FeII 7224.4 A, and CaII 8542.6 A, we examine the spectroscopic variation of sunspot penumbral and umbral structures at the heights of formation of these lines. These high resolution observations were acquired on 2004 July 30-31, of active region NOAA 10654, using the high order NSO adaptive optics system. We map the spatio-temporal variation of Doppler signatures in these spectral lines, from the photosphere to the chromosphere. From a 70-minute temporal average of individual 32-second cadence Doppler observations we find that the averaged velocities decrease with height, about 3.5 times larger in the deeper photosphere (FeII 7224.4 A; height-of-formation approx. 50 km) than in the upper photosphere FeI 6301.5 A; height-of-formation approx. 350 km), There is a remarkable coherence of Doppler signals over the height difference of 300 km. From a high-speed animation of the Doppler sequence we find evidence for what appears to be ejection of high speed gas concentrations from edges of penumbral filaments into the surrounding granular photosphere. The Evershed flow persists a few arcseconds beyond the traditionally demarcated penumbra-granulation boundary. We present these and other results and discuss the implications of these measurements for sunspot models.

Balasubramaniam, K. S.↗

Inhomogeneous conditions at open boundaries for wave propagation problems

Absorbing boundary conditions contain differential operators even for first-order systems. There is a fundamental difficulty with this, since the conditions are applied on the ingoing variables, and the approximations necessarily become weakly unstable. This difficulty is more pronounced for inhomogeneous boundary conditions which occur if there is a source outside the computational domain D, or if the initial data are nonzero outside D. In this paper this is further investigated and it is shown that reasonable estimates can still be obtained if the solution is smooth. However, it is demonstrated that the approximations are less robust. A cure for this is proposed and an implementation of high-order conditions for first-order systems is described.

Gustafsson, Bertil↗

Control of large flexible spacecraft by the independent modal-space control method

The problem of control of a large-order flexible structure in the form of a plate-like lattice by the Independent Modal-Space Control (IMSC) method is presented. The equations of motion are first transformed to the modal space, thus obtaining internal (plant) decoupling of the system. Then, the control laws are designed in the modal space for each mode separately, so that the modal equations of motion are rendered externally (controller) decoupled. This complete decoupling applies both to rigid-body modes and elastic modes. The application of linear optimal control, in conjunction with a quadratic performance index, is first reviewed. A solution for high-order systems is proposed here by the IMSC method, whereby the problem is reduced to a number of modal minimum-fuel problems for the controlled modes.

Meirovitch, L.↗

Advances in Rotor Performance and Turbulent Wake Simulation Using DES and Adaptive Mesh Refinement

Time-dependent Navier-Stokes simulations have been carried out for a rigid V22 rotor in hover, and a flexible UH-60A rotor in forward flight. Emphasis is placed on understanding and characterizing the effects of high-order spatial differencing, grid resolution, and Spalart-Allmaras (SA) detached eddy simulation (DES) in predicting the rotor figure of merit (FM) and resolving the turbulent rotor wake. The FM was accurately predicted within experimental error using SA-DES. Moreover, a new adaptive mesh refinement (AMR) procedure revealed a complex and more realistic turbulent rotor wake, including the formation of turbulent structures resembling vortical worms. Time-dependent flow visualization played a crucial role in understanding the physical mechanisms involved in these complex viscous flows. The predicted vortex core growth with wake age was in good agreement with experiment. High-resolution wakes for the UH-60A in forward flight exhibited complex turbulent interactions and turbulent worms, similar to the V22. The normal force and pitching moment coefficients were in good agreement with flight-test data.

Chaderjian, Neal M.↗

Tetrahedral-Mesh Simulation of Turbulent Flows with the Space-Time Conservative Schemes

Direct numerical simulations of turbulent flows are predominantly carried out using structured, hexahedral meshes despite decades of development in unstructured mesh methods. Tetrahedral meshes offer ease of mesh generation around complex geometries and the potential of an orientation free grid that would provide un-biased small-scale dissipation and more accurate intermediate scale solutions. However, due to the lack of consistent multi-dimensional numerical formulations in conventional schemes for triangular and tetrahedral meshes at the cell interfaces, numerical issues exist when flow discontinuities or stagnation regions are present. The space-time conservative conservation element solution element (CESE) method - due to its Riemann-solver-free shock capturing capabilities, non-dissipative baseline schemes, and flux conservation in time as well as space - has the potential to more accurately simulate turbulent flows using unstructured tetrahedral meshes. To pave the way towards accurate simulation of shock/turbulent boundary-layer interaction, a series of wave and shock interaction benchmark problems that increase in complexity, are computed in this paper with triangular/tetrahedral meshes. Preliminary computations for the normal shock/turbulence interactions are carried out with a relatively coarse mesh, by direct numerical simulations standards, in order to assess other effects such as boundary conditions and the necessity of a buffer domain. The results indicate that qualitative agreement with previous studies can be obtained for flows where, strong shocks co-exist along with unsteady waves that display a broad range of scales, with a relatively compact computational domain and less stringent requirements for grid clustering near the shock. With the space-time conservation properties, stable solutions without any spurious wave reflections can be obtained without a need for buffer domains near the outflow/farfield boundaries. Computational results for the isotropic turbulent flow decay, at a relatively high turbulent Mach number, show a nicely behaved spectral decay rate for medium to high wave numbers. The high-order CESE schemes offer very robust solutions even with the presence of strong shocks or widespread shocklets. The explicit formulation in conjunction with a close to unity theoretical upper Courant number bound has the potential to offer an efficient numerical framework for general compressible turbulent flow simulations with unstructured meshes.

Chang, Chau-Lyan↗

An analytic modeling and system identification study of rotor/fuselage dynamics at hover

A combination of analytic modeling and system identification methods have been used to develop an improved dynamic model describing the response of articulated rotor helicopters to control inputs. A high-order linearized model of coupled rotor/body dynamics including flap and lag degrees of freedom and inflow dynamics with literal coefficients is compared to flight test data from single rotor helicopters in the near hover trim condition. The identification problem was formulated using the maximum likelihood function in the time domain. The dynamic model with literal coefficients was used to generate the model states, and the model was parametrized in terms of physical constants of the aircraft rather than the stability derivatives resulting in a significant reduction in the number of quantities to be identified. The likelihood function was optimized using the genetic algorithm approach. This method proved highly effective in producing an estimated model from flight test data which included coupled fuselage/rotor dynamics. Using this approach it has been shown that blade flexibility is a significant contributing factor to the discrepancies between theory and experiment shown in previous studies. Addition of flexible modes, properly incorporating the constraint due to the lag dampers, results in excellent agreement between flight test and theory, especially in the high frequency range.

Hong, Steven W.↗