Engineering Papers⌕ Search

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.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

An engineering study of hybrid adaptation of wind tunnel walls for three-dimensional testing

Solid wall tunnels having only upper and lower walls flexing are described. An algorithm for selecting the wall contours for both 2 and 3 dimensional wall flexure is presented and numerical experiments are used to validate its applicability to the general test case of 3 dimensional lifting aircraft models in rectangular cross section wind tunnels. The method requires an initial approximate representation of the model flow field at a given lift with wallls absent. The numerical methods utilized are derived by use of Green's source solutions obtained using the method of images; first order linearized flow theory is employed with Prandtl-Glauert compressibility transformations. Equations are derived for the flexed shape of a simple constant thickness plate wall under the influence of a finite number of jacks in an axial row along the plate centerline. The Green's source methods are developed to provide estimations of residual flow distortion (interferences) with measured wall pressures and wall flow inclinations as inputs.

Brown, Clinton↗
NASA NTRS1988

Superbubbles in disk galaxies

Correlated supernovae from an OB association create a superbubble: a large, thin, shell of cold gas surrounding a hot pressurized interior. Because supernova blast waves usually become subsonic before reaching the walls of the shell or cooling radiatively, the energy input from supernovae may be reasonably approximated as a continuous luminosity. Using the Kompaneets (thin-shell) approximation, the growth of superbubbles in various stratified atmospheres is numerically modeled. A dimensionless quantity predicts whether a superbubble will blow out of the H I disk of a spiral galaxy (and begin to accelerate upward) or collapse. Superbubbles blow out of the H I layer when they have a radius in the plane between one and two scale heights. They blow out only one side of a disk galaxy if their centers are more than 50-60 p above the plane and the gas layer has density and scale height typical of the Milky Way. Fingers of warm interstellar gas intrude into the hot interior when the superbubble overtakes dense clouds.

Mac Low, Mordecai-Mark↗
NASA NTRS1988

Statistical analysis of effective singular values in matrix rank determination

A major problem in using SVD (singular-value decomposition) as a tool in determining the effective rank of a perturbed matrix is that of distinguishing between significantly small and significantly large singular values to the end, conference regions are derived for the perturbed singular values of matrices with noisy observation data. The analysis is based on the theories of perturbations of singular values and statistical significance test. Threshold bounds for perturbation due to finite-precision and i.i.d. random models are evaluated. In random models, the threshold bounds depend on the dimension of the matrix, the noisy variance, and predefined statistical level of significance. Results applied to the problem of determining the effective order of a linear autoregressive system from the approximate rank of a sample autocorrelation matrix are considered. Various numerical examples illustrating the usefulness of these bounds and comparisons to other previously known approaches are given.

Konstantinides, Konstantinos↗
NASA NTRS1989

Solitosynthesis: Cosmological evolution of non-topological solitons

The thermal creation, fusion, evaporation, and destruction of non-topological solitons (NTS) after a phase transition in the early universe is considered. By defining and following NTS statistical equilibrium and departures from it, and depending on particle physics parameters, one of three possible scenarios occurs. If reaction rates are high enough, a period of equilibrium occurs and relic abundances are determined by the freeze-out temperature. Equilibrium first drives most NTS's into their constituents (free phi particles) and then causes rapid fusion into large NTS's. If freeze-out occurs during the first phase, the NTS's are almost entirely destroyed, while if it occurs during the second phase, solitosynthesis occurs and NTS's may be cosmically relevant. For slow reaction rates the NTS's are born frozen out and have the abundance determined by the phase transition. Analytic approximations for determining the abundances are developed, and tested by numerically integrating a reaction network in an expanding universe. Unfortunately, for most of the parameter space considered, solito-destruction/evaporation occurs.

Griest, Kim↗
NASA NTRS1989

Protostellar formation in rotating interstellar clouds. VIII - Inner core formation

The results are presented of a variety of spherically symmetric one-dimensional (1D) calculations intended to determine the robustness of the dynamical hiccup phenomenon in protostellar cores. The 1D models show that the phenomenon is relatively insensitive to changes in the equations of state, numerical resolution, initial density and temperature, and the radiative transfer approximation. In 1D, the hiccup results in an explosive destruction of the entire inner protostellar core. Inner core formation is studied with a sequence of three-dimensional models which show that rapid inner core rotation stabilizes the hiccup instability. Instead, the inner core becomes quite flat and undergoes a cycle of binary fragmentation, binary decay into a single object surrounded by a bar, breakup of the bar into a binary, etc. When lesser amounts of rotation are involved, the inner core does hiccup somewhat, but mass is ejected in only a few directions, leading to several broad streams of ejecta.

Boss, Alan P.↗
NASA NTRS1991

Calculating Viscous/Inviscid Interactions

Report describes some experiences in numerical simulation of compressible, viscous flows by finite-difference approximations of Navier-Stokes equations. Focuses upon viscous/inviscid-interaction approach, in which flow mathematically modeled with interacting zones, within each of which different model of viscous or inviscid flow applies. Reduces task of computation by taking advantage of simpler models and coarser grids in some of zones.

Van Dalsem, W. R.↗
NASA NTRS1989

Unstructured mesh solution of the Euler and Navier-Stokes equations

Mesh generation procedures as well as solution algorithms for solving the Euler and Navier-Stokes equations on unstructured meshes are presented. The solution algorithms discussed utilize approximate Riemann solver, upwind differencing to achieve high spatial accuracy. Numerical results for Euler flow over single and multi-element airfoils are presented.

Barth, Timothy J.↗
NASA NTRS1991

The thermal conditions of Venus

Models of Venus' thermal evolution are examined. The following subject areas are covered: (1) modified approximation of parameterized convection; (2) description of the model; (3) numerical results and asymptotic solution of the MAPC equations; (4) magnetism and the thermal regime of the cores of Earth and Venus; and (5) the thermal regime of the Venusian crust.

Zharkov, Vladimir N.↗
NASA NTRS1991

Infrared radiation parameterizations in numerical climate models

This study presents various approaches to parameterizing the broadband transmission functions for utilization in numerical climate models. One-parameter scaling is applied to approximate a nonhomogeneous path with an equivalent homogeneous path, and the diffuse transmittances are either interpolated from precomputed tables or fit by analytical functions. Two-parameter scaling is applied to parameterizing the carbon dioxide and ozone transmission functions in both the lower and middle atmosphere. Parameterizations are given for the nitrous oxide and methane diffuse transmission functions.

Chou, Ming-Dah↗
NASA NTRS1991

Acoustic fatigue life prediction for non-linear structures

Using an approach based on a time domain analysis, a method of equivalent linearization is applied for an estimation of the strain response of complex nonlinear structures having nearly arbitrary complexity. Fatigue lives estimated for a nonlinear beam with random excitation using the approximate method were compared with results obtained using a conventional numerical simulation, yielding nearly identical results.

Sun, J. Q.↗
NASA NTRS1992

Analytical shape sensitivities and approximations of modal response of generally laminated tapered skew plates

In the present determination of the derivatives of natural frequencies and mode shapes of a generally laminated tapered skew plate, with respect to various shape parameters, springs are used to simulate the essential boundary conditions. The independent shape parameters are plate surface area, aspect ratio, taper ratio, and sweep angle. Eigenvalues and eigenvectors are approximated over the range of the variable using linear, exponential, and pseudoexponential approximation schemes, and are compared with reanalysis-obtained values. Numerical results are presented for symmetrically and unsymmetrically laminated plates.

Singhvi, Sarvesh↗
NASA NTRS1994

Splines and control theory

In this work, the relationship between splines and the control theory has been analyzed. We show that spline functions can be constructed naturally from the control theory. By establishing a framework based on control theory, we provide a simple and systematic way to construct splines. We have constructed the traditional spline functions including the polynomial splines and the classical exponential spline. We have also discovered some new spline functions such as trigonometric splines and the combination of polynomial, exponential and trigonometric splines. The method proposed in this paper is easy to implement. Some numerical experiments are performed to investigate properties of different spline approximations.

Zhang, Zhimin↗
NASA NTRS1995

The construction of high-accuracy schemes for acoustic equations

An accuracy analysis of various high order schemes is performed from an interpolation point of view. The analysis indicates that classical high order finite difference schemes, which use polynomial interpolation, hold high accuracy only at nodes and are therefore not suitable for time-dependent problems. Thus, some schemes improve their numerical accuracy within grid cells by the near-minimax approximation method, but their practical significance is degraded by maintaining the same stencil as classical schemes. One-step methods in space discretization, which use piecewise polynomial interpolation and involve data at only two points, can generate a uniform accuracy over the whole grid cell and avoid spurious roots. As a result, they are more accurate and efficient than multistep methods. In particular, the Cubic-Interpolated Psuedoparticle (CIP) scheme is recommended for computational acoustics.

Tang, Lei↗
NASA NTRS1994

Zero potential vorticity envelopes for the zonal-mean velocity of the Venus/Titan atmospheres

The diagnostic analysis of numerical simulations of the Venus/Titan wind regime reveals an overlooked constraint upon the latitudinal structure of their zonal-mean angular momentum. The numerical experiments, as well as the limited planetary observations, are approximately consistent with the hypothesis that within the latitudes bounded by the wind maxima the total Ertel potential vorticity associated with the zonal-mean motion is approximately well mixed with respect to the neutral equatorial value for a stable circulation. The implied latitudinal profile of angular momentum is of the form M equal to or less than M(sub e)(cos lambda)(exp 2/Ri), where lambda is the latitude and Ri the local Richardson number, generally intermediate between the two extremes of uniform angular momentum (Ri approaches infinity) and uniform angular velocity (Ri = 1). The full range of angular momentum profile variation appears to be realized within the observed meridional - vertical structure of the Venus atmosphere, at least crudely approaching the implied relationship between stratification and zonal velocity there. While not itself indicative of a particular eddy mechanism or specific to atmospheric superrotation, the zero potential vorticity (ZPV) constraint represents a limiting bound for the eddy - mean flow adjustment of a neutrally stable baroclinic circulation and may be usefully applied to the diagnostic analysis of future remote sounding and in situ measurements from planetary spacecraft.

Allison, Michael↗
NASA NTRS1995

Structural Damage Prediction and Analysis for Hypervelocity Impact: Properties of Largest Fragment Produced by Hypervelocity Impact of Aluminum Spheres with Thin Aluminum Sheets

Results of a series of hypervelocity impact tests are presented. In these tests, 1.275-g, 9.53-mm-diameter, 2017-T4 aluminum spheres were fired at normal incidence at eight thicknesses of 6061-T6 aluminum sheet. Bumper thickness to projectile diameter (t/D) ratio ranged from 0.026 to 0.424. Nominal impact velocity was 6.7 km/s. Results of five tests using 6.35, 9.53, and 12.70-mm-diameter aluminum spheres and other aluminum alloy bumpers are also given. A large chunky fragment of projectile was observed at the center of the debris clouds produced by the impacts. The equivalent diameter of this large fragment ranged from 5.5 mm for the lowest t/D ratio to a minimum of 0.6 mm for the case where maximum breakup of the projectile occurred (t/D approximately 0.2 to 0.3). When the t/D ratio was 0.42, numerous large flaky fragments were evenly distributed in the external bubble of bumper debris. Velocity of the large central fragments decreased continuously with increasing t/D ratio, ranging from about 99 percent to less than 80 percent of the impact velocity. The change in the velocity of small fragments spalling from the rear of the projectile was used to obtain a relationship showing a linear increase in the size of the central projectile fragment with decrease in the shock-induced stress in the projectile.

Piekutowski, Andrew J.↗
NASA NTRS1998

A Preconditioning Method for Shape Optimization Governed by the Euler Equations

We consider a classical aerodynamic shape optimization problem subject to the compressible Euler flow equations. The gradient of the cost functional with respect to the shape variables is derived with the adjoint method at the continuous level. The Hessian (second order derivative of the cost functional with respect to the shape variables) is approximated also at the continuous level, as first introduced by Arian and Ta'asan (1996). The approximation of the Hessian is used to approximate the Newton step which is essential to accelerate the numerical solution of the optimization problem. The design space is discretized in the maximum dimension, i.e., the location of each point on the intersection of the computational mesh with the airfoil is taken to be an independent design variable. We give numerical examples for 86 design variables in two different flow speeds and achieve an order of magnitude reduction in the cost functional at a computational effort of a full solution of the analysis partial differential equation (PDE).

Arian, Eyal↗
NASA NTRS2014

Effects of Nonsphericity on the Behavior of Lorenz-Mie Resonances in Scattering Characteristics of Liquid-Cloud Droplets

By using the results of highly accurate T-matrix computations for randomly oriented oblate and prolate spheroids and Chebyshev particles with varying degrees of asphericity, we analyze the effects of a deviation of water-droplet shapes from that of a perfect sphere on the behavior of Lorenz-Mie morphology-dependent resonances of various widths. We demonstrate that the positions and profiles of the resonances can change significantly with increasing asphericity. The absolute degree of asphericity required to suppress a Lorenz-Mie resonance is approximately proportional to the resonance width. Our results imply that numerical averaging of scattering characteristics of real cloud droplets over sizes may rely on a significantly coarser size-parameter resolution than that required for ideal, perfectly spherical particles.

matrices (mathematics)↗
NASA NTRSDate not supplied

Validating Drag and Heating Coefficients for Hollow Reentry Objects in Continuum Flow Using a Mach 7 Ludwieg Tube

Drag and heating coefficient databases and models are crucial to destructive reentry simulation. The NASA Orbital Debris Program Office (ODPO) develops, maintains, and performs analysis with the Object Reentry Survival Analysis Tool (ORSAT), which comprises drag and heating models for free molecular, transitional, and continuum flow regimes. These models have, in the past, only included solid, convex, blunt shapes (such as boxes, spheres, and cylinders). Previous work led by ODPO includes the extension of these models to hollow cylinders and square boxes in free molecular and transitional flow using the Direct Simulation Monte Carlo (DSMC) method. Since 2019, the ODPO has continued its program of DSMC simulations and extended the project to include analyses with the NASA Data Parallel Line Relaxation (DPLR) program on hollow cylinders and boxes (with varying wall thickness-diameter ratio). In fall 2022, the ODPO began a collaboration with the University of Texas San Antonio (UTSA) to use the Mach 7 Ludwieg Tube facility to validate the model built using numerical simulations. This facility can replicate (at a scale of approximately 100:1) the conditions seen by reentering objects near typical demise altitudes. We present here the drag and heating coefficients derived from the continued DSMC simulations, the new DPLR simulations, and the 26-test series at UTSA.

Chris Ostrom↗