SEARCH · Engineering Papers

Results for “numerical approximation & analysis”

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 505 records · Page 28

All sources NASA NTRS DOE OSTI

NASA NTRS1972

A numerical investigation of film cooling as an alternate thermal protection system for spacecraft during atmospheric reentry

Research was undertaken to determine the mass required to protect a surface during atmospheric reentry by injecting a cold gas tangentially along the surface. A consideration of this problem revealed that the boundary layer equations could be employed to approximate the most important flow configurations which occur in film cooling. However, even with this simplification, the equations to be solved are complex, nonlinear, coupled partial differential equations. Due to the relative weakness of classical methods of boundary layer analysis, finite difference methods were chosen as the best approach to the formidable problem under consideration. Instead of blindly adopting a proposed technique and trying to force it to the problem, a more systematic study was deemed to be not only desirable but necessary. A relatively basic study of the application of finite difference methods was undertaken with the objective of developing an algorithm capable of obtaining accurate solutions to film cooling problems.

Clausing, A. M.↗

NASA NTRS1987

Decoupled direct method for sensitivity analysis in combustion kinetics

An efficient, decoupled direct method for calculating the first order sensitivity coefficients of homogeneous, batch combustion kinetic rate equations is presented. In this method the ordinary differential equations for the sensitivity coefficients are solved separately from , but sequentially with, those describing the combustion chemistry. The ordinary differential equations for the thermochemical variables are solved using an efficient, implicit method (LSODE) that automatically selects the steplength and order for each solution step. The solution procedure for the sensitivity coefficients maintains accuracy and stability by using exactly the same steplengths and numerical approximations. The method computes sensitivity coefficients with respect to any combination of the initial values of the thermochemical variables and the three rate constant parameters for the chemical reactions. The method is illustrated by application to several simple problems and, where possible, comparisons are made with exact solutions and those obtained by other techniques.

Radhakrishnan, Krishnan↗

NASA NTRS1961

Method for Approximating the Vacuum Motions of Spinning Symmetrical Bodies with Nonconstant Spin Rates

The analysis includes non-constant spin rates and inertias and considers the effects of time-varying thrust misalignments, mass unbalance, and jet damping. The method was developed for bodies having small trans verse angular velocities. Results are presented in the form of equations for space-referenced Euler angles, flight-path angles, body-referenced attitude rates, and earth-referenced vehicle-trajectory coordinates. Also, equations for maximum wobble have been derived for certain input conditions. Comparisons with numerical solutions are included for two sample problems.

Martz, C. William↗

NASA NTRS1979

Performance optimization of an MHD generator with physical constraints

A technique has been described which optimizes the power out of a Faraday MHD generator operating under a prescribed set of electrical and magnetic constraints. The method does not rely on complicated numerical optimization techniques. Instead the magnetic field and the electrical loading are adjusted at each streamwise location such that the resultant generator design operates at the most limiting of the cited stress levels. The simplicity of the procedure makes it ideal for optimizing generator designs for system analysis studies of power plants. The resultant locally optimum channel designs are, however, not necessarily the global optimum designs. The results of generator performance calculations are presented for an approximately 2000 MWe size plant. The difference between the maximum power generator design and the optimal design which maximizes net MHD power are described. The sensitivity of the generator performance to the various operational parameters are also presented.

Pian, C. C. P.↗

NASA NTRS2016

The Quality Control Algorithms Used in the Process of Creating the NASA Kennedy Space Center Lightning Protection System Towers Meteorological Database

The methodology and the results of the quality control (QC) process of the meteorological data from the Lightning Protection System (LPS) towers located at Kennedy Space Center (KSC) launch complex 39B (LC-39B) are documented in this paper. Meteorological data are used to design a launch vehicle, determine operational constraints, and to apply defined constraints on day-of-launch (DOL). In order to properly accomplish these tasks, a representative climatological database of meteorological records is needed because the database needs to represent the climate the vehicle will encounter. Numerous meteorological measurement towers exist at KSC; however, the engineering tasks need measurements at specific heights, some of which can only be provided by a few towers. Other than the LPS towers, Tower 313 is the only tower that provides observations up to 150 m. This tower is located approximately 3.5 km from LC-39B. In addition, data need to be QC'ed to remove erroneous reports that could pollute the results of an engineering analysis, mislead the development of operational constraints, or provide a false image of the atmosphere at the tower's location.

Orcutt, John M.↗

NASA NTRS1999

Satellite Data Assimilation

Satellite altimetry provides the only routine observation of a dynamic variable of the global ocean. For instance, measurements from the TOPEX/POSEIDON altimeter have provided dramatic basin-wide images of the space-time evolution of sea level associated with the 1997-98 El Nino event. The nature of large-scale ocean circulation is being studied that underlies such sea level changes measured by satellite altimetry. To first approximation, properties of large-scale (order 1000 km and larger) sea level variability is latitudinally dependent. An analysis using a general circulation model shows that sea level changes are largely due to wind-driven baroclinic (depth dependent) circulation in the tropics (within 20 degrees of the equator), but are primarily due to the expansion and contraction of near surface water in temperate latitudes (between 20 deg. and 40 deg.) forced by seasonal heating and cooling. In contrast, wind-driven barotropic (depth independent) circulation dominates sea level variability in high latitudes (40 deg. and higher) characterized by periods that are as short as a few days. The presence of such inhomogeneity and the significance of high-frequency, large-scale sea level changes had not been fully recognized prior to this study, and are summarized in Fukumori et al. Numerical models provide theoretical relationships among properties that can be inverted using observations so as to estimate the entire state of the ocean, including properties that are otherwise difficult to measure remotely. The process is data assimilation. An approximate Kalman filter and smoother have been devised to assimilate three years of TOPEX/POSEIDON sea level data into a global ocean general circulation model. The figure below demonstrates the skill of the assimilation, and shows the altimeter assimilated estimate being in closer agreement than the simulation is with independent in situ measurements of subsur ace temperature and velocity, consistent with formal uncertainty estimates. Results from this study have been summarized and submitted for publication. The study demonstrates the feasibility of global ocean data assimilation and illustrates applications in monitoring and understanding of processes controlling the evolution of the ocean. Additional information is contained in the original.

Fukumori, Ichiro↗

NASA NTRS1979

An efficient algorithm for numerical airfoil optimization

A new optimization algorithm is presented. The method is based on sequential application of a second-order Taylor's series approximation to the airfoil characteristics. Compared to previous methods, design efficiency improvements of more than a factor of 2 are demonstrated. If multiple optimizations are performed, the efficiency improvements are more dramatic due to the ability of the technique to utilize existing data. The method is demonstrated by application to subsonic and transonic airfoil design but is a general optimization technique and is not limited to a particular application or aerodynamic analysis.

Vanderplaats, G. N.↗

NASA NTRS1985

A three-dimensional viscous flow analysis for the helicopter tip vortex generation problem

The tip vortex flow field occurring in the vicinity of the tip region of a a helicopter rotor blade is a very complicated three-dimensional, viscous flow phenomenon. The details of the flow in the tip region can have a major effect in determining the generated rotor noise and can significantly affect the performance and dynamic loading of the rotor blade. The three-dimensional viscous subsonic tip vortex generation processes is investigated by a numerical procedure which allows spatial forward-marching integration, utilizing flow approximations from the velocity-decomposition approach of Briley and McDonald. The approach has been applied to compute the laminar and turbulent tip vortex flows for a constant thickness slab airfoil with a square tip, a constant thickness slab airfoil with a half round tip and a NACA 0012 airfoil with a half round tip. The basic mechanism of the tip vortex generation process as well as the prediction of vortex appearance, strength and secondary flow shown by the calculations are in qualitative agreement with experimental results.

Lin, S. J.↗

NASA NTRS1995

Effect of various approximations on predicted progressive failure in plain weave composites

Three-dimensional finite element analysis was used to simulate progressive failure of a plain weave composite subjected to in-plane extension. The loading was parallel to one of the tow directions. The effects of various characteristics of the finite element model on predicted behavior were examined. The predicted behavior was found to be sensitive to quadrature order, mesh refinement, and the material degradation model. Also the sensitivity of the predictions to the tow waviness was studied. The predicted strength decreased considerably with increased waviness. More numerical studies and comparisons with experimental data are needed to establish reliable guidelines for accurate progressive failure prediction.

Whitcomb, John↗

NASA NTRS2004

Variable Step Integration Coupled with the Method of Characteristics Solution for Water-Hammer Analysis, A Case Study

One-dimensional water-hammer modeling involves the solution of two coupled non-linear hyperbolic partial differential equations (PDEs). These equations result from applying the principles of conservation of mass and momentum to flow through a pipe, and usually the assumption that the speed at which pressure waves propagate through the pipe is constant. In order to solve these equations for the interested quantities (i.e. pressures and flow rates), they must first be converted to a system of ordinary differential equations (ODEs) by either approximating the spatial derivative terms with numerical techniques or using the Method of Characteristics (MOC). The MOC approach is ideal in that no numerical approximation errors are introduced in converting the original system of PDEs into an equivalent system of ODEs. Unfortunately this resulting system of ODEs is bound by a time step constraint so that when integrating the equations the solution can only be obtained at fixed time intervals. If the fluid system to be modeled also contains dynamic components (i.e. components that are best modeled by a system of ODEs), it may be necessary to take extremely small time steps during certain points of the model simulation in order to achieve stability and/or accuracy in the solution. Coupled together, the fixed time step constraint invoked by the MOC, and the occasional need for extremely small time steps in order to obtain stability and/or accuracy, can greatly increase simulation run times. As one solution to this problem, a method for combining variable step integration (VSI) algorithms with the MOC was developed for modeling water-hammer in systems with highly dynamic components. A case study is presented in which reverse flow through a dual-flapper check valve introduces a water-hammer event. The predicted pressure responses upstream of the check-valve are compared with test data.

Turpin, Jason B.↗

NASA NTRS1973

Viscous effects in rapidly rotating stars with application to white-dwarf models. I, II.

A general approximate numerical technique is proposed for constructing evolutionary sequences of rapidly rotating axisymmetric barytropic equilibrium configurations, with allowance for angular momentum transfer by a nonconstant isotropic viscosity. The principal physical assumption involved is the constancy of the angular momentum per unit mass on cylinders about the axis of rotation. Rapidly rotating nonmagnetic white-dwarf models with a zero-temperature degenerate-electron equation of state are considered as a particular application. The viscosity used in the analysis is that of the degenerate electrons.

Durisen, R. H.↗

NASA NTRS1987

Near-symmetric instability for general Prandtl number

An analysis of the stability of a simple baroclinic flow to perturbations whose horizontal wave vector lies in and near the vertical plane containing the density gradient has been performed. The Ekman number (E) and 'azimuthal' wavenumber (alpha) were both assumed much less than 1, and expansions about these quantities were performed for the eigenmodes and Richardson number (R). Agreement with the previous analysis of Busse and Chen (1981) that the correction to the critical R was O(alpha) for Prandtl number (P) away from unity was obtained. An expression for the correction to the critical R for arbitrary P, and an approximate expression for P = 1 were obtained as new results. The correction for P = 1 has the same sign as that for P less than 1, and is O(alpha E exp 2/3). This result compares well with the numerical results of Miller and Antar (1986) for small E.

Reynolds, Nathaniel D.↗

NASA NTRS1992

Modal element method for scattering of sound by absorbing bodies

The modal element method for acoustic scattering from 2-D body is presented. The body may be acoustically soft (absorbing) or hard (reflecting). The infinite computational region is divided into two subdomains - the bounded finite element domain, which is characterized by complicated geometry and/or variable material properties, and the surrounding unbounded homogeneous domain. The acoustic pressure field is represented approximately in the finite element domain by a finite element solution, and is represented analytically by an eigenfunction expansion in the homogeneous domain. The two solutions are coupled by the continuity of pressure and velocity across the interface between the two subdomains. Also, for hard bodies, a compact modal ring grid system is introduced for which computing requirements are drastically reduced. Analysis for 2-D scattering from solid and coated (acoustically treated) bodies is presented, and several simple numerical examples are discussed. In addition, criteria are presented for determining the number of modes to accurately resolve the scattered pressure field from a solid cylinder as a function of the frequency of the incoming wave and the radius of the cylinder.

Baumeister, Kenneth J.↗

NASA NTRS1993

Modal element method for scattering and absorbing of sound by two-dimensional bodies

The modal element method for acoustic scattering from a 2-D body is presented. The body may be acoustically soft (absorbing) or hard (reflecting). The infinite computational region is divided into two subdomains - the bounded finite element domain, which is characterized by complicated geometry and/or variable material properties, and the surrounding unbounded homogeneous domain. The acoustic pressure field is represented approximately in the finite element domain by a finite element solution, and is represented analytically by an eigenfunction expansion in the homogeneous domain. The two solutions are coupled by the continuity of pressure and velocity across the interface between the two subdomains. Also, for hard bodies, a compact modal ring grid system is introduced for which computing requirements are drastically reduced. Analysis for 2-D scattering from solid and coated (acoustically treated) bodies is presented, and several simple numerical examples are discussed. In addition, criteria are presented for determining the number of modes to accurately resolve the scattered pressure field from a solid cylinder as a function of the frequency of the incoming wave and the radius of the cylinder.

Baumeister, K. J.↗

NASA NTRSDate not supplied

Quality Control Methods of Tower Data at Kennedy Space Center’s Launch Complex 39-B and the USAF Cape Canaveral 500 Foot Tower

The National Aeronautics and Space Administration (NASA) has long used meteorological data from weather towers located at Kennedy Space Center (KSC) and the United States Air Force’s Eastern Range(ER) in support of their various launch vehicles and climate studies. Some of the most valuable data is gathered at the Launch Complex 39-B (LC39-B) and the 500 foot tower located approximately two miles west of LC39-B and three miles north of the Vehicle Assembly Building (VAB). The data from LC39-B gives NASA Space Launch System (SLS) engineers valuable insight into the weather the vehicle can expect at the pad, while the 500 foot tower has a longer period of record and has been used by the Shuttle and SLS programs. However, numerous data quality (DQ) control issues have arisen when utilizing this data for analysis. In addition, instrumentation tends to deteriorate faster than normal due to the corrosive nature of the high salt content in the air. NASA’s Natural Environments branch has developed several QC databases of these towers, but these studies were for fixed periods of records. EV44 identified a need for a continually updated and QC’ed database of tower data not only to provide along term QC’ed database for vehicle and climate analyses, but also to provide the capability to investigate recent weather events such as downbursts and other high wind events. This study follows the methods developed in prior Natural Environment QC tower databases, which includes variable specific QC thresholds and checks used to generate QC flags. However, this study developed new techniques such as the development and design of the graphical user interface (GUI) for manual verification, and the creation of the final monthly QC files.

Quality Control↗

NASA NTRSDate not supplied

Computing the effective elasticity of anisotropic porous media from X-ray computed micro-tomography images

Development and optimization of composite materials designed for thermal protection of NASA’s spacecraft requires the understanding of their physical response to high-enthalpy environments. To predict their macro-scale properties and behavior, high-fidelity 3D simulations are performed at the micro-scale on realistic representations of these composites. The digital micro-structures are generated either synthetically or through X-ray micro-tomography reconstructions. One of the main challenges in the prediction of heatshield material structural response is the computation of the effective elasticity of the fibrous composite, as well as the understanding of the deformation and stresses generated at the micro-scale. These are driven by the fiber layout within the micro-structure and the distribution of the infused matrix. In this effort, the micro-mechanical linear elastic behavior of fibrous ablators is modeled through the use of a numerical method based on the Multi-Point Stress Approximation (MPSA) finite volume scheme. The MPSA, a generalization of the more commonly used Multi-Point Flux Approximation (MPFA), was discussed in a previous presentation at the 15th USCCN and in reference. To predict the behavior of fibrous and woven architectures, algorithms that compute the local fiber orientation are used. The implementation of the MPSA was verified using analytical solutions, engineering test cases and compared against legacy Finite Element Analysis (FEA) software. The stress analysis models were then applied to real geometries used by NASA in thermal protection systems such as fibrous preforms and woven materials and the results were compared to experimental data.

Federico Semeraro↗

NASA NTRS2001

Quasi 1-D Analysis of a Circular, Compressible, Turbulent Jet Laden with Water Droplets

Recent experimental studies indicate that presence of small amount of liquid droplets reduces the Overall Sound Pressure Level (OASPL) of a jet. Present study is aimed at numerically investigating the effect of liquid particles on the overall flow quantities of a heated, compressible round jet. The jet is assumed perfectly expanded. A quasi-1D model was developed for this purpose which uses area-averaged quantities that satisfy integral conservation equations. Special attention is given to represent the early development region since it is acoustically important. Approximate velocity and temperature profiles were assumed in this region to evaluate entrainment rate. Experimental correlations were used to obtain spreading rate of shear layer. The base flow thus obtained is then laden with water droplets at the exit of the nozzle. Mass, momentum and energy coupling between the two phases is represented using empirical relations. Droplet size and mass loading are varied to observe their effect on flow variables.

Source record↗

NASA NTRS1976

Technique for atmospheric rate chemistry calculations

The possibility that predictions of atmospheric photochemistry/transport models are sensitive to uncertainties in reaction rates and other inputs stresses the need for rapid numerical integration schemes in rate photochemistry problems. Reducing the computational burden has a major merit in facilitating sensitivity studies to assess the effect of uncertainties on predicted ozone diminutions from NOx (NO + NO2) in the exhaust plume of SST engines. The paper discusses the validity of an algorithmic approach to integration of rate chemistry problems in combustion, developed by Rubel and Baronti for an approximate calculation of the production rate of the i-th chemical species involved. An analysis of two projected SST engines confirms the validity of the proposed algorithm. Because of the relative arithmetical simplicity, it may be easier to treat diffusion rate chemistry calculations using the Rubel and Baronti approximation than would be possible by other approaches.

Matloff, G. L.↗