Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical schemes”

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 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
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 NTRS1991

On the application of ENO scheme with subcell resolution to conservation laws with stiff source terms

Two approaches are used to extend the essentially non-oscillatory (ENO) schemes to treat conservation laws with stiff source terms. One approach is the application of the Strang time-splitting method. Here the basic ENO scheme and the Harten modification using subcell resolution (SR), ENO/SR scheme, are extended this way. The other approach is a direct method and a modification of the ENO/SR. Here the technique of ENO reconstruction with subcell resolution is used to locate the discontinuity within a cell and the time evolution is then accomplished by solving the differential equation along characteristics locally and advancing in the characteristic direction. This scheme is denoted ENO/SRCD (subcell resolution - characteristic direction). All the schemes are tested on the equation of LeVeque and Yee (NASA-TM-100075, 1988) modeling reacting flow problems. Numerical results show that these schemes handle this intriguing model problem very well, especially with ENO/SRCD which produces perfect resolution at the discontinuity.

Chang, Shih-Hung↗
NASA NTRS1989

Efficient implementation of essentially non-oscillatory shock-capturing schemes. II

The present elaboration of essentially nonoscillatory (ENO) shock-capturing schemes proceeds with a novel, simplified expression for the ENO construction procedure having its basis in numerical fluxes rather than cell-averages. The ENO-local Lax-Friedrichs and ENO-Roe schemes introduced yield sharper shock transitions and greater overall accuracy than previous ENO-scheme implementations. Numerical fluxes and TVD Runge-Kutta time discretizations are used to apply Harten's subcell-resolution concept (1986) and Yang's artificial compression method to the current ENO schemes, in the interest of sharper contact discontinuities.

Shu, Chi-Wang↗
DOE OSTI2021

Medium-induced radiative kernel with the Improved Opacity Expansion

We calculate the fully differential medium-induced radiative spectrum at next-to-leading order (NLO) accuracy within the Improved Opacity Expansion (IOE) framework. This scheme allows us to gain analytical control of the radiative spectrum at low and high gluon frequencies simultaneously. The high frequency regime can be obtained in the standard opacity expansion framework in which the resulting power series diverges at the characteristic frequency ω c ~ q^L 2 . In the IOE, all orders in opacity are resumed systematically below ωc yielding an asymptotic series controlled by logarithmically suppressed remainders down to the thermal scale T « ω c , while matching the opacity expansion at high frequency. Furthermore, we demonstrate that the IOE at NLO accuracy reproduces the characteristic Coulomb tail of the single hard scattering contribution as well as the Gaussian distribution resulting from multiple soft momentum exchanges. Finally, we compare our analytic scheme with a recent numerical solution, that includes a full resummation of multiple scatterings, for LHC-inspired medium parameters. Furthermore, we find a very good agreement both at low and high frequencies showcasing the performance of the IOE which provides for the first time accurate analytic formulas for radiative energy loss in the relevant perturbative kinematic regimes for dense media.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Structural stability and artificial buckling modes in topology optimization

Abstract This paper demonstrates how a strain energy transition approach can be used to remove artificial buckling modes that often occur in stability constrained topology optimization problems. To simulate the structural response, a nonlinear large deformation hyperelastic simulation is performed, wherein the fundamental load path is traversed using Newton’s method and the critical buckling load levels are estimated by an eigenvalue analysis. The goal of the optimization is to minimize displacement, subject to constraints on the lowest critical buckling loads and maximum volume. The topology optimization problem is regularized via the Helmholtz PDE-filter and the method of moving asymptotes is used to update the design. The stability and sensitivity analyses are outlined in detail. The effectiveness of the energy transition scheme is demonstrated in numerical examples.

Dalklint, Anna (ORCID:0000000346195205)↗
DOE OSTI2021

Elastic Wave Propagation in Curvilinear Coordinates with Mesh Refinement Interfaces by a Fourth Order Finite Difference Method

In this work, we develop a fourth order accurate finite difference method for the three dimensional elastic wave equation in isotropic media with the piecewise smooth material property. In our model, the material property can be discontinuous at curved interfaces. The governing equations are discretized in second order form on curvilinear meshes by using a fourth order finite difference operator satisfying a summation-by-parts property. The method is energy stable and high order accurate. The highlight is that mesh sizes can be chosen according to the velocity structure of the material so that computational efficiency is improved. At the mesh refinement interfaces with hanging nodes, physical interface conditions are imposed by using ghost points and interpolation. With a fourth order predictor-corrector time integrator, the fully discrete scheme is energy conserving. Numerical experiments are presented to verify the fourth order convergence rate and the energy conserving property.

58 GEOSCIENCES↗
DOE OSTI2021

Sparse Approximate Multifrontal Factorization with Butterfly Compression for High-Frequency Wave Equations

In this work, we present a fast and approximate multifrontal solver for large-scale sparse linear systems arising from finite-difference, finite-volume or finite-element discretization of high-frequency wave equations. The proposed solver leverages the butterfly algorithm and its hierarchical matrix extension for compressing and factorizing large frontal matrices via graph-distance guided entry evaluation or randomized matrix-vector multiplication-based schemes. Complexity analysis and numerical experiments demonstrate $\mathcal{O}(N\log^2 N)$ computation and $\mathcal{O}(N)$ memory complexity when applied to an $N\times N$ sparse system arising from 3D high-frequency Helmholtz and Maxwell problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

On the Efficient Evaluation of the Exchange Correlation Potential on Graphics Processing Unit Clusters

The predominance of Kohn–Sham density functional theory (KS-DFT) for the theoretical treatment of large experimentally relevant systems in molecular chemistry and materials science relies primarily on the existence of efficient software implementations which are capable of leveraging the latest advances in modern high-performance computing (HPC). With recent trends in HPC leading toward increasing reliance on heterogeneous accelerator-based architectures such as graphics processing units (GPU), existing code bases must embrace these architectural advances to maintain the high levels of performance that have come to be expected for these methods. In this work, we purpose a three-level parallelism scheme for the distributed numerical integration of the exchange-correlation (XC) potential in the Gaussian basis set discretization of the Kohn–Sham equations on large computing clusters consisting of multiple GPUs per compute node. In addition, we purpose and demonstrate the efficacy of the use of batched kernels, including batched level-3 BLAS operations, in achieving high levels of performance on the GPU. We demonstrate the performance and scalability of the implementation of the purposed method in the NWChemEx software package by comparing to the existing scalable CPU XC integration in NWChem.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1972

Comparison of characteristic schemes for three-dimensional, steady, isentropic flow.

The accuracy and order of three bicharacteristic schemes based on the numerical method of characteristics for three-dimensional, steady, isentropic, supersonic flow were compared based on numerical calculations. Comparisons were made for source flows and for flows in contoured axisymmetric nozzles. The first two methods were tetrahedral line networks, differing only in the manner in which the cross derivatives were evaluated. The third method was the pentahedral line network proposed by Butler. All three methods were numerically stable and produced efficient and accurate solutions. A discussion of the relative merits and disadvantages of each method is presented.

Cline, M. C.↗
NASA NTRS1973

Nonlinear analysis of a shock-loaded membrane.

Results from a computer method for analyzing the unsteady interaction of a fluid stream and a flat circular elastic membrane are presented. The loading on the membrane is assumed to be caused by the firing of a shock tube. The fluid pressures and velocities are determined from a scheme based on the numerical method of characteristics, and the membrane is analyzed using exact relations for membrane strain. The interactive solution is found to give peak stresses 40% lower than a solution which assumes a pressure invariant in space and time.

Madden, R.↗
NASA NTRS1974

On the remote sensing of directional wave spectra and surface winds

The potential benefits of remote SEASAT-A satellite sensing of wind directional wave spectra by active microwave systems are elaborated. It is shown that a combined effort which uses supplemental satellite information in numerical forecasting and hindcasting schemes upgrades the accuracy of existing forecasting methods of an order of magnitude.

Shemdin, O. H.↗
NASA NTRS1977

Recent developments in the Green's function method

A recent computational development on the Green's function method (the method used in the computer program SOUSSA: Steady, Oscillatory and Unsteady Subsonic and Supersonic Aerodynamics) is presented. A scheme consisting of combined numerical (Gaussian quadrature) and analytical procedures for the evaluation of the source and doublet integrals used in the program is presented. This combination results in 80 to 90% reduction in computer time.

Tseng, K.↗
NASA NTRS1980

Boundary-fitted coordinates for regions with highly curved boundaries and reentrant boundaries

A procedure was developed using the differential equation approach for generating boundary fitted coordinates for regions with highly curved boundaries as well as reentrant boundaries, such as those encountered in breaking surface waves. The resulting coordinates are nearly orthogonal and can provide adequate resolution even in the reentrant region. Consistent treatment of end boundaries and the use of a systematic initialization scheme and advanced implicit numerical solution techniques make the procedure highly efficient. The method developed for implicit enforcement of the periodicity boundary condition should be beneficial in the analysis of turbomachinery flow applications.

Ghia, U.↗
NASA NTRS1984

A Shear Deformable Shell Element for Laminated Composites

A three-dimensional element based on the total Lagrangian description of the motion of a layered anisotropic composite medium is developed, validated, and used to analyze layered composite shells. The element contains the following features: geometric nonlinearity, dynamic (transient) behavior, and arbitrary lamination scheme and lamina properties. Numerical results of nonlinear bending, natural vibration, and transient response are presented to illustrate the capabilities of the element.

Chao, W. C.↗
NASA NTRS1983

Large deformation analysis of layered composite shells

A three-dimensional element based on the total Lagrangian description of large deformations of a layered anisotropic composite medium is developed, validated, and used to analyze layered composite shells. The element contains the following features: geometric nonlinearity, dynamic (transient) behavior, and arbitrary lamination scheme and lamina properties. Numerical results of nonlinear bending, natural vibration, and transient response are presented to illustrate the capabilities of the element.

Chao, W. C.↗
NASA NTRS1982

Combustion hot section technology

The Turbine Engine Hot Section Technology Combustion Program is briefly described. The overall objective of the project is to develop and verify improved and more accurate analysis methods for increasing the ability to design with confidence the combustion system for advanced aircraft turbine engines. The approach is to first assess and evaluate existing combustor aerothermal analysis models by means of a contracted effort initiated during FY-82. The program also includes both analytical and experimental research efforts in the areas of aerothermal modeling and liner cyclic life. It is expected that the combustor model development effort will generate improved understanding in the areas of high pressure flame radiation characteristics, model numerical methods and solution schemes, complex geometrical boundary conditions, fuel spray - flow field interactions, combustion kinetics, flow and mixing of dilution jets, turbulence and heat transfer, and soot and carbon formation.

Ercegovic, D. B.↗
NASA NTRS1984

Analysis of laminated composite shells using a degenerated 3-D element

A special three-dimensional element based on the total Lagrangian description of the motion of a layered anisotropic composite medium is developed, validated and employed to analyze laminated anisotropic composite shells. The element contains the following features: geometric nonlinearity, dynamic (transient) behavior and arbitrary lamination scheme and lamina properties. Numerical results of nonlinear bending, natural vibration, and transient response are presented to illustrate the capabilities of the element.

Chao, W. C.↗
NASA NTRS1986

Adaptive finite element methods for compressible flow problems

Some recent work on adaptive FEMs for solving transient Euler equations in two-dimensional domains is summarized. The formulation of an FEM model of the Euler equations is shown, and the application of the adaptive strategies to data management schemes is addressed. Sample numerical results from the application of the model and strategies to the flow over a step and to transient cases are given.

Oden, J. T.↗