SEARCH · Engineering Papers

Results for “Numerical approximations”

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 379 records · Page 21

All sources NASA NTRS DOE OSTI

NASA NTRS1994

A test of the adhesion approximation for gravitational clustering

We quantitatively compare a particle implementation of the adhesion approximation to fully nonlinear, numerical 'N-body' simulations. Our primary tool, cross-correlation of N-body simulations with the adhesion approximation, indicates good agreement, better than that found by the same test performed with the Zel'dovich approximation (hereafter ZA). However, the cross-correlation is not as good as that of the truncated Zel'dovich approximation (TZA), obtained by applying the Zel'dovich approximation after smoothing the initial density field with a Gaussian filter. We confirm that the adhesion approximation produces an excessively filamentary distribution. Relative to the N-body results, we also find that: (a) the power spectrum obtained from the adhesion approximation is more accurate that that from ZA to TZA, (b) the error in the phase angle of Fourier components is worse that that from TZA, and (c) the mass distribution function is more accurate than that from ZA or TZA. It appears that adhesion performs well statistically, but that TZA is more accurate dynamically, in the sense of moving mass to the right place.

Melott, Adrian L.↗

NASA NTRS2005

Reliability-Based Control Design for Uncertain Systems

This paper presents a robust control design methodology for systems with probabilistic parametric uncertainty. Control design is carried out by solving a reliability-based multi-objective optimization problem where the probability of violating design requirements is minimized. Simultaneously, failure domains are optimally enlarged to enable global improvements in the closed-loop performance. To enable an efficient numerical implementation, a hybrid approach for estimating reliability metrics is developed. This approach, which integrates deterministic sampling and asymptotic approximations, greatly reduces the numerical burden associated with complex probabilistic computations without compromising the accuracy of the results. Examples using output-feedback and full-state feedback with state estimation are used to demonstrate the ideas proposed.

Crespo, Luis G.↗

NASA NTRS1985

Factorization and reduction methods for optimal control of distributed parameter systems

A Chandrasekhar-type factorization method is applied to the linear-quadratic optimal control problem for distributed parameter systems. An aeroelastic control problem is used as a model example to demonstrate that if computationally efficient algorithms, such as those of Chandrasekhar-type, are combined with the special structure often available to a particular problem, then an abstract approximation theory developed for distributed parameter control theory becomes a viable method of solution. A numerical scheme based on averaging approximations is applied to hereditary control problems. Numerical examples are given.

Burns, J. A.↗

NASA NTRS1986

Convection in the vertical midplane of a horizontal cylinder Comparison of two-dimensional approximations with three-dimensional results

Experimental results and numerical data on the flow characteristics of core-driven and boundary-layer-driven regimes (CDRs) and (BLDRs) are compared to two-dimensional solutions for heated convective flow in the vertical midplane of a horizontal cylinder. The data were obtained from laser Doppler velocimetry (LDA) measurements and from three-dimensional numerical solutions. It is shown that the approximate solutions gave the correct Rayleigh-number-dependences and distributions for the velocities in the CDR and BLDR regimes. However, the transition between the CDR and the BLDR, as well as the magnitude of the distribution of velocity components, were found to be dependent on the type of two-dimensional approximation selected. In the case of CDR, core velocities were overestimated by as much as 30 percent by the two-dimensional approximations; and in the case of BLDR, the core velocities were overestimated by 40-50 percent as compared to the measured velocities. A schematic illustration of the flow patterns in the cylinder is provided.

Bontoux, P.↗

NASA NTRS2018

Tetrahedral-Mesh Simulations of Shock-Turbulence Interaction

Despite decades of development of unstructured mesh methods, direct numerical simulations (DNS) of turbulent flows are still predominantly performed on structured or unstructured hexahedral meshes with high-order finite-difference methods, weighted essentially nonoscillatory (WENO) schemes, or hybrid schemes formed by their combinations. Tetrahedral meshes offer easy mesh generation and adaptation around complex geometries and the potential of an orientation-free grid that would benefit the isotropic nature of small-scale dissipation, as well as the solution accuracy of intermediate scales. To advance the state of the art of unstructured-mesh simulation capabilities for shock/turbulence interaction, DNS using pure tetrahedral meshes are carried out with the space-time conservation element, solution element (CESE) method in this research. By its design, the CESE method is constructed based on a non-dissipative scheme and is a genuinely multidimensional numerical framework that is free from the use of an approximate Riemann-solver. The numerical framework also provides the ability to add numerical dissipation (the nondissipative scheme acts as the reference state like that of the reversible state in thermodynamics) when needed (with justification from mathematics/physics). The above-mentioned features along with the CESE method's consistent shock-capturing approach and strong enforcement of flux conservation in spacetime offers a novel method to accurately simulate turbulent flows and their interaction with shocks using tetrahedral meshes. Two canonical problems, namely, isotropic turbulence interaction with a normal shock and a Mach 2.9 turbulent boundary layer flow over a 24deg compression corner are investigated in this study. Computational results show reasonably good agreement with experimental data and results from structured-mesh, high-order simulations available in the literature. Successful validation of these canonical problems demonstrated here paves the way for future high-fidelity supersonic flow simulations involving complex-geometries.

Venkatachari, Balaji Shankar↗

NASA NTRS1983

Parameter identification in continuum models

Approximation techniques for use in numerical schemes for estimating spatially varying coefficients in continuum models such as those for Euler-Bernoulli beams are discussed. The techniques are based on quintic spline state approximations and cubic spline parameter approximations. Both theoretical and numerical results are presented.

Banks, H. T.↗

NASA NTRS1983

Parameter identification in continuum models

Banks, H. T.↗

DOE OSTI2021

Low-temperature T -linear resistivity due to umklapp scattering from a critical mode

In this study, we consider the transport properties of a model of fermions scattered by a critical bosonic mode. The mode is overdamped and scattering is mainly in the forward direction. Such a mode appears at the quantum critical point for an electronic nematic phase transition and in gauge theories for a U(1) spin liquid. It leads to a short fermion lifetime, violating Landau's criterion for a Fermi liquid. In spite of this, transport can be described by a Boltzmann equation. We include momentum relaxation by umklapp scattering, supplemented by weak impurity scattering. We find that above a very low temperature which scales with Δ q 3 , where Δ q is the minimum umklapp scattering vector, the resistivity is linear in T with a coefficient which is independent of the amount of disorder. This picture holds until the temperature dependent part of the resistivity exceeds that due to impurity scattering. We compare the relaxation time approximation with an exact numerical solution of the Boltzmann equation. Surprisingly we find that, unlike the resistivity, the Hall coefficient strongly deviates from the relaxation time approximation and shows a strong reduction with increasing temperature. We comment on possible comparisons with experiments on high T c cuprates.

36 MATERIALS SCIENCE↗

NASA NTRS1984

Stability and Convergence of Underintegrated Finite Element Approximations

The effects of underintegration on the numerical stability and convergence characteristics of certain classes of finite element approximations were analyzed. Particular attention is given to hourglassing instabilities that arise from underintegrating the stiffness matrix entries and checkerboard instabilities that arise from underintegrating constrain terms such as those arising from incompressibility conditions. A fundamental result reported here is the proof that the fully integrated stiffness is restored in some cases through a post-processing operation.

Oden, J. T.↗

DOE OSTI2021

Heat transport in liquid water from first-principles and deep neural network simulations

In this work, we compute the thermal conductivity of water within linear response theory from equilibrium molecular dynamics simulations, by adopting two different approaches. In one, the potential energy surface (PES) is derived on the fly from the electronic ground state of density functional theory (DFT) and the corresponding analytical expression is used for the energy flux. In the other, the PES is represented by a deep neural network (DNN) trained on DFT data, whereby the PES has an explicit local decomposition and the energy flux takes a particularly simple expression. By virtue of a gauge invariance principle, established by Marcolongo, Umari, and Baroni, the two approaches should be equivalent if the PES were reproduced accurately by the DNN model. We test this hypothesis by calculating the thermal conductivity, at the GGA (PBE) level of theory, using the direct formulation and its DNN proxy, finding that both approaches yield the same conductivity, in excess of the experimental value by approximately 60%. Besides being numerically much more efficient than its direct DFT counterpart, the DNN scheme has the advantage of being easily applicable to more sophisticated DFT approximations, such as meta-GGA and hybrid functionals, for which it would be hard to derive analytically the expression of the energy flux. We find in this way that a DNN model, trained on meta-GGA (SCAN) data, reduces the deviation from experiment of the predicted thermal conductivity by about 50%, leaving the question open as to whether the residual error is due to deficiencies of the functional, to a neglect of nuclear quantum effects in the atomic dynamics, or, likely, to a combination of the two.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1979

Population statistics of faint stellar and non-stellar objects

A disc and halo population model is constructed to fit star counts and color data down to V approximately 23 at absolute value of b = 90 deg. This model is used to predict star counts and colors down to V approximately 30. Deviations from these extrapolated relationships provide constraints on the number of faint quasars and black dwarf stars. It is shown that extra-galactic globular clusters start contributing significantly to star counts at V approximately 25 and are more numerous than stars for V 31. Morphological studies of galaxies with approximately 0.5, were made with the space telescope. Significant constraints on theoretical models that describe the evolution of clusters of galaxies are provided.

Vandenbergh, S.↗

DOE OSTI2021

Quasinormal modes and shadow of a Schwarzschild black hole with GUP

Highlights: • The physics of quantum corrected Schwarzschild Black Hole via GUP is investigated. • Quasinormal Modes (QNM) are obtained via WKB and numeric methods. • The relationship between shadow and QNM in the eikonal limit is established. • The quadratic and linear GUP parameters affect the results in a clear way. We consider quantum corrections for the Schwarzschild black hole metric by using the generalized uncertainty principle (GUP) to investigate quasinormal modes, shadow and their relationship in the eikonal limit. We calculate the quasinormal frequencies of the quantum-corrected Schwarzschild black hole by using the sixth-order Wentzel–Kramers–Brillouin (WKB) approximation, and also perform a numerical analysis that confirms the results obtained from this approach. We also find that the shadow radius is nonzero even at very small mass limit for finite GUP parameter.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1992

An implicit Navier-Stokes code for turbulent flow modeling

This paper presents a numerical approach to calculating turbulent flows employing advanced turbulence models. The main features include a line-by-line Gauss-Seidel algorithm using Roe's approximate Riemann solver, TVD numerical schemes, implicit boundary conditions and a decoupled turbulence-model solver. Based on the problems tested so far, the method has consistently demonstrated its ability in offering accuracy, boundedness and a fast rate of convergence to steady-state solution.

Huang, P. G.↗

NASA NTRS2003

MHD Turbulence at Moderate Magnetic Reynolds Number

In the present article we will consider the decay of MHD turbulence under the influence of a strong external magnetic field at moderate magnetic Reynolds numbers. Typical values of R(sub m) that are considered here range from R(sub m) approx. 0.1 to R(sub m) approx. 20. As a comparison, the initial kinetic Reynolds number common to all our simulations is Re(sub L) = 199. This means that the range of Prandtl numbers explored is 5 x 10(exp -4) to 10(exp -1). Our motivation is mainly to exhibit how the transition from the QS approximation to FMHD occurs. At the lowest values of R(sub m) studied here, the QS approximation is shown to model the flow faithfully. However, for the higher values of R(sub m) considered, it is clearly inadequate but can be replaced by another approximation which will be referred to as the Quasi-Linear (QL) approximation. Another objective of the present study is to describe how variations in the magnetic Reynolds number (while maintaining all other parameters constant) affect the dynamics of the flow. This complements past studies where variations in either the strength of the external magnetic field or the kinetic Reynolds number were considered. This article is organized as follows. In section 2 we recall the definition of the quasi-static approximation. Section 3 is devoted to the description of the numerical experiments performed using the quasi-static approximation and full MHD. In section 4 we describe the quasi-linear approximation and test it numerically against full MHD. A concluding summary is given in section 5.

Knaepen, B.↗

NASA NTRS1991

An efficient method for estimating neighboring steady-state numerical solutions to the Euler equations

The paper concentrates on a linear approximation method for predicting the changes occurring in steady-state numerical solutions of the Euler equations as a consequence of small changes in the independent variables which control the problem. The importance of proper boundary-condition treatment and other issues concerning the problem are covered along with the importance of proper algorithm selection for a fully supersonic inviscid flow. The method is applied to a subsonic nozzle involving variation of the pressure on the outflow boundary and to a supersonic inlet involving variation of the inflow Mach number. In the subsonic test case, the comparisons between the predicted and conventional numerical solutions are shown to be good, while in the supersonic test case, the agreement between the approximation method and conventional numerical solution starts out well but rapidly degenerates at some point in the flowfield as the perturbation of the boundary conditions is increased.

Taylor, Arthur C., III↗

NASA NTRS1973

Nonlinear vibrations of rectangular plates.

A finite-difference method is developed to determine the large amplitude dynamic responses of thin elastic plates subjected to uniform pressure pulse-type loads. Four different sets of boundary conditions are considered. Some specific problems are solved. The results are compared with approximate solutions obtained by Yamaki (1961). The numerical method presented provides an accurate and efficient approximate solution to the problem, and should be useful as a check on other approximate methods. The grid-size and the time-step necessary for obtaining numerical stability depend on the particular problem. For many cases the method converges rapidly and a rather large grid-size and time-step is adequate.

Bayles, D. J.↗

NASA NTRS1990

A numerical study of the separation wavenumber in the two-scale scattering approximation

The choice of the separation wavenumber Kd for the two-scale scattering approximation is examined by the comparison of two-scale calculations with exact numerical solutions using the moment method. It was found that the optimal choice of Kd is near k/2 when attempting to optimize for all incidence angles. When optimizing for specific incidence angles, Kd varied from k/5 at low incidence angles to just below k at higher incidence angles. These numerical results are in rough agreement with previous theoretical derivations of the separation wavenumber.

Durden, Stephen L.↗

NASA NTRS1992

Model reduction using new optimal Routh approximant technique

An optimal Routh approximant of a single-input single-output dynamic system is a reduced-order transfer function of which the denominator is obtained by the Routh approximation method while the numerator is determined by minimizing a time-response integral-squared-error (ISE) criterion. In this paper, a new elegant approach is presented for obtaining the optimal Routh approximants for linear time-invariant continuous-time systems. The approach is based on the Routh canonical expansion, which is a finite-term orthogonal series of rational basis functions, and minimization of the ISE criterion. A procedure for combining the above approach with the bilinear transformation is also presented in order to obtain the optimal bilinear Routh approximants of linear time-invariant discrete-time systems. The proposed technique is simple in formulation and is amenable to practical implementation.

Hwang, Chyi↗