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.

SourceResearch connectionFrom yearThrough year

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

At least 991 records · Page 55

All sources NASA NTRS DOE OSTI

NASA NTRS1977

Development of prediction techniques for multi-jet thermal ground flow fields and fountain formation

Avoiding detrimental ground interaction is important for practical V/STOL aircraft. This paper reports recent developments in a numerical method for estimating thermal ground footprints. Upwash and fountain formation for arbitrarily oriented jet arrangements is predicted. Flow asymmetry due to roll, pitch, differential thrust or ground inclination is included. The prediction methodology uses simple inviscid relations for energy and momentum conservation along with an empirical entrainment law, applied in independent sectors of the wall jet and upwash. Asymmetrical stagnation line prediction is compared with experiment. Detailed flow measurements for a three-jet interaction are also presented.

Siclari, M. J.↗

NASA NTRS1980

Computation of three-dimensional viscous supersonic flow in inlets

A new approach has been developed for the computation of the three-dimensional viscous supersonic flow with embedded subsonic regions adjacent to solid boundaries and is applied to a mixed-compression supersonic inlet typical of current designs. The approach uses a reduced form of the three-dimensional Navier-Stokes equations so that the resultant equations can be treated as an initial boundary value problem and thus be solved by non-iterative forward marching in space. The numerical procedure utilizes an efficient consistently-split linearized block implicit technique to solve the finite difference analogues to the set of governing partial differential equations.

Buggeln, R. C.↗

NASA NTRS1988

Applications of an exponential finite difference technique

An exponential finite difference scheme first presented by Bhattacharya for one dimensional unsteady heat conduction problems in Cartesian coordinates was extended. The finite difference algorithm developed was used to solve the unsteady diffusion equation in one dimensional cylindrical coordinates and was applied to two and three dimensional conduction problems in Cartesian coordinates. Heat conduction involving variable thermal conductivity was also investigated. The method was used to solve nonlinear partial differential equations in one and two dimensional Cartesian coordinates. Predicted results are compared to exact solutions where available or to results obtained by other numerical methods.

Handschuh, Robert F.↗

NASA NTRS1990

Fourier analysis of finite element preconditioned collocation schemes

The spectrum of the iteration operator of some finite element preconditioned Fourier collocation schemes is investigated. The first part of the paper analyses one-dimensional elliptic and hyperbolic model problems and the advection-diffusion equation. Analytical expressions of the eigenvalues are obtained with use of symbolic computation. The second part of the paper considers the set of one-dimensional differential equations resulting from Fourier analysis (in the tranverse direction) of the 2-D Stokes problem. All results agree with previous conclusions on the numerical efficiency of finite element preconditioning schemes.

Deville, Michel O.↗

NASA NTRS1993

Geology of Galileo Regio quadrangle, Ganymede

Galileo Regio quadrangle, including most of Galileo Regio, contains two basic geologic units: older dark, furrowed, and heavily cratered material and younger light, grooved, relatively less-cratered material. Dark material dominates, the light material occurring to the W and SW in about five percent of the quadrangle. Dominant structural features, of uncertain origin, are the NW-trending arcuate furrows and associated orthogonal and oblique furrows in the dark material; grooves in the light material are not as well-ordered nor as intensely developed as elsewhere on Ganymede. Fifteen palimpsests, located partly or totally in the quadrangle and representing the period when the crust became sufficiently strong to record impact events, are subdivided into three relative-age classes, as are the numerous impact craters, some of which are dome, moat, and rampart-like. Galileo Regio is believed to have formed when, because of cooling, differentiation, and expansion, the thin, weak, icy lithosphere broke into blocks and light material filled the space between, or overspread those blocks that subsided, or both. Grooves formed in most of the light material and both light and dark material were cratered by impact, the most recent of which created bright-ray craters, the youngest features visible in the quadrangle.

Underwood, James R., Jr.↗

NASA NTRS1993

Approximate optimal guidance for the advanced launch system

A real-time guidance scheme for the problem of maximizing the payload into orbit subject to the equations of motion for a rocket over a spherical, non-rotating earth is presented. An approximate optimal launch guidance law is developed based upon an asymptotic expansion of the Hamilton - Jacobi - Bellman or dynamic programming equation. The expansion is performed in terms of a small parameter, which is used to separate the dynamics of the problem into primary and perturbation dynamics. For the zeroth-order problem the small parameter is set to zero and a closed-form solution to the zeroth-order expansion term of Hamilton - Jacobi - Bellman equation is obtained. Higher-order terms of the expansion include the effects of the neglected perturbation dynamics. These higher-order terms are determined from the solution of first-order linear partial differential equations requiring only the evaluation of quadratures. This technique is preferred as a real-time, on-line guidance scheme to alternative numerical iterative optimization schemes because of the unreliable convergence properties of these iterative guidance schemes and because the quadratures needed for the approximate optimal guidance law can be performed rapidly and by parallel processing. Even if the approximate solution is not nearly optimal, when using this technique the zeroth-order solution always provides a path which satisfies the terminal constraints. Results for two-degree-of-freedom simulations are presented for the simplified problem of flight in the equatorial plane and compared to the guidance scheme generated by the shooting method which is an iterative second-order technique.

Feeley, T. S.↗

NASA NTRS1995

Effects of pressure gradients on turbulent premixed flames

The influence of a constant acceleration on a turbulent premixed flame is studied by direct numerical simulation. This acceleration induces a mean pressure gradient across the flame brush, leading to a modification of the turbulent flame structure due to differential buoyancy mechanisms between heavy cold fresh and light hot burnt gases. Such a pressure gradient may be encountered in practical applications in ducted flames. A favorable pressure gradient, i.e. the pressure decreases from unburnt to burnt gases, is found to decrease the flame wrinkling, the flame brush thickness, and the turbulent flame speed. A favorable pressure gradient also promotes counter-gradient turbulent transport. On the other hand, adverse pressure gradients tend to increase the flame brush thickness and turbulent flame speed, and promote classical gradient turbulent transport. The balance equation for the turbulent flux of the Favre averaged progress variable is also analyzed. The first results show that the fluctuating pressure term, cannot be neglected as generally assumed in models. Simple models assuming that a high mean pressure gradient may only be balanced by the cross-dissipation term seem too approximate. This analysis has to be continued to compare simulation data and closure schemes proposed for the transport equation. The analysis developed by Veynante et al.(1995) has been extended to imposed acceleration and mean pressure gradients. A simple model for the turbulent flux is proposed and validated from simulation data. Then, a modified criterion is derived to delineate between counter-gradient and gradient turbulent diffusion. In fact, counter-gradient diffusion may occur in most practical applications, especially for ducted flames.

Veynante, D.↗

NASA NTRS1996

Approach to Synchronization Control of Magnetic Bearings Using Fuzzy Logic

This paper presents a fuzzy-logic approach to the synthesis of synchronization control for magnetically suspended rotor system. The synchronization control enables a whirling rotor to undergo synchronous motion along the magnetic bearing axes; thereby avoiding the gyroscopic effect that degrade the stability of rotor systems when spinning at high speed. The control system features a fuzzy controller acting on the magnetic bearing device, in which the fuzzy inference system trained through fuzzy rules to minimize the differential errors between four bearing axes so that an error along one bearing axis can affect the overall control loop for the motion synchronization. Numerical simulations of synchronization control for the magnetically suspended rotor system are presented to show the effectiveness of the present approach.

Yang, Li-Farn↗

NASA NTRS2000

A Dream of a Mission: Stellar Imager and Seismic Probe

The Stellar Imager and Seismic Probe (SISP) is a mission to understand the various effects of magnetic fields of stars, the dynamos that generate them, and the internal structure and dynamics of the stars in which they exist. The ultimate goal is to achieve the best-possible forecasting of solar activity on times scales ranging up to decades, and an understanding of the impact of stellar magnetic activity on astrobiology and life in the Universe. The road to that goal will revolutionize our understanding of stars and stellar systems, the building blocks of the Universe. SISP will zoom in on what today - with few exceptions - we only know as point sources, revealing processes never before seen, thus providing a tool to astrophysics as fundamental as the microscope is to the study of life on Earth. SISP is an ultraviolet aperture-synthesis imager with 8-10 telescopes with meter-class apertures, and a central hub with focal-plane instrumentation that allows spectrophotometry in passbands as narrow as a few Angstroms up to hundreds of Angstroms. SISP will image stars and binaries with one hundred to one thousand resolution elements on their surface, and sound their interiors through asteroseismology to image internal structure, differential rotation, and large-scale circulations; this will provide accurate knowledge of stellar structure and evolution and complex transport processes, and will impact numerous branches of (astro)physics ranging from the Big Bang to the future of the Universe. Fitting naturally within the NASA long-term time line, SISP complements defined missions, and with them will show us entire other solar systems, from the central star to their orbiting planets.

Carpenter, Kenneth G.↗

NASA NTRS1956

On the Theory of Anisotropic Shallow Shells

A numerical analysis of thin-walled shallow shells is presented. Equations of equilibrium and relations between deformations and stresses are included along with fundamental differential equations.

Ambartsumyan, S. A.↗

NASA NTRS2001

Auxins and tropisms

Differential growth of plants in response to the changes in the light and gravity vectors requires a complex signal transduction cascade. Although many of the details of the mechanisms by which these differential growth responses are induced are as yet unknown, auxin has been implicated in both gravitropism and phototropism. Specifically, the redistribution of auxin across gravity or light-stimulated tissues has been detected and shown to be required for this process. The approaches by which auxin has been implicated in tropisms include isolation of mutants altered in auxin transport or response with altered gravitropic or phototropic response, identification of auxin gradients with radiolabeled auxin and auxin-inducible gene reporter systems, and by use of inhibitors of auxin transport that block gravitropism and phototropism. Proteins that transport auxin have been identified and the mechanisms which determine auxin transport polarity have been explored. In addition, recent evidence that reversible protein phosphorylation controls this process is summarized. Finally, the data in support of several hypotheses for mechanisms by which auxin transport could be differentially regulated during gravitropism are examined. Although many details of the mechanisms by which plants respond to gravity and light are not yet clear, numerous recent studies demonstrate the role of auxin in these processes.

NASA Discipline Plant Biology↗

NASA NTRS2009

Developing Lightning Prediction Tools for the CCAFS Dual-Polarimetric Radar

NASA Marshall Space Flight Center and the University of Alabama Huntsville are collaborating with the 45th Weather Squadron (45WS) to develop improved lightning prediction capabilities for the new C-band dual-polarimetric weather radar being acquired for use by 45WS and launch weather forecasters at Cape Canaveral Air Force Station (CCAFS). In particular, these algorithms will focus on lightning onset, cessation and combined lightning-radar applications for convective winds assessment. Research using radar reflectivity (Z) data for prediction of lightning onset has been extensively discussed in the literature and subsequently applied by launch weather forecasters as it pertains to lightning nowcasting. Currently the forecasters apply a relatively straight forward but effective temperature-Z threshold algorithm for assessing the likelihood of lightning onset in a given storm. In addition, a layered VIL above the freezing level product is used as automated guidance for the onset of lightning. Only limited research and field work has been conducted on lightning cessation using Z and vertically-integrated Z for determining cessation. Though not used operationally vertically-integrated Z (basis for VIL) has recently shown promise as a tool for use in nowcasting lightning cessation. The work discussed herein leverages and expands upon these and similar reflectivity-threshold approaches via the application/addition of over two decades of polarimetric radar research focused on distinct multi-parameter radar signatures of ice/mixed-phase initiation and ice-crystal orientation in highly electrified convective clouds. Specifically, our approach is based on numerous previous studies that have observed repeatable patterns in the behavior of the vertical hydrometeor column as it relates to the temporal evolution of differential reflectivity and depolarization (manifested in either LDR or p(sub hv)), development of in-situ mixed and ice phase microphysics, electric fields, and ensuing lightning in the sub-tropical/tropical convection typical of the southeastern U.S., Maritime Continent, and southwestern Amazon. The polarimetric signatures detected in this setting provide a basis for automated 3-D detection of hydrometeor types in fuzzy logic hydrometeor identification algorithms (HID). Our working hypothesis is that improvement in lightning onset warning lead time and specificity for a given storm, relative to application of a Z-threshold algorithm, should arise as a consequence of the ability of dual-polarimetric radar to unambiguously detect and identify (through HID algorithms) the updraft elevation of rain-water cores above the freezing level and subsequent onset of drop freezing, riming, and robust mixed phase processes leading to significant charge separation and lightning. This type of algorithm, though dependent on the quality of the polarimetric data should be less susceptible to variable Z-calibration that can impact a given Z-threshold approach. To facilitate development of the algorithm while the 45WS dual-pol radar is in its current test stages and to evaluate the impact of polarimetric data quality (e.g., modified scan parameters and sampling) on the ensuing algorithms, we are using the ARMOR C-band dual-pol radar in Huntsville combined with N. Alabama LMA data and ARMOR HID algorithms [NCAR algorithm modified for application at C-band] in a testbed fashion. For lightning cessation we are revisiting the application of differential propagation phase variables for the monitoring of ice crystal alignment driven by in-cloud electric fields combined with metrics of ice water path (i.e., vertically integrated reflectivity). Importantly it should be noted that this approach is still very much a research topic and as such, we will explore operational applications that involve radar frequencies other than C-Band by using the UAH MAX X-band dual-pol radar in slow staring modes.

Petersen, W. A.↗

NASA NTRS1988

Numerical modeling of multidimensional flow in seals and bearings used in rotating machinery

The rotordynamic behavior of turbomachinery is critically dependent on fluid dynamic rotor forces developed by various types of seals and bearings. The occurrence of self-excited vibrations often depends on the rotor speed and load. Misalignment and rotor wobbling motion associated with differential clearance were often attributed to stability problems. In general, the rotative character of the flowfield is a complex three dimensional system with secondary flow patterns that significantly alter the average fluid circumferential velocity. A multidimensional, nonorthogonal, body-fitted-grid fluid flow model is presented that describes the fluid dynamic forces and the secondary flow pattern development in seals and bearings. Several numerical experiments were carried out to demonstrate the characteristics of this complex flowfield. Analyses were performed by solving a conservation form of the three dimensional Navier-Stokes equations transformed to those for a rotating observer and using the general-purpose computer code PHOENICS with the assumptions that the rotor orbit is circular and that static eccentricity is zero. These assumptions have enabled a precise steady-state analysis to be used. Fluid injection from ports near the seal or bearing center increased fluid-film direct dynamic stiffness and, in some cases, significantly increased quadrature dynamic stiffness. Injection angle and velocity could be used for active rotordynamic control; for example, injection, when compared with no injection, increased direct dynamic stiffness, which is an important factor for hydrostatic bearings.

Hendricks, R. C.↗

NASA NTRS2022

A Dynamo Simulation Generating Saturn-Like Small Magnetic Dipole Tilts

Abstract Among planetary dynamos, the magnetic field of Saturn stands out in its exceptional level of axisymmetry. One of its peculiar features is that the magnetic dipole mode is tilted with respect to the planetary rotation axis by only ≈0.007° or less. Numerical dynamo simulations performed in this context have had great difficulty in producing such small dipole tilt angles without introducing ad hoc ingredients such as a latitudinally varying heat flux pattern in the outer layers or stably stratified layers (SSLs). Here we present a numerical dynamo simulation that generates a highly axisymmetric dynamo with a dipole tilt of about ≈0.0008° on average. The model consists of a deep dynamo layer and an overlying low-conductivity layer but without any SSLs. We highlight a novel mechanism where strong differential rotation generated in the atmospheric layer penetrates into the dynamo region, helping to maintain a very small magnetic dipole tilt. Plain Language Summary Saturn's dipole-dominant magnetic field exhibits a very peculiar feature: the dipole component of the planetary magnetic field is tilted by less than ≈0.007° with respect to the planetary spin axis. Numerical simulations performed in this context suggest that if a spatial heat-flux variation is imposed, along with a stably stratified region, on top of an active dynamo layer, then small dipole tilt values can be realized. Here we present a model where extremely small dipole tilt values can be achieved without these ad hoc ingredients. Our simulations demonstrate that dynamo theory allows extremely small dipole tilt values in a relatively simple model configuration.

Rakesh Kumar Yadav↗

NASA NTRS1985

A fully implicit scheme for global numerical weather prediction

A fast-slow factored scheme is presented for use with shallow-water primitive equation numerical weather prediction models. The technique was developed to reduce the rotational mode errors which arise when the fast and slow terms of the governing differential equations are treated simultaneously. The method factors out the fast and slow terms along the coordinate directions by means of a modified Crank-Nicolson scheme. A finite-difference spatial discretization is carried out in the zonal and meridional directions to reduce the factorization error to near-zero, and that time steps of 60-90 min can be used to obtain acceptably accurate results, even in the presence of fine spatial structures in the flow.

Augenbaum, J. M.↗

NASA NTRS1988

Three-dimensional adaptive grid generation for body-fitted coordinate system

This report describes a numerical method for generating 3-D grids for general configurations. The basic method involves the solution of a set of quasi-linear elliptic partial differential equations via pointwise relaxation with a local relaxation factor. It allows specification of the grid spacing off the boundary surfaces and the grid orthogonality at the boundary surfaces. It includes adaptive mechanisms to improve smoothness, orthogonality, and flow resolution in the grid interior.

Chen, S. C.↗

NASA NTRS1988

Three-dimensional adaptive grid generation for body-fitted coordinate system

Chen, S. C.↗

NASA NTRS1990

A fully Sinc-Galerkin method for Euler-Bernoulli beam models

A fully Sinc-Galerkin method in both space and time is presented for fourth-order time-dependent partial differential equations with fixed and cantilever boundary conditions. The Sinc discretizations for the second-order temporal problem and the fourth-order spatial problems are presented. Alternate formulations for variable parameter fourth-order problems are given which prove to be especially useful when applying the forward techniques to parameter recovery problems. The discrete system which corresponds to the time-dependent partial differential equations of interest are then formulated. Computational issues are discussed and a robust and efficient algorithm for solving the resulting matrix system is outlined. Numerical results which highlight the method are given for problems with both analytic and singular solutions as well as fixed and cantilever boundary conditions.

Smith, R. C.↗