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 1,459 records · Page 81
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Comparative study of computational efficiency of two LU schemes for non-equilibrium reacting flows

The systems of equations governing chemically nonequilibrium flows are solved numerically by two lower-upper(LU) algorithms, namely, LU-SSOR (Yoon and Jameson, 1987) and LU-SW (Steger and Warming, 1981). The methods of flux-Jacobian splitting differentiate the two LU algorithms. The LU-SSOR gives a simpler implicit operator with less temporal damping, while the LU-SW converges faster but requires more computation for each iteration. Relative efficiencies and mutual merits of these two LU algorithms are compared. The results show that LU-SW with upwind-differenced right-hand-side gives the best convergence.

Tsai, Y.-L. Peter↗
NASA NTRS1991

Dust particles from comets and asteroids collected at the Earth's orbit: Parent-daughter relationships

The relative contributions of comets and asteroids to the reservoir of dust in the interplanetary medium is not well known. There are direct observations of dust released from comets and there is evidence to associate the IRAS dust bands with possible collisions of Asteroids in the main belt. It is believed that one may combine lab analysis of the physics and chemistry of captured particles with orbital data in order to identify comet and asteroid parent bodies. It is possible to use the collected orbits of the dust to connect with its source in two ways. One is to consider the long time orbit evolution of the dust under Poynting-Robertson drag. The other is to look at the prompt orbit change of dust from comets onto trajectories that intersect the earth's orbit. In order to characterize the orbits of dust particles evolved over a long period of time, a study of its orbital evolution was undertaken. Various parameters associated with these dust orbits as they cross the Earth's orbit were considered in order to see if one may discriminate between particles evolved from comets and asteroids. The method was to calculate by a numerical procedure the orbits of dust particles after they left their parent bodies. It appears that as the particles pass the Earth's orbit, asteroidal grains and cometary grains can be differentiated on the basis of their measured orbital eccentricities even after much planetary perturbation. Broad parent daughter associations can be made on this basis from measurement of their trajectories intercepted in earth orbit.

Jackson, A. A.↗
NASA NTRS1990

Development of a model of space station solar array

Space structures, such as the space station solar arrays, must be extremely lightweight, flexible structures. Accurate prediction of the natural frequencies and mode shapes is essential for determining the structural adequacy of components, and designing a control system. The tension preload in the blanket of photovoltaic solar collectors, and the free/free boundary conditions of a structure in space, causes serious reservations on the use of standard finite element techniques of solution. In particular, a phenomena known as grounding, or false stiffening, of the stiffness matrix occurs during rigid body rotation. The grounding phenomena is examined in detail. Numerous stiffness matrices developed by others are examined for rigid body rotation capability, and found lacking. Various techniques are used for developing new stiffness matrices from the rigorous solutions of the differential equations, including the solution of the directed force problem. A new directed force stiffness matrix developed by the author provides all the rigid body capabilities for the beam in space.

Bosela, Paul A.↗
NASA NTRS1991

The analysis of control trajectories using symbolic and database computing

The research broadly concerned the symbolic computation, mixed numeric-symbolic computation, and data base computation of trajectories of dynamical systems, especially control systems. It was determined that trees can be used to compute symbolically series which approximate solutions to differential equations.

Grossman, Robert↗
NASA NTRS1993

Numerical Studies of Collisionless Current Layers

The purpose of this proposal was to investigate collisionless current layers using a variety of analytic and numerical tools. The first year of the contract was dedicated to analytical studies, to the porting and adaption of codes being used in this study, and to the numerical simulation of collisionless current layers. The second year entailed the development of multi-dimensional hybrid algorithms as well as the re-examination of the problem of integro-differential equations that occur in the linear stage of plasma instabilities.

Quest, Kevin B.↗
NASA NTRS2003

Fully Nonlinear Modeling and Analysis of Precision Membranes

High precision membranes are used in many current space applications. This paper presents a fully nonlinear membrane theory with forward and inverse analyses of high precision membrane structures. The fully nonlinear membrane theory is derived from Jaumann strains and stresses, exact coordinate transformations, the concept of local relative displacements, and orthogonal virtual rotations. In this theory, energy and Newtonian formulations are fully correlated, and every structural term can be interpreted in terms of vectors. Fully nonlinear ordinary differential equations (ODES) governing the large static deformations of known axisymmetric membranes under known axisymmetric loading (i.e., forward problems) are presented as first-order ODES, and a method for obtaining numerically exact solutions using the multiple shooting procedure is shown. A method for obtaining the undeformed geometry of any axisymmetric membrane with a known inflated geometry and a known internal pressure (i.e., inverse problems) is also derived. Numerical results from forward analysis are verified using results in the literature, and results from inverse analysis are verified using known exact solutions and solutions from the forward analysis. Results show that the membrane theory and the proposed numerical methods for solving nonlinear forward and inverse membrane problems are accurate.

Pai, P. Frank↗
NASA NTRS1949

Two-Dimensional Potential Flows

Contents include the following: Characteristic differential equations - initial and boundary conditions. Integration of the second characteristic differential equations. Direct application of Meyer's characteristic hodograph table for construction of two-dimensional potential flows. Prandtl-Busemann method. Development of the pressure variation for small deflection angles. Numerical table: relation between deflection, pressure, velocity, mach number and mach angle for isentropic changes of state according to Prandtl-Meyer for air (k = 1.405). References.

Schaefer, Manfred↗
NASA NTRS2006

Thermal Model of a Current-Carrying Wire in a Vacuum

A computer program implements a thermal model of an insulated wire carrying electric current and surrounded by a vacuum. The model includes the effects of Joule heating, conduction of heat along the wire, and radiation of heat from the outer surface of the insulation on the wire. The model takes account of the temperature dependences of the thermal and electrical properties of the wire, the emissivity of the insulation, and the possibility that not only can temperature vary along the wire but, in addition, the ends of the wire can be thermally grounded at different temperatures. The resulting second-order differential equation for the steady-state temperature as a function of position along the wire is highly nonlinear. The wire is discretized along its length, and the equation is solved numerically by use of an iterative algorithm that utilizes a multidimensional version of the Newton-Raphson method.

Border, James↗
NASA NTRS2012

A Posteriori Study of a DNS Database Describing Super critical Binary-Species Mixing

Currently, the modeling of supercritical-pressure flows through Large Eddy Simulation (LES) uses models derived for atmospheric-pressure flows. Those atmospheric-pressure flows do not exhibit the particularities of high densitygradient magnitude features observed both in experiments and simulations of supercritical-pressure flows in the case of two species mixing. To assess whether the current LES modeling is appropriate and if found not appropriate to propose higher-fidelity models, a LES a posteriori study has been conducted for a mixing layer that initially contains different species in the lower and upper streams, and where the initial pressure is larger than the critical pressure of either species. An initially-imposed vorticity perturbation promotes roll-up and a double pairing of four initial span-wise vortices into an ultimate vortex that reaches a transitional state. The LES equations consist of the differential conservation equations coupled with a real-gas equation of state, and the equation set uses transport properties depending on the thermodynamic variables. Unlike all LES models to date, the differential equations contain, additional to the subgrid scale (SGS) fluxes, a new SGS term that is a pressure correction in the momentum equation. This additional term results from filtering of Direct Numerical Simulation (DNS) equations, and represents the gradient of the difference between the filtered pressure and the pressure computed from the filtered flow field. A previous a priori analysis, using a DNS database for the same configuration, found this term to be of leading order in the momentum equation, a fact traced to the existence of high-densitygradient magnitude regions that populated the entire flow; in the study, models were proposed for the SGS fluxes as well as this new term. In the present study, the previously proposed constantcoefficient SGS-flux models of the a priori investigation are tested a posteriori in LES, devoid of or including, the SGS pressure correction term. The present pressure-correction model is different from, and more accurate as well as less computationally intensive than that of the a priori study. The constant-coefficient SGS-flux models encompass the Smagorinsky (SMC), in conjunction with the Yoshizawa (YO) model for the trace, the Gradient (GRC) and the Scale Similarity (SSC) models, all exercised with the a priori study constant coefficients calibrated at the transitional state. The LES comparison is performed with the filtered- and-coarsened (FC) DNS, which represents an ideal LES solution. Expectably, an LES model devoid of SGS terms is shown to be considerably inferior to models containing SGS effects. Among models containing SGS effects, those including the pressure-correction term are substantially superior to those devoid of it. The sensitivity of the predictions to the initial conditions and grid size are also investigated. Thus, it has been discovered that, additional to the atmospheric-pressure models currently used, a new model is necessary to simulate supercritical-pressure flows. This model depends on the thermodynamic characteristics of the chemical species involved.

Bellan, Josette↗
NASA NTRS2020

Evaluation of the Stratospheric and Tropospheric Bromine Burden over Fairbanks, Alaska Based on Column Retrievals of Bromine Monoxide

In spring 2011, columns of bromine monoxide (BrO) were retrieved over Fairbanks, Alaska using a ground-based multifunction differential optical absorption spectroscopy (MFDOAS) instrument. MFDOAS vertical column BrO is consistently lower than retrievals from the satellite-based Ozone Monitoring Instrument (OMI), with a relative bias of 20 ± 14%. Numerous tropical-based studies suggest that 5 ± 2 ppt of bromine from very short-lived substances (VSLS) reaches the stratosphere. We evaluate upper limits on the contribution of VSLS to stratospheric bromine by treating the column retrievals of BrO as purely stratospheric and modeling the ratio of BrO to total inorganic bromine. The OMI and MFDOAS retrievals respectively present 8 and 5 ppt upper limits on the stratospheric injection of VSLS, and kinetic uncertainties in the daytime partitioning of bromine species decrease both values by ~1.6 ppt. The OMI-based estimate is in agreement with the 5 ppt tropical-based value for stratospheric injection of VSLS if the tropospheric column of BrO is 1.5 × 10(exp 13) molecules/cu. cm over Fairbanks, which is within the range of uncertainty of a second ground-based instrument that monitored tropospheric BrO during the campaign. Because our ground-based instruments detected no BrO near the surface, this value for tropospheric BrO would originate from the free troposphere and is in agreement with previous retrievals of background tropospheric BrO. Our calculations of tropospheric BrO over Fairbanks are most sensitive to uncertainties in the stratospheric loading of VSLS, followed by the difference between the OMI and MFDOAS retrievals of BrO.

Pamela A. Wales↗
NASA NTRS2022

Continuum Covariance Propagation for Understanding Variance Loss in Advective Systems

Motivated by the spurious variance loss encountered during covariance propagation in atmospheric and other large-scale data assimilation systems, we consider the problem for state dynamics governed by the continuity and related hyperbolic partial differential equations. This loss of variance has been attributed to reduced-rank representations of the covariance matrix, as in ensemble methods for example, or else to the use of dissipative numerical methods. Through a combination of analytical work and numerical experiments, we demonstrate that significant variance loss, as well as gain, typically occurs during covariance propagation, even at full rank. The cause of this unusual behavior is a discontinuous change in the continuum covariance dynamics as correlation lengths become small, for instance in the vicinity of sharp gradients in the velocity field. This discontinuity in the covariance dynamics arises from hyperbolicity: the diagonal of the kernel of the covariance operator is a characteristic surface for advective dynamics. Our numerical experiments demonstrate that standard numerical methods for evolving the state are not adequate for propagating the covariance, because16they do not capture the discontinuity in the continuum covariance dynamics as correlations lengths tend to zero. Our analytical and numerical results show that this leads to significant, spurious variance loss in certain regions, and gain in others. The results suggest that developing local covariance propagation methods designed specifically to capture covariance evolution near the diagonal may prove a useful alternative to current methods of covariance propagation.

covariance propagation↗
NASA NTRS1985

Numerical simulation of convective response to mesoscale forcing

A two-dimensional mesoscale model was used to initialize a two-dimensional cloud model with both mesoscale thermodynamic and dynamic information utilized in the initialization. The Mesoscale-cloud mdoel linkage was used to examine differential convective response due to mesoscale variations for the April 24, 1982 Atmospheric Variability Experiment-Vertical Atmospheric Sounder (AVE-VAS), Case 4. On this day both data analyses and mesoscale model simulations indicated strong variations in thermodynamic and dynamic structure across the panhandle of Texas where a moderately strong convective line formed. Using the cloud model, the current research has shown a preferred area for strong convection to occur due to concomitant mesoscale convergence and mesoscale destabilization of the atmosphere.

Mcnider, R. T.↗
DOE OSTI2025

Two-Scale Neural Networks for Partial Differential Equations with Small Parameters

We propose a two-scale neural network method for solving partial differential equations (PDEs) with small parameters using physics-informed neural networks (PINNs). We directly incorporate the small parameters into the architecture of neural networks. The proposed method enables solving PDEs with small parameters in a simple fashion, without adding Fourier features or other computationally taxing searches of truncation parameters. Various numerical examples demonstrate reasonable accuracy in capturing features of large derivatives in the solutions caused by small parameters.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1974

The solution of three-variable duct-flow equations

This paper establishes a numerical method for the solution of three-variable problems and is applied here to rotational flows through ducts of various cross sections. An iterative scheme is developed, the main feature of which is the addition of a duplicate variable to the forward component of velocity. Two forward components of velocity result from integrating two sets of first order ordinary differential equations for the streamline curvatures, in intersecting directions across the duct. Two pseudo-continuity equations are introduced with source/sink terms, whose strengths are dependent on the difference between the forward components of velocity. When convergence is obtained, the two forward components of velocity are identical, the source/sink terms are zero, and the original equations are satisfied. A computer program solves the exact equations and boundary conditions numerically. The method is economical and compares successfully with experiments on bent ducts of circular and rectangular cross section where secondary flows are caused by gradients of total pressure upstream.

Stuart, A. R.↗
DOE OSTI2020

Assessing three closed-loop learning algorithms by searching for high-quality quantum control pulses

Designing a high-quality control is crucial for reliable quantum computation. Among the existing approaches, closed-loop leaning control is an effective choice. Its efficiency depends on the learning algorithm employed, thus deserving algorithmic comparisons for its practical applications. Here we assess three representative learning algorithms, including GRadient Ascent Pulse Engineering (GRAPE), improved Nelder-Mead (NMplus), and Differential Evolution (DE), by searching for high-quality control pulses to prepare the Bell state. We first implement each algorithm experimentally in a nuclear magnetic resonance system and then conduct a numerical study considering the impact of some possible significant experimental uncertainties. The experiments report the successful preparation of the high-fidelity target state by the three algorithms, while NMplus converges fastest, and these results coincide with the numerical simulations when potential uncertainties are negligible. However, under certain significant uncertainties, these algorithms possess distinct performance with respect to their resulting precision and efficiency, and DE shows the best robustness. Finally, this study provides insight to aid in the practical application of different closed-loop learning algorithms in realistic physical scenarios.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1971

Acquisition time for sideband modulated phase lock loops.

In this paper a definition for acquisition of a sideband modulated carrier is proposed for the phase lock loop, and the acquisition differential equation solved for a region of the modulation parameter space to determine the expected value of acquisition time. Principally the case of modulation by a single unmodulated subcarrier is considered. The numerical analysis is made for the noiseless case.

Gilbert, A. L.↗
NASA NTRS1973

Survey of computational methods for three-dimensional supersonic inviscid flows with shocks

Hyperbolic systems of partial differential equations governing supersonic inviscid flows are discussed and analyzed. Finite-difference analogues for integrating these systems in the interior of fluid domains are described from two points of view: a differential form approach and an integral form approach. The algorithms presented are analyzed for stability and accuracy. The concept of time splitting is discussed and supplied to these methods to achieve increased numerical efficiency. Techniques for treating conditions at the boundaries of the fluid domain and shock-wave discontinuities at surfaces within the domain are described.

Maccormack, R. W.↗