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,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

Convenient total variation diminishing conditions for nonlinear difference schemes

Convenient conditions for nonlinear difference schemes to be total-variation diminishing (TVD) are reviewed. It is shown that such schemes share the TVD property, provided their numerical fluxes meet a certain positivity condition at extrema values but can be arbitrary otherwise. The conditions are invariant under different incremental representations of the nonlinear schemes, and thus provide a simplified generalization of the TVD conditions due to Harten and others.

Tadmor, Eitan↗
NASA NTRS1990

Fourier analysis of finite element preconditioned collocation schemes

The spectrum of the iteration operator of some finite element preconditioned Fourier collocation schemes is investigated. The first part of the paper analyses one-dimensional elliptic and hyperbolic model problems and the advection-diffusion equation. Analytical expressions of the eigenvalues are obtained with use of symbolic computation. The second part of the paper considers the set of one-dimensional differential equations resulting from Fourier analysis (in the tranverse direction) of the 2-D Stokes problem. All results agree with previous conclusions on the numerical efficiency of finite element preconditioning schemes.

Deville, Michel O.↗
NASA NTRS1991

Modeling flow at the nozzle of a solid rocket motor

The mechanical behavior of a rocket motor internal flow field results in a system of nonlinear partial differential equations which can be solved numerically. The accuracy and the convergence of the solution of the system of equations depends largely on how precisely the sharp gradients can be resolved. An adaptive grid generation scheme is incorporated into the computer algorithm to enhance the capability of numerical modeling. With this scheme, the grid is refined as the solution evolves. This scheme significantly improves the methodology of solving flow problems in rocket nozzle by putting the refinement part of grid generation into the computer algorithm.

Chow, Alan S.↗
NASA NTRS1992

Nonlinearly stable compact schemes for shock calculations

The applications of high-order, compact finite difference methods in shock calculations are discussed. The main concern is to define a local mean which will serve as a reference for introducing a local nonlinear limiting to control spurious numerical oscillations while maintaining the formal accuracy of the scheme. For scalar conservation laws, the resulting schemes can be proven total-variation stable in one space dimension and maximum-norm stable in multiple space dimensions. Numerical examples are shown to verify accuracy and stability of such schemes for problems containing shocks. These ideas can also be applied to other implicit schemes such as the continuous Galerkin finite element methods.

Cockburn, Bernardo↗
NASA NTRS1984

A numerical study of the two- and three-dimensional unsteady Navier-Stokes equations in velocity-vorticity variables using compact difference schemes

A compact finite-difference approximation to the unsteady Navier-Stokes equations in velocity-vorticity variables is used to numerically simulate a number of flows. These include two-dimensional laminar flow of a vortex evolving over a flat plate with an embedded cavity, the unsteady flow over an elliptic cylinder, and aspects of the transient dynamics of the flow over a rearward facing step. The methodology required to extend the two-dimensional formulation to three-dimensions is presented.

Gatski, T. B.↗
DOE OSTI2024

Precision spectroscopy of fast, hot, exotic isotopes using machine-learning-assisted event-by-event Doppler correction

We propose an experimental scheme for performing sensitive, high-precision laser spectroscopy studies on fast exotic isotopes. By inducing a stepwise resonant ionization of the atoms traveling inside an electric field and subsequently detecting the ion and the corresponding electron, time-, and position-sensitive measurements of the resulting particles can be performed. Using a mixture density network, we can leverage this information to predict the initial energy of individual atoms and thus apply a Doppler correction of the observed transition frequencies on an event-by-event basis. We conduct numerical simulations of the proposed experimental scheme and show that kHz-level uncertainties can be achieved for ion beams produced at extreme temperatures ( > 10 8 K), with energy spreads as large as 10 keV and nonuniform velocity distributions. The ability to perform in-flight spectroscopy, directly on highly energetic beams, offers unique opportunities to study short-lived isotopes with lifetimes in the millisecond range and below, produced in low quantities, in hot and highly contaminated environments, without the need for cooling techniques. Such species are of marked interest for nuclear structure, astrophysics, and new physics searches. Published by the American Physical Society 2024

Udrescu, S. M. (ORCID:000000021989576X)↗
DOE OSTI2021

A Bimodal Diagnostic Cloud Fraction Parameterization. Part I: Motivating Analysis and Scheme Description

Cloud fraction parameterizations are beneficial to regional, convection-permitting numerical weather prediction. For its operational regional midlatitude forecasts, the Met Office uses a diagnostic cloud fraction scheme that relies on a unimodal, symmetric subgrid saturation-departure distribution. This scheme has been shown before to underestimate cloud cover and hence an empirically based bias correction is used operationally to improve performance. This first of a series of two papers proposes a new diagnostic cloud scheme as a more physically based alternative to the operational bias correction. The new cloud scheme identifies entrainment zones associated with strong temperature inversions. Additionally, for model grid boxes located in this entrainment zone, collocated moist and dry Gaussian modes are used to represent the subgrid conditions. The mean and width of the Gaussian modes, inferred from the turbulent characteristics, are then used to diagnose cloud water content and cloud fraction. It is shown that the new scheme diagnoses enhanced cloud cover for a given gridbox mean humidity, similar to the current operational approach. It does so, however, in a physically meaningful way. Using observed aircraft data and ground-based retrievals over the southern Great Plains in the United States, it is shown that the new scheme improves the relation between cloud fraction, relative humidity, and liquid water content. An emergent property of the scheme is its ability to infer skewed and bimodal distributions from the large-scale state that qualitatively compare well against observations. A detailed evaluation and resolution sensitivity study will follow in Part II.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1974

Propagation of sound through a real jet flowfield

The sound propagation of harmonic disturbances through a real jet flowfield, which contains mean flow gradients in pressure and velocities, has been studied. A finite-difference approximation for the equations which govern the acoustic disturbance is obtained for a subsonic axisymmetric jet in terms of acoustic pressure and acoustic energy flux. Experimental results are used to describe the mean flowfield of a model jet. The equations are solved numerically using a Newton-type iterative scheme. The directivity at various distances from a point source located at two jet diameters downstream of the jet exit on its center line compared qualitatively well with the numerical result of Schubert (1972) and with the experimental result of Grande (1966).

Liu, C. H.↗
NASA NTRS1987

A mapped finite difference study of noise propagation in nonuniform ducts with mean flow

The primary objective of this work is to study noise propagation in acoustically lined variable area ducts with mean fluid flow. The method of study is numerical in nature and involves a body-fitted grid mapping procedure in conjunction with a factored-implicit finite difference technique. The mean fluid flow model used is two-dimensional, inviscid, irrotational, incompressible, and nonheat conducting. Fully-coupled solutions of the linearized gasdynamic equations are obtained for both positive and negative Mach numbers as well as for hard and soft wall conditions. The factored-implicit finite difference technique used did give rise to short wavelength perturbations, but these were dampened by the introduction of higher order artificial dissipation terms into the scheme. Results compared favorably with available numerical and experimental data.

Raad, Peter E.↗
NASA NTRS1988

Computation of hypersonic flow through a narrow expansion slot

The compressible Navier-Stokes equations are numerically solved for hypersonic flow over a three-dimensional ramp with a narrow expansion slot. In a two-dimensional test case, it is shown that a higher order scheme is needed to avoid too much numerical dissipation and hence too much total-pressure loss. The calculation demonstrates the role of viscosity in the expansion process of hypersonic flow through the narrow slot. Cases with various wall temperatures and slot widths are studied. Calculations show that wall cooling reduces the thickness of boundary layer, and hence increases the flow expansion substantially. In the lower portion of the slot, inviscidly, the flow is dominated by a highly expanded low density fluid, and viscously, by a viscous layer. As a direct consequence, as soon as the wedge angle is large enough, the mass flux through the slot is almost constant.

Hung, Ching-Mao↗
NASA NTRS1987

Time-dependent approximation schemes for some problems of parameter estimation in distributed systems

A parameter estimation method that can be used to estimate functional parameters in delay differential equations and moving boundary problems is discussed. In either problem, the original model equation (which is infinite-dimensional) is approximated by a system of ordinary differential equations that can be solved numerically in an efficient way. The approximation scheme is based on time-dependent spline elements. For the delay equation with time-varying delay, convergence results are presented that indicate the estimates obtained using the approximating system. Numerical test examples converge in some sense to a best-fit parameter for the original system are included by means of which time-varying and state-dependent delays and a time-varying diffusion coefficient in a one-phase, one-dimensional Stefan problem are estimated.

Murphy, K. A.↗
NASA NTRS2018

Predicting Global Ground Geoelectric Field With Coupled Geospace and Three-Dimensional Geomagnetic Induction Models

We forecast the global effects of space weather on the geoelectric and geomagnetic fields using a novel combination of methods. We use a realistic three-dimensional (3-D) model of Earth's electrical conductivity and a realistic representation of magnetospheric and ionospheric current systems. Our scheme involves the following steps: (1) We run a global magnetohydrodynamic model of the magnetosphere coupled to an electrostatic model of the ionosphere. (2) We calculate a global time series of the ground magnetic field resulting from the ionospheric, field-aligned, and magnetospheric currents of the global magnetohydrodynamic model. (3) We approximate this external field by an equivalent source current flowing in a thin shell above Earth. (4) We calculate a global time series of geoelectric and geomagnetic fields from the equivalent current and a 3-D conductivity model of Earth that also takes into account the coast effect due to large horizontal conductivity gradient. We verify our implementation by comparing the results against known analytic and numeric solutions, and then apply our scheme to the geomagnetic storm of 14 and 15 December 2006. In particular, we show that accounting for 3-D structure of Earth's conductivity results in significantly enhanced geoelectric field at large lateral gradients of conductivity, especially in coastal regions, both at middle and high latitudes. In the studied geomagnetic storm the largest values of 3-D geoelectric field are detected at high latitudes reaching 2.5 volts per kilometer and the 3-D effect extends inland by a few hundred kilometers.

Honkonen, I.↗
NASA NTRS1995

Application of the space-time conservation element and solution element method to two-dimensional advection-diffusion problems

The existing 2-D alpha-mu scheme and alpha-epsilon scheme based on the method of space-time conservation element and solution element, which were constructed for solving the linear 2-D unsteady advection-diffusion equation and unsteady advection equation, respectively, are tested. Also, the alpha-epsilon scheme is modified to become the V-E scheme for solving the nonlinear 2-D inviscid Burgers equation. Numerical solutions of six test problems are presented in comparison with their exact solutions or numerical solutions obtained by traditional finite-difference or finite-element methods. It is demonstrated that the 2-D alpha-mu, alpha-epsilon, and nu-epsilon schemes can be used to obtain numerical results which are more accurate than those based on some of the traditional methods but without using any artificial tuning in the computation. Similar to the previous 1-D test problems, the high accuracy and simplicity features of the space-time conservation element and solution element method have been revealed again in the present 2-D test results.

Wang, Xiao-Yen↗
NASA NTRS1995

Efficient implementation of weighted ENO schemes

In this paper, we further analyze, test, modify and improve the high order WENO (weighted essentially non-oscillatory) finite difference schemes of Liu, Osher and Chan. It was shown by Liu et al. that WENO schemes constructed from the r-th order (in L1 norm) ENO schemes are (r+1)-th order accurate. We propose a new way of measuring the smoothness of a numerical solution, emulating the idea of minimizing the total variation of the approximation, which results in a 5-th order WENO scheme for the case r = 3, instead of the 4-th order with the original smoothness measurement by Liu et al. This 5-th order WENO scheme is as fast as the 4-th order WENO scheme of Liu et al., and both schemes are about twice as fast as the 4-th order ENO schemes on vector supercomputers and as fast on serial and parallel computers. For Euler systems of gas dynamics, we suggest computing the weights from pressure and entropy instead of the characteristic values to simplify the costly characteristic procedure. The resulting WENO schemes are about twice as fast as the WENO schemes using the characteristic decompositions to compute weights, and work well for problems which do not contain strong shocks or strong reflected waves. We also prove that, for conservation laws with smooth solutions, all WENO schemes are convergent. Many numerical tests, including the 1D steady state nozzle flow problem and 2D shock entropy wave interaction problem, are presented to demonstrate the remarkable capability of the WENO schemes, especially the WENO scheme using the new smoothness measurement, in resolving complicated shock and flow structures. We have also applied Yang's artificial compression method to the WENO schemes to sharpen contact discontinuities.

Jiang, Guang-Shan↗
DOE OSTI2023

An introduction to financial option pricing on a qudit-based quantum computer

The financial sector is anticipated to be one of the first industries to benefit from the increased computational power of quantum computers, in areas such as portfolio optimisation and risk management to financial derivative pricing. Financial mathematics, and derivative pricing in particular, are not areas quantum physicists are traditionally trained in despite the fact that they often have the raw technical skills needed to understand such topics. On the other hand, most quantum algorithms have largely focused on qubits, which are comprised of two discrete states, as the information carriers. However, discrete higher-dimensional qudits, in addition to possibly possessing increased noise robustness and allowing for novel error correction protocols in certain hardware implementations, also have logarithmically greater information storage and processing capacity. In the current NISQ era of quantum computing, a wide array of hardware paradigms are still being studied and any potential advantage a platform offers is worth exploring. Here we introduce the basic concepts behind financial derivatives for the unfamiliar enthusiast as well as outline in great detail the quantum algorithm routines needed to price a European option, the simplest derivative. This is done within the context of a quantum computer comprised of qudits and employing the natural higher-dimensional analogue of a qubit-based pricing algorithm with its various subroutines. From these pieces, one should relatively easily be able to tailor the scheme to more complex, realistic financial derivatives. Finally, the entire stack is numerically simulated with the results demonstrating how the qudit-based scheme's payoff quickly approaches that of both a similarly-resourced classical computer as well as the true payoff, within error, for a modest increase in qudit dimension.

Bornman, Nicholas↗
NASA NTRS1993

A numerical study of the effect of geometry variation, turbulence models, and dissipation on the flow past control surfaces

The steady 3D viscous flow past the ONERA M6 wing and a slender delta wing-body with trailing edge control surfaces has been computed. A cell-centered finite-volume Navier-Stokes patched zonal method has been used for the numerical simulation. Both diagonalized and LUSGS schemes have been implemented. Besides the standard nonplanar zonal interfacing techniques, a new virtual zone capability has been employed. For code validation, the transonic flow past the ONERA M5 wing is calculated for angles-of-attack of 3.06 deg and 5.06 deg and compared with the available experiments. The wing-body computational results are compared with experimental data for both trailing-edge flaps deflected. The experimental flow conditions are M subinfinity = 0.4, a turbulent Reynolds number of 5.41 million based on a mean aerodynamic chord of 25.959 inches, adiabatic wall, and angles-of-attack varying from 0 deg to 23.85 deg. The computational results are presented for the 23.85 deg angle-of-attack case. The effects of the base flow due to a model sting, the varying second and fourth order numerical dissipation, and the turbulence model are all considered.

Chaussee, Denny S.↗
NASA NTRS1991

Calculations of separated 3-D flows with a pressure-staggered Navier-Stokes equations solver

A Navier-Stokes equations solver based on a pressure correction method with a pressure-staggered mesh and calculations of separated three-dimensional flows are presented. It is shown that the velocity pressure decoupling, which occurs when various pressure correction algorithms are used for pressure-staggered meshes, is caused by the ill-conditioned discrete pressure correction equation. The use of a partial differential equation for the incremental pressure eliminates the velocity pressure decoupling mechanism by itself and yields accurate numerical results. Example flows considered are a three-dimensional lid driven cavity flow and a laminar flow through a 90 degree bend square duct. For the lid driven cavity flow, the present numerical results compare more favorably with the measured data than those obtained using a formally third order accurate quadratic upwind interpolation scheme. For the curved duct flow, the present numerical method yields a grid independent solution with a very small number of grid points. The calculated velocity profiles are in good agreement with the measured data.

Kim, S.-W.↗