Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical differentiation”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
NASA NTRS1980

The Finite Analytic Method for steady and unsteady heat transfer problems

A new numerical method called the Finite Analytical Method for solving partial differential equations is introduced. The basic idea of the finite analytic method is the incorporation of the local analytic solution in obtaining the numerical solution of the problem. The finite analytical method first divides the total region of the problem into small subregions in which local analytic solutions are obtained. Then an algebraic equation is derived from the local analytic solution for each subregion relating an interior nodal value at a point P in the subregion to its neighboring nodal values. The assembly of all the local analytic solutions thus provides the finite-analytic numerical solution of the problem. In this paper the finite analytic method is illustrated in solving steady and unsteady heat transfer problems.

Chen, C.-J.↗
DOE OSTI2025

A compartmentalized model of multiphase chemical kinetics

There are significant challenges in predicting multiphase chemical kinetics due to the complex coupling of reaction and mass transport across a phase boundary (i.e., interface). Here, we describe a framework for predicting multiphase kinetics that embeds the elementary kinetic steps of reaction, solvation, and diffusion into a coarse grain spatial description of two phases. The model is constructed to bridge the short-timescale interfacial dynamics observed in molecular simulations with the longer timescales observed in kinetic experiments. A simple set of governing differential equations is derived, which, when solved numerically or analytically, yield accurate predictions of multiphase kinetics in microdroplets. Although the equations are formulated for gas-liquid reactions, the underlying conceptual framework is general and can be applied to transformations in other two-phase systems (solid-liquid, liquid-liquid, etc.).

Chemical kinetics and dynamics↗
DOE OSTI2024

Microscopic optical potentials for medium-mass isotopes derived at the first order of Watson multiple-scattering theory

We perform a first-principles calculation of optical potentials for nucleon elastic scattering off medium-mass isotopes. Fully based on a saturating chiral Hamiltonian, the optical potentials are derived by folding nuclear density distributions computed with self-consistent Green's function theory with a nucleon-nucleon t matrix computed with a consistent chiral interaction. The dependence on the folding interaction as well as the convergence of the target densities are investigated. Numerical results are presented and discussed for differential cross sections and analyzing powers, with focus on elastic proton scattering off calcium and nickel isotopes. Our optical potentials generally show a remarkable agreement with the available experimental data for laboratory energies in the range 65–200 MeV. We study the evolution of the scattering observables with increasing proton-neutron asymmetry by computing theoretical predictions of the cross section and analyzing power over the calcium and nickel isotopic chains. Published by the American Physical Society 2024

Physics↗
DOE OSTI2022

Adaptive Space-Time Methods for Large Scale Optimal Design

When modeling complex physical systems with advanced dynamics, such as shocks and singularities, many classic methods for solving partial differential equations can return inaccurate or unusable results. One way to resolve these complex dynamics is through r-adaptive refinement methods, in which a fixed number of mesh points are shifted to areas of high interest. The mesh refinement map can be found through the solution of the Monge-Ampére equation, a highly nonlinear partial differential equation. Due to its nonlinearity, the numerical solution of the Monge-Ampére equation is nontrivial and has previously required computationally expensive methods. In this report, we detail our novel optimization-based, multigrid-enabled solver for a low-order finite element approximation of the Monge-Ampére equation. This fast and scalable solver makes r-adaptive meshing more readily available for problems related to large-scale optimal design. Beyond mesh adaptivity, our report discusses additional applications where our fast solver for the Monge-Ampére equation could be easily applied.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1972

Assessment of a transitional boundary layer theory at low hypersonic Mach numbers

An investigation was carried out to assess the accuracy of a transitional boundary layer theory in the low hypersonic Mach number regime. The theory is based upon the simultaneous numerical solution of the boundary layer partial differential equations for the mean motion and an integral form of the turbulence kinetic energy equation which controls the magnitude and development of the Reynolds stress. Comparisions with experimental data show the theory is capable of accurately predicting heat transfer and velocity profiles through the transitional regime and correctly predicts the effects of Mach number and wall cooling on transition Reynolds number. The procedure shows promise of predicting the initiation of transition for given free stream disturbance levels. The effects on transition predictions of the pressure dilitation term and of direct absorption of acoustic energy by the boundary layer were evaluated.

Shamroth, S. J.↗
NASA NTRS1973

Computer simulation of space station computer steered high gain antenna

The mathematical modeling and programming of a complete simulation program for a space station computer-steered high gain antenna are described. The program provides for reading input data cards, numerically integrating up to 50 first order differential equations, and monitoring up to 48 variables on printed output and on plots. The program system consists of a high gain antenna, an antenna gimbal control system, an on board computer, and the environment in which all are to operate.

Beach, S. W.↗
NASA NTRS1974

Inelastic, nonlinear analysis of stiffened shells of revolution by numerical integration

This paper describes the latest addition to the STARS system of computer programs, STARS-2P, for the plastic, large deflection analysis of axisymmetrically loaded shells of revolution. The STARS system uses a numerical integration scheme to solve the governing differential equations. Several unique features for shell of revolution programs that are included in the STARS-2P program are described. These include orthotropic nonlinear kinematic hardening theory, a variety of shell wall cross sections and discrete ring stiffeners, cyclic and nonproportional mechanical and thermal loading capability, the coupled axisymmetric large deflection elasto-plastic torsion problem, an extensive restart option, arbitrary branching capability, and the provision for the inelastic treatment of smeared stiffeners, isogrid, and waffle wall constructions. To affirm the validity of the results, comparisons with available theoretical and experimental data are presented.

Levine, H. S.↗
NASA NTRS1976

Stabilization by modification of the Lagrangian

In order to reduce the error growth during a numerical integration, a method for stabilizing the differential equations of Keplerian motion is offered. It is characterized by the use of the eccentric anomaly as independent variable in such a way that the time transformation is given by a generalized Lagrange formalism. The control terms in the equations of motion obtained by this modified Lagrangian give immediately a completely Lyapunov-stable set of differential equations. The equation of time integration is modified by a control term which leads to an integral which defined the time element for the perturbed Keplerian motion.

Baumgarte, J. W.↗
NASA NTRS1981

Multigrid Methods

Developments in numerical solution of certain types of partial differential equations by rapidly converging sequences of operations on supporting grids that range from very fine to very coarse are presented.

Source record↗
NASA NTRS1982

The optimization of convergence for Chebyshev polynomial methods in an unbounded domain

Grosch and Orszag (1977) have performed a numerical analysis of the problem of solving differential equations in a semiinfinite or infinite domain using Chebyshev polynomials. The principal limitation of the conducted study was that it was entirely empirical. Various differential equations were solved in different ways and the numbers were compared. The present investigation has the objective to extend the studies conducted by Grosch and Orszag by deriving asymptotic approximations to the Chebyshev coefficients of simple model functions. This approach makes it possible to conduct more systematic comparisons of different methods, extend the range of comparisons, and, perhaps most important, give simple analytic formulas for choosing the optimum domain size or mapping parameter L for various situations.

Boyd, J. P.↗
NASA NTRS1982

Symmetric baroclinic instability of a Hadley cell

A symmetric baroclinic instability is examined in terms of a Boussinesq fluid contained between two horizontal plates to determine the effects of the Ekman and thermal layers. Governing equations are written for a rotating reference frame, taking into account the Rossby, Ekman, and Prandtl numbers. Equations are defined for the perturbation functions, treated as an eigenvalue problem, and a numerical integration of the full eighth order differential system is performed by a shooting technique. An instability is found to occur in the Hadley cell containing both Ekman and thermal boundary layers when the Richardson number is close to unity. If the Prandtl number is fixed the critical Richardson number decreases with an increasing Ekman number until the Ekman number reaches a certain value, at which time the fluid is stable.

Antar, B. N.↗
NASA NTRS1983

Computation of shock wave/target interaction

Computational results of shock waves impinging on targets and the ensuing diffraction flowfield are presented. A number of two-dimensional cases are computed with finite difference techniques. The classical case of a shock wave/cylinder interaction is compared with shock tube data and shows the quality of the computations on a pressure-time plot. Similar results are obtained for a shock wave/rectangular body interaction. Here resolution becomes important and the use of grid clustering techniques tend to show good agreement with experimental data. Computational results are also compared with pressure data resulting from shock impingement experiments for a complicated truck-like geometry. Here of significance are the grid generation and clustering techniques used. For these very complicated bodies, grids are generated by numerically solving a set of elliptic partial differential equations.

Mark, A.↗
NASA NTRS1982

Elliptic grid generation

Various types of generating systems for boundary-conforming coordinate systems based on the numerical solution of systems of elliptic partial differential equations are discussed. Particular emphasis is given to the determination of functions in these equations which control the distribution of the curvilinear coordinate lines in the field.

Thompson, J. F.↗
NASA NTRS1988

The identification of a distributed parameter model for a flexible structure

A computational method is developed for the estimation of parameters in a distributed model for a flexible structure. The structure we consider (part of the RPL experiment) consists of a cantilevered beam with a thruster and linear accelerometer at the free end. The thruster is fed by a pressurized hose whose horizontal motion effects the transverse vibration of the beam. The Euler-Bernoulli theory is used to model the vibration of the beam and treat the hose thruster assembly as a lumped or point mass dashpot spring system at the tip. Measurements of linear acceleration at the tip are used to estimate the hose parameters (mass, stiffness, damping) and a Voigt-Kelvin viscoelastic structural damping parameter for the beam using a least squares fit to the data. Spline based approximations are considered to the hybrid (coupled ordinary and partial differential equations) systems; theoretical convergence results and numerical studies with both simulation and actual experimental data obtained from the structure are presented and discussed.

Banks, H. T.↗
NASA NTRS1989

Nonlinear evolution of interacting oblique waves on two-dimensional shear layers

The effects of critical layer nonlinearity are considered on spatially growing oblique instability waves on nominally two-dimensional shear layers between parallel streams. The analysis shows that three-dimensional effects cause nonlinearity to occur at much smaller amplitudes than it does in two-dimensional flows. The nonlinear instability wave amplitude is determined by an integro-differential equation with cubic type nonlinearity. The numerical solutions to this equation are worked out and discussed in some detail. The numerical solutions always end in a singularity at a finite downstream distance.

Goldstein, M. E.↗
NASA NTRS1989

Cirrus crystal nucleation by homogeneous freezing of solution droplets

A numerical model consisting of a system of differential equations is used to study cirrus crystal nucleation in a rising parcel containing a distribution of cloud condensation nuclei. The evolution of the particle population and the thermodynamic variables in the parcel are examined. The results suggest that, if homogeneous freezing is not considered, liquid water should be detected below -40 C. If homogeneous freezing is considered, the rapid growth of ice crystals and vapor depletion prevent water saturation from being reached. It is shown that the likelihood of a droplet being frozen is increased by lower temperatures, larger droplet diameter, or lower solution density.

Heymsfield, Andrew J.↗
NASA NTRS1989

Nonlinear evolution of interacting oblique waves on two-dimensional shear layers

The effects of critical layer nonlinearity are considered on spatially growing oblique instability waves on nominally two-dimensional shear layers between parallel streams. The analysis shows that three-dimensional effects cause nonlinearity to occur at much smaller amplitudes than it does in two-dimensional flows. The nonlinear instability wave amplitude is determined by an integro-differential equation with cubic type nonlinearity. The numerical solutions to this equation are worked out and discussed in some detail. The numerical solutions always end in a singularity at a finite downstream distance.

Goldstein, M. E.↗