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,171 records · Page 65

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Nonlinear evolution of resistive tearing mode instability with shear flow and viscosity

The effect of shear flow on the nonlinear evolution of the tearing mode is investigated via numerical solutions of the resistive MHD equations in slab geometry, using a finite-difference alternative-direction implicit method. It was found that, when the shear flow is small (V less than 0.3), the tearing mode saturates within one resistive time, whereas for larger flows the nonlinear saturation develops on longer time scales. The magnetic energy release decreases and the saturation time increases with increasing values of V for both small and large resistivity. Shear flow was found to decrease the saturated magnetic island width and to generate currents far from the tearing layer. Results suggest that equilibrium shear flow may improve the confinement of tokamak plasma.

Ofman, L.↗

NASA NTRS1994

Low-dissipation and -dispersion Runge-Kutta schemes for computational acoustics

In this paper, we investigate accurate and efficient time advancing methods for computational acoustics, where non-dissipative and non-dispersive properties are of critical importance. Our analysis pertains to the application of Runge-Kutta methods to high-order finite difference discretization. In many CFD applications multi-stage Runge-Kutta schemes have often been favored for their low storage requirements and relatively large stability limits. For computing acoustic waves, however, the stability consideration alone is not sufficient, since the Runge-Kutta schemes entail both dissipation and dispersion errors. The time step is now limited by the tolerable dissipation and dispersion errors in the computation. In the present paper, it is shown that if the traditional Runge-Kutta schemes are used for time advancing in acoustic problems, time steps greatly smaller than that allowed by the stability limit are necessary. Low-Dissipation and -Dispersion Runge-Kutta (LDDRE) schemes are proposed, based on an optimization that minimizes the dissipation and dispersion errors for wave propagation. Order optimizations of both single-step and two-step alternating schemes are considered. The proposed LDDRK schemes are remarkably more efficient than the classical Runge-Kutta schemes for acoustic computations. Moreover, low storage implementations of the optimized schemes are discussed. Special issues of implementing numerical boundary conditions in the LDDRK schemes are also addressed.

Hu, F. Q.↗

NASA NTRS1995

Application of low dissipation and dispersion Runge-Kutta schemes to benchmark problems in computational aeroacoustics

We investigate accurate and efficient time advancing methods for computational aeroacoustics, where non-dissipative and non-dispersive properties are of critical importance. Our analysis pertains to the application of Runge-Kutta methods to high-order finite difference discretization. In many CFD applications, multi-stage Runge-Kutta schemes have often been favored for their low storage requirements and relatively large stability limits. For computing acoustic waves, however, the stability consideration alone is not sufficient, since the Runge-Kutta schemes entail both dissipation and dispersion errors. The time step is now limited by the tolerable dissipation and dispersion errors in the computation. In the present paper, it is shown that if the traditional Runge-Kutta schemes are used for time advancing in acoustic problems, time steps greatly smaller than that allowed by the stability limit are necessary. Low Dissipation and Dispersion Runge-Kutta (LDDRK) schemes are proposed, based on an optimization that minimizes the dissipation and dispersion errors for wave propagation. Optimizations of both single-step and two-step alternating schemes are considered. The proposed LDDRK schemes are remarkably more efficient than the classical Runge-Kutta schemes for acoustic computations. Numerical results of each Category of the Benchmark Problems are presented. Moreover, low storage implementations of the optimized schemes are discussed. Special issues of implementing numerical boundary conditions in the LDDRK schemes are also addressed.

Hu, F. Q.↗

NASA NTRS1996

Nonstandard Finite Difference Schemes: Relations Between Time and Space Step-Sizes in Numerical Schemes for PDE's That Follow from Positivity Condition

A large class of physical phenomena can be modeled by evolution and wave type Partial Differential Equations (PDE). Few of these equations have known explicit exact solutions. Finite-difference techniques are a popular method for constructing discrete representations of these equations for the purpose of numerical integration. However, the solutions to the difference equations often contain so called numerical instabilities; these are solutions to the difference equations that do not correspond to any solution of the PDE's. For explicit schemes, the elimination of this behavior requires functional relations to exist between the time and space steps-sizes. We show that such functional relations can be obtained for certain PDE's by use of a positivity condition. The PDE's studied are the Burgers, Fisher, and linearized Euler equations.

Mickens, Ronald E.↗

NASA NTRS1998

A Virtual World of Visualization

In 1990, Sterling Software, Inc., developed the Flow Analysis Software Toolkit (FAST) for NASA Ames on contract. FAST is a workstation based modular analysis and visualization tool. It is used to visualize and animate grids and grid oriented data, typically generated by finite difference, finite element and other analytical methods. FAST is now available through COSMIC, NASA's software storehouse.

Source record↗

NASA NTRS2018

Does Orbital Angular Momentum Have Effect on Laser’s Scattering by Molecular Atmosphere?

Lasers with orbital angular momentum (OAM) have potential applications in communication technology, manipulation of particles, and remote sensing. Because of its unusual light-scattering properties, the OAM laser’s interaction with a molecular atmosphere must be studied to ensure that it is not lossy for communication or remote-sensing applications that involve its transmission through an atmospheric environment. In this study, the finite-difference time-domain (FDTD) method [21] is applied to calculate the light scattering of the purely azimuthal (the radial mode number is assumed to be zero) Laguerre-Gaussian (LG) beams with OAM by very small dielectric particles. Not like Lorentz-Mie solutions, the FDTD method can calculate for particles off the central axis of the LG beam. It is found that when the particles are very small, and the topological charge number of the OAM of a laser is not extremely large, the laser’s OAM has little effect on the scattering phase function. This suggests that Rayleigh theory can be applied directly to calculate the light scattering by atmospheric molecules. The transmission of a laser beam with OAM in a molecular atmosphere is not different from that of a regular Gaussian beam.

Molecular atmosphere↗

NASA NTRS1976

Advances in the application of fast semidirect computational methods in transonic flow

The paper uses finite-difference algorithms called 'fast direct elliptic solvers' within an iteration scheme for the rapid solution of the equations of inviscid transonic aerodynamics. The methods are called 'direct' (or semidirect) because the entire computational field is solved at once rather than in successive traverses over the field. These semidirect iterative methods have been limited here to the investigation of two-dimensional steady inviscid flow over airfoils in subsonic free stream.

Martin, E. D.↗

NASA NTRS1984

Structural sensitivity analysis: Methods, applications, and needs

Some innovative techniques applicable to sensitivity analysis of discretized structural systems are reviewed. These techniques include a finite-difference step-size selection algorithm, a method for derivatives of iterative solutions, a Green's function technique for derivatives of transient response, a simultaneous calculation of temperatures and their derivatives, derivatives with respect to shape, and derivatives of optimum designs with respect to problem parameters. Computerized implementations of sensitivity analysis and applications of sensitivity derivatives are also discussed. Finally, some of the critical needs in the structural sensitivity area are indicated along with Langley plans for dealing with some of these needs.

Adelman, H. M.↗

NASA NTRS1996

Algebraic Multigrid by Smoothed Aggregation for Second and Fourth Order Elliptic Problems

Multigrid methods are very efficient iterative solvers for system of algebraic equations arising from finite element and finite difference discretization of elliptic boundary value problems. The main principle of multigrid methods is to complement the local exchange of information in point-wise iterative methods by a global one utilizing several related systems, called coarse levels, with a smaller number of variables. The coarse levels are often obtained as a hierarchy of discretizations with different characteristic meshsizes, but this requires that the discretization is controlled by the iterative method. To solve linear systems produced by existing finite element software, one needs to create an artificial hierarchy of coarse problems. The principal issue is then to obtain computational complexity and approximation properties similar to those for nested meshes, using only information in the matrix of the system and as little extra information as possible. Such algebraic multigrid method that uses the system matrix only was developed by Ruge. The prolongations were based on the matrix of the system by partial solution from given values at selected coarse points. The coarse grid points were selected so that each point would be interpolated to via so-called strong connections. Our approach is based on smoothed aggregation introduced recently by Vanek. First the set of nodes is decomposed into small mutually disjoint subsets. A tentative piecewise constant interpolation (in the discrete sense) is then defined on those subsets as piecewise constant for second order problems, and piecewise linear for fourth order problems. The prolongation operator is then obtained by smoothing the output of the tentative prolongation and coarse level operators are defined variationally.

Vanek, Petr↗

NASA NTRS2015

Impedance Eduction in Sound Fields With Peripherally Varying Liners and Flow

A two-dimensional impedance eduction theory is extended to three-dimensional sound fields and peripherally varying duct liners. The approach is to first measure the acoustic pressure field at a series of flush-mounted wall microphones located around the periphery of the flow duct. The numerical solution for the acoustic pressure field at these microphones is also obtained by solving the three-dimensional convected Helmholtz equation using the finite element method. A quadratic objective function based on the difference between the measured and finite element solution is constructed and the unknown impedance function is obtained by minimizing this objective function. Impedance spectra educed for two uniform-structure liners (a wire-mesh and a conventional liner) and a hard-soft-hard peripherally varying liner (for which the soft segment is that of the conventional liner) are presented. Results are presented at three mean flow Mach numbers and fourteen sound source frequencies. The impedance spectra of the uniform-structure liners are also computed using a two-dimensional impedance eduction theory. The primary conclusions of the study are: 1) when measured data is used with the uniform-structure liners, the three-dimensional theory reproduces the same impedance spectra as the two-dimensional theory except for frequencies corresponding to very low or very high liner attenuation; and 2) good agreement between the educed impedance spectra of the uniform structure conventional liner and the soft segment of the peripherally varying liner is obtained.

Watson, W. R.↗

NASA NTRS1971

Numerical calculation of turbulent boundary layers including suction or injection with binary diffusion.

The modelling of turbulent transport processes with an eddy viscosity representation of the inner, wall region for cases with suction or injection is considered in detail. It is shown how the van Driest solid wall model can be logically extended for such cases. Further, a new inductive model is developed. Finite difference calculations using the new inductive model are shown to be in good agreement with experimental data over a wide range of conditions. Lastly, a new method for calculating the wall shear within a finite difference solution is developed. This procedure greatly relaxes the requirement for fine grid spacing in the vertical direction.

Schetz, J. A.↗

NASA NTRS1991

Finite difference time domain calculation of transients in antennas with nonlinear loads

Determining transient electromagnetic fields in antennas with nonlinear loads is a challenging problem. Typical methods used involve calculating frequency domain parameters at a large number of different frequencies, then applying Fourier transform methods plus nonlinear equation solution techniques. If the antenna is simple enough so that the open circuit time domain voltage can be determined independently of the effects of the nonlinear load on the antennas current, time stepping methods can be applied in a straightforward way. Here, transient fields for antennas with more general geometries are calculated directly using Finite Difference Time Domain (FDTD) methods. In each FDTD cell which contains a nonlinear load, a nonlinear equation is solved at each time step. As a test case, the transient current in a long dipole antenna with a nonlinear load excited by a pulsed plane wave is computed using this approach. The results agree well with both calculated and measured results previously published. The approach given here extends the applicability of the FDTD method to problems involving scattering from targets, including nonlinear loads and materials, and to coupling between antennas containing nonlinear loads. It may also be extended to propagation through nonlinear materials.

Luebbers, Raymond J.↗

NASA NTRS1981

The use of adaptive grids in conjunction with shock capturing methods

The use of shock capturing finite-difference techniques in computing flow fields containing shocks results in a smeared or oscillatory solution in the vicinity of the shocks. This smearing or oscillatory behavior is due to the discretized form of the governing differential equations used to compute the solution. The discretization error can be reduced by a proper clustering of mesh points in the region of the shock and by using shock aligned grids. This paper uses a simple method that was developed earlier to cluster points near the shocks and serves to introduce a new method of generating a shock aligned mesh. Applications to the one-dimensional inviscid Burgers' equation, supersonic flow over a wedge with the associated straight oblique shock, one- and two-dimensional inviscid flows through an expanding duct and the problem of a straight oblique shock in a uniform supersonic freestream are presented. Significant reduction in the oscillatory behavior of the solution is demonstrated.

Rai, M. M.↗

NASA NTRS1989

Direct simulations of turbulent flow using finite-difference schemes

A high-order accurate finite-difference approach is presented for calculating incompressible turbulent flow. The methods used include a kinetic energy conserving central difference scheme and an upwind difference scheme. The methods are evaluated in test cases for the evolution of small-amplitude disturbances and fully developed turbulent channel flow. It is suggested that the finite-difference approach can be applied to complex geometries more easilty than highly accurate spectral methods. It is concluded that the upwind scheme is a good candidate for direct simulations of turbulent flows over complex geometries.

Rai, Man Mohan↗

NASA NTRS1980

Generations of orthogonal surface coordinates

Two generation methods were developed for three dimensional flows where the computational domain normal to the surface is small. With this restriction the coordinate system requires orthogonality only at the body surface. The first method uses the orthogonal condition in finite-difference form to determine the surface coordinates with the metric coefficients and curvature of the coordinate lines calculated numerically. The second method obtains analytical expressions for the metric coefficients and for the curvature of the coordinate lines.

Blottner, F. G.↗

NASA NTRS1984

A method for modeling finite-core vortices in wake-flow calculations

A numerical method for computing nonplanar vortex wakes represented by finite-core vortices is presented. The approach solves for the velocity on an Eulerian grid, using standard finite-difference techniques; the vortex wake is tracked by Lagrangian methods. In this method, the distribution of continuous vorticity in the wake is replaced by a group of discrete vortices. An axially symmetric distribution of vorticity about the center of each discrete vortex is used to represent the finite-core model. Two distributions of vorticity, or core models, are investigated: a finite distribution of vorticity represented by a third-order polynomial, and a continuous distribution of vorticity throughout the wake. The method provides for a vortex-core model that is insensitive to the mesh spacing. Results for a simplified case are presented. Computed results for the roll-up of a vortex wake generated by wings with different spanwise load distributions are presented; contour plots of the flow-field velocities are included; and comparisons are made of the computed flow-field velocities with experimentally measured velocities.

Stremel, P. M.↗

DOE OSTI2021

Evaluating dynamic thermal performance of building envelope components using small-scale calibrated hot box tests

The hot box test method has been applied to evaluate both the steady-state (U-value) and dynamic thermal properties of building envelopes. However, the high construction cost of full-scale hot box apparatus and the testing time required (usually several days) may prevent its wider adoption. To overcome the limitations of full-scale hot box tests, here we propose a novel method to evaluate the dynamic thermal performance of building envelope components using a small-scale calibrated hot box and scaled-down specimen. In this paper, the scaling relationships of thermal properties evaluated using a full-size specimen and a scaled-down specimen are established based on the Laplace transform of the heat transfer equations. In addition, dynamic thermal properties obtained from scaled-down experimental tests are compared to the values calculated by the EN ISO13768 (ISO) method. A small-scale hot box with a 355 mm × 355 mm metering area was constructed and calibrated to validate the correlations. Three scaled-down concrete sandwich wall panels were then tested and the scaling relationship was cross-validated using the experimental results, finite difference (FD) simulations, and the ISO method. The results indicate that the dynamic thermal properties obtained from a scaled-down hot box test can be correlated to its full-size counterpart when certain conditions are met. The scaled-down hot box test is demonstrated to be an effective yet economical alternative to a full-scale test with significantly reduced experimentation cost and turn-around time.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2017

Treatment of the polar coordinate singularity in axisymmetric wave propagation using high-order summation-by-parts operators on a staggered grid

In this paper, we develop a high-order finite difference scheme for axisymmetric wave propagation in a cylindrical conduit filled with a viscous fluid. The scheme is provably stable, and overcomes the difficulty of the polar coordinate singularity in the radial component of the diffusion operator. The finite difference approximation satisfies the principle of summation-by-parts (SBP), which is used to establish stability using the energy method. To treat the coordinate singularity without losing the SBP property of the scheme, a staggered grid is introduced and quadrature rules with weights set to zero at the endpoints are considered. Finally, the accuracy of the scheme is studied both for a model problem with periodic boundary conditions at the ends of the conduit and its practical utility is demonstrated by modeling acoustic-gravity waves in a magmatic conduit.

42 ENGINEERING↗