SEARCH · Search NASA

Results for “finite difference method”

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,225 records · Page 68

All sources NASA NTRS DOE OSTI

NASA NTRS1987

Steady flow past sudden expansions at large Reynolds number. II - Navier-Stokes solutions for the cascade expansion

The equations of motion in the steady laminar flow past a sudden expansion at large Reynolds number, R, reduce to the boundary layer equations as R tends to the freestream valve, when the longitudinal length scale of the separated eddy increases linearly and indefinitely with R. A global Newton method is presently used to obtain finite difference solutions to the steady Navier-Stokes equations up to R of 1000 for a uniform inflow past a cascade of sudden expansions. For large expansion ratio values, the eddy length increases linearly with R; for smaller values, however, where boundary layer equation solutions could not be found, the steady solutions to the Navier-Stokes equations approach the limit of an inviscid eddy in length with increasing R.

Milos, Frank S.↗

NASA NTRS1987

A new PNS code for three-dimensional chemically reacting flows

A new parabolized Navier-Stokes (PNS) code has been developed to compute the hypersonic, viscous, chemically reacting flow fields around three-dimensional bodies. The flow medium is assumed to be a multicomponent mixture of thermally perfect but calorically imperfect gases. The new PNS code solves the gasdynamic and species conservation equations in a coupled manner using a noniterative, implicit, approximately-factored, finite-difference algorithm. The space-marching method is made well-posed by special treatment of the streamwise pressure gradient term. The code has been used to compute hypersonic laminar flow of chemically reacting air over cones at angles of attack. The results of the computations are compared with the results of reacting boundary-layer computations and show excellent agreement.

Prabhu, D. K.↗

NASA NTRS1988

A semianalytical technique for sensitivity analysis of unsteady aerodynamic computations

A semianalytical approach is developed for the sensitivity analysis of linear unsteady aerodynamic loads. The semianalytical approach is easier to implement than the analytical approach. It is also computationally less expensive than the finite difference approach when used with panel methods, which require a large number of panels. The semianalytical approach is applied to an isolated airfoil in a 2-D flow and rotating propfan blades in 3-D flow. Sensitivity coefficients with respect to non-shape-dependent variables are shown for some cases. It is expected that the semianalytical approach will be useful in aeroelastic design procedures particularly when mistuning is present, and that it is potentially useful for shape sensitivity analysis of linear unsteady aerodynamics.

Murthy, Durbha V.↗

NASA NTRS1988

Numerical computation of space shuttle orbiter flow field

A new parabolized Navier-Stokes (PNS) code has been developed to compute the hypersonic, viscous chemically reacting flow fields around 3-D bodies. The flow medium is assumed to be a multicomponent mixture of thermally perfect but calorically imperfect gases. The new PNS code solves the gas dynamic and species conservation equations in a coupled manner using a noniterative, implicit, approximately factored, finite difference algorithm. The space-marching method is made well-posed by special treatment of the streamwise pressure gradient term. The code has been used to compute hypersonic laminar flow of chemically reacting air over cones at angle of attack. The results of the computations are compared with the results of reacting boundary-layer computations and show excellent agreement.

Tannehill, John C.↗

NASA NTRS1988

Application of a semianalytical technique for sensitivity analysis of unsteady aerodynamic computations

A semianalytical approach is developed for the sensitivity analysis of linear unsteady aerodynamic loads. The semianalytical approach is easier to implement than the analytical approach. It is also computationally less expensive than the finite difference approach when used with panel methods which require a large number of panels. The semianalytical approach is applied to an isolated airfoil in a two-dimensional flow and rotating propfan blades in three-dimensional flow. Sensitivity coefficients with respect to nonshape-dependent variables are shown for some cases.

Murthy, Durbha V.↗

NASA NTRS1988

Navier-Stokes simulation of unsteady three-dimensional blade-vortex interactions

The unsteady, viscous, three-dimensional flow field of a helicopter rotor blade encountering a passing vortex is calculated by solving the thin layer Navier-Stokes equations by a finite-difference numerical procedure. A prescribed vortex method is adopted to preserve the structure of the interacting concentrated vortex. The test cases considered correspond to the experimental model rotor test conditions of Caradonna et al.

Srinivasan, G. R.↗

NASA NTRS1989

Flow of rarefied gases over two-dimensional bodies

A kinetic-theory analysis is made of the flow of rarefied gases over two-dimensional bodies of arbitrary curvature. The Boltzmann equation simplified by a model collision integral is written in an arbitrary orthogonal curvilinear coordinate system, and solved by means of finite-difference approximation with the discrete ordinate method. A numerical code is developed which can be applied to any two-dimensional submerged body of arbitrary curvature for the flow regimes from free-molecular to slip at transonic Mach numbers. Predictions are made for the case of a right circular cylinder.

Jeng, Duen-Ren↗

NASA NTRS1989

Unsteady interaction of a rotor with a vortex

The unsteady, three-dimensional flow field of a helicopter rotor blade encountering a passing vortex is calculated by solving the Euler/thin layer Navier-Stokes equations by a finite-difference numerical procedure. A prescribed vortex method is adopted to preserve the structure of the interacting vortex. The cases considered for computation correspond to the experimental model rotor test conditions of Caradonna, et al. and consist of parallel and oblique interactions. Comparison of the numerical results with test data show good agreement for both parallel and oblique interactions at subsonic and transonic tip speeds.

Srinivasan, G. R.↗

NASA NTRS1989

Rarefied gas flow through two-dimensional nozzles

A kinetic theory analysis is made of the flow of a rarefied gas from one reservoir to another through two-dimensional nozzles with arbitrary curvature. The Boltzmann equation simplified by a model collision integral is solved by means of finite-difference approximations with the discrete ordinate method. The physical space is transformed by a general grid generation technique and the velocity space is transformed to a polar coordinate system. A numerical code is developed which can be applied to any two-dimensional passage of complicated geometry for the flow regimes from free-molecular to slip. Numerical values of flow quantities can be calculated for the entire physical space including both inside the nozzle and in the outside plume. Predictions are made for the case of parallel slots and compared with existing literature data. Also, results for the cases of convergent or divergent slots and two-dimensional nozzles with arbitrary curvature at arbitrary knudsen number are presented.

De Witt, Kenneth J.↗

NASA NTRS1991

Evaluating lightning hazards to building environments using explicit numerical solutions of Maxwell's equations

The objective here is to describe the lightning hazards to buildings and their internal environments using advanced formulations of Maxwell's Equations. The method described is the Three Dimensional Finite Difference Time Domain Solution. It can be used to solve for the lightning interaction with such structures in three dimensions with the inclusion of a considerable amount of detail. Special techniques were developed for including wire, plumbing, and rebar into the model. Some buildings have provisions for lightning protection in the form of air terminals connected to a ground counterpoise system. It is shown that fields and currents within these structures can be significantly high during a lightning strike. Time lapse video presentations were made showing the electric and magnetic field distributions on selected cross sections of the buildings during a simulated lightning strike.

Collier, Richard S.↗

NASA NTRS1991

A parametric sensitivity and optimization study for the active flexible wing wind-tunnel model flutter characteristics

In this paper an effort is made to improve the analytical open-loop flutter predictions for the Active Flexible Wing wind-tunnel model using a sensitivity based optimization approach. The sensitivity derivatives of the flutter frequency and dynamic pressure of the model with respect to the lag terms appearing in the Roger's unsteady aerodynamics approximations are evaluated both analytical and by finite differences. Then, the Levenberg-Marquardt method is used to find the optimum values for these lag-terms. The results obtained here agree much better with the experimental (wind tunnel) results than those found in the previous studies.

Rais-Rohani, Masoud↗

NASA NTRS1993

Numerical modeling for an electric-field hyperthermia applicator

Hyperthermia, in conjunction with radiation and chemotherapy for treatment of cancers, is an area of current concern. Experiments have shown that hyperthermia can increase the potency of many chemotherapy drugs and the effectiveness of radiation for treating cancer. A combination of whole body or regional hyperthermia with chemotherapy or radiation should improve treatment results. Conventional methods for inducing whole body hyperthermia, such as exposing a patient in a radiant cabinet or under a hot water blanket, conduct heat very slowly from the skin to the body core. Thus a more efficient system, such as the three-plate electric-field hyperthermia applicator (EHA), is developed. This three-plate EHA has one top plate over and two lower plates beneath the patient. It is driven at 27.12 MHz with 500 Watts through a matching circuit. Using this applicator, a 50 kg pig was successfully heated to 42 C within 45 minutes. However, phantom and animal studies have indicated non-uniform heating near the side of the body. In addition, changes in the size and distance between the electrode plates can affect the heating (or electromagnetic field) pattern. Therefore, numerical models using the method of moments (MOM) or the finite difference time domain (FDTD) technique are developed to optimize the heating pattern of this EHA before it is used for human trials. The accuracy of the numerical modeling has been achieved by the good agreement between the MOM and FDTD results for the three-plate EHA without a biological body. The versatile FDTD technique is then applied to optimize the EHA design with a human body. Both the numerical and measured data in phantom blocks will be presented. The results of this study will be used to design an optimized system for whole body or regional hyperthermia.

Wu, Te-Kao↗

NASA NTRS2017

On-Orbit Structural Dynamics Performance of the Roll-Out Solar Array

The Roll-Out Solar Array (ROSA) flight experiment was launched to the International Space Station (ISS) on June 3rd, 2017. ROSA is an innovative, lightweight solar array with a flexible substrate that makes use of the stored strain energy in its composite structural members to provide deployment without the use of motors. This paper will discuss the results of various structural dynamics experiments conducted on the ISS during the weeks following launch. Data gathered from instrumentation on the solar array wing during the experiments are compared with pre-flight predictions from two different finite element modeling efforts. Two distinct methods were used to reconstruct the modal characteristics of ROSA from the data collected on orbit. Of particular interest in this effort are the first few system modes and mode shapes of the array, the amount of structural damping present, and degree of structural thermal interaction seen during eclipse exit. Discrepancies between the behavior predicted by the models and that observed on orbit are identified and discussed. The goal in this effort was to better understand the performance of ROSA and to improve modeling efforts for future designs of similar solar arrays.

Chamberlain, Matthew K.↗

NASA NTRS1994

3D Electromagnetic Plasma Particle Simulations on the Intel Delta Parallel Computer

A three-dimensional electromagnetic PIC code has been developed on the 512 node Intel Touchstone Delta MIMD parallel computer. This code is based on the General Concurrent PIC algorithm which uses a domain decomposition to divide the computation among the processors. The 3D simulation domain can be partitioned into 1-, 2-, or 3-dimensional subdomains. Particles must be exchanged between processors as they move among the subdomains.

PIC↗

NASA NTRS1977

Calculation methods for compressible turbulent boundary layers, 1976

Equations and closure methods for compressible turbulent boundary layers are discussed. Flow phenomena peculiar to calculation of these boundary layers were considered, along with calculations of three dimensional compressible turbulent boundary layers. Procedures for ascertaining nonsimilar two and three dimensional compressible turbulent boundary layers were appended, including finite difference, finite element, and mass-weighted residual methods.

Bushnell, D. M.↗

DOE OSTI2018

Exploiting Geometric Partitioning in Task Mapping for Parallel Computes

We present a new method for mapping applications' MPI tasks to cores of a parallel computer such that applications' communication time is reduced. We address the case of sparse node allocation, where the nodes assigned to a job are not necessarily located in a contiguous block nor within close proximity to each other in the network, although our methods generalize to contiguous allocations as well. The goal is to assign tasks to cores so that interdependent tasks are performed by "nearby' cores, thus lowering the distance messages must travel, the amount of congestion in the network, and the overall cost of communication. Our new method applies a geometric partitioning algorithm to both the tasks and the processors, and assigns task parts to the corresponding processor parts. We also present a number of algorithmic optimizations that exploit specific features of the network or application. We show that, for the structured finite difference mini-application MiniGhost, our mapping methods reduced communication time up to 75% relative to MiniGhost's default mapping on 128K cores of a Cray XK7 with sparse allocation. For the atmospheric modeling code E3SM/HOMME, our methods reduced communication time up to 31% on 32K cores of an IBM BlueGene/Q with contiguous allocation.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2019

GPR Imaging for Deeply Buried Objects: A Comparative Study Based on Compositing of Scanning Frequencies and a Chirp Excitation Function

Compositing of ground penetrating radar (GPR) scans of differing frequencies have been found to produce cleaner images at depth using the Gaussian mixture model (GMM) feature of the expectation-maximization (EM) algorithm. GPR scans at various heights (“Stand Off”), as well as ground-based scans, have been studied. In this paper, we compare the GPR response from a chirp excitation function-based radar with the response from the EM GMM algorithm compositing process, using the same mix of frequencies. A chirp excitation pulse was found to be effective in delineating the defined buried object, but the resulting image is less sharp than the GMM EM method.

58 GEOSCIENCES↗

DOE OSTI2018

Numerical similarity solution for a variable coefficient K(m, n) equation

A technique for finding numerical similarity solutions to an initial boundary value problem (IBVP) for generalized K(m, n) equations is described. The equation under consideration is nonlinear and has variable coefficients. The original problem is transformed with the aid of Lie symmetries to an initial value problem (IVP) for a nonlinear third-order ordinary differential equation. The existence and uniqueness of the solution are examined, and the problem is consequently solved with the aid of a finite-difference scheme for various values of the governing parameters. In lieu of an exact symbolic solution, the scheme is validated by comparing the numerical solutions with the approximate analytic solutions obtained with the aid of the method of successive approximations in their region of validity. The accuracy, efficiency, and consistency of the scheme are demonstrated. Numerical solutions to the original initial boundary value problem are constructed for selected parameter values with the aid of the transforms. The qualitative behavior of the solutions as a function of the governing parameters is analyzed, and it is found that the examined IBVPs for generalized K(m, n) equations with variable coefficients that are functions of time, do not admit solitary wave or compacton solutions.

97 MATHEMATICS AND COMPUTING↗