SEARCH · Engineering Papers

Results for “numerical solution”

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 541 records · Page 30

All sources NASA NTRS DOE OSTI

NASA NTRS1991

Navier-Stokes computation of compressible turbulent flows with a second order closure

The objective was the development of a complete second order closure for wall bounded flows, including all components of the dissipation rate tensor and a numerical solution procedure for the resulting system of equations. The main topics discussed are the closure of the pressure correlations and the viscous destruction terms in the dissipation rate equations and the numerical solution scheme based on a block-tridiagonal solver for the nine equations required for the prediction of plane or axisymmetric flows.

Dingus, C.↗

NASA NTRS1980

Flow oscillations of spike-tipped bodies

The oscillating flow field around a spike-tipped body has been simulated by the numerical solution of the 3-D compressible Navier-Stokes equations at a Mach number of 3.0 and a nominal Reynolds number of 7.87 x 10 to the 6th/m. Computations were performed using a vectorized 3-D Navier-Stokes program on the STAR 100 computer. Numerical solutions confirmed the experimental result that the self-sustained oscillation occurs within a limited range of the protruded spike length to shoulder height ratio. The numerical result predicted correctly the discrete frequency range as well as the rms pressure intensity. The detailed flow structure is also presented and discussed.

Shang, J. S.↗

NASA NTRS1979

A split finite element algorithm for the compressible Navier-Stokes equations

An accurate and efficient numerical solution algorithm is established for solution of the high Reynolds number limit of the Navier-Stokes equations governing the multidimensional flow of a compressible essentially inviscid fluid. Finite element interpolation theory is used within a dissipative formulation established using Galerkin criteria within the Method of Weighted Residuals. An implicit iterative solution algorithm is developed, employing tensor product bases within a fractional steps integration procedure, that significantly enhances solution economy concurrent with sharply reduced computer hardware demands. The algorithm is evaluated for resolution of steep field gradients and coarse grid accuracy using both linear and quadratic tensor product interpolation bases. Numerical solutions for linear and nonlinear, one, two and three dimensional examples confirm and extend the linearized theoretical analyses, and results are compared to competitive finite difference derived algorithms.

Baker, A. J.↗

NASA NTRS1990

A numerical method for phase-change problems

A highly accurate and efficient finite-difference method for phase-change problems with multiple moving boundaries of irregular shape is developed by employing a coordinate transformation that immobilizes moving boundaries and preserves the conservative forms of the original governing equations. The numerical method is first presented for one-dimensional phase-change problems (involving large density variation between phases, heat generation, and multiple moving boundaries) and then extended to solve two-dimensional problems (without change of densities between phases). Numerical solutions are obtained non-iteratively using an explicit treatment of the interfacial mass and energy balances and an implicit treatment of the temperature field equations. The accuracy and flexibility of the present numerical method are verified by solving some phase-change problems and comparing the results with existing analytical, semi-analytical and numerical solutions. Results indicate that one- and two-dimensional phase-change problems can be handled easily with excellent accuracies.

Kim, Charn-Jung↗

NASA NTRS2011

Closed-Form and Numerically-Stable Solutions to Problems Related to the Optimal Two-Impulse Transfer Between Specified Terminal States of Keplerian Orbits

The first part of the paper presents some closed-form solutions to the optimal two-impulse transfer between fixed position and velocity vectors on Keplerian orbits when some constraints are imposed on the magnitude of the initial and final impulses. Additionally, a numerically-stable gradient-free algorithm with guaranteed convergence is presented for the minimum delta-v two-impulse transfer. In the second part of the paper, cooperative bargaining theory is used to solve some two-impulse transfer problems when the initial and final impulses are carried by different vehicles or when the goal is to minimize the delta-v and the time-of-flight at the same time.

Senent, Juan↗

NASA NTRS1983

Finite element or Galerkin type semidiscrete schemes

A finite element of Galerkin type semidiscrete method is proposed for numerical solution of a linear hyperbolic partial differential equation. The question of stability is reduced to the stability of a system of ordinary differential equations for which Dahlquist theory applied. Results of separating the part of numerical solution which causes the spurious oscillation near shock-like response of semidiscrete scheme to a step function initial condition are presented. In general all methods produce such oscillatory overshoots on either side of shocks. This overshoot pathology, which displays a behavior similar to Gibb's phenomena of Fourier series, is explained on the basis of dispersion of separated Fourier components which relies on linearized theory to be satisfactory. Expository results represented.

Durgun, K.↗

NASA NTRS1981

On a fourth order accurate implicit finite difference scheme for hyperbolic conservation laws. I - Nonstiff strongly dynamic problems

An implicit finite difference method of fourth order accuracy in space and time is introduced for the numerical solution of one-dimensional systems of hyperbolic conservation laws. The basic form of the method is a two-level scheme which is unconditionally stable and nondissipative. The scheme uses only three mesh points at level t and three mesh points at level t + delta t. The dissipative version of the basic method given is conditionally stable under the CFL (Courant-Friedrichs-Lewy) condition. This version is particularly useful for the numerical solution of problems with strong but nonstiff dynamic features, where the CFL restriction is reasonable on accuracy grounds. Numerical results are provided to illustrate properties of the proposed method.

Harten, A.↗

NASA NTRS1998

Preconditioning the Helmholtz Equation for Rigid Ducts

An innovative hyperbolic preconditioning technique is developed for the numerical solution of the Helmholtz equation which governs acoustic propagation in ducts. Two pseudo-time parameters are used to produce an explicit iterative finite difference scheme. This scheme eliminates the large matrix storage requirements normally associated with numerical solutions to the Helmholtz equation. The solution procedure is very fast when compared to other transient and steady methods. Optimization and an error analysis of the preconditioning factors are present. For validation, the method is applied to sound propagation in a 2D semi-infinite hard wall duct.

Baumeister, Kenneth J.↗

DOE OSTI2023

A fast Fourier transform-based solver for elastic micropolar composites

This work presents a spectral micromechanical formulation for obtaining the full-field and homogenized response of elastic micropolar composites. The algorithm relies on a coupled set of convolution integral equations for the micropolar strains, where periodic Green’s operators associated with a linear homogeneous reference medium are convolved with functions of the Cauchy and couple stress fields that encode the material’s heterogeneity, as well as any potential material nonlinearity. Such convolution integral equations take an algebraic form in the reciprocal Fourier space that can be solved iteratively. In this vein, the fast Fourier transform (FFT) algorithm is leveraged to accelerate the numerical solution, resulting in a mesh-free formulation in which the periodic unit cell representing the heterogeneous material can be discretized by a regular grid of pixels in two dimensions (or voxels in three dimensions). For verification, the numerical solutions obtained with the micropolar FFT solver are compared with analytical solutions for a matrix with a dilute circular inclusion subjected to plane strain loading. The developed computational framework is then used to study length-scale effects and effective (micropolar) moduli of composites with various topological configurations.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1991

Improved visualization of flow field measurements

A capability was developed that makes it possible to apply to measured flow field data the visualization tools developed to display numerical solutions for computational fluid dynamic problems. The measurement monitor surface (MMS) procedure was applied to the analysis of flow field measurements within a low aspect ratio transonic axial flow fan rotor obtained with 2-D laser anemometry. The procedure generates input for the visualization tools developed to display numerical solutions for computational fluid dynamics problems. The relative Mach number contour plots obtained by this method resemble the conventional contour plots obtained by more traditional methods. The results show that the MMS procedure can be used to generate input for the multidimensional processing and analysis tools developed for data from numerical flow field simulations. They show that an experimenter can apply the MMS procedure to his data and then use an interactive graphics program to display scalar quantities like the Mach number by profiles, carpet plots, contour lines, and surfaces using various colors. Also, flow directionality can be shown by display of vector fields and particle traces.

Miles, Jeffrey Hilton↗

NASA NTRS1973

Unsteady, combined radiation and conduction in an absorbing, scattering, and emitting medium.

The transient cooldown of a gray, absorbing, isotropic scattering, emitting, and conducting medium bounded by gray, diffusely emitting and reflecting parallel plates is considered. Numerical solutions are obtained for the initial boundary-value problem with a discontinuous decrease in temperature at one boundary. The quasi-steady equation of radiative transfer is solved using Gaussian quadrature and a matrix eigenvector technique together with explicit numerical solution of the unsteady energy equation. Temperature and energy flux distributions are presented for variations of optical thickness, boundary emissivity, albedo, and conduction-radiation parameter.

Weston, K. C.↗

NASA NTRS1983

Real gas flow fields about three dimensional configurations

Real gas, inviscid supersonic flow fields over a three-dimensional configuration are determined using a factored implicit algorithm. Air in chemical equilibrium is considered and its local thermodynamic properties are computed by an equilibrium composition method. Numerical solutions are presented for both real and ideal gases at three different Mach numbers and at two different altitudes. Selected results are illustrated by contour plots and are also tabulated for future reference. Results obtained compare well with existing tabulated numerical solutions and hence validate the solution technique.

Balakrishnan, A.↗

NASA NTRS2005

Benchmark Problems Used to Assess Computational Aeroacoustics Codes

The field of computational aeroacoustics (CAA) encompasses numerical techniques for calculating all aspects of sound generation and propagation in air directly from fundamental governing equations. Aeroacoustic problems typically involve flow-generated noise, with and without the presence of a solid surface, and the propagation of the sound to a receiver far away from the noise source. It is a challenge to obtain accurate numerical solutions to these problems. The NASA Glenn Research Center has been at the forefront in developing and promoting the development of CAA techniques and methodologies for computing the noise generated by aircraft propulsion systems. To assess the technological advancement of CAA, Glenn, in cooperation with the Ohio Aerospace Institute and the AeroAcoustics Research Consortium, organized and hosted the Fourth CAA Workshop on Benchmark Problems. Participants from industry and academia from both the United States and abroad joined to present and discuss solutions to benchmark problems. These demonstrated technical progress ranging from the basic challenges to accurate CAA calculations to the solution of CAA problems of increasing complexity and difficulty. The results are documented in the proceedings of the workshop. Problems were solved in five categories. In three of the five categories, exact solutions were available for comparison with CAA results. A fourth category of problems representing sound generation from either a single airfoil or a blade row interacting with a gust (i.e., problems relevant to fan noise) had approximate analytical or completely numerical solutions. The fifth category of problems involved sound generation in a viscous flow. In this case, the CAA results were compared with experimental data.

Dahl, Milo D.↗

NASA NTRS2019

Comparative Study of Active Flow Control Strategies for Lift Enhancement of a Simplified High-Lift Configuration

Numerical simulations have been performed for a simplified high-lift (SHL) version of the Common Research Model (CRM) configuration, where the Fowler flaps of the conventional high-lift (CRM-HL) configuration are replaced by a set of simple hinged flaps. These hinged flaps are equipped with integrated modular active flow control (AFC) cartridges on the suction surface, and the resulting geometry is known as the CRM-SHL-AFC configuration. The main objective is to make use of AFC devices on the CRM-SHL-AFC configuration to recover the aerodynamic performance (lift) of the CRM-HL configuration. In the current paper, a Lattice Boltzmann method-based computational fluid dynamics (CFD) code, known as PowerFLOWQ® is used to simulate the entire flow field associated with the CRM-SHL-AFC configuration equipped with several different types of AFC devices. The transonic version of the PowerFLOWQ® code that has been validated for high speed flows is used to accurately simulate the flow field generated by the high-momentum actuators required to mitigate reversed flow regions on the suction surfaces of the main wing and the flap. The numerical solutions predict the expected trends in aerodynamic forces as the actuation levels are increased. More efficient AFC systems and actuator arrangements emerged based on the parametric studies performed prior to a Fall 2018 wind tunnel test. Preliminary comparisons of the numerical solutions for lift and surface pressures are presented here with the experimental data, demonstrating the usefulness of CFD for predicting the flow field and lift characteristics of AFC-enabled high-lift configurations.

Vatsa, Veer N.↗

NASA NTRS2000

Slow Crack Growth Analysis of Advanced Structural Ceramics Under Combined Loading Conditions: Damage Assessment in Life Prediction Testing

Slow crack growth analysis was performed with three different loading histories including constant stress-rate/constant stress-rate testing (Case I loading), constant stress/constant stress-rate testing (Case II loading), and cyclic stress/constant stress-rate testing (Case III loading). Strength degradation due to slow crack growth arid/or damage accumulation was determined numerically as a Function of percentage of interruption time between the two loading sequences for a given loading history. The numerical solutions were examined with the experimental data determined at elevated temperatures using four different advanced ceramic materials, two silicon nitrides, one silicon carbide and one alumina for the Case I loading history, and alumina for the Case II loading history. The numerical solutions were in reasonable agreement with the experimental data, indicating that notwithstanding some degree of creep deformation presented for some test materials slow crack growth was a governing mechanism associated with failure for all the test material&

Choi, Sung R.↗

NASA NTRS2000

Slow Crack Growth Analysis of Advanced Structural Ceramics Under Combined Loading Conditions: Damage Assessment in Life Prediction Testing

Slow crack growth analysis was performed with three different loading histories including constant stress-rate/constant stress-rate testing (Case 1 loading), constant stress/constant stress-rate testing (Case 2 loading), and cyclic stress/constant stress-rate testing (Case 2 loading). Strength degradation due to slow crack growth and/or damage accumulation was determined numerically as a function of percentage of interruption time between the two loading sequences for a given loading history. The numerical solutions were examined with the experimental data determined at elevated temperatures using four different advanced ceramic materials, two silicon nitrides, one silicon carbide and one alumina for the Case 1 loading history, and alumina for the Case 3 loading history. The numerical solutions were in reasonable agreement with the experimental data, indicating that notwithstanding some degree of creep deformation presented for some test materials slow crack growth was a governing mechanism associated with failure for all the test materials.

Choi, Sung R.↗

NASA NTRS2001

Slow Crack Growth Analysis of Advanced Structural Ceramics Under Combined Loading Conditions: Damage Assessment in Life Prediction Testing

Slow crack growth analysis was performed with three different loading histories including constant stress- rate/constant stress-rate testing (Case I loading), constant stress/constant stress-rate testing (Case II loading), and cyclic stress/constant stress-rate testing (Case III loading). Strength degradation due to slow crack growth and/or damage accumulation was determined numerically as a function of percentage of interruption time between the two loading sequences for a given loading history. The numerical solutions were examined with the experimental data determined at elevated temperatures using four different advanced ceramic materials, two silicon nitrides, one silicon carbide and one alumina for the Case I loading history, and alumina for the Case II loading history. The numerical solutions were in reasonable agreement with the experimental data, indicating that notwithstanding some degree of creep deformation presented for some test materials slow crack growth was a governing mechanism associated with failure for all the rest materials.

Choi, S. R.↗

NASA NTRS2018

Numerical Simulation of a Simplified High-Lift CRM Configuration Embedded with Fluidic Actuators

Numerical simulations have been performed for a simplified high-lift configuration that is representative of a modern transport airplane. This configuration includes a leading-edge slat, fuselage, wing, nacelle-pylon and a simple hinged flap. The suction surface of the flap is embedded with multiple rows of fluidic actuators to reduce the extent of reversed flow regions and improve the aerodynamic performance of the configuration with flap in a deployed state. In the current paper, a Lattice Boltzmann Method based high-fidelity computational fluid dynamics (CFD) code, known as PowerFLOW® is used to simulate the entire flow field associated with this configuration, including the flow inside the actuators. A fully compressible version of the PowerFLOW® code that has been validated for high speed flows is used for the present simulations to accurately represent the transonic flow regimes that are encountered in the flow field generated by the actuators operating at higher mass flow (momentum) rates required to mitigate reverse flow regions on the suction surfaces of the main wing and the flap. The numerical solutions predict the expected trends in aerodynamic forces as the actuation levels are increased. More efficient active flow control (AFC) systems and actuator arrangement for lift augmentation are emerging based on the parametric studies conducted here prior to wind tunnel tests. These numerical solutions will be compared with experimental data, once such data becomes available.

Vatsa, Veer N.↗