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,423 records · Page 79
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Dissipation measures in weakly collisional plasmas

ABSTRACT The physical foundations of the dissipation of energy and the associated heating in weakly collisional plasmas are poorly understood. Here, we compare and contrast several measures that have been used to characterize energy dissipation and kinetic-scale conversion in plasmas by means of a suite of kinetic numerical simulations describing both magnetic reconnection and decaying plasma turbulence. We adopt three different numerical codes that can also include interparticle collisions: the fully kinetic particle-in-cell vpic, the fully kinetic continuum Gkeyll, and the Eulerian Hybrid Vlasov–Maxwell (HVM) code. We differentiate between (i) four energy-based parameters, whose definition is related to energy transfer in a fluid description of a plasma, and (ii) four distribution function-based parameters, requiring knowledge of the particle velocity distribution function. There is an overall agreement between the dissipation measures obtained in the PIC and continuum reconnection simulations, with slight differences due to the presence/absence of secondary islands in the two simulations. There are also many qualitative similarities between the signatures in the reconnection simulations and the self-consistent current sheets that form in turbulence, although the latter exhibits significant variations compared to the reconnection results. All the parameters confirm that dissipation occurs close to regions of intense magnetic stresses, thus exhibiting local correlation. The distribution function-based measures show a broader width compared to energy-based proxies, suggesting that energy transfer is co-localized at coherent structures, but can affect the particle distribution function in wider regions. The effect of interparticle collisions on these parameters is finally discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1972

Transient, compressible heat and mass transfer in porous media using the strongly implicit iteration procedure.

Coupled nonlinear partial differential equations describing heat and mass transfer in a porous matrix are solved in finite difference form with the aid of a new iterative technique (the strongly implicit procedure). Example numerical results demonstrate the characteristics of heat and mass transport in a porous matrix such as a charring ablator. It is emphasized that multidimensional flow must be considered when predicting the thermal response of a porous material subjected to nonuniform boundary conditions.

Curry, D. M.↗
NASA NTRS1983

A noncanonical analytic solution to the J2 perturbed two-body problem

The motion of a satellite subject to an inverse-square gravitational force of attraction and a perturbation due to the earth's oblateness as the J2 term is analyzed, and a uniform, analytic solution correct to first-order in J2, is obtained using a noncanonical approach. The basis for the solution is the transformation and uncoupling of the differential equations for the model. The resulting solution is expressed in terms of elementary functions of the independent variable (the 'true anomaly'), and is of a compact and simple form. Numerical results are comparable to existing solutions.

Jezewski, D. J.↗
NASA NTRS1985

Mixed models and reduction method for dynamic analysis of anisotropic shells

A time-domain computational procedure is presented for predicting the dynamic response of laminated anisotropic shells. The two key elements of the procedure are: (1) use of mixed finite element models having independent interpolation (shape) functions for stress resultants and generalized displacements for the spatial discretization of the shell, with the stress resultants allowed to be discontinuous at interelement boundaries; and (2) use of a dynamic reduction method, with the global approximation vectors consisting of the static solution and an orthogonal set of Lanczos vectors. The dynamic reduction is accomplished by means of successive application of the finite element method and the classical Rayleigh-Ritz technique. The finite element method is first used to generate the global approximation vectors. Then the Rayleigh-Ritz technique is used to generate a reduced system of ordinary differential equations in the amplitudes of these modes. The temporal integration of the reduced differential equations is performed by using an explicit half-station central difference scheme (Leap-frog method). The effectiveness of the proposed procedure is demonstrated by means of a numerical example and its advantages over reduction methods used with the displacement formulation are discussed.

Noor, A. K.↗
NASA NTRS1987

An Exponential Finite Difference Technique for Solving Partial Differential Equations

An exponential finite difference algorithm, as first presented by Bhattacharya for one-dimensianal steady-state, heat conduction in Cartesian coordinates, has been extended. The finite difference algorithm developed was used to solve the diffusion equation in one-dimensional cylindrical coordinates and applied to two- and three-dimensional problems in Cartesian coordinates. The method was also used to solve nonlinear partial differential equations in one (Burger's equation) and two (Boundary Layer equations) dimensional Cartesian coordinates. Predicted results were compared to exact solutions where available, or to results obtained by other numerical methods. It was found that the exponential finite difference method produced results that were more accurate than those obtained by other numerical methods, especially during the initial transient portion of the solution. Other applications made using the exponential finite difference technique included unsteady one-dimensional heat transfer with temperature varying thermal conductivity and the development of the temperature field in a laminar Couette flow.

Handschuh, Robert F.↗
NASA NTRS1987

Degree of cure and viscosity of Hercules HBRF-55 resin

The rate of cure and viscosity were measured for Hercules HBRF-55 resin. The rate of cure was measured by differential scanning calorimetry, while the viscosity was measured by a parallel disk and plate type apparatus. The data were fitted to analytic expressions to make them suitable for use in numerical calculations.

Bhi, S. T.↗
NASA NTRS1989

Asymptotic structure and similarity solutions for three-dimensional turbulent boundary layers

The asymptotic structure of the three-dimensional turbulent boundary layer is investigated in the limit of large Reynolds numbers. A self-consistent, but relatively complex, two-layer structure exists and the simplest situation, corresponding to a plane of symmetry, is considered in this paper as a first step. The adjustment of the streamwise velocity to relative rest, through an outer defect layer and then an inner wall layer, is similar to that in two-dimensional flow. The adjustment of the cross-streamwise velocity is more complicated and it is shown that two terms in the expansion are required to obtain useful results, and in particular to obtain the velocity skew angle at the wall near the symmetry plane. The conditions under which self-similarity is achieved near a plane of symmetry are investigated. A set of ordinary differential equations is developed which describe the streamwise and cross-streamwise velocities near a plane of symmetry in a self-similar flow through two orders of magnitude. Calculated numerical solutions of these equations yield trends which are consistent with experimental observations.

Degani, A. T.↗
NASA NTRS1990

Radiative configuration factors from cylinders to coaxial axisymmetric bodies

Exact solutions are obtained for the radiative configuration factor between differential elements of arbitrary orientation and cylinders. A general method is then proposed for calculating the view factor from a cylinder to a coaxial axisymmetric body using only a single numerical integration. The method is illustrated for axisymmetric bodies with function generators described by a power law equation. The method may be useful in calculating radiative heat transfer between cylindrical bodies and high-density exhaust gases and between annular radiative fins and their bases.

Saltiel, C.↗
NASA NTRS1990

SIAM Conference on Parallel Processing for Scientific Computing, 4th, Chicago, IL, Dec. 11-13, 1989, Proceedings

Attention is given to such topics as an evaluation of block algorithm variants in LAPACK and presents a large-grain parallel sparse system solver, a multiprocessor method for the solution of the generalized Eigenvalue problem on an interval, and a parallel QR algorithm for iterative subspace methods on the CM2. A discussion of numerical methods includes the topics of asynchronous numerical solutions of PDEs on parallel computers, parallel homotopy curve tracking on a hypercube, and solving Navier-Stokes equations on the Cedar Multi-Cluster system. A section on differential equations includes a discussion of a six-color procedure for the parallel solution of elliptic systems using the finite quadtree structure, data parallel algorithms for the finite element method, and domain decomposition methods in aerodynamics. Topics dealing with massively parallel computing include hypercube vs. 2-dimensional meshes and massively parallel computation of conservation laws. Performance and tools are also discussed.

Dongarra, Jack↗
NASA NTRS1995

Algorithm For Computation Of Chemically Reacting Flow

Efficient algorithm developed for use in solving differential equations of transient, chemically reacting flows at all speeds from zero to high mach numbers. Mathematical models represent coupling of thermal, chemical, and dynamical effects. Original intended application is numerical simulation of flows in rocket engines; also applicable to other complex flows affected by finite-rate chemistry - for example, flows in turbines and in internal-combustion engines.

Chen, Yen-Sen↗
NASA NTRS1995

Computational Study of Flow Establishment in a Ram Accelerator

The temporal evolution of the combustion process established during projectile transition from the launch tube into the ram accelerator section containing an explosive hydrogen-oxygen-argon gas mixture is studied. The Navier-Stokes equations for chemically reacting flow are solved in a fully coupled manner, using an implicit, time accurate algorithm. The solution procedure is based on a spatially second order total variation diminishing scheme and a temporally second order, variable-step, backward differentiation formula method. The hydrogen-oxygen chemistry is modeled with a 9-species, 19-step mechanism. The accuracy of the solution method is first demonstrated by several benchmark calculations. Numerical simulations of two ram accelerator configurations are then presented. In particular, the temporal developments of shock-induced combustion and thrust forces are followed. Positive thrust is established in both cases; however, in one of the ram accelerator configurations studied, combustion in the boundary layer enhances its separation, ultimately resulting in unstart.

Yungster, S.↗
NASA NTRS1995

Computational Study of Flow Establishment in Hypersonic Pulse Facilities

This paper presents a study of the temporal evolution of the combustion flowfield established by the interaction of ram-accelerator-type projectiles with an explosive gas mixture accelerated to hypersonic speeds in an expansion tube. The Navier-Stokes equations for a chemically reacting gas are solved in a fully coupled manner using an implicit, time accurate algorithm. The solution procedure is based on a spatially second order, total variation diminishing (TVD) scheme and a temporally second order, variable-step, backward differentiation formula method. The hydrogen-oxygen chemistry is modeled with a 9-species, 19-step mechanism. The accuracy of the solution method is first demonstrated by several benchmark calculations. Numerical simulations of expansion tube flowfields are then presented for two different configurations. In particular, the development of the shock-induced combustion process is followed. In one case, designed to ensure ignition only in the boundary layer, the lateral extent of the combustion front during the initial transient phase was surprisingly large. The time histories of the calculated thrust and drag forces on the ram accelerator projectile are also presented.

Yungster, S.↗
NASA NTRS1996

Regular and Chaotic Quantum Dynamics of Two-Level Atoms in a Selfconsistent Radiation Field

Dynamics of two-level atoms interacting with their own radiation field in a single-mode high-quality resonator is considered. The dynamical system consists of two second-order differential equations, one for the atomic SU(2) dynamical-group parameter and another for the field strength. With the help of the maximal Lyapunov exponent for this set, we numerically investigate transitions from regularity to deterministic quantum chaos in such a simple model. Increasing the collective coupling constant b is identical with 8(pi)N(sub 0)(d(exp 2))/hw, we observed for initially unexcited atoms a usual sharp transition to chaos at b(sub c) approx. equal to 1. If we take the dimensionless individual Rabi frequency a = Omega/2w as a control parameter, then a sequence of order-to-chaos transitions has been observed starting with the critical value a(sub c) approx. equal to 0.25 at the same initial conditions.

Konkov, L. E.↗
NASA NTRS1996

A Multilevel Algorithm for the Solution of Second Order Elliptic Differential Equations on Sparse Grids

A multilevel algorithm is presented that solves general second order elliptic partial differential equations on adaptive sparse grids. The multilevel algorithm consists of several V-cycles. Suitable discretizations provide that the discrete equation system can be solved in an efficient way. Numerical experiments show a convergence rate of order Omicron(1) for the multilevel algorithm.

Pflaum, Christoph↗
NASA NTRS2016

Analysis of Summer-Time Ozone and Precursor Species in the Southeast United States

Ozone (O3) is a greenhouse gas and toxic pollutant which plays a major role in air quality and atmospheric chemistry. The understanding and ability to model the horizontal and vertical structure of O3 mixing ratios is difficult due to the complex formation/destruction processes and transport pathways that cause large variability of O3. The Environmental Protection Agency has National Ambient Air Quality Standards for O3 set at 75 ppb with future standards proposed to be as low as 65 ppb. These lower values emphasize the need to better understand/simulate the transport processes, emission sources, and chemical processes controlling precursor species (e.g., NOx, VOCs, and CO) which influence O3 mixing ratios. The uncertainty of these controlling variables is particularly large in the southeast United States (US) which is a region impacted by multiple different emission sources of precursor species (anthropogenic and biogenic) and transport processes resulting in complex spatio-temporal O3 patterns. During this work we will evaluate O3 and precursor species in the southeast US applying models, ground-based and airborne in situ data, and lidar observations. In the summer of 2013, the UAH O3 Differential Absorption Lidar (DIAL) (part of the Tropospheric Ozone Lidar Network (TOLNet)) measured vertical O3 profiles from the surface up to approximately 12 km. During this period, the lidar observed numerous periods of dynamic temporal and vertical O3 structures. In order to determine the sources/processes impacting these O3 mixing ratios we will apply the CTM GEOS-Chem (v9-02) at a 0.25 deg x 0.3125 deg resolution. Using in situ ground-based (e.g., SEARCH Network, CASTNET), airborne (e.g., NOAA WP-3D - SENEX 2013, DC-8 - SEAC4RS), and TOLNet lidar data we will first evaluate the model to determine the capability of GEOS-Chem to simulate the spatio-temporal variability of O3 in the southeast US. Secondly, we will perform model sensitivity studies in order to quantify which emission sources (e.g., anthropogenic, biogenic, lighting, wildfire) and transport processes (e.g., stratospheric, long-range, local scale) are contributing to these TOLNet-observed dynamic O3 patterns. Results from the evaluation of the model and the study of sources/processes impacting observed O3 mixing ratios will be presented.

Optical Radar↗
NASA NTRSDate not supplied

Continuum Covariance Propagation for Understanding Variance Loss in Advective Systems

We demonstrate for state dynamics governed by the continuity equation and related hyperbolic partial differential equations that significant, spurious variance loss occurs during covariance propagation by traditional methods used in data assimilation, even at full rank. This inaccurate variance evolution is caused not by numerical dissipation, but rather by a discontinuous change in the continuum covariance dynamics as correlation lengths tend to zero.

Shay Gilpin↗
NASA NTRS1980

Thermal Pollution Mathematical Model. Volume 3: User's Manual for One-Dimensional Numerical Model for the Seasonal Thermocline

A user's manual for a one dimensional thermal model to predict the temperature profiles of a deep body of water for any number of annual cycles is presented. The model is essentially a set of partial differential equations which are solved by finite difference methods using a high speed digital computer. The model features the effects of area change with depth, nonlinear interaction of wind generated turbulence and buoyancy, adsorption of radiative heat flux below the surface, thermal discharges, and the effects of vertical convection caused by discharge. The main assumption in the formulation is horizontal homogeneity. The environmental impact of thermal discharges from power plants is emphasized. Although the model is applicable to most lakes, a specific site (Lake Keowee, S.C.) application is described in detail. The programs are written in FORTRAN 5.

Lee, S. S.↗
NASA NTRS1986

Stability with large step sizes for multistep discretizations of stiff ordinary differential equations

One-leg and multistep discretizations of variable-coefficient linear systems of ODEs having both slow and fast time scales are investigated analytically. The stability properties of these discretizations are obtained independent of ODE stiffness and compared. The results of numerical computations are presented in tables, and it is shown that for large step sizes the stability of one-leg methods is better than that of the corresponding linear multistep methods.

Majda, George↗