SEARCH · Engineering Papers

Results for “finite difference methods”

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 1,603 records · Page 89

All sources NASA NTRS DOE OSTI

NASA NTRS2017

A Requirements-Driven Optimization Method for Acoustic Liners Using Analytic Derivatives

More than ever, there is flexibility and freedom in acoustic liner design. Subject to practical considerations, liner design variables may be manipulated to achieve a target attenuation spectrum. But characteristics of the ideal attenuation spectrum can be difficult to know. Many multidisciplinary system effects govern how engine noise sources contribute to community noise. Given a hardwall fan noise source to be suppressed, and using an analytical certification noise model to compute a community noise measure of merit, the optimal attenuation spectrum can be derived using multidisciplinary systems analysis methods. In a previous paper on this subject, a method deriving the ideal target attenuation spectrum that minimizes noise perceived by observers on the ground was described. A simple code-wrapping approach was used to evaluate a community noise objective function for an external optimizer. Gradients were evaluated using a finite difference formula. The subject of this paper is an application of analytic derivatives that supply precise gradients to an optimization process. Analytic derivatives improve the efficiency and accuracy of gradient-based optimization methods and allow consideration of more design variables. In addition, the benefit of variable impedance liners is explored using a multi-objective optimization.

Aircraft Noise↗

NASA NTRS1971

Mixing of supersonic jets including the effects of supersonic jets including the effects of transverse pressure gradient using difference methods

The usual boundary-layer equations describing the steady-state mixing of parallel jets are supplemented by the momentum equation in the direction normal to the flow. This allows detailed computation of the flow field in the mixing region and simultaneous computation of the outer inviscid flow. An explicit and an implicit finite difference scheme have been developed and applied in several illustrative examples. The examples include mixing of planar and axisymmetric supersonic jets of different composition with both matched and unmatched static pressures. Numerical results were compared with available experimental data obtained for the unmatched pressure case.

Kurkov, A. P.↗

NASA NTRS1995

A finite element propagation model for extracting normal incidence impedance in nonprogressive acoustic wave fields

A propagation model method for extracting the normal incidence impedance of an acoustic material installed as a finite length segment in a wall of a duct carrying a nonprogressive wave field is presented. The method recasts the determination of the unknown impedance as the minimization of the normalized wall pressure error function. A finite element propagation model is combined with a coarse/fine grid impedance plane search technique to extract the impedance of the material. Results are presented for three different materials for which the impedance is known. For each material, the input data required for the prediction scheme was computed from modal theory and then contaminated by random error. The finite element method reproduces the known impedance of each material almost exactly for random errors typical of those found in many measurement environments. Thus, the method developed here provides a means for determining the impedance of materials in a nonprogressirve wave environment such as that usually encountered in a commercial aircraft engine and most laboratory settings.

Watson, Willie R.↗

NASA NTRS1976

Transonic solutions of the Euler equations by the finite volume method

An investigation is conducted of the time-accurate convergence of representative transonic flows to a steady state under given constraints of time-accuracy. Factored explicit and implicit difference operators are used to accelerate the calculations. Attention is given to flow at a Mach number of 1.35 past a circular cylinder and supersonic flow past a NACA 0012 airfoil for three different Mach numbers. Questions of transonic wave behavior are considered along with the equations of motion and the characteristics of the mesh network.

Rizzi, A.↗

NASA NTRS1997

Hierarchical Parallelism in Finite Difference Analysis of Heat Conduction

Based on the concept of hierarchical parallelism, this research effort resulted in highly efficient parallel solution strategies for very large scale heat conduction problems. Overall, the method of hierarchical parallelism involves the partitioning of thermal models into several substructured levels wherein an optimal balance into various associated bandwidths is achieved. The details are described in this report. Overall, the report is organized into two parts. Part 1 describes the parallel modelling methodology and associated multilevel direct, iterative and mixed solution schemes. Part 2 establishes both the formal and computational properties of the scheme.

Padovan, Joseph↗

NASA NTRS1973

Decay of Far-Flowfield in Trailing Vortices

Methods for reduction of velocities in trailing vortices of large aircraft are of current interest for the purpose of shortening the waiting time between landings at central airports. We have made finite-difference calculations of the flow in turbulent wake vortices as an aid to interpretation of wind-tunnel and flight experiments directed toward that end. Finite-difference solutions are capable of adding flexibility to such investigations if they are based on an adequate model of turbulence. Interesting developments have been taking place in the knowledge of turbulence that may lead to a complete theory in the future. In the meantime, approximate methods that yield reasonable agreement with experiment are appropriate. The simplified turbulence model we have selected contains features that account for the major effects disclosed by more sophisticated models in which the parameters are not yet established. Several puzzles are thereby resolved that arose in previous theoretical investigations of wake vortices.

Baldwin, B. S.↗

NASA NTRSDate not supplied

Discretization Error Estimation and Control for Farfield Acoustic Signatures

We investigate the utility of adjoint-based error estimates for sonic boom farfield simulations governed by solutions of the augmented Burgers’ equation. Solution of this nonlinear system uses operator splitting with a second-order finite volume discretization in space and second-order Runge-Kutta time marching, while the absorption and molecular relaxation are solved using second-order central differencing. The discretization error in selected ground sonic boom cost functionals is estimated using the method of adjoint-weighted residuals. Key elements of the implementation process are emphasized with details provided on the practical aspects as appliedto the sonic boom farfield propagation. We establish the accuracy of the adjoint solutions usingcomplex step and finite difference approaches, and examine the accuracy of the error estimates using analytical N-wave solutions. We then apply it to a pressure waveform corresponding to the X-59 research aircraft. The investigations demonstrate that the method of adjoint-weighted residuals accurately predicts the level of discretization error present in sonic boom farfield simulations while offering insight into which features of the near field signal are the primary drivers of ground noise metrics. The numerical results indicate that at sampling frequencies as low as50kHz, discretization error in the propagation is under 0.01 dB[A] for realistically complex examples.

CST↗

NASA NTRS1988

An explicit Runge-Kutta method for 3D turbulent incompressible flows

A computer code has been developed to solve for the steady-state solution of the 3D incompressible Reynolds-averaged Navier-Stokes equations. The approach is based on the cell-center, central-difference, finite-volume formulation and an explicit one-step, multistage Runge-Kutta time-stepping scheme. The Baldwin-Lomax turbulence model is used. Techniques to accelerate the rate of convergence to a steady-state solution include the preconditioned method, the local time stepping, and the implicit residual smoothing. Improvements in computational efficiency have been demonstrated in several areas. This numerical procedure has been used to simulate the turbulent horseshoe vortex flow around an airfoil/flat-plate juncture.

Sung, Chao-Ho↗

NASA NTRS2024

An Experimental and Predictive Study of Bent-Chamber Acoustic Liners

Acoustic liner samples with rectangular, sharp chamber bends are 3D-printed and tested in the NASA Normal Incidence Tube. The angle of the bend is varied from 0 to 180 degrees. Four different methods are used to predict the impedance behavior of these samples: (1) midline-length assumption, (2) Cummings correction for curved duct bends, (3) Cummings prediction for 180 degrees sharp duct bends, and (4) Helmholtz equation via finite element method. Trends in the experimental data are discussed, predictions are compared to the experimental data, and recommendations are made on prediction method usage. Follow-on work is suggested, including capturing additional test points at higher frequencies and with different chamber geometries. Future modeling strategies are also briefly described.

chamber↗

NASA NTRS1972

Computation of space shuttle flow fields using noncentered finite-difference schemes.

Second- and third-order, noncentered finite-difference schemes are described for the numerical solution of the hyperbolic equations of fluid dynamics. The advantages of noncentered methods over the more conventional centered schemes are: simpler programming logic, nonhomogeneous terms are easily included, and generalization to multidimensional problems is direct. Second- and third-order methods are compared with regard to dissipative and dispersive errors and shock-capturing ability. These schemes are then used in a shock-capturing technique to determine the inviscid, supersonic flow field surrounding space shuttle vehicles (SSV). Resulting flow fields about typical pointed and blunted, delta-winged SSVs at angle of attack are presented and compared with experiment.

Kutler, P.↗

NASA NTRS1974

Efficient stream distributions in radiative transfer theory

This paper discusses a new, computationally-efficient method for approximating the integro-differential equation of radiative transfer with a finite set of coupled differential equations for discrete streams. The method uses recommended spatial distributions of streams that are quite different from those typically used in that they are based on the symmetry of several regular Platonic solids. To facilitate the use of such distributions, an explicit, one-parameter relationship between the physical radiance and the abstract stream is formulated. The parameter is used to determine the minimum number of streams required in the radiative transfer model. Accuracy and computational efficiency are shown to be served best by choosing a stream distribution that is invariant to a large number of three space rotations. For various values of the above-mentioned parameter, the resulting recommended stream distribution is shown to be more computationally efficient than more conventional stream distributions. Finally, the incorporation of polarization in the stream definition is described.

Whitney, C.↗

NASA NTRS1979

A review and assessment of the stress-intensity factors for surface cracks

The stress-intensity factor solutions proposed for a surface crack in a finite plate subjected to uniform tension are reviewed. Fourteen different solutions obtained over the past 16 years using approximate analytical methods, experimental methods, and engineering estimates are compared. The accuracy of the various solutions is assessed by correlating fracture data on surface-cracked tension specimens made of a brittle epoxy material. Fracture of the epoxy material was characterized by a constant value of stress-intensity factor as failure. Thus, the correctness of various solutions is judged by the variations in the stress-intensity factors at failure. The solutions were ranked in order of minimum standard deviation. The highest ranking solutions correlated 95% of data analyzed within + or - 10% whereas the lowest ranking solutions correlated 95% of data analyzed within + or - 20%. Some solutions could be applied to all data considered, whereas others were limited with respect to crack shapes and crack sizes that could be analyzed.

Newman, J. C., Jr.↗

NASA NTRS1982

Transonic perturbation analysis of wing-fuselage-nacelle-pylon configurations with powered jet exhausts

A transonic small perturbation method has been developed for the analysis of general wing-fuselage-nacelle-pylon configurations with powered jet exhausts. Finite difference successive line relaxation algorithm is used to solve the small disturbance potential equation in conservative form. The nacelle tangency condition and the jet exhaust plume contact conditions are fulfilled in a quasi-cylindrical fashion on a surface fitting the Cartesian grid. The pylon tangency condition is treated in a quasi-planar manner as for the wing. Viscous displacement effects on the wing are modeled by suitable shape changes including the placement of a viscous ramp at the base of the shock. Computed results of a transport configuration show satisfactory correlation with test data.

Wai, J. C.↗

NASA NTRS1986

Polynomial elimination theory and non-linear stability analysis for the Euler equations

Numerical methods are presented that exploit the polynomial properties of discretizations of the Euler equations. It is noted that most finite difference or finite volume discretizations of the steady-state Euler equations produce a polynomial system of equations to be solved. These equations are solved using classical polynomial elimination theory, with some innovative modifications. This paper also presents some preliminary results of a new non-linear stability analysis technique. This technique is applicable to determining the stability of polynomial iterative schemes. Results are presented for applying the elimination technique to a one-dimensional test case. For this test case, the exact solution is computed in three iterations. The non-linear stability analysis is applied to determine the optimal time step for solving Burgers' equation using the MacCormack scheme. The estimated optimal time step is very close to the time step that arises from a linear stability analysis.

Kennon, S. R.↗

NASA NTRS1988

Minimum weight design of rotorcraft blades with multiple frequency and stress constraints

Minimum weight designs of helicopter rotor blades with constraints on multiple coupled flap-lag natural frequencies are studied. Constraints are imposed on the minimum value of the blade autorotational inertia to ensure sufficient rotary inertia to autorotate in case of engine failure and on stresses to guard against structural failure due to blade centrifugal forces. Design variables include blade taper ratio, dimensions of the box beam located inside the airfoil and magnitudes of nonstructural weights. The program CAMRAD is used for the blade modal analysis; the program CONMIN is used for the optimization. A linear approximation involving Taylor series expansion is used to reduce the analysis effort. The procedure contains a sensitivity analysis consisting of analytical derivatives for objective function and constraints on autorotational inertia and stresses. Central finite difference derivatives are used for frequency constraints. Optimal designs are obtained for both rectangular and tapered blades. Using this method, it is possible to design a rotor blade with reduced weight, when compared to a baseline blade, while satisfying all the imposed design requirements.

Chattopadhyay, Aditi↗

NASA NTRS1992

Incompressible viscous flow computations for the pump components and the artificial heart

A finite difference, three dimensional incompressible Navier-Stokes formulation to calculate the flow through turbopump components is utilized. The solution method is based on the pseudo compressibility approach and uses an implicit upwind differencing scheme together with the Gauss-Seidel line relaxation method. Both steady and unsteady flow calculations can be performed using the current algorithm. Here, equations are solved in steadily rotating reference frames by using the steady state formulation in order to simulate the flow through a turbopump inducer. Eddy viscosity is computed by using an algebraic mixing-length turbulence model. Numerical results are compared with experimental measurements and a good agreement is found between the two.

Kiris, Cetin↗

NASA NTRS1990

Comparisons of thermochemical nonequilibrium viscous flowfield predictions for AFE vehicle

Numerical solutions for the Aeroassist Flight Experiment vehicle were obtained from three methods at a trajectory point corresponding to the maximum aerodynamic heating. The flow regime and vehicle's speed require a viscous model of laminar flow and finite-rate chemical and thermal modeling. The computational domain covers both forebody and base such that the shock layer and near wake flowfield are included. Because of differences in computational grids, methods of solution, and models of rate equations, the results are generally in poor agreement. Temperature and species concentrations are strongly affected by the physical model equations and associated parameters. The chemistry model based on 11 species is found to yield lower translational temperature profile near stagnation than those from a seven-species model. The vibrational temperature varies according to the modeling details. All solutions indicate that strong neutral and/or molecular dissociation and weak ionization take place at the forebody and vibrational freezing is present in the afterbody expansion region where the vibrational temperature is higher than the translational temperature. Some forms of shear layer emanating from the aerobrake skirt coalesce in the region of reversed flow behind the vehicle.

Tam, Luen T.↗

NASA NTRS1988

Transonic stability and control of aircraft using CFD methods

Implementation of a capability to calculate longitudinal short-period response in the CAP-TSD (Computational Aeroelasticity Program - Transonic Small Disturbance) finite-difference code is described. The code, developed recently at the NASA Langley Research Center, is capable of solving steady and unsteady flows about complete aircraft configurations and is used primarily for aeroelastic calculations in the critical transonic speed range. The longitudinal short-period equations of motion in state-space form have been coupled to the time-accurate lift and moment calculated by the program. Transient responses to an elevator pulse for free-flying aircraft demonstrate the new capability. A trim routine is also added to the code to obtain trim automatically during steady-state flow field convergence. Stability and control derivatives are estimated from the calculated transient response by a maximum likelihood estimation program. Results for a fighter configuration and a general aviation configuration are presented to assess the capability.

Vinh, Lam-Son↗