Engineering Papers⌕ Search

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.

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

At least 433 records · Page 24

Heat generation in aircraft tires

A method was developed for calculating the internal temperature distribution in an aircraft tire while free rolling under load. The method uses an approximate stress analysis of each point in the tire as it rolls through the contact patch, and from this stress change the mechanical work done on each volume element may be obtained and converted into a heat release rate through a knowledge of material characteristics. The tire cross-section is then considered as a body with internal heat generation, and the diffusion equation is solved numerically with appropriate boundary conditions of the wheel and runway surface. Comparison with data obtained with buried thermocouples in tires shows good agreement.

Clark, S. K.↗

Preferred latitudes of the intertropical convergence zone

The latitude preference of the intertropical convergence zone (ITCZ) is examined on the basis of observations, theory, and a modeling analysis. Observations show that convection is enhanced at latitudes of about 4 deg to 10 deg relative to the equator, even in regions where the sea surface temperature (SST) is maximum on the equator. Both linear shallow-water theory and a moist primitive equation model suggest a new explanation for the off-equatorial latitude preference of the ITCZ that requires neither the existence of zonally propagating disturbances nor an off-equatorial maximum in SST. The shallow-water theory indicates that a finite-width, zonally oriented, midtropospheric heat source (i.e., an ITCZ) produces the greatest local low-level convergence when placed a finite distance away from the equator. This result suggests that an ITCZ is most likely to be supported via low-level convergence of moist energy when located at these "preferred" latitudes away from the equator. For a plausible range of heating widths and damping parameters, the theoretically predicted latitude is approximately equal to the observed position (s) of the ITCZ (s). Analysis with an axially symmetric, moist, primitive equation model indicates that when the latent heating field is allowed to be determined internally, a positive feedback develops between the midtropospheric latent heating and the low-level convergence, with the effect of enhancing the organization of convection at latitudes of about 4 deg to 12 deg. Numerical experiments show that (1) two peaks in convective precipitation develop straddling the equator when the SST maximum is located on the equator; (2) steady ITCZ-like structures form only when the SST maximum is located away from the equator; and (3) peaks in convection can develop away from the maximum in SST, with a particular preference for latitudes of about 4 deg to 12 deg, even in the ('cold') hemisphere without the SST maximum. The relationship between this mechanism and earlier theories is discussed, as are implications for the coupled ocean-atmosphere system and the roles played by midlevel latent heating and SST gradients in forcing the low-level atmospheric circulation in the tropics.

Waliser, Duane E.↗

Impact ejecta on the moon

The response of a lunar-sized object to the impact of meteoroids no more than about 100 km in radius is studied by means of a numerical model. The partitioning of impact energy into the kinetic and internal energy of the ejecta is obtained by using the conservation of mass, momentum, and energy conservation equations in finite-difference form within an Eulerian framework with approximate equations of state. The calculations are performed for a 15 km/sec impact of an iron object 5 cm in radius on a gabbroic anorthosite surface. Ejecta ballistic analysis is then performed. Most of the material lost escaping the moon is lunar crust material. Only 0.2% of the meteoroid escapes, all in the vapor phase.

Okeefe, J. D.↗

Numerical simulation of natural convection in a spherical container due to cooling at the center (idealization of the Lal/Kroes experiment)

Natural convection in a spherical container with cooling at the center was numerically simulated using a numerical fluid dynamics computer program. The numerical analysis was simplified by assuming axisymmetric flow in the spherical container, with the symmetry axis being a sphere diagonal parallel to the gravity vector. This axisymmetric spherical geometry was intended as an idealization of the proposed Lal/Kroes crystal growing experiment to be performed on Spacelab. Results were obtained for a range of Rayleigh numbers from 25 to 10,000. The computed velocities were found to be approximately proportional to the Rayleigh number over the range of Rayleigh numbers investigated.

Robertson, S. J.↗

An evaluation of higher-order model methods for calculating transient structural response

A higher-order modal method proposed by Leung for transient structural analysis entitled the force-derivative method is evaluated. This method repeatedly integrates by parts with respect to time the convolution-integral form of the structural response to produce successively better approximations to the contribution of the higher modes which are neglected in the modal summation. Comparisons are made of the force-derivative, the mode-displacement, and the mode-acceleration methods for several numerical example problems for various times, levels of damping, and forcing functions. The example problems include a tip-loaded cantilevered beam and a simply-supported multispan beam. The force-derivative method is shown to converge to an accurate solution in fewer modes than either the mode-displacement or the mode-acceleration methods. In addition, for problems in which there are a large number of closely-spaced frequencies whose mode shapes have a negligible contribution to the response, the force-derivative method is very effective in representing the effect of the important, but otherwise neglected, higher modes.

Camarda, Charles J.↗

Probabilistic structural analysis algorithm development for computational efficiency

The PSAM (Probabilistic Structural Analysis Methods) program is developing a probabilistic structural risk assessment capability for the SSME components. An advanced probabilistic structural analysis software system, NESSUS (Numerical Evaluation of Stochastic Structures Under Stress), is being developed as part of the PSAM effort to accurately simulate stochastic structures operating under severe random loading conditions. One of the challenges in developing the NESSUS system is the development of the probabilistic algorithms that provide both efficiency and accuracy. The main probability algorithms developed and implemented in the NESSUS system are efficient, but approximate in nature. In the last six years, the algorithms have improved very significantly.

Wu, Y.-T.↗

The Evolution of Disks and Winds in Dwarf Nova Outbursts - FUSE

This project was a project to study the FUV spectra of two proto-typical dwarf novae, U Gem and SS Cygni, through an outburst cycle. The luminosity of the boundary layer in the two systems, as evidenced by earlier EUVE observations, is different in the two systems. Our intensive study of the two systems was intended to (1) probe the ionization and kinematic structure of the wind as a function of system brightness, (2) isolate the contributions of the disk to the FUV spectra, and (3) examine physical conditions and abundances of material just being accreted onto the disk from the secondary. The U Gem and SS Cyg observations took place in March and October 2000, respectively. The data obtained with FUSE was of excellent quality. Analysis of the both observations is now essentially complete, although some modeling of the SS Cyg spectra is ongoing, as we complete an ApJ manuscript on this object. Our main results for U Gem are as follows: The plateau spectra have continuum shapes and fluxes that are approximated by steady state accretion disk model spectra with an accretion rate 7x10(exp 9) Msolar/yr. The spectra also show numerous absorption lines of H I, He II, and 2-5 times ionized transitions of C, N, O, P, S, and Si. There are no emission features in the spectra, with the possible exception of a weak feature on the red wing of the 0 VI doublet. The absorption lines are narrow (FWHM approx. 50 km/s), too narrow to arise from the disk photosphere, and at low velocities (less than or equal to 700 km/s). The S VI and O VI doublets are optically thick. The absorption lines in the plateau spectra show orbital variability: in spectra obtained at orbital phases between 0.53 and 0.79, low-ionization absorption lines appear and the central depths of the preexisting lines increase. The increase in line absorption occurs at the same orbital phases as previously observed EUV and X-ray light-curve dips. If the absorbing material is in (near-) Keplerian rotation around the disk, it must be located at large disk radii. The final observation occurred when U Gem was about 2 mag from optical quiescence. The spectra are dominated by emission from an approx. 43,000 K, metal-enriched white dwarf (WD). The inferred radius of the WD is 4.95x10(exp 8) cm, close to that observed in quiescence. Allowing for a hot heated region on the surface of the WD improves the fit to the spectrum at short (less than 960 A) wavelengths. Our main results for SS Cyg are as follows: The first two of four observations of SS Cyg show disk dominated spectra with accretion rates of order 10(exp -8) Msolar/yr. Except for narrow interstellar features (atomic and molecular H), the lines are all broad consistent with a disk or wind origin. The O VI line in the spectra is mostly of wind origin as detailed modeling with our Monte Carlo code (developed in part using funds from this project) show. The continua from spectra in observations 3 and 4, observed during the decline phase, are not well fit with steady-state disks, and show considerable resemblance to quiescent spectra obtained with HUT. The most probable interpretation for the emission features seen in the spectrum in the last two observations is that they arise from a photo-illuminated choronosphere above the disk, rather than a wind.

Long, Knox↗

Single Degree-of-Freedom Modeling of SLS Liquid Hydrogen Pre-Valve Flow Guide to Enable Rapid Transient Analysis

A unique single degree-of-freedom approximation technique has been developed to enable rapid application of a temporally-defined multi-spectral semi-narrow-band loading for generation of realistic stress/cycle values compared to a resonant analysis. The technique uses the harmonic analysis at resonance of a high-fidelity finite element model to produce a transfer function, which is then used to calibrate the response of the SDOF model. A standard numerical ordinary differential equation solver is then used to obtain the temporal response, and its histogram is used in a fatigue/fracture model. This technique is related to other SDOF methods used widely in industry, such as Miles' Equation and the Shock Response Spectra, but it is unique in that it produces a realistic time history of the response. The most obvious error in the process, which is the effect of closely-spaced modes, was also assessed using the parallel application of several SDOF models, and the error is shown to be small. The application of this unique and tractable reduced-order methodology has enabled the SLS program to avoid substantial cost and schedule penalties if a redesign or change of material were required. It has also enabled quick analysis of a number of other structures undergoing the same or similar excitation fields, and quick assessment when the excitation and structural configuration has been altered due to design changes in the system.

Brown, Andrew M.↗

Propulsion-Airframe Integration for Conceptual Redesign of a Low-Boom Supersonic Transport

A low-boom supersonic transport was generated in a previous low-boom multidisciplinary optimization (MDO) study that used computational fluid dynamics (CFD) off-body pressure to compute the undertrack sonic boom ground signature and calibrated low-fidelity aerodynamics analyses to compute the mission performance metrics. This low-boom aircraft, referred to as the Mach 1.7 40-PAX concept, can carry 40 passengers for a low-boom overland mission with cruise Mach 1.7 for the airport pairs over the continental US and has the potential to achieve an undertrack sonic boom ground noise level below 70 perceived level of decibels at the start of overland cruise (SOC). The engine for the Mach 1.7 40-PAX concept was designed using the Numerical Propulsion System Simulation (NPSS) for the mission analysis and modeled as a flow-through nacelle for CFD-based sonic boom analysis. To understand how the engine plume affects the undertrack ground signature, the Mach 1.7 40-PAX concept is redesigned in this paper after replacing the flow-through nacelles by CFD engines for sonic boom analysis using aeropropulsive CFD simulation. A process is developed for approximation of the NPSS engine at SOC by a CFD engine for aeropropulsive CFD simulation. The generated CFD engine has the identical nozzle boundary conditions and approximately the same mass flow and thrust as those of the NPSS engine at SOC. Then, the outer mold line of the configuration with the CFD engines is redesigned to approximately restore the low-boom characteristics of the Mach 1.7 40-PAX concept at SOC. Finally, the CFD simulation data for the redesigned concept with the CFD engines is used to calibrate the low-fidelity aerodynamic analyses for mission analysis of this concept. The cyclic dependency of the involved disciplinary analyses is resolved using an iteration method for a consistent coupling of the mission analysis, NPSS engine analysis, and low-boom redesign using the aeropropulsive CFD simulation. This low-boom redesign study is used as an example to demonstrate how the propulsion-airframe integration could be implemented for conceptual design of supersonic transports that satisfy both the low-boom and mission performance requirements.

Low-boom supersonic transport↗

Dynamics analysis of electrodynamic satellite tethers. Equations of motion and numerical solution algorithms for the tether

The equations of motion are developed for a perfectly flexible, inelastic tether with a satellite at its extremity. The tether is attached to a space vehicle in orbit. The tether is allowed to possess electrical conductivity. A numerical solution algorithm to provide the motion of the tether and satellite system is presented. The resulting differential equations can be solved by various existing standard numerical integration computer programs. The resulting differential equations allow the introduction of approximations that can lead to analytical, approximate general solutions. The differential equations allow more dynamical insight of the motion.

Nacozy, P. E.↗

Continuation of tailored composite structures of ordered staple thermoplastic material

The search for the cost effective composite structure has motivated the investigation of several approaches to develop composite structure from innovative material forms. Among the promising approaches is the conversion of a planar sheet to components of complex curvature through sheet forming or stretch forming. In both cases, the potential for material stretch in the fiber direction appears to offer a clear advantage in formability over continuous fiber systems. A framework was established which allows the simulation of the anisotropic mechanisms of deformation of long discontinuous fiber laminates wherein the matrix phase is a viscous fluid. Predictions for the effective viscosities of a hyper-anisotropic medium consisting of collimated, discontinuous fibers suspended in viscous matrix were extended to capture the characteristics of typical polymers including non-Newtonian behavior and temperature dependence. In addition, the influence of fiber misorientation was also modeled by compliance averaging to determine ensemble properties for a given orientation distribution. A design tool is presented for predicting the effect of material heterogeneity on the performance of curved composite beams such as those used in aircraft fuselage structures. Material heterogeneity can be induced during manufacturing processes such as sheet forming and stretch forming of thermoplastic composites. This heterogeneity can be introduced in the form of fiber realignment and spreading during the manufacturing process causing radial and tangential gradients in material properties. Two analysis procedures are used to solve the beam problems. The first method uses separate two-dimensional elasticity solutions for the stresses in the flange and web sections of the beam. The separate solutions are coupled by requiring that forces and displacements match section boundaries. The second method uses an approximate Rayleigh-Ritz technique to find the solutions for more complex beams. Analyses are performed for curved beams of various cross-sections loaded in pure bending and with a uniform distributed load. Preliminary results show that the geometry of the beam dictates the effect of heterogeneity on performance. The role of heterogeneity is larger in beams with a small average radius-to-depth ration, R/t, where R is the average radius of the beam and t is the difference between the inside and outside radii. Results of the anlysis are in the form of stresses and displacements and are compared to both mechanics of materials and numerical solutions obtained using finite element analysis.

Santare, Michael H.↗

Transition Delay via Vortex Generators in a Hypersonic Boundary Layer at Flight Conditions

The potential of realizable, stationary streaks undergoing non-modal growth to stabilize a hypersonic boundary-layer flow and, subsequently, delay the laminar-turbulent transition onset, is studied via numerical computations. The geometry and flow conditions are selected to match a relevant trajectory location from the ascent phase of the HIFiRE-1 flight experiment, namely, a 7-degree half-angle cone with 2.5 mm nose radius, freestream Mach number of 5.30, freestream unit Reynolds number equal to 13.42 x 10(exp 6)/m, and wall-to-adiabatic temperature ratio of approximately 0.35 over most of the test article. This paper investigates flow modifications induced by wall-mounted vortex generators (VGs), followed by an analysis of the modal instability of the perturbed, streaky boundary-layer flow. Results are presented both for a single array of VGs that was designed on the basis of optimal growth theory and for a VG configuration involving two separate arrays with opposite orientations that ware designed to provide staged control of flow instabilities while simultaneously reducing the amplification of streak instabilities resulting from the control devices. Earlier research had shown that the onset of transition during the HIFiRE-1 flight experiment, which did not include any control devices, correlated with an amplification factor of N = 14.7 for the planar Mack modes. If one assumes that the transition N -factor is not affected by the introduction of the VGs, then the control configurations based on a single array of VGs and two separate arrays would result in a transition delay of 17% and 40%, respectively. These findings suggest a passive flow control s to induce streaks that would delay transition in hypersonic boundary dominated by Mack-mode instabilities.

Incompressible boundary layer↗

Design, fabrication, and evaluation of charge-coupled devices with aluminum-anodized-aluminum gates

A 4-phase, 49 1/2 bit CCD shift register was designed and fabricated using two levels of aluminum metallization with anodic Al2O3 insulation separating the layers. Test circuitry was also designed and constructed. A numerical analysis of an MOS-RC transmission line was made and results are given to characterize performance for various conductivities. The electrical design of the CCD included a low-noise dual-gate input and a balanced floating diffusion output circuit. Metallization was accomplished both by low voltage DC sputtering and thermal evaporation. The audization was according to published procedures using a buffered tartaric acid bath. Approximately 20 wafers were processed with 50 complete chips per wafer. All devices failed by shorting between the metal levels at some point. Experimental procedures eliminated temperature effects from sintering and drying, anodic oxide thickness, edge effects, photoresist stripping procedures, and metallization techniques as the primary causes of failure. It was believed from a study of SEM images that protuberances (hillocks) grow up from the first level metal through the oxide either causing a direct short or producing a weak, highly stressed insulation point which fails at low voltage. The cause of these hillocks is unknown; however, they have been observed to grow during temperature excursions to 470 C.

Gassaway, J. D.↗

Calculating Effects of Reflector-Antenna Distortions

Mathematical numerical model accurately predicts effects of systematic distortions on radiation pattern of reflector antenna. Techniques for predicting radiation patterns developed previously for reflectors with perfect surfaces and surfaces with random deviations from perfection. New model treats reflectors distorted by thermal, gravitational, and dynamic effects. Diffraction analysis of reflector antenna depends on interpolation of surface description for farfield radiation integral. Local interpolation scheme uses bicubic or higher order polynomial approximation of surface to express desired surface point in terms of nearby known surface points.

Rahmat-Samii, Y.↗

Postbuckling of laminated anisotropic panels

A two-part study of the buckling and postbuckling of laminated anisotropic plates with bending-extensional coupling is presented. The first part involves the development and application of a modified Rayleigh-Ritz analysis technique. Modifications made to the classical technique can be grouped into three areas. First, known symmetries of anisotropic panels are exploited in the selection of approximation functions. Second, a reduced basis technique based on these same symmetries is applied in the linear range. Finally, geometric boundary conditions are enforced via an exterior penalty function approach, rather than relying on choice of approximation functions to satisfy these boundary conditions. Numerical results are presented for both the linear and nonlinear range, with additional studies made to determine the effect of variation in penalty parameter and number of basis vectors. In the second part, six panels possessing anisotropy and bending-extensional coupling are tested. Detailed comparisons are made between experiment and finite element results in order to gain insight into the postbuckling and failure characteristics of such panels. The panels are constructed using two different lamination sequences, and panels with three different aspect ratios were constructed for each lamination sequence.

Jeffrey, Glenda L.↗

User's guide for NETS/PROSSS

Expensive analysis programs are often combined with optimization procedures to solve engineering problems. To obtain an optimal solution requires numerous iterations between the analysis program and the optimizer. This often becomes prohibitive due to the cost and amount of computer time needed to converge to an optimal solution. NETS/PROSSS was developed to address this problem. The purpose of this paper is to serve as a user's guide for NETS/PROSSS. The key features include the neural network, determining the training pairs for the neural network, and the approximated analysis/optimization process. A small problem is given to serve as an example of how to apply the system.

Rogers, James L.↗

A determination of the external forces required to move the benchmark active controls testing model in pure plunge and pure pitch

In view of the strong need for a well-documented set of experimental data which is suitable for the validation and/or calibration of modern Computational Fluid Dynamics codes, the Benchmark Models Program was initiated by the Structural Dynamics Division of the NASA Langley Research Center. One of the models in the program, the Benchmark Active Controls Testing Model, consists of a rigid wing of rectangular planform with a NACA 0012 profile and three control surfaces (a trailing-edge control surface, a lower-surface spoiler, and an upper-surface spoiler). The model is affixed to a flexible mount system which allows only plunging and/or pitching motion. An approximate analytical determination of the forces required to move this model, with its control surfaces fixed, in pure plunge and pure pitch at a number of test conditions is included. This provides a good indication of the type of actuator system required to generate the aerodynamic data resulting from pure plunging and pure pitching motion, in which much interest was expressed. The analysis makes use of previously obtained numerical results.

Dcruz, Jonathan↗

H2SO4 photolysis: A souce of sulfur dioxide in the upper stratosphere

Numerous absoption lines of stratospheric sulfer dioxide (SO2) were identified in solar occulation spectra recorded by the Atmospheric Trace Molecule Spectroscopy (ATMOS) Fourier transform spectrometer during the Atmospheric Laboratory for Applications and Science (ATLAS)-1 shuttle mission (March 24-April 2, 1992). based on their analysis, a volume mixing ratio profile of SO2 increasing from (13 +/- 4) p.p.t.v. (parts per 10(exp -12) by volume) at 16 mbar (approximately 28 km) to 455 +/- 90 p.p.t.v. at 0.63 mbar (approximately 52 km) was measured with no significant profile differences between 20 deg N and 60 deg S latitude. The increase in the SO2 mixing ratios with altitude indicates the presence of a source of SO2 in the upper stratosphere. Profiles retrieved from ATMOS spectra recorded during shuttle flights in April-May 1985 and April 1993 show similar vertical distributions but lower concentrations. Two-dimensional model calculations with SO2 assumed as the end product of H2SO4 photolysis produce SO2 profiles consistent with the ATMOS measuremnts to within about a factor 2.

Rinsland, C. P.↗