Search NASA⌕ Search

SEARCH · Search NASA

Results for “generalized finite differences”

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.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

Efficient grid generation

Because the governing equations in fluid dynamics contain partial differentials and are too difficult in most cases to solve analytically, these differentials are generally replaced by finite difference terms. These terms contain terms in the solution at nearby states. This procedure discretizes the field into a finite number of states. These states, when plotted, form a grid, or mesh, of points. It is at these states, or field points, that the solution is found. The optimum choice of states, the x, y, z coordinate values, minimizes error and computational time. But the process of finding these states is made more difficult by complex boundaries, and by the need to control step size differences between the states, that is, the need to control the spacing of field points. One solution technique uses a different set of state variables, which define a different coordinate system, to generate the grid more easily. A new method, developed by Dr. Joseph Steger, combines elliptic and hyperbolic partial differential equations into a mapping function between the physical and computational coordinate systems. This system of equations offers more control than either equation provides alone. The Steger algorithm was modified in order to allow bodies with stronger concavities to be used, offering the possibility of generating a single grid about multiple bodies. Work was also done on identifying areas where grid breakdown occurs.

Seki, Rycichi↗

An energy-conserving method for stochastic Maxwell equations with multiplicative noise

In this paper, it is shown that three-dimensional stochastic Maxwell equations with multiplicative noise are stochastic Hamiltonian partial differential equations possessing a geometric structure (i.e. stochastic multi-symplectic conservation law), and the energy of system is a conservative quantity almost surely. We propose a stochastic multi-symplectic energy-conserving method for the equations by using the wavelet collocation method in space and stochastic symplectic method in time. Numerical experiments are performed to verify the excellent abilities of the proposed method in providing accurate solution and preserving energy. The mean square convergence result of the method in temporal direction is tested numerically, and numerical comparisons with finite difference method are also investigated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A computational method for the determination of the response of a linear system

A number of issues are discussed relating to the analysis of a linear damped oscillator equation with a forcing term where the right side of the equation, g(t), is known only at fixed, equal time intervals. (A particular class of such equations has application to civil earthquake engineering.) It is noted that the equation is not a differential equation because the condition of uniqueness does not obtain. A general computational method is presented for calculating x(t), based on an unconditionally stable finite difference technique. The work generalizes the results given in the recent paper of Ly (1984).

Mickens, R. E.↗

Numerical analysis and parametric studies of the buckling of composite orthotropic compression and shear panels

A computer program is presented which was developed for the combined compression and shear of stiffened variable thickness orthotropic composite panels on discrete springs: boundary conditions are general and include elastic boundary restraints. Buckling solutions are obtained by using a newly developed trigonometric finite difference procedure which improves the solution convergence rate over conventional finite difference methods. The classical general shear buckling results which exist only for simply supported panels over a limited range of orthotropic properties, were extended to the complete range of these properties for simply supported panels and, in addition, to the complete range of orthotropic properties for clamped panels. The program was also applied to parametric studies which examine the effect of filament orientation upon the buckling of graphite-epoxy panels. These studies included an examination of the filament orientations which yield maximum shear or compressive buckling strength for panels having all four edges simply supported or clamped over a wide range of aspect ratios. Panels with such orientations had higher buckling loads than comparable, equal weight, thin skinned aluminum panels. Also included among the parameter studies were examinations of combined axial compression and shear buckling and examinations of panels with rotational elastic edge restraints.

Housner, J. M.↗

Uncertainty quantification for generalized Langevin dynamics

Here, we present efficient finite difference estimators for goal-oriented sensitivity indices with applications to the generalized Langevin equation (GLE). In particular, we apply these estimators to analyze an extended variable formulation of the GLE where other well known sensitivity analysis techniques such as the likelihood ratio method are not applicable to key parameters of interest. These easily implemented estimators are formed by coupling the nominal and perturbed dynamics appearing in the finite difference through a common driving noise or common random path. After developing a general framework for variance reduction via coupling, we demonstrate the optimality of the common random path coupling in the sense that it produces a minimal variance surrogate for the difference estimator relative to sampling dynamics driven by independent paths. In order to build intuition for the common random path coupling, we evaluate the efficiency of the proposed estimators for a comprehensive set of examples of interest in particle dynamics. These reduced variance difference estimators are also a useful tool for performing global sensitivity analysis and for investigating non-local perturbations of parameters, such as increasing the number of Prony modes active in an extended variable GLE.

97 MATHEMATICS AND COMPUTING↗

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

In this paper transient fields for antennas with more general geometries are calculated directly using Finite Difference Time Domain 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.↗

The Relation of Finite Element and Finite Difference Methods

Finite element and finite difference methods are examined in order to bring out their relationship. It is shown that both methods use two types of discrete representations of continuous functions. They differ in that finite difference methods emphasize the discretization of independent variable, while finite element methods emphasize the discretization of dependent variable (referred to as functional approximations). An important point is that finite element methods use global piecewise functional approximations, while finite difference methods normally use local functional approximations. A general conclusion is that finite element methods are best designed to handle complex boundaries, while finite difference methods are superior for complex equations. It is also shown that finite volume difference methods possess many of the advantages attributed to finite element methods.

Vinokur, M.↗

L{sup ∞} error bound of conservative compact difference scheme for the generalized symmetric regularized long-wave (GSRLW) equations

In this paper, we design a compact finite difference scheme which preserves the original conservative properties to solve the generalized symmetric regularized long-wave equations. The existence of the difference solution is proved by the Brouwer fixed-point theorem. Applying the discrete energy method, the convergence and stability of the difference scheme is obtained, and its numerical convergence order is O(τ{sup 2}+h{sup 4}) in the L{sup ∞}-norm for u and ρ. For computing the nonlinear algebraic system generated by the compact scheme, a decoupled iterative algorithm is constructed and proved to be convergent. Numerical experiment results show that the theory is accurate and the method is efficient and reliable.

97 MATHEMATICS AND COMPUTING↗

Solution of three-dimensional time-dependent viscous flows

A procedure for solving three-dimensional, time-dependent turbulent flows is presented. The consistently split Linearized Block Implicit (LBI) scheme is used in conjunction with the QR Operator scheme to solve an approximate form of the Navier-Stokes equations in generalized nonorthogonal coordinates employing physical velocity components. Results of computations for both second order finite differences and the fourth order generalized Operator Compact Implicit (OCI) schemes are presented.

Weinberg, B. C.↗

Computation of three-dimensional nozzle-exhaust flow fields with the GIM code

A methodology is introduced for constructing numerical analogs of the partial differential equations of continuum mechanics. A general formulation is provided which permits classical finite element and many of the finite difference methods to be derived directly. The approach, termed the General Interpolants Method (GIM), can combined the best features of finite element and finite difference methods. A quasi-variational procedure is used to formulate the element equations, to introduce boundary conditions into the method and to provide a natural assembly sequence. A derivation is given in terms of general interpolation functions from this procedure. Example computations for transonic and supersonic flows in two and three dimensions are given to illustrate the utility of GIM. A three-dimensional nozzle-exhaust flow field is solved including interaction with the freestream and a coupled treatment of the shear layer. Potential applications of the GIM code to a variety of computational fluid dynamics problems is then discussed in terms of existing capability or by extension of the methodology.

Spradley, L. W.↗

Assessment of advances in structural analysis for space shuttle

The latest technology for structural analysis in relation to the design tasks that lie ahead for the space shuttle is reviewed. For shell-of-revolution structures, the analysis can be formulated as a one-dimensional problem which is readily solved by using finite-difference or numerical-integration techniques. For more general asymmetric shells, a two-dimensional formulation is required. However, the governing equations are readily formulated and are amenable to solution by finite-difference techniques. For a completely general structural arrangement, such as structural frameworks, recourse is usually made to discretized formulations using finite elements. Of course, the finite-element programs could be used for shell structures, but at a loss in accuracy and increase in computer time compared with the special purpose programs.

Dixon, S. C.↗

An implicit finite-difference code for a two-equation turbulence model for three-dimensional flows

An implicit finite difference code was developed which solves the transport equations for the turbulence kinetic energy and its dissipation rate in generalized coordinates in three dimensions. The finite difference equations are solved using the Beam-Warming algorithm. The kinetic energy-dissipation code, KEM, provides the closure; i.e., the turbulent viscosity for calculation of either compressible or incompressible flows. Turbulent internal flow over a backward-facing step has been calculated using the present code in conjunction with the Incompressible Navier-Stokes Code, INS3D. The results are in good agreement with experiments and two dimensional computations of other researchers.

Kaul, U. K.↗

Higher-order numerical methods derived from three-point polynomial interpolation

Higher-order collocation procedures resulting in tridiagonal matrix systems are derived from polynomial spline interpolation and Hermitian finite-difference discretization. The equations generally apply for both uniform and variable meshes. Hybrid schemes resulting from different polynomial approximations for first and second derivatives lead to the nonuniform mesh extension of the so-called compact or Pade difference techniques. A variety of fourth-order methods are described and this concept is extended to sixth-order. Solutions with these procedures are presented for the similar and non-similar boundary layer equations with and without mass transfer, the Burgers equation, and the incompressible viscous flow in a driven cavity. Finally, the interpolation procedure is used to derive higher-order temporal integration schemes and results are shown for the diffusion equation.

Rubin, S. G.↗

A three-dimensional incompressible flow simulation method and its application to the Space Shuttle main engine. II Turbulent flow

An implicit finite difference code cast in general curvilinear coordinates is further developed for three-dimensional incompressible turbulent flows. The code is based on the method of pseudocompressibility and utilizes the Beam and Warming implicit approximate factorization algorithm to achieve computational efficiency. A multiple-zone method is further extended to include composite-grids to overcome the excessive computer memory required for solving turbulent flows in complex three-dimensional geometries. A simple turbulence model is proposed for internal flows. The code is being used for the Space Shuttle Main Engine (SSME) internal flow analyses.

Chang, J. L. C.↗

Computational and experimental study of turbulent flows in 180 degree bends

A theoretical and experimental investigation of turbulent flows inside a 180-deg strongly curved channel is included in the present study. A 10 by 100 centimeter water channel with a 10 centimeter radius turn was employed for the experimental study. Data measured are velocity profiles and distributions of longitudinal turbulence intensity. Numerical investigations are performed by solving the incompressible Navier-Stokes equations and turbulence models using a finite difference flow solver in generalized curvilinear coordinate systems. A widely used two-equation turbulence model and its modified version, which includes the effect of strong curvature on the turbulence length scale through a relation of algebraic Reynolds stress model, are employed to predict the turbulence quantities inside the flow field. Improvements on the predictions are obtained using the present approach.

Chen, Y. S.↗

Numerical study of laminar and turbulent flows inside a turnaround duct with and without guide vanes

The purpose of this study is to examine in detail incompressible laminar and turbulent flows inside a turnaround duct with and without guide vanes and to investigate the effects of vanes on the flow characteristics. To perform this study, an implicit finite difference code cast in general curvilinear coordinates is further developed. The code is based on the method of pseudo-compressibility and utilize ADI or implicit approximate factorization algorithm to achieve computational efficiency. Method of segmental sweeping is developed to overcome the multiple-zone problem due to imposition of guide vanes. In the present study, several test cases have been computed. These test cases include laminar and turbulent flows inside a turnaround duct without and with two or three guide vanes. The study reveals that: (1) there exists large recirculation zones inside the duct if no vanes are present; (2) properly shaped and positional guide vanes are effective in eliminating flow separation; and (3) laminar and turbulent flows have similar flow features. But turbulent flow has less total pressure drop.

Lin, S.-J.↗