SEARCH · Engineering Papers

Results for “Numerical Methods”

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,279 records · Page 71

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Vibration and noise analysis of a gear transmission system

This paper presents a comprehensive procedure to predict both the vibration and noise generated by a gear transmission system under normal operating conditions. The gearbox vibrations were obtained from both numerical simulation and experimental studies using a gear noise test rig. In addition, the noise generated by the gearbox vibrations was recorded during the experimental testing. A numerical method was used to develop linear relationships between the gearbox vibration and the generated noise. The hypercoherence function is introduced to correlate the nonlinear relationship between the fundamental noise frequency and its harmonics. A numerical procedure was developed using both the linear and nonlinear relationships generated from the experimental data to predict noise resulting from the gearbox vibrations. The application of this methodology is demonstrated by comparing the numerical and experimental results from the gear noise test rig.

Choy, F. K.↗

NASA NTRS1993

Vibration and noise analysis of a gear transmission system

Choy, F. K.↗

NASA NTRS1999

Some Continuum Aspects of Data Assimilation

A long-sought goal in data assimilation is to build more fully the time dimension into the data assimilation process. Truly four-dimensional data assimilation requires evolving second-moment information, namely the estimation error covariance, along with the state estimate itself. The continuum evolution equations for the estimation error covariance constitute a system of partial differential equations in six space dimensions, forced by a model error covariance and by the observations themselves. The high dimensionality of this system poses a difficult computational problem. Numerous methods have been proposed for approximate, discrete solution of this system of equations. Four-dimensional variational schemes solve these equations indirectly, while schemes based on Kalman filtering solve the equations more directly, usually on a small subspace of the full six-dimensional space. In both cases, most proposed solution methods are not derived from the continuum covariance evolution equations themselves. Instead, most methods simply inherit the discretization used for evolving the state estimate, either in the form of a discrete tangent linear model or, in the case of ensemble schemes, in the form of a fully nonlinear discrete model. In this lecture we show that solutions of the continuum covariance evolution equations possess simple properties that are not readily satisfied upon discretization of the equations. we give numerical examples illustrating that "inherited" discretizations sometimes fail to produce meaningfully accurate solutions. Finally, we suggest alternative discretization methods that may yield more faithful approximate solutions of the continuum problem.

Cohn, Stephen E.↗

NASA NTRS1982

Numerical techniques in radiative heat transfer for general, scattering, plane-parallel media

The study of radiative heat transfer with scattering usually leads to the solution of singular Fredholm integral equations. The present paper presents an accurate and efficient numerical method to solve certain integral equations that govern radiative equilibrium problems in plane-parallel geometry for both grey and nongrey, anisotropically scattering media. In particular, the nongrey problem is represented by a spectral integral of a system of nonlinear integral equations in space, which has not been solved previously. The numerical technique is constructed to handle this unique nongrey governing equation as well as the difficulties caused by singular kernels. Example problems are solved and the method's accuracy and computational speed are analyzed.

Sharma, A.↗

NASA NTRS1980

Validation Methods Research for Fault-Tolerant Avionics and Control Systems Sub-Working Group Meeting. CARE 3 peer review

A computer aided reliability estimation procedure (CARE 3), developed to model the behavior of ultrareliable systems required by flight-critical avionics and control systems, is evaluated. The mathematical models, numerical method, and fault-tolerant architecture modeling requirements are examined, and the testing and characterization procedures are discussed. Recommendations aimed at enhancing CARE 3 are presented; in particular, the need for a better exposition of the method and the user interface is emphasized.

Trivedi, K. S.↗

DOE OSTI2015

On control of continuous dynamical polysystems in discrete times

This paper considers control systems with a finite number of control parameters, i.e. dynamical polysystems. It is assumed that control switchings can occur only at certain discrete times. Statistical methods in number theory are used to investigate the controllability of these systems. Existence of control switching times is established through analysis of Diophantine equations. The values of control switching times are found through suitable numerical methods. Certain assertions on the controllability of the polysystems of the class under consideration are proved. Some examples are considered.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2017

Brownian dynamics of confined suspensions of active microrollers

In this study, we develop efficient numerical methods for performing many-body Brownian dynamics simulations of a recently observed fingering instability in an active suspension of colloidal rollers sedimented above a wall [M. Driscoll, B. Delmotte, M. Youssef, S. Sacanna, A. Donev, and P. Chaikin, Nat. Phys. (2016), preprint arXiv:1609.08673. We present a stochastic Adams-Bashforth integrator for the equations of Brownian dynamics, which has the same cost but is more accurate than the widely used Euler-Maruyama scheme, and use a random finite difference to capture the stochastic drift proportional to the divergence of the configuration-dependent mobility matrix. We generate the Brownian increments using a Krylov method and show that for particles confined to remain in the vicinity of a no-slip wall by gravity or active flows, the number of iterations is independent of the number of particles. Our numerical experiments with active rollers show that the thermal fluctuations set the characteristic height of the colloids above the wall, both in the initial condition and the subsequent evolution dominated by active flows. Finally, the characteristic height in turn controls the time scale and wavelength for the development of the fingering instability

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1986

Semiempirical methods for computing turbulent flows

Two semiempirical theories which provide a basis for determining the turbulent friction and heat exchange near a wall are presented: (1) the Prandtl-Karman theory, and (2) the theory utilizing an equation for the energy of turbulent pulsations. A comparison is made between exact numerical methods and approximate integral methods for computing the turbulent boundary layers in the presence of pressure, blowing, or suction gradients. Using the turbulent flow around a plate as an example, it is shown that, when computing turbulent flows with external turbulence, it is preferable to construct a turbulence model based on the equation for energy of turbulent pulsations.

Belov, I. A.↗

NASA NTRS1990

Numerical simulation of transition in boundary layers

Different implications of employment of the temporal or spatial model in direct numerical simulations of transition based on the complete Navier-Stokes equations are outlined, and it is noted that in situations where approximations are not justified, the temporal model cannot be employed. The discussion is then concentrated on developing numerical methods for solving the Navier-Stokes equations based on the spatial model that are applicable for realistic situations for which the temporal model would not be valid. It is concluded that in the future numerical simulations based on the spatial approach will replace, to a large extent, the temporal simulations.

Fasel, H.↗

DOE OSTI2014

Analytical calculation for the solid angle subtended at an arbitrarily positioned point source by an ellipse - Paper 38

The geometrical solid angle can be used to calibrate radiation detectors or to determine the activity of radioactive sources. In the present work, we introduce a direct analytical method for calculating the solid angle subtended at an arbitrarily positioned point source by a closed contour (ellipse). The validity of the derived analytical expressions was successfully confirmed by the comparison with some published data (Numerical Method). Direct mathematical expressions to calculate the geometrical solid angle subtended by an elliptical disk detector with respect to arbitrarily positioned radiating point sources have been derived. The agreement between the calculated and the published solid angle values is excellent; the discrepancies were (0.0%). This means that the present approach is efficient and sufficiently powerful to evaluate the geometrical efficiency of the elliptical disk detectors. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗

DOE OSTI2018

Evaluation of Neutron-induced Cross Sections and their Related Covariances with Physical Constraints

Nuclear data, along with numerical methods and the associated calculation schemes, continue to play a key role in reactor design, reactor core operating parameters calculations, fuel cycle management and criticality safety calculations. Due to the intensive use of Monte-Carlo calculations reducing numerical biases, the final accuracy of neutronic calculations increasingly depends on the quality of nuclear data used. This paper gives a broad picture of all ingredients treated by nuclear data evaluators during their analyses. After giving an introduction to nuclear data evaluation, we present implications of using the Bayesian inference to obtain evaluated cross sections and related uncertainties. In particular, a focus is made on systematic uncertainties appearing in the analysis of differential measurements as well as advantages and drawbacks one may encounter by analyzing integral experiments. The evaluation work is in general done independently in the resonance and in the continuum energy ranges giving rise to inconsistencies in evaluated files. For future evaluations on the whole energy range, we call attention to two innovative methods used to analyze several nuclear reaction models and impose constraints. Finally, we discuss suggestions for possible improvements in the evaluation process to master the quantification of uncertainties. These are associated with experiments (microscopic and integral), nuclear reaction theories and the Bayesian inference.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Revealing Inconsistencies in X-ray Width Methods for Nanomaterials.

Abstract not provided.

Kunka, Cody↗

DOE OSTI2013

The backward phase flow method for the Eulerian finite time Lyapunov exponent computations

We propose a simple Eulerian approach to compute the moderate to long time flow map for approximating the Lyapunov exponent of a (periodic or aperiodic) dynamical system. The idea is to generalize a recently proposed backward phase flow method which is specially designed for long time level set propagation. Unlike the original phase flow method or the backward phase flow method, which is applicable only to autonomous systems, the current approach can also be applied to any time-dependent (periodic or aperiodic) flow. We will discuss the stability of the proposed method. Numerical examples will be given to demonstrate the effectiveness of the algorithm.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1992

Aerodynamics of loaded cascades in subsonic flows subject to unsteady three-dimensional vortical disturbances

A highly efficient numerical method is developed for three-dimensional, periodic, vortical flows around a cascade of loaded airfoils. The method linearizes Euler's equations about the mean flow of the cascade and thus fully accounts for the effects of distortion of the vortical disturbances as they propogate and interact with the cascade mean flow. The numerical scheme is based on splitting the unsteady velocity into vortical and potential parts. The latter is governed by a non-constant coefficient inhomogeneous convective wave equation. A new and computationally suitable out-flow conditions are derived and avoid the difficulties associated with the singular velocity downstream. Solutions were obtained in the frequency domain by using a body-fitted coordinate system. Results are presented to demonstrate the effects of the out-flow boundary conditions, cascade spacing, mean blade loading and gust upstream conditions on the aerodynamics response and unsteady pressure field of a cascade.

Fang, J.↗

NASA NTRS1993

Non-oscillatory and non-diffusive solution of convection problems by the iteratively reweighted least-squares finite element method

A comparative description is presented for the least-squares FEM (LSFEM) for 2D steady-state pure convection problems. In addition to exhibiting better control of the streamline derivative than the streamline upwinding Petrov-Galerkin method, numerical convergence rates are obtained which show the LSFEM to be virtually optimal. The LSFEM is used as a framework for an iteratively reweighted LSFEM yielding nonoscillatory and nondiffusive solutions for problems with contact discontinuities; this method is shown to convect contact discontinuities without error when using triangular and bilinear elements.

Jiang, Bo-Nan↗

NASA NTRS1999

A Parallel, Finite-Volume Algorithm for Large-Eddy Simulation of Turbulent Flows

A parallel, finite-volume algorithm has been developed for large-eddy simulation (LES) of compressible turbulent flows. This algorithm includes piecewise linear least-square reconstruction, trilinear finite-element interpolation, Roe flux-difference splitting, and second-order MacCormack time marching. Parallel implementation is done using the message-passing programming model. In this paper, the numerical algorithm is described. To validate the numerical method for turbulence simulation, LES of fully developed turbulent flow in a square duct is performed for a Reynolds number of 320 based on the average friction velocity and the hydraulic diameter of the duct. Direct numerical simulation (DNS) results are available for this test case, and the accuracy of this algorithm for turbulence simulations can be ascertained by comparing the LES solutions with the DNS results. The effects of grid resolution, upwind numerical dissipation, and subgrid-scale dissipation on the accuracy of the LES are examined. Comparison with DNS results shows that the standard Roe flux-difference splitting dissipation adversely affects the accuracy of the turbulence simulation. For accurate turbulence simulations, only 3-5 percent of the standard Roe flux-difference splitting dissipation is needed.

Bui, Trong T.↗

NASA NTRS1983

Symposium on Numerical and Physical Aspects of Aerodynamic Flows, 2nd, California State University, Long Beach, CA, January 17-20, 1983, Proceedings

The present conference covers topics concerning the measurement and calculation of interactive flows, together with problems posed by subsonic and transonic wings, missiles, and ships. Discussions are presented on the time-dependent finite difference simulation of unsteady interactive flows, Navier-Stokes equation methods, numerical solutions for spatially periodic boundary layers, the application of unsteady laminar tripple deck theory to viscous-inviscid interaction, the coupling of boundary layer and Euler equation solutions, and viscous-inviscid flow interactions. Also discussed are leading and trailing edge flows, three-dimensional wing flows, small disturbance calculations including entropy corrections, an inviscid computational method for tactical missiles, and boundary layer and flow separation characteristics of bodies of revolution at incidence.

Source record↗

DOE OSTI2020

Closed-Form Solutions for the Equations of Motion of the Heavy Symmetrical Top with One Point Fixed

The Equations of Motion (EOM) for the Heavy Symmetrical Top with One Point Fixed are highly non-linear. The literature describes the numerical methods that are used to resolve this Classical system including modern tools i.e. the Runge-Kutta Fourth Order method. It is more difficult to derivate closed-form solutions for the EOM and as mentioned in the literature it is not always possible to find the close-form solution for all the EOM. Fortunately, there are a few examples available that will serve as a guide to move further on in this topic. It is the purpose of this paper to find a methodology that will produce the solutions for a given subset of EOM that fulfill certain requisites. The report is organized as follows: it starts with a very short summary of the literature available on this topic and quickly follows into the derivation of the EOM using the Euler-Lagrange method. The Routhian will be used to reduce the size of the expression. It continues with the formulation of the Classical cubic function (f(u)) through a novel process. The roots of f(u) are of the outmost importance to be able to find the EOM closed-form solution, and once the final roots are selected the general method that will produce the closed-form solutions is presented. Two sets of examples are included to show the validity of the process and comparisons of the results from the closed-form solutions vs. the numerical results for these examples are shown.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗