SEARCH · Search NASA

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 235 records · Page 13

All sources NASA NTRS DOE OSTI

NASA NTRS1990

Counterrotating prop-fan simulations which feature a relative-motion multiblock grid decomposition enabling arbitrary time-steps

Improvements are presented of a computer algorithm developed for the time-accurate flow analysis of rotating machines. The flow model is a finite volume method utilizing a high-resolution approximate Riemann solver for interface flux definitions. The numerical scheme is a block LU implicit iterative-refinement method which possesses apparent unconditional stability. Multiblock composite gridding is used to orderly partition the field into a specified arrangement of blocks exhibiting varying degrees of similarity. Block-block relative motion is achieved using local grid distortion to reduce grid skewness and accommodate arbitrary time step selection. A general high-order numerical scheme is applied to satisfy the geometric conservation law. An even-blade-count counterrotating unducted fan configuration is chosen for a computational study comparing solutions resulting from altering parameters such as time step size and iteration count. The solutions are compared with measured data.

Janus, J. Mark↗

NASA NTRS1989

Spectral element simulations of laminar and turbulent flows in complex geometries

Spectral element methods are high-order weighted residual techniques based on spectral expansions of variables and geometry for the Navier-Stokes (NS) and transport equations. Here, practical aspects of these methods and their efficient implementation are examined, and several examples of flows in truly complex geometries are presented. The spectral element discretization for NS equations is introduced, and the convergence of the method is addressed. An efficient data management scheme is discussed in the context of parallel processing computations. The method is validated by comparing the spectral element solutions with the exact eigensolutions for the Orr-Sommerfeld equations in two and three dimensions. Computer-aided flow visualizations are presented for an impulsive flow past a sharp edge wedge. Three-dimensional states of channel flow disrupted by an array of cylindrical eddy promoters are studied, and the results of a direct simulation of the turbulent flow in a plane channel are presented.

Karniadakis, George EM↗

NASA NTRS1990

Multiloop Balanced Bridge Feedback in application to precision pointing

The Balanced Bridge Feedback technique is applied to the multiloop pointing control problem. Using colocated torque and angular velocity sensors and high-order compensators, the motor loop is decoupled from the flexible plant, the feedback bandwidth is increased in the motor and plant loops, and the accuracy improved by orders of magnitude compared to the results achievable with conventional control techniques.

Lurie, B. J.↗

NASA NTRS1990

Estimating short-period dynamics using an extended Kalman filter

An extended Kalman filter is used to estimate the parameters of a low-order model from aircraft transient response data. The low-order model is a state-space model derived from the short-period approximation of the longitudinal aircraft dynamics. The model corresponds to the pitch rate to stick force transfer function currently used in flying qualities analysis. The parameters are estimated from flight data as well as from a 6-DOF nonlinear simulation of the aircraft. These two estimates are then compared, and the discrepancies noted. The low-order model is able to satisfactorily match both flight data and simulation data from a high-order computer simulation.

Bauer, Jeffrey E.↗

NASA NTRS1989

Spectral element methods for the incompressible Navier-Stokes equations

Spectral element methods are high-order weighted-residual techniques for partial differential equations that combine the geometric flexibility of finite element techniques with the rapid convergence rate of spectral schemes. The theoretical foundations and numerical implementation of spectral element methods for the incompressible Navier-Stokes equations are presented, considering the construction and analysis of optimal-order spectral element discretizations for elliptic and saddle (Stokes) problems, as well as the efficient solution of the resulting discrete equations by rapidly convergent tensor-product-based iterative procedures. Several examples of spectral element simulation of moderate Reynolds number unsteady flow in complex geometry are presented.

Maday, Yvon↗

NASA NTRS1990

Supersonic aerodynamic characteristics of a Mach 3 high-speed civil transport configuration

The NASA-Langley AST31 high-speed civil transport (HSCT) configuration, a blended wing/fuselage concept scaled for 250-passenger carriage over 6500 nmi, has been subjected to both CFD and wind tunnel tests; the wind tunnel tests were conducted in the Mach 1.6-3.6 range. The inability of the linear theory method to accurately predict stability levels as well as nonlinear pitching moment characteristics demonstrates the need to undertake high-order analyses as early in the design process as possible. An HSCT laminar flow capability experimental study was conducted in order to simulate 50-percent laminar flow conditions, and thereby ascertain the difference in drag relative to fully turbulent conditions.

Covell, Peter F.↗

NASA NTRS1990

A closer look at the induced drag of crescent-shaped wings

Recent interest in the induced drag characteristics of crescent-shaped wings has led to a closer look at the methods used for determination of induced drag from computational aerodynamic methods. Induced drag may be computed by integration of surface pressure, or by evaluation of a contour integral in the Trefftz plane. A high-order panel method was used to study the induced drag of crescent and elliptical wings using both techniques. Induced drag computations using surface-pressure integration were strongly affected by panel density and angle of attack. Accurate drag computations for the crescent wing were obtained only when the spanwise as well as chordwise panel density was extremely high. Trefftz-plane results for the two wing planforms do not show this sensitivity to panel density or angle of attack. Span efficiencies of 0.994 for the crescent wing and 0.987 for the elliptical wing were computed by the Trefftz-plane technique. Substitution of a force-free, rolled-up wake geometry on the crescent wing did not change the pressure-integrated drag significantly. The slightly higher span efficiency of the crescent wing is attributed to a more nearly elliptical spanwise lift distribution. The chord distribution of the elliptical wing was modified to produce an elliptical span-loading on a wing with an unswept quarter-chord line. This wing demonstrated a span efficiency equal to that of the crescent wing.

Smith, Stephen C.↗

NASA NTRS1989

The determination of third order linear models from a seventh order nonlinear jet engine model

Results are presented that demonstrate how good reduced-order models can be obtained directly by recursive parameter identification using input/output (I/O) data of high-order nonlinear systems. Three different methods of obtaining a third-order linear model from a seventh-order nonlinear turbojet engine model are compared. The first method is to obtain a linear model from the original model and then reduce the linear model by standard reduction techniques such as residualization and balancing. The second method is to identify directly a third-order linear model by recursive least-squares parameter estimation using I/O data of the original model. The third method is to obtain a reduced-order model from the original model and then linearize the reduced model. Frequency responses are used as the performance measure to evaluate the reduced models. The reduced-order models along with their Bode plots are presented for comparison purposes.

Lalonde, Rick J.↗

NASA NTRS1989

Understanding transition and turbulence through direct simulations

Direct simulations consist in solving the full Navier-Stokes equations, without any turbulence model, and describing all the detailed features of the flow. Usually the flows are three-dimensional and time-dependent and contain both coarse and fine structures, which makes the numerical task very challenging in terms of both the algorithm and the computational effort. Most of the work until now has involved spectral methods, which are highly accurate but not very flexible in terms of geometry or complex equations. For that reason, future work will also rely on high-order finite-difference or other methods. Direct simulations complement experimental work, and both contribute to the theory and the empirical knowledge of turbulence. Once such a simulation has been shown to be accurate, the flow field is completely known in three dimensions and time, including the pressure, the vorticity and any other quantity. On the other hand, most simulations to date solved the incompressible equations in rather simple geometries, and direct simulations will always be limited to moderate Reynolds numbers. Extensive simulations have been conducted in homogeneous turbulence, channel flows, boundary layers, and mixing layers. Much effort is devoted to addressing flows with compressibility and chemical reactions, and to new geometries such as a backward-facing step.

Spalart, P. R.↗

NASA NTRS1990

Numerical solution of the incompressible Navier-Stokes equations

The current work is initiated in an effort to obtain an efficient, accurate, and robust algorithm for the numerical solution of the incompressible Navier-Stokes equations in two- and three-dimensional generalized curvilinear coordinates for both steady-state and time-dependent flow problems. This is accomplished with the use of the method of artificial compressibility and a high-order flux-difference splitting technique for the differencing of the convective terms. Time accuracy is obtained in the numerical solutions by subiterating the equations in psuedo-time for each physical time step. The system of equations is solved with a line-relaxation scheme which allows the use of very large pseudo-time steps leading to fast convergence for steady-state problems as well as for the subiterations of time-dependent problems. Numerous laminar test flow problems are computed and presented with a comparison against analytically known solutions or experimental results. These include the flow in a driven cavity, the flow over a backward-facing step, the steady and unsteady flow over a circular cylinder, flow over an oscillating plate, flow through a one-dimensional inviscid channel with oscillating back pressure, the steady-state flow through a square duct with a 90 degree bend, and the flow through an artificial heart configuration with moving boundaries. An adequate comparison with the analytical or experimental results is obtained in all cases. Numerical comparisons of the upwind differencing with central differencing plus artificial dissipation indicates that the upwind differencing provides a much more robust algorithm, which requires significantly less computing time. The time-dependent problems require on the order of 10 to 20 subiterations, indicating that the elliptical nature of the problem does require a substantial amount of computing effort.

Rogers, Stuart E.↗

NASA NTRS1990

A numerical method for parameterization of atmospheric chemistry - Computation of tropospheric OH

An efficient and stable computational scheme for parameterization of atmospheric chemistry is described. The 24-hour-average concentration of OH is represented as a set of high-order polynomials in variables such as temperature, densities of H2O, CO, O3, and NO(t) (defined as NO + NO2 + NO3 + 2N2O5 + HNO2 + HNO4) as well as variables determining solar irradiance: cloud cover, density of the overhead ozone column, surface albedo, latitude, and solar declination. This parameterization of OH chemistry was used in the three-dimensional study of global distribution of CH3CCl3. The proposed computational scheme can be used for parameterization of rates of chemical production and loss or of any other output of a full chemical model.

Spivakovsky, C. M.↗

NASA NTRS1990

Modeling the pilot in visually controlled flight

The simplest model for a human operator is a gain with a time delay. However, there have been no comprehensive studies evaluating human control strategies in visually controlled flight. The results of preliminary studies on this topic are described. Human visually guided flight control is important both in low-level flight, where it predominates, and in higher-altitude flights, where instrument failure is always a potential danger. Two general approaches to this problem, one founded on high-order perceptual psychophysics and the other on control systems engineering, are described. Initial results show that the use of control engineering modeling techniques, together with a psychophysical analysis of information in the perspective scene, holds promise for capturing the manual control strategies used during visual flight.

Johnson, Walter W.↗

NASA NTRS1990

Flux-difference split scheme for turbulent transport equations

This paper describes a high accuracy upwind method to solve the mean compressible flow equations closed with two-equation turbulence models. A flux-difference splitting scheme is developed to compute the flowfield without introducing oscillations into the solution. A high-order accuracy scheme is used to predict both the mean flow variables and the turbulence variables. Results comparing two different two-equation turbulence models are presented for supersonic flow over a flat plate and a compressible, free shear-layer. The non-oscillatory behavior of the upwind scheme is demonstrated on the free shear-layer.

Morrison, J. H.↗

NASA NTRS1990

An approach to optimal guidance of an advanced launch vehicle concept

An approximate solution for the maximum payload trajectory of a two-stage launch vehicle using a regular perturbation technique is presented. A zero-order solution for a two-stage vehicle based on a flat-earth approximation and negligible atmospheric effects is obtained in closed form. High-order correction terms are obtained from the solution of nonhomogeneous, first-order linear differential equations by quadrature. This promises the capability for an onboard optimal guidance law implementation.

Leung, Martin S.↗

NASA NTRS1990

Numerical experiments on the accuracy of ENO and modified ENO schemes

Numerical experiments have been performed using different ENO schemes and different time discretizations, in order to assess accuracy-degeneracy phenomena of the type described by Rogerson and Meiburg (1990). A modified ENO scheme recovers the correct order of accuracy for all the test problems with smooth initial conditions and gives results comparable to the original ENO schemes for discontinuous problems. It is concluded that ENO schemes may lose the full high-order accuracy predicted by local truncation-error analysis for scalar linear conservation laws with smooth initial conditions. A modified ENO scheme can overcome this accuracy-degeneracy problem for the test problems without increasing the computational cost.

Shu, Chi-Wang↗

NASA NTRS1988

Maximum likelihood identification for large space structures

This paper examines the use of on-orbit identification based on Maximum Likelihood Estimation (MLE) to provide these high-order, high-accuracy control design models for large space structures (LSS's). First, it outlines a general MLE identification algorithm, together with a covariance-analysis procedure to assess algorithm performance in terms of systematic and stochastic errors. Next, it examines various simplifications appropriate for the LSS identification application. Simplified analytical performance results are presented, as are numerical results to support these analyses. Finally, a graphical interpretation of these results is given.

Barrett, Michael F.↗

DOE OSTI2025

Coupling Approaches with Non-matching Grids for Classical Linear Elasticity and Bond-based Peridynamic Models in 1D

Local-nonlocal coupling approaches provide a means to combine the computational efficiency of local models and the accuracy of nonlocal models. To facilitate the coupling of the two models, non-matching grids are often desirable as nonlocal grids usually require a finer resolution than local grids. In that case, it is often convenient to resort to interpolation operators so that models can exchange information in the overlap regions when nodes from the two grids do not coincide. This paper studies three existing coupling approaches, namely 1) a method that enforces matching displacements in an overlap region, 2) a variant that enforces a constraint on the stresses instead, and 3) a method that considers a variable horizon in the vicinity of the interfaces. Further, the effect of the interpolation order and of the grid ratio on the performance of the three coupling methods with non-matching grids is carefully studied on one-dimensional examples using polynomial manufactured solutions. The numerical results show that the degree of the interpolants should be chosen with care to avoid introducing additional modeling errors, or simply minimize these errors, in the coupling approach.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

An extension of the localized artificial diffusivity method for immiscible and high density ratio flows

The localized artificial diffusivity (LAD) method is widely regarded as the preferred multi-material regularization scheme for the compact finite difference method, because it is conservative, easy to implement, and generally robust for a wide range of multi-material problems. However, traditional LAD methods face significant challenges when applied to flows with large density ratios and when maintaining thermodynamic equilibrium across material interfaces. These limitations arise from the formulation of the artificial diffusivity flux and the reliance on enthalpy diffusion for interface regularization. Additionally, traditional LAD methods struggle to ensure stability under large density ratio conditions, fail to maintain a finite interface thickness, and are therefore unsuitable for modeling immiscible interfaces. Here, in this work, we discuss the origins of these issues in traditional LAD methods and propose modifications which enable the simulation of large density ratio and immiscible flows. The proposed method targets the artificial diffusion fluxes at gradients and ringing in the volume fraction, rather than the mass fraction in traditional methods, to consistently regularize large density ratio interfaces. Furthermore, the proposed method introduces an artificial bulk density diffusion term to enforce equilibrium conditions across interfaces. To address the challenge of modeling immiscible flows, a conservative diffuse interface term is incorporated into the formulation to ensure a finite interface thickness. Specific consideration is taken in the design of the method to ensure that these crucial properties are maintained for N -material flows. The effectiveness of the proposed method is demonstrated through a series of canonical test cases, and its accuracy is validated by comparison with experimental data on micro-bubble collapse in water. These results highlight the method’s robustness and its ability to overcome the limitations of traditional LAD approaches.

Artificial diffusivity↗