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 1,027 records · Page 57

Aeroelastic Modeling of a Nozzle Startup Transient

Lateral nozzle forces are known to cause severe structural damage to any new rocket engine in development during test. While three-dimensional, transient, turbulent, chemically reacting computational fluid dynamics methodology has been demonstrated to capture major side load physics with rigid nozzles, hot-fire tests often show nozzle structure deformation during major side load events, leading to structural damages if structural strengthening measures were not taken. The modeling picture is incomplete without the capability to address the two-way responses between the structure and fluid. The objective of this study is to develop a tightly coupled aeroelastic modeling algorithm by implementing the necessary structural dynamics component into an anchored computational fluid dynamics methodology. The computational fluid dynamics component is based on an unstructured-grid, pressure-based computational fluid dynamics formulation, while the computational structural dynamics component is developed under the framework of modal analysis. Transient aeroelastic nozzle startup analyses at sea level were performed, and the computed transient nozzle fluid-structure interaction physics presented,

Wang, Ten-See↗

Parasitic current collection by PASP Plus solar arrays

Solar cells at potentials positive with respect to a surrounding plasma collect electrons. Current is collected by the exposed high voltage surfaces: the interconnects and the sides of the solar cells. This current is a drain on the array power that can be significant for high-power arrays. In addition, this current influences the current balance that determines the floating potential of the spacecraft. One of the objectives of the Air Force (PL/GPS) PASP Plus (Photovoltaic Array Space Power Plus Diagnostics) experiment is an improved understanding fo parasitic current collection. We have done computer modeling of parasitic current collection and have examined current collection flight data from the first year of operations. Prior to the flight we did computer modeling to improve our understanding of the physical processes that control parasitic current collection. At high potentials, the current rapidly rises due to a phenomenon called snapover. Under snapover conditions, the equilibrium potential distribution across the dielectric surface is such that part of the area is at potentials greater than the first crossover of the secondary yield curve. Therefore, each incident electron generates more than one secondary electron. The net effect is that the high potential area and the collecting area increase. We did two-dimensional calculations for the various geometries to be flown. The calculations span the space of anticipated plasma conditions, applied potential, and material parameters. We used the calculations and early flight data to develop an analytic formula for the dependence of the current on the primary problem variables. The analytic formula was incorporated into the EPSAT computer code. EPSAT allows us to easily extend the results to other conditions. PASP Plus is the principal experiment integrated onto the Advanced Photovoltaic and Electronics Experiments (APEX) satellite bus. The experiment is testing twelve different solar array designs. Parasitic current collection is being measured for eight of the designs under various operational and environment conditions. We examined the current collected as a function of the various parameters for the six non-concentrator designs. The results are similar to those obtained in previous experiments and predicted by the calculations. We are using the flight data to validate the analytic formula developed. The formula can be used to quantify the parasitic current collected. Anticipating the parasitic current value allows the spacecraft designer to include this interaction when developing the design.

Davis, Victoria Ann↗

High speed inlet calculations with real gas effects

A 2-D steady-state Navier-Stokes solver has been upgraded to include the effects of frozen and equilibrium air chemistry for applications to high speed flight vehicles. To provide a computationally economical first order approximation to the high temperature physics, variable thermodynamic data is used for the chemically frozen mode to allow for a variation with temperature of the air specific heats and enthalpy. For calculations involving air in chemical equilibrium, a specially modified version of the NASA Lewis Chemical Equilibrium Code, CEC, is used to compute the chemical composition and resultant thermochemical properties. The upgraded solver is demonstrated by comparing results from calorically perfect (C sub p=constant), thermally perfect (frozen) and equilibrium air calculations for a variety of geometries, and flight Mach numbers.

Coirier, William J.↗

High speed inlet calculations with real gas effects

A 2-D steady-state Navier-Stokes solver has been upgraded to include the effects of frozen and equilibrium air chemistry for applications to high speed flight vehicles. To provide a computationally economical first order approximation to the high temperature physics, variable thermodynamic data is used for the chemically frozen mode to allow for a variation with temperature of the air specific heats and enthalpy. For calculations involving air in chemical equilibrium, a specially modified version of the NASA Lewis Chemical Equilibrium Code, CEC, is used to compute the chemical composition and resultant thermochemical properties. The upgraded solver is demonstrated by comparing results from calorically perfect (C sub p=constant), thermally perfect (frozen) and equilibrium air calculations for a variety of geometries, and flight Mach numbers.

Coirier, William J.↗

The Laboratory for Oceans Computing Facility

The first VAX computer in the Laboratory for Oceans Computing Facility (LOCF) was installed and the facility was largely expanded. The growth is not only in hardware and software, but also in the number of users and in supporting research and development projects. The LOCF serves as a general purpose computing facility for: ocean color research projects, sea ice research projects, processing of the Nimbus-7 Coastal Zone Color Scanner data set, real time ingest and analysis of TIROS-N satellite data, study of the Synthetic Aperture Radar data, study of LANDSAT data, and many others. The physical space and the electrical power layout of the computing room were modified to accommodate all the equipment. The LOCF has several image processing stations which include two International Imaging Systems (IIS) model 75 processors and one Adage processor. The facility has the capability of ingesting the TIROS-N HRPT satellite data on a real time basis. More than 30 software packages were installed on the systems. System software packages, network software, FORTRAN and C compilers, database management software, image processing software, graphics, mathematics and statistics packages, TAE, Catalog Manager, GEMPAK, LAS and many other software developed on the LOCF computers such as SEAPAK have greatly advanced the capability of the LOCF.

Kao, R.↗

A review of high-speed, convective, heat-transfer computation methods

The objective is to provide useful engineering formulations and to instill a modest degree of physical understanding of the phenomena governing convective aerodynamic heating at high flight speeds. Some physical insight is not only essential to the application of the information presented here, but also to the effective use of computer codes which may be available to the reader. Given first is a discussion of cold-wall, laminar boundary layer heating. A brief presentation of the complex boundary layer transition phenomenon follows. Next, cold-wall turbulent boundary layer heating is discussed. This topic is followed by a brief coverage of separated flow-region and shock-interaction heating. A review of heat protection methods follows, including the influence of mass addition on laminar and turbulent boundary layers. Next is a discussion of finite-difference computer codes and a comparison of some results from these codes. An extensive list of references is also provided from sources such as the various AIAA journals and NASA reports which are available in the open literature.

Tauber, Michael E.↗

Multiphase Transport in Porous Media: Gas-Liquid Separation Using Capillary Pressure Gradients International Space Station (ISS) Flight Experiment Development

Advances in the understanding of multiphase flow characteristics under variable gravity conditions will ultimately lead to improved and as of yet unknown process designs for advanced space missions. Such novel processes will be of paramount importance to the success of future manned space exploration as we venture into our solar system and beyond. In addition, because of the ubiquitous nature and vital importance of biological and environmental processes involving airwater mixtures, knowledge gained about fundamental interactions and the governing properties of these mixtures will clearly benefit the quality of life here on our home planet. The techniques addressed in the current research involving multiphase transport in porous media and gas-liquid phase separation using capillary pressure gradients are also a logical candidate for a future International Space Station (ISS) flight experiment. Importantly, the novel and potentially very accurate Lattice-Boltzmann (LB) modeling of multiphase transport in porous media developed in this work offers significantly improved predictions of real world fluid physics phenomena, thereby promoting advanced process designs for both space and terrestrial applications.This 3-year research effort has culminated in the design and testing of a zero-g demonstration prototype. Both the hydrophilic (glass) and hydrophobic (Teflon) media Capillary Pressure Gradient (CPG) cartridges prepared during the second years work were evaluated. Results obtained from ground testing at 1-g were compared to those obtained at reduced gravities spanning Martian (13-g), Lunar (16-g) and zero-g. These comparisons clearly demonstrate the relative strength of the CPG phenomena and the efficacy of its application to meet NASAs unique gas-liquid separation (GLS) requirements in non-terrestrial environments.LB modeling software, developed concurrently with the zero-g test effort, was shown to accurately reproduce observed CPG driven gas-liquid separation phenomena. The design and fabrication of a micropost plate-lamina Hele-Shaw (HS) cell was performed which served as a computationally attainable geometric structure facilitating direct comparison between physical phenomena observed in our laboratory and the LB software predictions.

Wheeler, Richard R., Jr.↗

Pre-Test CFD for the Design and Execution of the Enhanced Injection and Mixing Project at NASA Langley Research Center

With the increasing costs of physics experiments and simultaneous increase in availability and maturity of computational tools it is not surprising that computational fluid dynamics (CFD) is playing an increasingly important role, not only in post-test investigations, but also in the early stages of experimental planning. This paper describes a CFD-based effort executed in close collaboration between computational fluid dynamicists and experimentalists to develop a virtual experiment during the early planning stages of the Enhanced Injection and Mixing project at NASA Langley Research Center. This projects aims to investigate supersonic combustion ramjet (scramjet) fuel injection and mixing physics, improve the understanding of underlying physical processes, and develop enhancement strategies and functional relationships relevant to flight Mach numbers greater than 8. The purpose of the virtual experiment was to provide flow field data to aid in the design of the experimental apparatus and the in-stream rake probes, to verify the nonintrusive measurements based on NO-PLIF, and to perform pre-test analysis of quantities obtainable from the experiment and CFD. The approach also allowed for the joint team to develop common data processing and analysis tools, and to test research ideas. The virtual experiment consisted of a series of Reynolds-averaged simulations (RAS). These simulations included the facility nozzle, the experimental apparatus with a baseline strut injector, and the test cabin. Pure helium and helium-air mixtures were used to determine the efficacy of different inert gases to model hydrogen injection. The results of the simulations were analyzed by computing mixing efficiency, total pressure recovery, and stream thrust potential. As the experimental effort progresses, the simulation results will be compared with the experimental data to calibrate the modeling constants present in the CFD and validate simulation fidelity. CFD will also be used to investigate different injector concepts, improve understanding of the flow structure and flow physics, and develop functional relationships. Both RAS and large eddy simulations (LES) are planned for post-test analysis of the experimental data.

Drozda, Tomasz G.↗

Dissertation Defense Computational Fluid Dynamics Uncertainty Analysis for Payload Fairing Spacecraft Environmental Control Systems

Spacecraft thermal protection systems are at risk of being damaged due to airflow produced from Environmental Control Systems. There are inherent uncertainties and errors associated with using Computational Fluid Dynamics to predict the airflow field around a spacecraft from the Environmental Control System. This paper describes an approach to quantify the uncertainty in using Computational Fluid Dynamics to predict airflow speeds around an encapsulated spacecraft without the use of test data. Quantifying the uncertainty in analytical predictions is imperative to the success of any simulation-based product. The method could provide an alternative to traditional "validation by test only" mentality. This method could be extended to other disciplines and has potential to provide uncertainty for any numerical simulation, thus lowering the cost of performing these verifications while increasing the confidence in those predictions. Spacecraft requirements can include a maximum airflow speed to protect delicate instruments during ground processing. Computational Fluid Dynamics can be used to verify these requirements; however, the model must be validated by test data. This research includes the following three objectives and methods. Objective one is develop, model, and perform a Computational Fluid Dynamics analysis of three (3) generic, non-proprietary, environmental control systems and spacecraft configurations. Several commercially available and open source solvers have the capability to model the turbulent, highly three-dimensional, incompressible flow regime. The proposed method uses FLUENT, STARCCM+, and OPENFOAM. Objective two is to perform an uncertainty analysis of the Computational Fluid Dynamics model using the methodology found in "Comprehensive Approach to Verification and Validation of Computational Fluid Dynamics Simulations". This method requires three separate grids and solutions, which quantify the error bars around Computational Fluid Dynamics predictions. The method accounts for all uncertainty terms from both numerical and input variables. Objective three is to compile a table of uncertainty parameters that could be used to estimate the error in a Computational Fluid Dynamics model of the Environmental Control System /spacecraft system. Previous studies have looked at the uncertainty in a Computational Fluid Dynamics model for a single output variable at a single point, for example the re-attachment length of a backward facing step. For the flow regime being analyzed (turbulent, three-dimensional, incompressible), the error at a single point can propagate into the solution both via flow physics and numerical methods. Calculating the uncertainty in using Computational Fluid Dynamics to accurately predict airflow speeds around encapsulated spacecraft in is imperative to the success of future missions.

CFD↗

Dissertation Defense: Computational Fluid Dynamics Uncertainty Analysis for Payload Fairing Spacecraft Environmental Control Systems

Spacecraft thermal protection systems are at risk of being damaged due to airflow produced from Environmental Control Systems. There are inherent uncertainties and errors associated with using Computational Fluid Dynamics to predict the airflow field around a spacecraft from the Environmental Control System. This paper describes an approach to quantify the uncertainty in using Computational Fluid Dynamics to predict airflow speeds around an encapsulated spacecraft without the use of test data. Quantifying the uncertainty in analytical predictions is imperative to the success of any simulation-based product. The method could provide an alternative to traditional validation by test only mentality. This method could be extended to other disciplines and has potential to provide uncertainty for any numerical simulation, thus lowering the cost of performing these verifications while increasing the confidence in those predictions.Spacecraft requirements can include a maximum airflow speed to protect delicate instruments during ground processing. Computational Fluid Dynamics can be used to verify these requirements; however, the model must be validated by test data. This research includes the following three objectives and methods. Objective one is develop, model, and perform a Computational Fluid Dynamics analysis of three (3) generic, non-proprietary, environmental control systems and spacecraft configurations. Several commercially available and open source solvers have the capability to model the turbulent, highly three-dimensional, incompressible flow regime. The proposed method uses FLUENT, STARCCM+, and OPENFOAM. Objective two is to perform an uncertainty analysis of the Computational Fluid Dynamics model using the methodology found in Comprehensive Approach to Verification and Validation of Computational Fluid Dynamics Simulations. This method requires three separate grids and solutions, which quantify the error bars around Computational Fluid Dynamics predictions. The method accounts for all uncertainty terms from both numerical and input variables. Objective three is to compile a table of uncertainty parameters that could be used to estimate the error in a Computational Fluid Dynamics model of the Environmental Control System spacecraft system.Previous studies have looked at the uncertainty in a Computational Fluid Dynamics model for a single output variable at a single point, for example the re-attachment length of a backward facing step. For the flow regime being analyzed (turbulent, three-dimensional, incompressible), the error at a single point can propagate into the solution both via flow physics and numerical methods. Calculating the uncertainty in using Computational Fluid Dynamics to accurately predict airflow speeds around encapsulated spacecraft in is imperative to the success of future missions.

Fluid Dynamics↗

Utilization of Streamtubes to Analyze the Physical Interaction of a Dispersed Cloud with the CRM65 Hybrid Midspan Model

Recently there have been numerous efforts to identify the relevant icing physics related to the formation of complex three dimensional features, sometimes referred to as ’scallops’, on swept wings. However, much of the physics is still not well understood. This paper computationally investigates the interaction of the icing cloud with the 65 percent Common Research Model (CRM65). Both the interaction with an uniced model, and the interaction with a representative simulated three dimensional ice accretion are analyzed. Preliminary results suggest that the liquid water content increases near the aerodynamic body. For small droplets on the uniced geometry, the particle velocity vector becomes nearly parallel with the aerodynamic model resulting in a low value of collection efficiency. When a three dimensional feature with a length scale much smaller than the leading edge of the airfoil is introduced into the flow, the impingement of small particles can become significantly more perpendicular to the particle velocity vector in the region of these features. Since the liquid water content near the body has increased due to the interaction with the larger scale features of the aerodynamic model, i.e. the leading edge of the swept wing, extremely high collection efficiency is observed. These results suggest that three dimensional features are likely a significant physical driver that should be modeled in some capacity when simulating the impingement of a cloud on an aerodynamic model.

Icing↗

Utilization of Streamtubes to Analyze the Physical Interaction of a Dispersed Cloud with the CRM65 Hybrid Midspan Model

Recently there have been numerous efforts to identify the relevant icing physics related to the formation of complex three dimensional features, sometimes referred to as ’scallops’, on swept wings. However, much of the physics is still not well understood. This paper computationally investigates the interaction of the icing cloud with the 65 percent Common Research Model (CRM65). Both the interaction with an uniced model, and the interaction with a representative simulated three dimensional ice accretion are analyzed. Preliminary results suggest that the liquid water content increases near the aerodynamic body. For small droplets on the uniced geometry, the particle velocity vector becomes nearly parallel with the aerodynamic model resulting in a low value of collection efficiency. When a three dimensional feature with a length scale much smaller than the leading edge of the airfoil is introduced into the flow, the impingement of small particles can become significantly more perpendicular to the particle velocity vector in the region of these features. Since the liquid water content near the body has increased due to the interaction with the larger scale features of the aerodynamic model, i.e. the leading edge of the swept wing, extremely high collection efficiency is observed. These results suggest that three dimensional features are likely a significant physical driver that should be modeled in some capacity when simulating the impingement of a cloud on an aerodynamic model.

Icing↗

Partitioning and packing mathematical simulation models for calculation on parallel computers

The development of multiprocessor simulations from a serial set of ordinary differential equations describing a physical system is described. Degrees of parallelism (i.e., coupling between the equations) and their impact on parallel processing are discussed. The problem of identifying computational parallelism within sets of closely coupled equations that require the exchange of current values of variables is described. A technique is presented for identifying this parallelism and for partitioning the equations for parallel solution on a multiprocessor. An algorithm which packs the equations into a minimum number of processors is also described. The results of the packing algorithm when applied to a turbojet engine model are presented in terms of processor utilization.

Arpasi, D. J.↗

Mars Global Climate Modeling.

Scientists use Global Climate Models (GCMs) to better under the current and past climate states of terrestrial (solid surface) bodies in our solar system and beyond, and the physical processes that control them. GCMs are complex, multi-dimensional computer codes that can generally be divided into two parts: the geophysical fluid dynamics (GFD) framework, which represents accelerations and spatially resolved processes, and the physics routines, which provide the forcing functions for the circulation. Producing a GCM that is appropriate for a particular body—Mars, for example—requires implementing the appropriate physics routines (e.g., radiative transfer, planetary boundary layer physics, dust lifting physics to generate dust storms, etc.) for that body. In this talk, I will give an overview of the components of the NASA Ames Mars GCM and discuss some of the scientific questions we address with this state-of-the-art numerical model.

Melinda A. Kahre↗

Direct numerical simulation of sheared turbulent flow

The summer assignment to study sheared turbulent flow was divided into three phases which were: (1) literature survey, (2) computational familiarization, and (3) pilot computational studies. The governing equations of fluid dynamics or Navier-Stokes equations describe the velocity, pressure, and density as functions of position and time. In principle, when combined with conservation equations for mass, energy, and thermodynamic state of the fluid a determinate system could be obtained. In practice the Navier-Stokes equations have not been solved due to the nonlinear nature and complexity of these equations. Consequently, the importance of experiments in gaining insight for understanding the physics of the problem has been an ongoing process. Reasonable computer simulations of the problem have occured as the computational speed and storage of computers has evolved. The importance of the microstructure of the turbulence dictates the need for high resolution grids in extracting solutions which contain the physical mechanisms which are essential to a successful simulation. The recognized breakthrough occurred as a result of the pioneering work of Orzag and Patterson in which the Navier-Stokes equations were solved numerically utilizing a time saving toggling technique between physical and wave space, known as a spectral method. An equally analytically unsolvable problem, containing the same quasi-chaotic nature as turbulence, is known as the three body problem which was studied computationally as a first step this summer. This study was followed by computations of a two dimensional (2D) free shear layer.

Harris, Vascar G.↗

A computational study of multiple jet and wall interaction

In the fuel preburner of the Space Shuttle Main Engine, face plate, injector, and baffle erosion have been observed. The observed patterns of erosion suggest that flame attachment to the walls is a contributing factor. To better understand the physical phenomena involved, a portion of the preburner was modeled computationally. The simulated 'preburner' had three two-dimensional jets entering a cavity adjacent to a baffle. The computational model employed the Patankar Spalding algorithm with upwind differencing. The turbulence model was a standard k-epsilon model with wall functions. The effects of incoming boundary conditions on turbulent kinetic energy and dissipation, k and epsilon, was studied. The results indicate a very strong sensitivity to these boundary conditions over certain ranges of values.

Campbell, W.↗

An accurate radiative transfer model for use in the direct physical inversion of HIRS2 and MSU temperature sounding data

The direct computation of atmospheric transmittance and clear column radiances for the channels of HIRS2 and MSU as a function of atmospheric and surface conditions is described in detail. A comparison is made between the observations and the calculated radiances derived from colocated oceanic radiosondes. It is found that under clear conditions, calculated brightness temperatures for the HIRS2 have a standard deviation of the order of 0.7 C compared with observations, whereas MSU channels have a standard deviation of approximately 1 C. In some channels, small biases are found that can be removed by an empirical 'tuning' with coefficients that can be successfully transferred from one season to another. Less satisfactory agreement is obtained from a comparison of calculations with 'reconstructed' clear radiances, which are used in analyzing sounding data under partially cloudy conditions.

Susskind, J.↗

The study of the physics of cometary nuclei

The development and utilization of an optimized computer program to analyze orbital stabilization by repeated calculations is presented. The stability of comets in the Opik-Oort Cloud about the Sun against perturbations by the Galactic center involve the same basic type of calculation. The supposed persistence of these bodies in orbits over the life of the solar system, depends upon the stability of bodies of negligible mass in orbits around a body whose mass is small compared to the central mass about which they revolve. The question remains of preferential orientation of extremely eccentric comet orbits, possibly to explain the asymmetry observed among new comet motions. A third application of the computing programs is suited to meteoroids that may exist in orbits about asteroids and that may endanger science spacecraft making flybys too near to asteroids. As in the double-comet case, solar activity and solar gravitational perturbations limit the attendance to an asteroid by small meteroids in their orbits. It is found that the mass distances planned for asteroid fly-bys are adequate.

Whipple, F. L.↗