Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical modelling”

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,207 records · Page 67

Tidal disruption of inviscid protoplanets

Roche showed that equilibrium is impossible for a small fluid body synchronously orbiting a primary within a critical radius now termed the Roche limit. Tidal disruption of orbitally unbound bodies is a potentially important process for planetary formation through collisional accumulation, because the area of the Roche limit is considerably larger then the physical cross section of a protoplanet. Several previous studies were made of dynamical tidal disruption and different models of disruption were proposed. Because of the limitation of these analytical models, we have used a smoothed particle hydrodynamics (SPH) code to model the tidal disruption process. The code is basically the same as the one used to model giant impacts; we simply choose impact parameters large enough to avoid collisions. The primary and secondary both have iron cores and silicate mantles, and are initially isothermal at a molten temperature. The conclusions based on the analytical and numerical models are summarized.

Boss, Alan P.↗

Flowfield analysis of modern helicopter rotors in hover by Navier-Stokes method

The viscous, three-dimensional, flowfields of UH60 and BERP rotors are calculated for lifting hover configurations using a Navier-Stokes computational fluid dynamics method with a view to understand the importance of planform effects on the airloads. In this method, the induced effects of the wake, including the interaction of tip vortices with successive blades, are captured as a part of the overall flowfield solution without prescribing any wake models. Numerical results in the form of surface pressures, hover performance parameters, surface skin friction and tip vortex patterns, and vortex wake trajectory are presented at two thrust conditions for UH60 and BERP rotors. Comparison of results for the UH60 model rotor show good agreement with experiments at moderate thrust conditions. Comparison of results with equivalent rectangular UH60 blade and BERP blade indicates that the BERP blade, with an unconventional planform, gives more thrust at the cost of more power and a reduced figure of merit. The high thrust conditions considered produce severe shock-induced flow separation for UH60 blade, while the BERP blade develops more thrust and minimal separation. The BERP blade produces a tighter tip vortex structure compared with the UH60 blade. These results and the discussion presented bring out the similarities and differences between the two rotors.

Srinivasan, G. R.↗

Modeled Seasonal Variations of Firn Density Induced by Steady State Surface Air Temperature Cycle

Seasonal variations of firn density in ice-sheet firn layers have been attributed to variations in deposition processes or other processes within the upper firn. A recent high-resolution (mm scale) density profile, measured along a 181 m core from Antarctica, showed small-scale density variations with a clear seasonal cycle that apparently was not-related to seasonal variations in deposition or known near-surface processes (Gerland and others 1999). A recent model of surface elevation changes (Zwally and Li, submitted) produced a seasonal variation in firn densification, and explained the seasonal surface elevation changes observed by satellite radar altimeters. In this study, we apply our 1-D time-dependent numerical model of firn densification that includes a temperature-dependent formulation of firn densification based on laboratory measurements of grain growth. The model is driven by a steady-state seasonal surface temperature and a constant accumulation rate appropriate for the measured Antarctic ice core. The modeled seasonal variations in firn density show that the layers of snow deposited during spring to mid-summer with the highest temperature history compress to the highest density, and the layers deposited during later summer to autumn with the lowest temperature history compress to the lowest density. The initial amplitude of the seasonal difference of about 0.13 reduces to about 0.09 in five years and asymptotically to 0.92 at depth, which is consistent with the core measurements.

Jun, Li↗

Aircraft Engine Systems

The objective is to develop the capability to numerically model the performance of gas turbine engines used for aircraft propulsion. This capability will provide turbine engine designers with a means of accurately predicting the performance of new engines in a system environment prior to building and testing. The 'numerical test cell' developed under this project will reduce the number of component and engine tests required during development. As a result, the project will help to reduce the design cycle time and cost of gas turbine engines. This capability will be distributed to U.S. turbine engine manufacturers and air framers. This project focuses on goals of maintaining U.S. superiority in commercial gas turbine engine development for the aeronautics industry.

Veres, Joseph P.↗

Radiated Sound Power from a Curved Honeycomb Panel

The validation of finite element and boundary element model for the vibro-acoustic response of a curved honeycomb core composite aircraft panel is completed. The finite element and boundary element models were previously validated separately. This validation process was hampered significantly by the method in which the panel was installed in the test facility. The fixture used was made primarily of fiberboard and the panel was held in a groove in the fiberboard by a compression fitting made of plastic tubing. The validated model is intended to be used to evaluate noise reduction concepts from both an experimental and analytic basis simultaneously. An initial parametric study of the influence of core thickness on the radiated sound power from this panel, using this numerical model was subsequently conducted. This study was significantly influenced by the presence of strong boundary condition effects but indicated that the radiated sound power from this panel was insensitive to core thickness primarily due to the offsetting effects of added mass and added stiffness in the frequency range investigated.

Robinson, Jay H.↗

Energetic Ion Interactions with the Galilean Satellites

The principal research tasks of this investigation are: (1) specification of the energetic (keV to MeV) ion environments upstream of the four Galilean satellites and (2) data analysis and numerical modeling of observed ion interactions with the satellites. Differential flux spectra are being compiled for the most abundant ions (protons, oxygen, and sulfur) from measurements at 20 keV to 100 MeV total energy by the Energetic Particle Detector (EPD) experiment and at higher ion energies by the Heavy Ion Counter (HIC) experiment. Runge-Kutta and other numerical techniques are used to propagate test particles sampled from the measured upstream spectra to the satellite surface or spacecraft through the local magnetic and corotational electric field environment of each satellite. Modeling of spatial variations in directional flux anisotropies measured during each close flyby provides limits on atomic charge states for heavy (O, S) magnetospheric ions and on internal or induced magnetic fields of the satellites. Validation of models for magnetic and electric field configurations then allows computation of rates for ion implantation, sputtering, and energy deposition into the satellite surfaces for further modeling of observable chemical changes induced by irradiation. Our ongoing work on production of oxidants and other secondary species by ice irradiation on Europa's surface has significant applications, already acknowledged in current literature, to astrobiological evolution. Finally, the work will improve understanding of energetic ion sources and sinks at the satellite orbits for improved modeling of magnetospheric transport processes. The scope of the research effort mainly includes data from the primary Galileo mission (1995-1997) but may also include some later data where directly relevant (e.g., comparison of J0 and I27 data for Io) to the primary mission objectives. Funding for this contract also includes partial support for our related education and public outreach activity, 'Interaction of Magnetospheric Particles Applied to Classroom Teaching (IMPACT).'

Cooper, John F.↗

A numerically efficient finite element hydroelastic analysis. Volume 2: Implementation in NASTRAN, part 1

Details are presented of the implementation of the new formulation into NASTRAN including descriptions of the DMAP statements required for conversion of the program and details pertaining to problem definition and bulk data considerations. Details of the current 1/8-scale space shuttle external tank mathematical model, numerical results and analysis/test comparisons are also presented. The appendices include a description and listing of a FORTRAN program used to develop harmonic transformation bulk data (multipoint constraint statements) and sample bulk data information for a number of hydroelastic problems.

Coppolino, R. N.↗

Equatorial waves in the NCAR stratospheric general circulation model

Equatorially trapped wave modes are very important in the tropical stratospheric momentum balance. Kelvin waves and mixed Rossby-gravity waves are believed to be responsible for the quasi-biennial oscillation of the zonal winds in the equatorial lower stratosphere. Both Kelvin and mixed Rossby-gravity waves have been identified in observations and in numerical models. Kelvin and mixed Rossby-gravity waves are identified in a general circulation model extending from the surface into the mesosphere and looks at the effect on the waves of lowering the top of the model.

Boville, B. A.↗

Orbiter Boundary Layer Transition Stability Modeling at Flight Entry Conditions

State of the art boundary layer stability modeling capabilities are increasingly seeing application to entry flight vehicles. With the advent of user friendly and robust implementations of two-dimensional chemical nonequilibrium stability modeling with the STABL/PSE-CHEM software, the need for flight data to calibrate such analyses capabilities becomes more critical. Recent efforts to perform entry flight testing with the Orbiter geometry related to entry aerothermodynamics and boundary layer transition is allowing for a heightened focus on the Orbiter configuration. A significant advancement in the state of the art can likely be achieved by establishing a basis of understanding for the occurrence of boundary layer transition on the Orbiter due to discrete protruding gap fillers and the nominal distributed roughness of the actual thermal protection system. Recent success in demonstrating centerline two-dimensional stability modeling on the centerline of the Orbiter at flight entry conditions provides a starting point for additional investigations. The more detailed paper will include smooth Orbiter configuration boundary layer stability results for several typical orbiter entry conditions. In addition, the numerical modeling approach for establishing the mean laminar flow will be reviewed and the method for determining boundary layer disturbance growth will be overviewed. In addition, if actual Orbiter TPS surface data obtained via digital surface scans become available, it may be possible to investigate the effects of an as-flown flight configuration on boundary layer transition compared to a smooth CAD reference.

Bartkowicz, Matt↗

Improvements in the Scalability of the NASA Goddard Multiscale Modeling Framework for Hurricane Climate Studies

Improving our understanding of hurricane inter-annual variability and the impact of climate change (e.g., doubling CO2 and/or global warming) on hurricanes brings both scientific and computational challenges to researchers. As hurricane dynamics involves multiscale interactions among synoptic-scale flows, mesoscale vortices, and small-scale cloud motions, an ideal numerical model suitable for hurricane studies should demonstrate its capabilities in simulating these interactions. The newly-developed multiscale modeling framework (MMF, Tao et al., 2007) and the substantial computing power by the NASA Columbia supercomputer show promise in pursuing the related studies, as the MMF inherits the advantages of two NASA state-of-the-art modeling components: the GEOS4/fvGCM and 2D GCEs. This article focuses on the computational issues and proposes a revised methodology to improve the MMF's performance and scalability. It is shown that this prototype implementation enables 12-fold performance improvements with 364 CPUs, thereby making it more feasible to study hurricane climate.

Shen, Bo-Wen↗

An experimental study of the effects of rapid rotation on turbulence

Experiments, large eddy simulations and direct numerical simulations all show that rapid rotation has a striking effect on homogeneous nearly isotropic turbulence. The cascade process is effectively inhibited by rotation, and thus dissipation is greatly reduced. Some attempts have been made to incorporate this effect in turbulence models. Numerical simulation further showed the somewhat surprising result that anisotropic homogeneous turbulence subjected to rotation tended towards an isotropic state; however, the residual anisotropy was not zero. Reynolds (1989) performed a detailed analysis using Rapid Distortion Theory (RDT) and showed that a reduction in the anisotropy is indeed expected and if the anisotropy is produced by irrotational strain then the anisotropy tensor b(sub ij) is asymptotically driven to half its initial value. Our objective is to extend the work of the experiments mentioned above to lower turbulent Rossby numbers and to confirm some of the results obtained by Reynolds (1989) for anisotropic turbulence.

Veeravalli, S. V.↗

Bubble Dynamics on a Heated Surface

In this work, we study steady and oscillatory thermocapillary and natural convective flows generated by a bubble on a heated solid surface. The interaction between gas and vapor bubbles with the surrounding fluid is of interest for both space and ground-based processing. A combined numerical-experimental approach is adopted here. The temperature field is visualized using Mach-Zehnder and/or Wollaston Prism Interferometry and the flow field is observed by a laser sheet flow visualization technique. A finite element numerical model is developed which solves the transient two-dimensional continuity, momentum, and energy equations and includes the effects of temperature-dependent surface tension and bubble surface deformation. Below the critical Marangoni number, the steady state low-g and 1-g temperature and velocity fields predicted by the finite element model are in excellent agreement with both the visualization experiments in our laboratory and recently published experimental results in the literature. Above the critical Marangoni number, the model predicts an oscillatory flow which is also closely confirmed by experiments. It is shown that the dynamics of the oscillatory flow are directly controlled by the thermal and hydrodynamic interactions brought about by combined natural and thermocapillary convection. Therefore, as numerical simulations show, there are considerable differences between the 1-g and low-g temperature and flow fields at both low and high Marangoni numbers. This has serious implications for both materials processing and fluid management in space.

Kassemi, M.↗

Development and Testing of Coupled Land-surface, PBL and Shallow/Deep Convective Parameterizations within the MM5

The objective of this investigation was to study the role of shallow convection on the regional water cycle of the Mississippi and Little Washita Basins using a 3-D mesoscale model, the PSUINCAR MM5. The underlying premise of the project was that current modeling of regional-scale climate and moisture cycles over the continents is deficient without adequate treatment of shallow convection. It was hypothesized that an improved treatment of the regional water cycle can be achieved by using a 3-D mesoscale numerical model having a detailed land-surface parameterization, an advanced boundary-layer parameterization, and a more complete shallow convection parameterization than are available in most current models. The methodology was based on the application in the MM5 of new or recently improved parameterizations covering these three physical processes. Therefore, the work plan focused on integrating, improving, and testing these parameterizations in the MM5 and applying them to study water-cycle processes over the Southern Great Plains (SGP): (1) the Parameterization for Land-Atmosphere-Cloud Exchange (PLACE) described by Wetzel and Boone; (2) the 1.5-order turbulent kinetic energy (TKE)-predicting scheme of Shafran et al.; and (3) the hybrid-closure sub-grid shallow convection parameterization of Deng. Each of these schemes has been tested extensively through this study and the latter two have been improved significantly to extend their capabilities.

Stauffer, David R.↗

Theoretical overview and modeling of the sodium and potassium atmospheres of the moon

A general theoretical overview for the sources, sinks, gas-surface interactions, and transport dynamics of sodium and potassium in the exospheric atmosphere of the Moon is given. These four factors, which control the spatial distribution of these two alkali-group gases about the Moon, are incorporated in numerical models. The spatial nature and relative importance of the initial source atoms atmosphere (which must be nonthermal to explain observational data) and the ambient (ballistic hopping) atom atmosphere are examined. The transport dynamics, atmospheric structure, and lunar escape of the nonthermal source atoms are time variable with season of the year and lunar phase because of their dependence on the radiation acceleration experienced by sodium and potassium atoms as they resonantly scatter solar photons. The dynamic transport time of fully thermally accomodated ambient atoms along the surface because of solar radiation acceleration (only several percent of surface gravity) is larger than the photoionization lifetimes and hence unimportant in determining the local density, although for potassium the situation is borderline. The sodium model was applied to analyze sodium observations of the sunward brightness profiles acquired near last quarter by Potter & Morgan (1988b) extending from the surface to an altitude of 1200 km, and near first quarter by Mendillo, Baumgardner, & Flynn (1991), extending in altitude from approximately 1430 to approximately 7000 km. The observations at larger altitudes could be fitted only for source atoms having a velocity distribution with a tail that is mildly nonthermal (like an approximately 1000 K Maxwell-Boltzmann distribution). Solar wind sputtering appears to a be a viable source atom mechanism for the sodium observations with photon-simulated desorption also possible but highly uncertain, although micrometeoroid impact vaporization appears to have a source that is too small and too hot, with likely an incorrect angular distribution about the Moon.

Smyth, William H.↗

Dynamic modeling of coronal and interplanetary responses to solar events

Recent progress in the dynamic modeling of responses of the corona and interplanetary medium to solar events (such as surges, eruptive prominences, flares, etc.) is reviewed. In particular, coronal transients and wave phenomena are discussed in some detail including pertinent mathematical requirements. Within the context of hydrodynamics and magnetohydrodynamics, a summary of both one- and two-dimensional time-dependent models is presented. A comparison of theoretical results with ground-based optical, radio, ATM Skylab observations, Pioneer 9 and Pioneer 10 solar wind observations is also presented. It is illustrated that: (1) substantial progress has been made in the theoretical (i.e., numerical) modeling of coronal and interplanetary responses within the last few years; (2) two-dimensional, time-dependent modeling is needed to examine the details of nonlinear wave coupling; and (3) theoretical results of modeling appear to reproduce physical consequences successfully.

Wu, S. T.↗

The Magnetic Field in the Outer Heliosphere

One of the great achievements of Parker was the prediction that the solar magnetic field would be drawn into Archimedian spirals as it is carried away from the Sun by the solar wind. This prediction has been amply confirmed by many in situ measurements in the intervening four decades. But, Parker made his prediction for a solar wind that expands into infinite space while we now know that the local interstellar medium (LISM) is far from empty and, in fact, confines the solar wind to a finite volume, known as the heliosphere, that extends to approximately 100 AU in the upstream direction (the solar system is moving through the LISM). Voyagers 1/2, presently at -80 AU, are approaching the upstream boundaries of the heliosphere and returning data on the properties of the magnetic field. This is important for understanding how galactic cosmic rays (GCRs) reach the Earth. Voyagers show that the IMF at 10-80 AU behaves much as Parker predicted - with two important exceptions. This is not surprising since the field is essentially passively advected by the solar wind out to 80 AU. But, new models say that nearer the heliosphere boundaries the field plays a major role in the solar wind-LISM interaction. However, of the many physical ingredients that constitute the outer heliosphere, the magnetic field poses some of the most interesting and difficult numerical modeling problems. Presently, only a few results have been published and much remains to be done. Here I will summarize the expected and measured behavior of the magnetic field at 80 AU. Then I will describe modeling predictions beyond 80 AU: magnetic "tornadoes", polarity envelopes, the Axford-Cranfill effect, inner and outer magnetic walls and more. I will also list what I believe to be important new modeling objectives. Finally, I will speculate on what is happening with the magnetic field near the nose of the heliosphere. My conclusion is that models of GCR modulation rarely incorporate even crudely realistic magnetic fields so it is a wonder that they are as successful as they are and no surprise that there are still important discrepancies between GCR modulation observations and the models.

Suess, S. T.↗

The NASA NEESPI Data Portal: Products, Information, and Services

Studies have indicated that land cover and use changes in Northern Eurasia influence global climate system. However, the procedures are not fully understood and it is challenging to understand the interactions between the land changes in this region and the global climate. Having integrated data collections form multiple disciplines are important for studies of climate and environmental changes. Remote sensed and model data are particularly important die to sparse in situ measurements in many Eurasia regions especially in Siberia. The NASA GES DISC (Goddard Earth Sciences Data and Information Services Center) NEESPI data portal has generated infrastructure to provide satellite remote sensing and numerical model data for atmospheric, land surface, and cryosphere. Data searching, subsetting, and downloading functions are available. ONe useful tool is the Web-based online data analysis and visualization system, Giovanni (Goddard Interactive Online Visualization ANd aNalysis Infrastructure), which allows scientists to assess easily the state and dynamics of terrestrial ecosystems in Northern Eurasia and their interactions with global climate system. Recently, we have created a metadata database prototype to expand the NASA NEESPI data portal for providing a venue for NEESPI scientists fo find the desired data easily and leveraging data sharing within NEESPI projects. The database provides product level information. The desired data can be found through navigation and free text search and narrowed down by filtering with a number of constraints. In addition, we have developed a Web Map Service (WMS) prototype to allow access data and images from difference data resources.

Shen, Suhung↗

Aerothermodynamics of Blunt Body Entry Vehicles

In this chapter, the aerothermodynamic phenomena of blunt body entry vehicles are discussed. Four topics will be considered that present challenges to current computational modeling techniques for blunt body environments: turbulent flow, non-equilibrium flow, rarefied flow, and radiation transport. Examples of comparisons between computational tools to ground and flight-test data will be presented in order to illustrate the challenges existing in the numerical modeling of each of these phenomena and to provide test cases for evaluation of Computational Fluid Dynamics (CFD) code predictions.

Hollis, Brian R.↗