Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic boundary”

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 253 records · Page 14

Boundary-Layer Flow Simulations Over Ablating Woven Thermal Protection System Material

Spallation is the mechanical removal of small chunks of material gets removed typically due to high shear conditions of the flow field. This reduces the ability of the thermal protection system (TPS) material to protect the spacecraft as well as cause turbulence in the flow causing higher heating rates. In this work, we focus on the material removal through ablation and high shear flow within the boundary layer region of woven TPS material. Woven TPS (WTPS) material is the latest class of material developed by NASA, to be used within the next generation of space flights. They are complex interlocked weaves designed to create a rigid structure that is highly resistant to heat and can be easily designed and tailored for a wide variety of entry environments. Due to material removal resulting from chemical degradation, the structural integrity of TPS material is affected. Spallation occurs when this structurally compromised material is exposed to the high shear flow conditions within the boundary layer. In order to understand the spallation mechanism within WTPS material, we first perform the material removal simulations which occur primarily through oxidation to obtain the microstructure at various stages of degradation. These simulations are performed using the Porous Microstructure Analysis (PuMA) software developed at NASA Ames. The micro-structure geometry used within these simulations were generated artificially to be similar to the 3D weave architecture of MSR-EEV (Mars Sample Return - Earth Entry Vehicle). The various eroded TPS micro-structures are then subjected to the boundary layer flow conditions to obtain critical surface quantities which contribute to the structural failure mechanism such as heat flux, pressure, and shear stress. The direct simulation Monte Carlo (DSMC) methodology is used to perform these simulations in order to accurately capture the strong gradients within the high-temperature boundary layer flow over the intricate geometry of WTPS material. The boundary layer profile is directly taken from the Computational Fluid Dynamics (CFD) simulation and provided as boundary conditions to the DSMC inlet and outlet. Further, the variation of these properties as the microstructure undergoes changes due to oxidation is also investigated. Finally, these quantities are used as input in PuMA to understand the material expansion/compression and strain within the woven TPS geometry and help in developing a comprehensive spallation and structure failure model.

microstructure↗

Boundary-Layer Flow Simulations Over Ablating Woven Thermal Protection System Material

Spallation is the mechanical removal of small chunks of material gets removed typically due to high shear conditions of the flow field. This reduces the ability of the thermal protection system (TPS) material to protect the spacecraft as well as cause turbulence in the flow causing higher heating rates. In this work, we focus on the material removal through ablation and high shear flow within the boundary layer region of woven TPS material. Woven TPS (WTPS) material is the latest class of material developed by NASA, to be used within the next generation of space flights. They are complex interlocked weaves designed to create a rigid structure that is highly resistant to heat and can be easily designed and tailored for a wide variety of entry environments. Due to material removal resulting from chemical degradation, the structural integrity of TPS material is affected. Spallation occurs when this structurally compromised material is exposed to the high shear flow conditions within the boundary layer. In order to understand the spallation mechanism within WTPS material, we first perform the material removal simulations which occur primarily through oxidation to obtain the microstructure at various stages of degradation. These simulations are performed using the Porous Microstructure Analysis (PuMA) software developed at NASA Ames. The micro-structure geometry used within these simulations were generated artificially to be similar to the 3D weave architecture of MSR-EEV (Mars Sample Return - Earth Entry Vehicle). The various eroded TPS micro-structures are then subjected to the boundary layer flow conditions to obtain critical surface quantities which contribute to the structural failure mechanism such as heat flux, pressure, and shear stress. The direct simulation Monte Carlo (DSMC) methodology is used to perform these simulations in order to accurately capture the strong gradients within the high-temperature boundary layer flow over the intricate geometry of WTPS material. The boundary layer profile is directly taken from the Computational Fluid Dynamics (CFD) simulation and provided as boundary conditions to the DSMC inlet and outlet. Further, the variation of these properties as the microstructure undergoes changes due to oxidation is also investigated. Finally, these quantities are used as input in PuMA to understand the material expansion/compression and strain within the woven TPS geometry and help in developing a comprehensive spallation and structure failure model.

microstructure↗

Finding the magnetopause location using soft X-ray observations and a statistical inverse method

Variability in the location and shape of the dayside magnetopause is attributed to magnetic reconnection, a fundamental process that enables the transfer of mass, energy, and momentum from the solar wind into the magnetosphere. The spatial and temporal properties of the magnetopause, under varying solar and magnetospheric conditions, remain largely unknown because empirical studies using in-situ observations are challenging to interpret. Global wide field-of-view (FOV) imaging is the only means to simultaneously observe the spatial distribution of the plasma properties over the vast dayside magnetospheric region and, subsequently, quantify the energy transport from the interplanetary medium into the terrestrial magnetosphere. Two upcoming missions, ESA/CAS SMILE and NASA’s LEXI will provide wide-field imagery of the dayside magnetosheath in soft X-rays, an emission generated by charge exchange interactions between high charge-state heavy ions of solar wind origin and exospheric neutral atoms. High-cadence two-dimensional observations of the magnetosheath will allow the estimation of dynamic properties of its inner boundary, the magnetopause, and enable studies of its response to changes in the solar wind dynamic pressure and interplanetary magnetic field orientation. This work introduces a statistically-based estimation approach based on inverse theory to estimate the spatial distribution of magnetosheath soft X-ray emissivities and, with this, identify the location of the magnetopause over the Sun−Earth line. To do so, we simulate the magnetosheath structure using the MHD-based OpenGGCM model and generate synthetic soft X-ray images using LEXI’s orbit and attitude information. Our results show that 3-D estimations using the described statistically-based technique are robust against Poisson-distributed shot noise inherent to soft X-ray images. Also, our proposed methodology shows that the accuracy of both three-dimensional (3-D) estimation and the magnetopause standoff distance calculation highly depends on the observational point.

Gonzalo Cucho-Padin↗

Modal Testing of a Flexible Wing on a Dynamically Active Test Fixture Using the Fixed Base Correction Method

In modal testing and finite element model correlation, analysts desire modal results using free-free or rigid boundary conditions to ease comparisons of test versus analytical data. It is often expensive both in cost and schedule to build and test with boundary conditions that replicate the free-free or rigid boundaries. Static test fixtures for load testing are often large, heavy, and unyielding, and not provide adequate boundaries for modal tests because they are dynamically too flexible and often contain natural frequencies within the test article frequency range of interest. The dynamic coupling between the test article and test fixture complicates the model updating process because significant effort needs to be spent on modeling the test fixture and boundary conditions in addition to the test article. If the modal results could be corrected for fixture coupling, then setups used for other structural testing could be adequate for modal testing and would allow significant schedule and cost savings by eliminating a unique setup for only modal testing. To simplify future modal tests, this report describes a Fixed Base Correction method that was investigated during modal testing of a full-scale, half-span, flexible wing cantilevered from a static test fixture. The results of this Fixed Base Correction approach look very promising. The method aided in producing similar wing modal characteristics for two different physical boundary configurations of a dynamically active test fixture.

Natalie D Spivey↗

Removing spurious reflections from CFD solutions by using the Complex Cepstrum

The Complex Cepstrum is shown to remove spurious reflections from artificial boundaries in computational fluid dynamic (CFD) solutions. First, the Complex Cepstrum theory is presented. A model time sequence consisting of a direct signal and reflections is analyzed theoretically with the Complex Cepstrum, and it is shown that the direct signal uncontaminated by reflections may be recovered in the time domain. Next, the Complex Cepstrum is applied to one- and three-dimensional CFD solutions, and spurious reflections from the boundary conditions are removed. By eliminating spurious reflections introduced by artificial boundary conditions, the applicability of CFD methods to aeroacoustic problems is greatly enhanced.

Meadows, Kristine R.↗

Modal Testing of a Flexible Wing on a Dynamically Active Test Fixture Using Fixed Base Correction Method - IFASD 2019

In modal testing and finite element model correlation, analysts desire modal results using free-free or rigid boundary conditions to ease comparisons of test versus analytical data. It is often expensive both in cost and schedule to build and test with boundary conditions that replicate the free-free or rigid boundaries. Static test fixtures for load testing are often large, heavy, and unyielding, but do not provide adequate boundaries for modal tests because they are dynamically too flexible and often contain natural frequencies within the test article frequency range of interest. The dynamic coupling between the test article and test fixture complicates the model updating process because significant effort needs to be spent on modeling the test fixture and boundary conditions in addition to the test article. If the modal results could be corrected for fixture coupling, then setups used for other structural testing could be adequate for modal testing and would allow significant schedule and cost savings by eliminating a unique setup for only modal testing. To simplify future modal tests, this paper reports on a fixed base correction method that was investigated during modal testing of a full-scale, half-span, flexible wing cantilevered from a static test fixture. The results of this fixed base correction approach look very promising. The method aided in producing similar wing modal characteristics for two different physical boundary configurations of a dynamically active test fixture.

modal testing↗

Prediction of axial-flow instabilities in a turbojet engine by use of a multistage compressor simulation on the digital computer

A method of estimating the undistorted stall line for an axial-flow compressor by using the digital computer is presented. The method involves linearization of nonlinear dynamic equations about an operating point on a speed line, and then application of the first method of Lyapunov to determine the stability of the nonlinear system from the stability of the linear system. The method is applied to a simulation of the J85 compressor, which utilizes stage stacking and lumped volume techniques for the interstage regions to simulate steady-state and dynamic compressor performance. The stability boundary predicted by the digital simulation compares quite well with the stall line predicted by a dynamic simulation of the J85 compressor programmed on the analog computer. Since previous studies have shown that the analog-predicted stall line agrees well with the stall line of the compressor, the digital method presented is also a good means of estimating the stall line.

Daniele, C. J.↗

Mixed convective-dynamic roll vortices and their effects on initial wind and temperature profiles

A new nonlinear 14-coefficient spectral model of two-dimensional shallow Boussinesq flow is developed and used to investigate the onset and development of both dynamically and convectively forced boundary-layer rolls. The model is developed to accept arbitrary basic-state wind profiles as dynamic forcing, using an Ekman profile to provide a means for easy comparison with other studies. The results are qualitatively compared with those of previous theoretical and observational investigations. The rolls are shown to significantly alter the initial wind profile in the sense found by Faller and Kaylor (1967) and Brown (1970), but via a mechanism independent of the Coriolis force.

Haack, Tracy↗

Assessing the Tundra-taiga Boundary with Multi-Sensor Satellite Data

Monitoring the dynamics of the circumpolar boreal forest (taiga) and Arctic tundra boundary is important for understanding the causes and consequences of changes observed in these areas. This ecotone, the world's largest, stretches for over 13,400 km and marks the transition between the northern limits of forests and the southern margin of the tundra. Because of the inaccessibility and large extent of this zone, remote sensing data can play an important role for mapping the characteristics and monitoring the dynamics. Basic understanding of the capabilities of existing space borne instruments for these purposes is required. In this study we examined the use of several remote sensing techniques for identifying the existing tundra- taiga ecotone. These include Landsat-7, MISR, MODIS and RADARSAT data. Historical cover maps, recent forest stand measurements and high-resolution IKONOS images were used for local ground truth. It was found that a tundra-taiga transitional area can be characterized using multi- spectral Landsat ETM+ summer images, multi-angle MISR red band reflectance images, RADARSAT images with larger incidence angle, or multi-temporal and multi-spectral MODIS data. Because of different resolutions and spectral regions covered, the transition zone maps derived from different data types were not identical, but the general patterns were consistent.

Ranson, K. J.↗

Impact of Assimilating AIRS Cloud-Cleared Radiances on Atmospheric Dynamics and Polar Low Representation at High Latitudes

This study explores the sensitivity of planetary boundary layer height and related atmospheric dynamics to the assimilation of cloud-cleared AIRS (Atmospheric Infrared Sounder) radiances in the Goddard Earth Observing System (GEOS, version 5) data assimilation and forecast system during the boreal fall 2014 season using observing system experiments (OSEs). Examined here are comparisons between the current, operational approach of assimilating AIRS clear-sky radiances against the assimilation of cloud-cleared radiances (CCR). In polar regions, assimilation of AIRS CCRs is particularly beneficial because of the sparsity of conventional observations and the prevalence of extended low-level stratus cloud cover, which limit the ingestion of clear-sky data. Assimilation of hyperspectral infrared information from AIRS over the Arctic region slightly modifies the lower midtropospheric temperature structure, which in turn contributes to adjustments in geopotential height, affecting the baroclinic instability properties over the entire hemisphere and explaining the overall improvement in global forecast skill. In addition, it is shown that the assimilation of CCRs benefits the representation of convectively-driven small-scale cyclones at high latitudes in the same way as previously noted for tropical cyclones. Specifically, assimilation of CCRs create a temperature dipole over the top of meteorologically active and strongly convective systems such as polar, arctic, and antarctic lows, which helps constrain the analyzed representation of their scale and vertical structure.

Mcgrath-Spangler, E. L.↗

Embedded cavity drag in steady and unsteady flows

The numerical solution of the laminar boundary-layer flow over an embedded cavity is studied. The purpose is to examine the relevant drag characteristics of laminar cavity flow. The solution field is obtained in terms of velocity and vorticity variables, with the stream function and pressure derivable from the directly computed variables. An analysis and comparison is made among four square cavities, ranging in size from 0.25 to 1.00 boundary-layer thicknesses deep. The dominant flow features are examined in the vicinity of the cavity by means of the stream function and iso-vorticity contours. The dominant physics in the overall drag characteristics of the flow is examined by an analysis of the pressure and wall shear stress distributions in the cavity, and upstream and downstream of the cavity. Pressure forces and frictional forces in, and in the vicinity of, the cavity are determined. Stress relaxation distances, both upstream and downstream of the cavity, are calculated and analyzed. The flow dynamics of the boundary-layer flow over an embedded cavity is summarized. Finally, the relevance of the results to the control of flow separation in such flows is discussed.

Gatski, T. B.↗

Embedded cavity drag in steady and unsteady flows

The numerical solution of the laminar boundary-layer flow over an embedded cavity is studied. The purpose of the study is to examine the relevant drag characteristics of laminar cavity flow. The solution field is obtained in terms of velocity and vorticity variables, with the stream function and pressure derivable from the directly computed variables. An analysis and comparison is made among four square cavities, ranging in size from 0.25 to 1.00 boundary-layer thicknesses deep. The dominant flow features are examined in the vicinity of the cavity by means of the stream function and iso-vorticity contours. The dominant physics in the overall drag characteristics of the flow is examined by an analysis of the pressure and wall shear stress distributions in the cavity, and upstream and downstream of the cavity. Pressure forces and frictional forces in, and in the vicinity of, the cavity are determined. Stress relaxation distances, both upstream and downstream of the cavity, are calculated and analyzed. The flow dynamics of the boundary-layer flow over an embedded cavity is summarized. Finally, the relevance of the present results to the control of flow separation in such flows is discussed.

Gatski, T. B.↗

Modeling the Extended Neutral Atmosphere and Plasma Environment near Saturn

In the three years of this study we have published five papers in refereed journals. We have first examined satellite sources and their contribution to the observed neutral cloud. Based on the total calculated satellite sources and the spatial distribution of sputtered H20, we concluded that they cannot produce observed OH cloud. The E-ring contribution has been also studied in detail.In order to produce observed OH cloud we suggested that the E-ring might be the dominant source in inner Saturnian magnetosphere. We proposed a possible resupply mechanism which is needed to keep both E-ring and OH cloud in the present state: collisions between E-ring grains and remains of a disrupted satellite near Enceladus Lagrangian points. In this scenario a large amount of material, ranging from a few micrometers to hundred of meters, which is optically invisible at present, is likely to exist there. The fourth paper compares the magnetosheaths of the outer planets. A surprising result is that the hot proton component comprises about 40% of the total density, much larger than predicted by shock theory. Gas dynamic models of the boundaries show that the magnetospheres of Jupiter and Saturn are flattened at the poles. The last paper was published in GRL and is the first based of the model of neutrals developed as a main goal of this grant and which is now operational.This Monte Carlo collision code self- consistently determines the neutral distribution. from the rings and satellites until they are lost by ionization, by collisions with rings, moons, or Saturn, or by escape from Saturn. Our model is unique in that it includes the effects of plasma chemistry and both plasma-neutral and neutral neutral collisions to determine the dynamical evolution of the water group neutrals in Saturn's magnetosphere. The lifetimes of the neutrals against loss to photoionization, charge exchange, electron dissociation and electron impact dissociation are based on the model given by (Richardson et al. 1998) and vary with position in the magnetosphere. The dominant neutral dissociation channels H20->OH+H, H20-> O+H2, and OH->O+H2 are considered.

Richardson, John D.↗

Analysis of Flowfields over Four-Engine DC-X Rockets

The objective of this study is to validate a computational methodology for the aerodynamic performance of an advanced conical launch vehicle configuration. The computational methodology is based on a three-dimensional, viscous flow, pressure-based computational fluid dynamics formulation. Both wind-tunnel and ascent flight-test data are used for validation. Emphasis is placed on multiple-engine power-on effects. Computational characterization of the base drag in the critical subsonic regime is the focus of the validation effort; until recently, almost no multiple-engine data existed for a conical launch vehicle configuration. Parametric studies using high-order difference schemes are performed for the cold-flow tests, whereas grid studies are conducted for the flight tests. The computed vehicle axial force coefficients, forebody, aftbody, and base surface pressures compare favorably with those of tests. The results demonstrate that with adequate grid density and proper distribution, a high-order difference scheme, finite rate afterburning kinetics to model the plume chemistry, and a suitable turbulence model to describe separated flows, plume/air mixing, and boundary layers, computational fluid dynamics is a tool that can be used to predict the low-speed aerodynamic performance for rocket design and operations.

Wang, Ten-See↗

Vibration Response Predictions for Heavy Panel Mounted Components from Panel Acreage Environment Specifications

The development of new launch vehicles in the Aerospace industry often relies on response measurements taken from previously developed vehicles during various stages of liftoff and ascent, and from wind tunnel models. These measurements include sound pressure levels, dynamic pressures in turbulent boundary layers and accelerations. Rigorous statistical scaling methods are applied to the data to derive new environments and estimate the performance of new skin panel structures. Scaling methods have proven to be reliable, particularly for designs similar to the vehicles used as the basis for scaling, and especially in regions of smooth acreage without exterior protuberances or heavy components mounted to the panel. To account for response attenuation of a panel-mounted component due to its apparent mass at higher frequencies, the vibroacoustics engineer often reduces the acreage vibration according to a weight ratio first suggested by Barrett. The accuracy of the reduction is reduced with increased weight of the panel-mounted component, and does not account for low-frequency amplification of the component/panel response as a system. A method is proposed that combines acreage vibration from scaling methods with finite element analysis to account for the frequency-dependent dynamics of heavy panel-mounted components. Since the acreage and mass-loaded skins respond to the same dynamic input pressure, such pressure may be eliminated in favor of a frequency-dependent scaling function applied to the acreage vibration to predict the mass-loaded panel response. The scaling function replaces the Barrett weight ratio, and contains all of the dynamic character of the loaded and unloaded skin panels. The solution simplifies for spatially uncorrelated and fully correlated input pressures. Since the prediction uses finite element models of the loaded and unloaded skins, a rich suite of response data are available to the design engineer, including interface forces, stress and strain, as well as acceleration and displacement. An extension of the method is also developed to incorporate the effect of a local protuberance near a heavy component. Acreage environments from traditional scaling methods with and without protuberance effects serve as the basis for the extension. Authors:

Harrison, Phillip↗

The effect of latent heat release on synoptic-to-planetary wave interactions and its implication for satellite observations: Theoretical modeling

The project objectives are to develop process models to investigate the interaction of planetary and synoptic-scale waves including the effects of latent heat release (precipitation), nonlinear dynamics, physical and boundary-layer processes, and large-scale topography; to determine the importance of latent heat release for temporal variability and time-mean behavior of planetary and synoptic-scale waves; to compare the model results with available observations of planetary and synoptic wave variability; and to assess the implications of the results for monitoring precipitation in oceanic-storm tracks by satellite observing systems. Researchers have utilized two different models for this project: a two-level quasi-geostrophic model to study intraseasonal variability, anomalous circulations and the seasonal cycle, and a 10-level, multi-wave primitive equation model to validate the two-level Q-G model and examine effects of convection, surface processes, and spherical geometry. It explicitly resolves several planetary and synoptic waves and includes specific humidity (as a predicted variable), moist convection, and large-scale precipitation. In the past year researchers have concentrated on experiments with the multi-level primitive equation model. The dynamical part of that model is similar to the spectral model used by the National Meteorological Center for medium-range forecasts. The model includes parameterizations of large-scale condensation and moist convection. To test the validity of results regarding the influence of convective precipitation, researchers can use either one of two different convective schemes in the model, a Kuo convective scheme or a modified Arakawa-Schubert scheme which includes downdrafts. By choosing one or the other scheme, they can evaluate the impact of the convective parameterization on the circulation. In the past year researchers performed a variety of initial-value experiments with the primitive-equation model. Using initial conditions typical of climatological winter conditions, they examined the behavior of synoptic and planetary waves growing in moist and dry environments. Surface conditions were representative of a zonally averaged ocean. They found that moist convection associated with baroclinic wave development was confined to the subtropics.

Branscome, Lee E.↗

Transition Analysis for the Pazy Wing

The Pazy wing, designed as a benchmark for highly flexible wings that experience large deformations, presents intriguing nonlinear phenomena due to its unique geometry with spanwise varying static deformation and low Reynolds number flow conditions. These phenomena include limit cycle oscillation related to laminar separation bubbles, boundary-layer transition, and dynamic stall. The primary focus of this study is to investigate transition behavior on the sagged Pazy wing under static loading by using the shear-stress-transport (SST)-based Langtry-Menter 𝜸 − 𝑹𝒆𝜽𝒕 transition model in the NASA OVERFLOW and FUN3D solvers. Based on the steady-state results obtained via different variants of the Langtry-Menter model, we investigate the likely transition behavior over a range of incidence angles. The predicted flowfield on the sagged Pazy wing is strongly three dimensional and the application of turbulence index as a viable indicator of the transition front in such flows is also explored. The study demonstrates the limitations of the Langtry-Menter model when applied to low Reynolds number flows. We introduce a spanwise periodic configuration to simulate a narrow-span section of the Pazy wing, specifically centered on one of the 14 nearly equidistant ribs spanning the wing. The computational results demonstrate that the spanwise periodic configuration provides a straightforward environment for studying transition characteristics as well as demonstrating systematic grid convergence and successful code-to-code comparisons. The research offers useful insights into the transition characteristics of flexible wings under various loading conditions, contributing to the field of aeroelastic simulations.

CFD Modeling↗

Transition Analysis for the Pazy Wing

The Pazy wing, designed as a benchmark for highly flexible wings that experience large deformations, presents intriguing nonlinear phenomena due to its unique geometry with spanwise varying static deformation and low Reynolds number flow conditions. These phenomena include limit cycle oscillation related to laminar separation bubbles, boundary-layer transition, and dynamic stall. The primary focus of this study is to investigate transition behavior on the sagged Pazy wing under static loading by using the shear-stress-transport (SST)-based Langtry-Menter 𝜸 − 𝑹𝒆𝜽𝒕 transition model in the NASA OVERFLOW and FUN3D solvers. Based on the steady-state results obtained via different variants of the Langtry-Menter model, we investigate the likely transition behavior over a range of incidence angles. The predicted flowfield on the sagged Pazy wing is strongly three dimensional and the application of turbulence index as a viable indicator of the transition front in such flows is also explored. The study demonstrates the limitations of the Langtry-Menter model when applied to low Reynolds number flows. We introduce a spanwise periodic configuration to simulate a narrow-span section of the Pazy wing, specifically centered on one of the 14 nearly equidistant ribs spanning the wing. The computational results demonstrate that the spanwise periodic configuration provides a straightforward environment for studying transition characteristics as well as demonstrating systematic grid convergence and successful code-to-code comparisons. The research offers useful insights into the transition characteristics of flexible wings under various loading conditions, contributing to the field of aeroelastic simulations.

CFD Modeling↗