SEARCH · Search NASA

Results for “Finite difference scheme”

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,063 records · Page 59

All sources NASA NTRS DOE OSTI

NASA NTRS1984

Analysis of airfoil transitional separation bubbles

A previously developed local inviscid-viscous interaction technique for the analysis of airfoil transitional separation bubbles, ALESEP (Airfoil Leading Edge Separation), has been modified to utilize a more accurate windward finite difference procedure in the reversed flow region, and a natural transition turbulence model has been incorporated for the prediction of transition within the separation bubble. Numerous calculations and experimental comparisons are presented to demonstrate the effects of the windward differencing scheme and the natural transition turbulence model. Grid sensitivity and convergence capabilities of this inviscid-viscous interaction technique are briefly addressed. A major conclusion of this paper is that a second, counter-rotating eddy has been found to exist in the wall layer of the primary separation bubble with the use of windward differencing.

Davis, R. L.↗

NASA NTRS1985

Full two-dimensional transient solutions of electrothermal aircraft blade deicing

Two finite difference methods are presented for the analysis of transient, two-dimensional responses of an electrothermal de-icer pad of an aircraft wing or blade with attached variable ice layer thickness. Both models employ a Crank-Nicholson iterative scheme, and use an enthalpy formulation to handle the phase change in the ice layer. The first technique makes use of a 'staircase' approach, fitting the irregular ice boundary with square computational cells. The second technique uses a body fitted coordinate transform, and maps the exact shape of the irregular boundary into a rectangular body, with uniformally square computational cells. The numerical solution takes place in the transformed plane. Initial results accounting for variable ice layer thickness are presented. Details of planned de-icing tests at NASA-Lewis, which will provide empirical verification for the above two methods, are also presented.

Masiulaniec, K. C.↗

NASA NTRS1989

Numerical analysis of flow about a total temperature sensor

The unsteady flowfield about an airfoil-shaped inlet temperature sensor has been investigated using the thin-layer and full Navier-Stokes equations. A finite-volume formulation of the governing equations was used in conjunction with a Runge-Kutta time stepping scheme to analyze the flow about the sensor. Flow characteristics for this configuration were established at Mach numbers of 0.5 and 0.8 for different Reynolds numbers. The results were obtained for configurations of increasing complexity; important physical phenomena such as shock formation, boundary-layer separation, and unsteady wake formation were noted. Based on the computational results, recommendations for further study and refinement of the inlet temperature sensor were made.

Von Lavante, Ernst↗

NASA NTRS1993

Comparison of Numerical Model Results with Diffusion Flames in Microgravity

The effects of gravity on methane diffusion flames were studied using a finite difference numerical model. Variables under consideration included velocity, pressure, temperature, enthalpy, reactive species, and inert carrier gas. Results obtained showed that flames in zero gravity were wider and taller than those in normal gravity. The use of a six reaction scheme produced similar predictions of flame width and height compared to those of the single reaction case. The numerical models predicted a considerable difference in the flame base region between normal gravity and low gravity. This region is considered to be beyond the boundaries of current analytical models.

Hyer, Paul V.↗

NASA NTRS1996

Simulation of Spiral Slot Antennas on Composite Platforms

The project goals, plan and accomplishments up to this point are summarized in the viewgraphs. Among the various accomplishments, the most important have been: the development of the prismatic finite element code for doubly curved platforms and its validation with many different antenna configurations; the design and fabrication of a new slot spiral antennas suitable for automobile cellular, GPS and PCs communications; the investigation and development of various mesh truncation schemes, including the perfectly matched absorber and various fast integral equation methods; and the introduction of a frequency domain extrapolation technique (AWE) for predicting broadband responses using only a few samples of the response. This report contains several individual reports most of which have been submitted for publication to referred journals. For a report on the frequency extrapolation technique, the reader is referred to the UM Radiation Laboratory report A total of 14 papers have been published or accepted for publication with the full or partial support of this grant. Several more papers are in preparation.

Volakis, John L.↗

NASA NTRS2015

Structured Overlapping Grid Simulations of Contra-rotating Open Rotor Noise

Computational simulations using structured overlapping grids with the Launch Ascent and Vehicle Aerodynamics (LAVA) solver framework are presented for predicting tonal noise generated by a contra-rotating open rotor (CROR) propulsion system. A coupled Computational Fluid Dynamics (CFD) and Computational AeroAcoustics (CAA) numerical approach is applied. Three-dimensional time-accurate hybrid Reynolds Averaged Navier-Stokes/Large Eddy Simulation (RANS/LES) CFD simulations are performed in the inertial frame, including dynamic moving grids, using a higher-order accurate finite difference discretization on structured overlapping grids. A higher-order accurate free-stream preserving metric discretization with discrete enforcement of the Geometric Conservation Law (GCL) on moving curvilinear grids is used to create an accurate, efficient, and stable numerical scheme. The aeroacoustic analysis is based on a permeable surface Ffowcs Williams-Hawkings (FW-H) approach, evaluated in the frequency domain. A time-step sensitivity study was performed using only the forward row of blades to determine an adequate time-step. The numerical approach is validated against existing wind tunnel measurements.

Grid↗

DOE OSTI2022

A hybrid nodal-staggered pseudo-spectral electromagnetic particle-in-cell method with finite-order centering

Electromagnetic particle-in-cell (PIC) codes are widely used to perform computer simulations of a variety of physical systems, including fusion plasmas, astrophysical plasmas, plasma wakefield particle accelerators, and secondary photon sources driven by ultra-intense lasers. In a PIC code, Maxwell's equations are solved on a grid with a numerical method of choice. This article focuses on pseudo-spectral analytical time-domain (PSATD) algorithms and presents a novel hybrid PSATD PIC scheme that combines the respective advantages of standard nodal and staggered methods. The novelty of the hybrid scheme consists in using finite-order centering of grid quantities between nodal and staggered grids, in order to combine the solution of Maxwell's equations on a staggered grid with the deposition of charges and currents and the gathering of electromagnetic forces on a nodal grid. The correctness and performance of the novel hybrid scheme are assessed by means of numerical tests that employ different classes of PSATD equations in a variety of physical scenarios, ranging from the modeling of electron-positron pair creation in vacuum to the simulation of laser-driven and particle beam-driven plasma wakefield acceleration. It is shown that the novel hybrid scheme offers significant numerical and computational advantages, compared to purely nodal or staggered methods, for all the test cases presented.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Inverse design of plasma metamaterial devices with realistic elements

In an expansion of a previous study, we apply inverse design methods to produce two-dimensional plasma metamaterial devices with realistic plasma elements which incorporate quartz envelopes, collisionality (loss), non-uniform density profiles, and resistance to experimental error/perturbation. Finite difference frequency domain simulations are used along with forward-mode differentiation to design waveguides and demultiplexers operating under the transverse magnetic polarization. Additionally, optimal devices with realistic elements are compared to previous devices with idealized elements, and several parameter initialization schemes for the optimization algorithm are explored, yielding a robust procedure for producing such devices. Demultiplexing and waveguiding are demonstrated for microwave-regime devices composed of plasma elements with reasonable space-averaged plasma frequencies ~10GHz and a collision frequency ~1 GHz, allowing for future in-situ training and experimental realization of these designs.

42 ENGINEERING↗

DOE OSTI2017

Two-color field enhancement at an STM junction for spatiotemporally resolved photoemission

Here, we report measurements and numerical simulations of ultrafast laser-excited carrier flow across a scanning tunneling microscope (STM) junction. The current from a nanoscopic tungsten tip across a ~1 nm vacuum gap to a silver surface is driven by a two-color excitation scheme that uses an optical delay-modulation technique to extract the two-color signal from background contributions. The role of optical field enhancements in driving the current is investigated using density functional theory and full three-dimensional finite-difference time-domain computations. We find that simulated field-enhanced two-photon photoemission (2PPE) currents are in excellent agreement with the observed exponential decay of the two-color photoexcited current with increasing tip–surface separation, as well as its optical-delay dependence. The results suggest an approach to 2PPE with simultaneous subpicosecond temporal and nanometer spatial resolution.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1985

Thermomechanically induced pre- and postbuckling of general structure

An algorithmic solution strategy which enables handling the positive/indefinite stiffness characteristics associated with the pre and postbuckling of structures subject to complex thermomechanical loading fields was developed. The flexibility of the procedure is such that it can be applied to both finite difference and element type simulations. Due to the generality of the algorithmic approach developed, both kinematic and thermal/mechanical type material nonlinearity including inelastic effects can be treated. This includes the possibility of handling completely general thermomechanical boundary conditions. To demonstrate the scheme, the results of several benchmark problems are presented.

Padovan, J.↗

NASA NTRS1990

High accuracy solutions of incompressible Navier-Stokes equations

In recent years, high accuracy finite difference approximations were developed for partial differential equations of elliptic type, with particular emphasis on the convection-diffusion equation. These approximations are of compact type, have a local truncation error of fourth order, and allow the use of standard iterative schemes to solve the resulting systems of algebraic equations. These high accuracy approximations are extended to the solution of Navier-Stokes equations. Solutions are obtained for the model problem of driven cavity and are compared with solutions obtained using other approximations and those obtained by other authors. It is discovered that the high order approximations do indeed produce high accuracy solutions and have a potential for use in solving important problems of viscous fluid flows.

Gupta, Murli M.↗

NASA NTRS1993

Conservative high-order-accurate finite-difference methods for curvilinear grids

Two fourth-order-accurate finite-difference methods for numerically solving hyperbolic systems of conservation equations on smooth curvilinear grids are presented. The first method uses the differential form of the conservation equations; the second method uses the integral form of the conservation equations. Modifications to these schemes, which are required near boundaries to maintain overall high-order accuracy, are discussed. An analysis that demonstrates the stability of the modified schemes is also provided. Modifications to one of the schemes to make it total variation diminishing (TVD) are also discussed. Results that demonstrate the high-order accuracy of both schemes are included in the paper. In particular, a Ringleb-flow computation demonstrates the high-order accuracy and the stability of the boundary and near-boundary procedures. A second computation of supersonic flow over a cylinder demonstrates the shock-capturing capability of the TVD methodology. An important contribution of this paper is the dear demonstration that higher order accuracy leads to increased computational efficiency.

Rai, Man M.↗

NASA NTRS1993

An upwind, kinetic flux-vector splitting method for flows in chemical and thermal non-equilibrium

We have developed new upwind kinetic difference schemes for flows with non-equilibrium thermodynamics and chemistry. These schemes are derived from the Boltzmann equation with the resulting Euler schemes developed as moments of the discretized Boltzmann scheme with a locally Maxwellian velocity distribution. Splitting the velocity distribution at the Boltzmann level is seen to result in a flux-split Euler scheme and is called Kinetic Flux Vector Splitting (KFVS). Extensions to flows with finite-rate chemistry and vibrational relaxation is accomplished utilizing nonequilibrium kinetic theory. Computational examples are presented comparing KFVS with the schemes of Van Leer and Roe for a quasi-one-dimensional flow through a supersonic diffuser, inviscid flow through two-dimensional inlet, and viscous flow over a cone at zero angle-of-attack. Calculations are also shown for the transonic flow over a bump in a channel and the transonic flow over an NACA 0012 airfoil. The results show that even though the KFVS scheme is a Riemann solver at the kinetic level, its behavior at the Euler level is more similar to the existing flux-vector splitting algorithms than to the flux-difference splitting scheme of Roe.

Eppard, W. M.↗

DOE OSTI2024

The Transient Multi-Level method for Monte Carlo reactor statics calculations

The Transient Multi-Level (TML) method is applied to a time-dependent Monte Carlo transport solver to offload some of the computational burden of the expensive Monte Carlo solve to lower-order Coarse Mesh Finite Difference (CMFD) and Exact Point Kinetics Equations (EPKE) solvers via factorization of the neutron flux at the transport and CMFD levels using the Predictor Corrector Quasi-Static Method (PCQM). The Monte Carlo transient is solved by a modified fission source iteration scheme that introduces a single transient source bank. The method is implemented in the production-level Monte Carlo code, Shift, and verified with prescribed reactivity ramps from the two-dimensional version of the C5G7-TD reactor benchmark. The results show that, as compared to other quasi-static methods, the TML reduces the stochastic noise inherent to the transient Monte Carlo solver by factors of ~2 to 6 for various norm comparisons of the reactor power amplitude. Finally, the TML additionally reduces the number of Monte Carlo evaluations needed to simulate the transient, leading to roughly an order of magnitude improvement in CPU time relative to the standard PCQM for the problems tested.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2018

Efficient three-dimensional Laplace-Fourier domain acoustic-wave simulations using discontinuous finite-difference meshes with embedded boundaries

© 2018 SEG. We develop embedded boundary and discontinuous mesh methods to handle arbirarily shaped topography and accurately simulate acoustic seismic wave propagation in Laplace-Fourier domain. The purpose of the embedded boundary method is to enhance accurate wave simulation near the surface and the discontinuous mesh method is used to achieve considerable savings in both computation time and memory savings relative to fixed mesh schemes.

AlSalem, Hussain↗

NASA NTRS1978

The development of a three-dimensional partially elliptic flow computer program for combustor research

A three dimensional, partially elliptic, computer program was developed. Without requiring three dimensional computer storage locations for all flow variables, the partially elliptic program is capable of predicting three dimensional combustor flow fields with large downstream effects. The program requires only slight increase of computer storage over the parabolic flow program from which it was developed. A finite difference formulation for a three dimensional, fully elliptic, turbulent, reacting, flow field was derived. Because of the negligible diffusion effects in the main flow direction in a supersonic combustor, the set of finite-difference equations can be reduced to a partially elliptic form. Only the pressure field was governed by an elliptic equation and requires three dimensional storage; all other dependent variables are governed by parabolic equations. A numerical procedure which combines a marching integration scheme with an iterative scheme for solving the elliptic pressure was adopted.

Pan, Y. S.↗

NASA NTRS1980

Numerical solution techniques for unsteady transonic aerodynamics problems

Basic concepts of finite difference solution techniques for unsteady transonic flows are presented. The hierarchy of mathematical forumulations that approximate the Navier-Stokes equations are reviewed. The basic concepts involved in constructing numerical algorthms to solve these formulations are given. Semi-implicit and implicit schemes are constructed and analyzed. The discussion focuses primarily on techniques for solving the low frequency transonic small disturbance equation. This is the simplest formulation that contains the essence of inviscid unsteady transonic flow physics. The low frequency formulation is emphasized here because codes based on this theory can be run in minutes of processor time on currently available computers. Furthermore, numerical techniques involved in solving this simple formulation also apply to the more complicated formulations. Extensions to these formulations are briefly described. An indication of the present capability for solving unsteady transonic flows is provided. Important areas of future research for the advancement of computational unsteady transonic aerodynamics are described.

Ballhaus, W. F.↗

NASA NTRS1984

A full-implicit-continuous-Eulerian (FICE) scheme for multidimensional transient magnetohydrodynamic (MHD) flows

A full implicit continuous Eulerian (FICE) scheme is developed for solving multidimensional transient MHD flow problems. The physical system under consideration is a general case of a transient MHD flow in which an initial steady state is subject to a finite amplitude disturbance. The governing equations are described, their finite difference formulation is presented and the FICE algorithm is given. The boundary conditions are treated by classifying them into physical and computational ones. The usefulness of the FICE algorithm is demonstrated using a physical example concerning the dynamical response of the static solar atmosphere due to a representative photospheric disturbance.

Wu, S. T.↗