Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computational Physics”

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 721 records · Page 40

On the numerical solution of the dynamically loaded hydrodynamic lubrication of the point contact problem

The transient analysis of hydrodynamic lubrication of a point-contact is presented. A body-fitted coordinate system is introduced to transform the physical domain to a rectangular computational domain, enabling the use of the Newton-Raphson method for determining pressures and locating the cavitation boundary, where the Reynolds boundary condition is specified. In order to obtain the transient solution, an explicit Euler method is used to effect a time march. The transient dynamic load is a sinusoidal function of time with frequency, fractional loading, and mean load as parameters. Results include the variation of the minimum film thickness and phase-lag with time as functions of excitation frequency. The results are compared with the analytic solution to the transient step bearing problem with the same dynamic loading function. The similarities of the results suggest an approximate model of the point contact minimum film thickness solution.

Lim, Sang G.↗

Inflection point caustic problems and solutions for high-gain dual-shaped reflectors

The singular nature of the uniform geometrical theory of diffraction (UTD) subreflector scattered field at the vicinity of the main reflector edge (for a high-gain antenna design) is investigated. It is shown that the singularity in the UTD edge-diffracted and slope-diffracted fields is due to the reflection distance parameter approaching infinity in the transition functions. While the geometrical optics (GO) and UTD edge-diffracted fields exhibit singularities of the same order, the edge slope-diffracted field singularity is more significant and is substantial for greater subreflector edge tapers. The diffraction analysis of such a subreflector in the vicinity of the main reflector edge has been carried out efficiently and accurately by a stationary phase evaluation of the phi-integral, whereas the theta-integral is carried out numerically. Computational results from UTD and physical optics (PO) analysis of a 34-m ground station dual-shaped reflector confirm the analytical formulations for both circularly symmetric and offset asymmetric subreflectors. It is concluded that the proposed PO(theta)GO(phi) technique can be used to study the spillover or noise temperature characteristics of a high-gain reflector antenna efficiently and accurately.

Galindo-Israel, Victor↗

Models of Wolf-Rayet stars

The current status of knowledge about formation, structure and evolution of Wolf-Rayet stars is reviewed, with emphasis on a discussion of corresponding stellar models. The relevance of the LBV-scenario for WR star formation is outlined. Hydrogenless WR stars are shown to closely follow simple relations for the dependence of luminosity, radius, and surface temperature as a function of their mass. The use of these relations for simplified WR evolution calculations is demonstrated. Surface abundance predictions for the different WR types are discussed, with special emphasis to the WN + WC spectral type. Details are presented concerning the WR phase of a recent 60 solar mass evolutionary calculation, which was computed with the same input physics which reproduced the progenitor evolution of SN 1987 A in a 20 solar mass case, and which may be a representative case concerning WR stars in many respects.

Langer, Norbert↗

Potential of mixed formulations for advanced analysis systems

Two recent applications of the combined use of mixed formulations with splitting methodologies and reduction methods are presented. Discussion focuses on how the combination can significantly improve the efficiency of the computational process and enhance the physical understanding of the response. Numerical examples are presented which demonstrate the unique features and potential of mixed formulations in advanced analysis systems.

Noor, A. K.↗

Heat pipe dynamic behavior

The vapor flow in a heat pipe was mathematically modeled and the equations governing the transient behavior of the core were solved numerically. The modeled vapor flow is transient, axisymmetric (or two-dimensional) compressible viscous flow in a closed chamber. The two methods of solution are described. The more promising method failed (a mixed Galerkin finite difference method) whereas a more common finite difference method was successful. Preliminary results are presented showing that multi-dimensional flows need to be treated. A model of the liquid phase of a high temperature heat pipe was developed. The model is intended to be coupled to a vapor phase model for the complete solution of the heat pipe problem. The mathematical equations are formulated consistent with physical processes while allowing a computationally efficient solution. The model simulates time dependent characteristics of concern to the liquid phase including input phase change, output heat fluxes, liquid temperatures, container temperatures, liquid velocities, and liquid pressure. Preliminary results were obtained for two heat pipe startup cases. The heat pipe studied used lithium as the working fluid and an annular wick configuration. Recommendations for implementation based on the results obtained are presented. Experimental studies were initiated using a rectangular heat pipe. Both twin beam laser holography and laser Doppler anemometry were investigated. Preliminary experiments were completed and results are reported.

Issacci, F.↗

Towards a methodology for robust parameter identification

Consideration is given to the problem of estimating, from experimental data, real parameters for a model with uncertainty in the form of both additive noise and norm-bounded perturbations. Such models frequently arise in robust control theory, and a framework is introduced for the consideration of experimental data in robust control analysis problems. If the analysis tools applied include robust stability tests for real parameter variations (real mu), the framework can be used to address the problem of robust parameter identification. While the techniques discussed can quickly become computationally overwhelming when applied to physical systems and real data, the approach introduces a novel way of looking at the identification problem and may be helpful in arriving at a more tractable methodology.

Smith, Roy S.↗

The 3D rocket combustor acoustics model

The theory and procedures for determining the characteristics of pressure oscillations in rocket engines with prescribed burning rate oscillations are presented. Analyses including radial and hub baffles and absorbers can be performed in one, two, and three dimensions. Pressure and velocity oscillations calculated using this procedure are presented for the SSME to show the influence of baffles and absorbers on the burning rate oscillations required to achieve neutral stability. Comparisons are made between the results obtained utilizing 1-D, 2-D, and 3-D assumptions with regards to capturing the physical phenomena of interest and computational requirements.

Priem, Richard J.↗

Diagnostic evaluation of numerical model simulations using the tendency equation

A procedure is proposed to expand the diagnostic capabilities of the pressure tendency equation of a primitive equation NWP model by computing the pressure tendency in physical coordinates. The advantages of isolating the density advection as a diagnostic tool to understand pressure changes is shown. By simple thermodynamic arguments it is demonstrated that in areas of synoptic-scale cyclonic development, the vertically integrated density advection is more than sufficient to explain the depletion of mass over a growing depression. Consequently, the joint contribution of the net divergence and vertical motion opposes the pressure fall. This is illustrated for a case of rapid cyclogenesis in southern South America.

Jusem, Juan C.↗

3D rocket combustor acoustics model

The theory and procedures for determining the characteristics of pressure oscillations in rocket engines with prescribed burning rate oscillations are presented. Analyses including radial and hub baffles and absorbers can be performed in one, two, or three dimensions. Pressure and velocity oscillations calculated using this procedure are presented for the SSME to show the influence of baffles and absorbers on the burning rate oscillations required to achieve neutral stability. Comparisons are made between the results obtained utilizing 1D, 2D, and 3D assumption with regards to capturing the physical phenomena of interest and computational requirements.

Priem, Richard J.↗

A method for modeling contact dynamics for automated capture mechanisms

Logicon Control Dynamics develops contact dynamics models for space-based docking and berthing vehicles. The models compute contact forces for the physical contact between mating capture mechanism surfaces. Realistic simulation requires proportionality constants, for calculating contact forces, to approximate surface stiffness of contacting bodies. Proportionality for rigid metallic bodies becomes quite large. Small penetrations of surface boundaries can produce large contact forces.

Williams, Philip J.↗

Computation of unsteady shock-induced combustion using logarithmic species conservation equations

Numerical simulations are used to investigate periodic combustion instabilities observed in ballistic-range experiments of blunt bodies flying at supersonic speeds through hydrogen-air mixtures. The computations are validated by comparing experimental shadowgraphs with shadowgraphs created from the computed flowfields and by comparing the experimentally measured instability frequencies with computed frequencies. The numerical simulations use a logarithmic transformation of the species conservation equations as a way to reduce the grid requirements for computing shock-induced combustion. The transformation is applied to the Euler equations coupled to a detailed hydrogen-air chemical reaction mechanism with 13 species and 33 reactions. The resulting differential equations are solved using a finite volume formulation and a two-step predictor-corrector scheme to advance the solution in time. Results are presented and compared for both a flux-vector splitting scheme and an upwind TVD scheme. The computations add insight to the physical processes observed in the experiments and the numerical methods needed to simulate them. The usefulness of the ballistic-range experiments for the validation of numerical techniques and chemical kinetic models is also demonstrated.

Wilson, Gregory J.↗

Clouding tracing: Visualization of the mixing of fluid elements in convection-diffusion systems

This paper describes a highly interactive method for computer visualization of the basic physical process of dispersion and mixing of fluid elements in convection-diffusion systems. It is based on transforming the vector field from a traditionally Eulerian reference frame into a Lagrangian reference frame. Fluid elements are traced through the vector field for the mean path as well as the statistical dispersion of the fluid elements about the mean position by using added scalar information about the root mean square value of the vector field and its Lagrangian time scale. In this way, clouds of fluid elements are traced and are not just mean paths. We have used this method to visualize the simulation of an industrial incinerator to help identify mechanisms for poor mixing.

Ma, Kwan-Liu↗

Application of the target decomposition theorem to a polarimetric random media model

With advances in polarimetric radar measurements of land surfaces, the need for understanding the underlying scattering mechanisms and dominant target features has become the focus of many studies. In particular, the maximum use of the polarimetric information to identify and/or separate parameters related to the surface features such as vegetation thickness, structure, water content, and soil surface characteristics will enhance the possibility of using polarimetric radars for monitoring the earth's surface from space. In this paper, Cloude's decomposition theorem is applied to a polarimetric random media model to simulate the radar measurements of vegetated canopies. The vegetated canopies are modeled as a three layer discrete random medium with leaves and branches in the first layer, tree trunks in the second layer, and a half space of homogeneous ground with rough interface as the bottom layer. The distorted born approximation (DBA) has been used to compute full Mueller matrix of the canopy, using canonical dielectric objects such as thin discs and cylinders as leaves, branches, and trunks, respectively. The Mueller matrix and the derived covariance matrix contain information on the second order statistics of radar signals at various polarizations from the canopy. To decompose the covariance matrix to its constituent targets, the eigenvalues and eigenvectors of the covariance matrix are computed in terms of the physical parameters of the canopy. In addition, each eigenvector explicitly shows the scattering mechanisms such as odd and even reflections in the canopy. Cloude's decomposition theorem is applied using the Pauli spin matrices as a basis and an expression for the degree of disorder or the entropy for the vegetated surface is found. Then, the physical parameters estimated from in situ measurements are used in the random media to obtain realistic covariance matrices. As a result, the sensitivity of the eigenvalue spectrums and the coefficients resulting from the target decomposition theorem to the physical parameters of the canopy are examined and the possible use of Cloude's theorem to estimate vegetation and soil parameters is discussed.

Saatchi, Sasan S.↗

Radiation and scattering from cylindrically conformal printed antennas

Microstrip patch antennas offer considerable advantages in terms of weight, aerodynamic drag, cost, flexibility, and observables over more conventional protruding antennas. These flat patch antennas were first proposed over thirty years ago by Deschamps in the United States and Gutton and Baisinot in France. Such antennas have been analyzed and developed for planar as well as curved platforms. However, the methods used in these designs employ gross approximations, suffer from extreme computational burden, or require expensive physical experiments. The goal of this thesis is to develop accurate and efficient numerical modeling techniques which represent actual antenna structures mounted on curved surfaces with a high degree of fidelity. In this thesis, the finite element method is extended to cavity-backed conformal antenna arrays embedded in a circular, metallic, infinite cylinder. Both the boundary integral and absorbing boundary mesh closure conditions will be used for terminating the mesh. These two approaches will be contrasted and used to study the scattering and radiation behavior of several useful antenna configurations. An important feature of this study will be to examine the effect of curvature and cavity size on the scattering and radiation properties of wraparound conformal antenna arrays.

Kempel, Leo C.↗

Quantitative reflectance spectroscopy of buddingtonite from the Cuprite mining district, Nevada

Buddingtonite, an ammonium-bearing feldspar diagnostic of volcanic-hosted alteration, can be identified and, in some cases, quantitatively measured using short-wave infrared (SWIR) reflectance spectroscopy. In this study over 200 samples from Cuprite, Nevada, were evaluated by X ray diffraction, chemical analysis, scanning electron microscopy, and SWIR reflectance spectroscopy with the objective of developing a quantitative remote-sensing technique for rapid determination of the amount of ammonium or buddingtonite present, and its distribution across the site. Based upon the Hapke theory of radiative transfer from particulate surfaces, spectra from quantitative, physical mixtures were compared with computed mixture spectra. We hypothesized that the concentration of ammonium in each sample is related to the size and shape of the ammonium absorption bands and tested this hypothesis for samples of relatively pure buddingtonite. We found that the band depth of the 2.12-micron NH4 feature is linearly related to the NH4 concentration for the Cuprite buddingtonite, and that the relationship is approximately exponential for a larger range of NH4 concentrations. Associated minerals such as smectite and jarosite suppress the depth of the 2.12-micron NH4 absorption band. Quantitative reflectance spectroscopy is possible when the effects of these associated minerals are also considered.

Felzer, Benjamin↗

X-ray induced ionization in the winds of near-main-sequence O and B stars

We investigate the effects of X-rays on the ionization state of stellar winds for O and early-B stars along the main sequence. In our investigation, detailed statistical equilibrium, radiative transfer, and atomic physics models are used to compute ionization distributions for H, He, C, N, O, and Si. X-rays are modeled as a spatially distributed source within the wind, with parameters constrained by ROSAT and Einstein observations. Our results indicate that the ionization balance in the winds of stars with spectral type B0 V and later is significantly altered by the X-ray radiation field. Unlike the case of denser O star winds, where the X-rays tend to perturb the level of ionization, the ionization state of the bulk wind of early-B stars can be significantly increased by soft X-rays. We examine in detail the case of tau Sco (B0 V), which has been well-studied at UV and X-ray wavelengths. Comparisons are made between calculated ionization fractions and those deduced from UV observations. In addition, we address the sensitivity of our results to the X-ray source characteristics, the wind temperature, and the photospheric extreme ultraviolet (EUV) flux. Our results suggest the possibility that for early-B stars X-rays play a critical role in both influencing the radiation line driving force, as well as ionizing and heating the wind all the way down to the top of the photosphere.

Macfarlane, J. J.↗

Efficient GO2/GH2 Injector Design: A NASA, Industry and University Cooperative Effort

Developing new propulsion components in the face of shrinking budgets presents a significant challenge. The technical, schedule and funding issues common to any design/development program are complicated by the ramifications of the continuing decrease in funding for the aerospace industry. As a result, new working arrangements are evolving in the rocket industry. This paper documents a successful NASA, industry, and university cooperative effort to design efficient high performance GO2/GH2 rocket injector elements in the current budget environment. The NASA Reusable Launch Vehicle (RLV) Program initially consisted of three vehicle/engine concepts targeted at achieving single stage to orbit. One of the Rocketdyne propulsion concepts, the RS 2100 engine, used a full-flow staged-combustion cycle. Therefore, the RS 2100 main injector would combust GO2/GH 2 propellants. Early in the design phase, but after budget levels and contractual arrangements had been set the limitations of the current gas/gas injector database were identified. Most of the relevant information was at least twenty years old. Designing high performance injectors to meet the RS 2100 requirements would require the database to be updated and significantly enhanced. However, there was no funding available to address the need for more data. NASA proposed a teaming arrangement to acquire the updated information without additional funds from the RLV Program. A determination of the types and amounts of data needed was made along with test facilities with capabilities to meet the data requirements, budget constraints, and schedule. After several iterations a program was finalized and a team established to satisfy the program goals. The Gas/Gas Injector Technology (GGIT) Program had the overall goal of increasing the ability of the rocket engine community to design efficient high-performance, durable gas/gas injectors relevant to RLV requirements. First, the program would provide Rocketdyne with data on preliminary gas/gas injector designs which would enable discrimination among candidate injector designs. Secondly, the program would enhance the national gas/gas database by obtaining high-quality data that increases the understanding of gas/gas injector physics and is suitable for computational fluid dynamics (CFD) code validation. The third program objective was to validate CFD codes for future gas/gas injector design in the RLV program.

Tucker, P. K.↗

Photochemical Phenomenology Model for the New Millennium

The "Photochemical Phenomenology Model for the New Millennium" project tackles the issue of reengineering and extension of validated physics-based modeling capabilities ("legacy" computer codes) to application-oriented software for use in science and science-support activities. While the design and architecture layouts are in terms of general particle distributions involved in scattering, impact, and reactive interactions, initial Photochemical Phenomenology Modeling Tool (PPMT) implementations are aimed at construction and evaluation of photochemical transport models with rapid execution for use in remote sensing data analysis activities in distributed systems. Current focus is on the Composite Infrared Spectrometer (CIRS) data acquired during the CASSINI flyby of Jupiter. Overall, the project has stayed on the development track outlined in the Year 1 annual report and most Year 2 goals have been met. The issues that have required the most attention are: implementation of the core photochemistry algorithms; implementation of a functional Java Graphical User Interface; completion of a functional CORBA Component Model framework; and assessment of performance issues. Specific accomplishments and the difficulties encountered are summarized in this report. Work to be carried out in the next year center on: completion of testing of the initial operational implementation; its application to analysis of the CASSINI/CIRS Jovian flyby data; extension of the PPMT to incorporate additional phenomenology algorithms; and delivery of a mature operational implementation.

Bishop, James↗