SEARCH · Engineering Papers
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.
Experiments are reported which show that plasticity-induced crack closure can account for crack growth retardation following an overload. The finite element method is used to provide evidence which supports the experimental observations of crack closure. Finally, a simple model is presented which predicts with limited success the retardation transient following an overload.
Three-dimensional effects related to the analysis of stitched textile composites are discussed. The method of calculation is based on the virtual crack closure technique (VCCT), and models that model the upper and lower surface of the delamination or debond with two-dimensional (2D) plate elements rather than three-dimensional (3D) solid elements. The major advantages of the plate element modeling technique are a smaller model size and simpler geometric modeling. Details of the modeling of the laminated plate and the stitching are discussed.
Report consists of: 1. List of accomplishments 2. List of publications 3. Abstracts of published or submitted papers and 4. Subject invention disclosure. The accomplishments of the grant listed are: 1. Improved the third-order turbulence closure in cloud resolving models to remove the liquid water oscillation. 2. Used the University of California-Los Angeles (UCLA) large-eddy simulation (LES) model to provide data for radiation transfer testing. 3. Revised shortwave k-distribution models based on HITRAN 2000. 4. Developed a gamma-weighted two-stream radiative transfer model for radiation budget estimate applications. 5. Estimated the effect of spherical geometry to the earth radiation budget. 6. Estimated top-of-atmosphere irradiance over snow and sea ice surfaces. 7. Estimated the aerosol direct radiative effect at the top of the atmosphere. 8. Estimated the top-of-atmosphere reflectance of the clear-sky molecular atmosphere over ocean. 9. Developed and validated new set of Angular Distribution Models for the CERES TRMM satellite instrument (tropical) 10. Developed and validated new set of Angular Distribution Models for the CERES Terra satellite instrument (global) 11. Quantified the top-of-atmosphere direct radiative effect of aerosols over global oceans from merged CERES and MODIS observations 12 Clarified the definition of TOA flux reference level for radiation budget studies 13. Developed new algorithm for unfaltering CERES measured radiances 14. Used multiangle POLDER measurements to produce narrowband angular distribution models and examine the effect of scene identification errors on TOA albedo estimates 15. Developed and validated a novel algorithm called the Multidirectional Reflectance Matching (MRM) model for inferring TOA albedos from ice clouds using multi-angle satellite measurements. 16. Developed and validated a novel algorithm called the Multidirectional Polarized Reflectance Matching (MPRM) model for inferring particle shapes from ice clouds using multi-angle polarized satellite measurements. 17. Developed 4 advanced light scattering models including the three-dimensional (3D) uniaxial perfectly matched layer (UPML) finite-difference time-domain (FDTD) model. 18. Develop sunglint in situ measurement and study reflectance distribution in the sunglint area. 19. Lead a balloon-borne radiometer TOA albedo validation effort. 20. Developed a CERES surface UVB, UVA, and UV index product.
The radiative heating increase due to increased CO2 concentration is the primary source for the rapid adjustment of atmospheric circulation and clouds. In this study, we investigate the rapid adjustment resulting from doubling of CO2 and its physical mechanism using a multiscale modeling framework (MMF). The MMF includes an advanced higher-order turbulence closure in its cloud-resolving model component and simulates realistic shallow and deep cloud climatology and boundary layer turbulence. The rapid adjustment over the tropics is characterized by 1) reduced ascent and descent strengths over the ocean, 2) increased lower tropospheric stability (LTS) over the subsidence region, 3) shoaling of planetary boundary layers over the ocean, 4) increased deep convection over lands and shift of cloud coverage from the ocean to lands, and 5) reduced sensible (SH) and latent heat (LH) fluxes over the oceanic deep convective regions. Unlike conventional general circulation models and another MMF, a reduction in the global-mean shortwave cloud radiative cooling is not simulated, due to the increase in low clouds at lower altitudes over the ocean, resulting from reduced cloud-top entrainment due to strengthened inversion. Changes in regional circulation play a key role in cloud changes and shift of cloud coverage to lands. Weaker energy transport resulting from water vapor and cloud CO2 masking effects reduces the upward motion and convective clouds in the oceanic regions. The ocean-land transports are linked to the partitioning of surface SH and LH fluxes that increases humidity over lands and enhances deep convection over the tropical lands.
We combined detailed bio-optical measurements and radiative transfer (RT) modeling to perform an optical closure experiment for optically complex and biologically productive Chesapeake Bay waters. We used this experiment to evaluate certain assumptions commonly used when modeling bio-optical processes, and to investigate the relative importance of several optical characteristics needed to accurately model and interpret remote sensing ocean-color observations in these Case 2 waters. Direct measurements were made of the magnitude, variability, and spectral characteristics of backscattering and absorption that are critical for accurate parameterizations in satellite bio-optical algorithms and underwater RT simulations. We found that the ratio of backscattering to total scattering in the mid-mesohaline Chesapeake Bay varied considerably depending on particulate loading, distance from land, and mixing processes, and had an average value of 0.0128 at 530 nm. Incorporating information on the magnitude, variability, and spectral characteristics of particulate backscattering into the RT model, rather than using a volume scattering function commonly assumed for turbid waters, was critical to obtaining agreement between RT calculations and measured radiometric quantities. In situ measurements of absorption coefficients need to be corrected for systematic overestimation due to scattering errors, and this correction commonly employs the assumption that absorption by particulate matter at near infrared wavelengths is zero.
A methodology is presented for a modal suppression control law design using flight test data instead of mathematical models to obtain the required gain and phase information about the flexible airplane. This approach is referred to as BODEDIRECT. The purpose of the BODEDIRECT program is to provide a method of analyzing the modal phase relationships measured directly from the airplane. These measurements can be achieved with a frequency sweep at the control surface input while measuring the outputs of interest. The measured Bode-models can be used directly for analysis in the frequency domain, and for control law design. Besides providing a more accurate representation for the system inputs and outputs of interest, this method is quick and relatively inexpensive. To date, the BODEDIRECT program has been tested and verified for computational integrity. Its capabilities include calculation of series, parallel and loop closure connections between Bode-model representations. System PSD, together with gain and phase margins of stability may be calculated for successive loop closures of multi-input/multi-output systems. Current plans include extensive flight testing to obtain a Bode-model representation of a commercial aircraft for design of a structural stability augmentation system.
This paper develops improved stochastic models for the description of a large variety of fatigue crack growth phenomena that occur in components of considerable importance to the functionality and reliability of complex engineering structures. In essence, the models are based on the McGill-Markov and Closure-Lognormal stochastic processes. Not only do these models have the capability of predicting the statistical dispersion of crack growth rates, they also, by incorporating the concept of crack closure, have the capability of transferring stochastic crack growth properties measured under ideal laboratory conditions to situations of industrial significance, such as those occurring under adverse loading and/or environmental conditions. The primary data required in order to be in a position to estimate the pertinent parameters of these stochastic models are obtained from a statistically significant number of replicate tests. In this paper, both the theory and the experimental technique are illustrated using a Ti-6Al-4V alloy. Finally, important structural integrity, reliability, availability and maintainability concepts are developed and illustrated.
This viewgraph presentation discusses (1) turbulence modeling: challenges in turbulence modeling, desirable attributes of turbulence models, turbulence models in FLUENT, and examples using FLUENT; and (2) combustion modeling: turbulence-chemistry interaction and FLUENT equilibrium model. As of now, three turbulence models are provided: the conventional k-epsilon model, the renormalization group model, and the Reynolds-stress model. The renormalization group k-epsilon model has broadened the range of applicability of two-equation turbulence models. The Reynolds-stress model has proved useful for strongly anisotropic flows such as those encountered in cyclones, swirlers, and combustors. Issues remain, such as near-wall closure, with all classes of models.
The aerosol-ice formation closure pilot study (AEROICESTUDY) was a field campaign conducted at the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility at the Southern Great Plains (SGP) observatory from October 7 to 28, 2019. The purpose of this campaign was to test a field observational approach for conducting an aerosol-ice formation closure study. In other words, the goal was to predict the number of ice nucleating particles (INPs) from the measurements of physical and chemical aerosol properties and compare those to measured INP number concentrations. This project was motivated by the fact that the aerosol community has widely conducted aerosol radiative closure, aerosol-cloud condensation nuclei (CCN), and CCN-droplet closure studies to test the physical models and parameterizations that cloud-resolving and climate models rely on to perform reliable simulations of the Earth system and energy budget. However, very few closure studies related to INPs have been conducted, and to our knowledge, none using robust, size-resolved, ambient aerosol composition measurements as model inputs
We estimate the impact of North Atlantic aerosols on the net short-wave flux at the tropopause by combining satellite-derived aerosol optical depth (AOD) maps with model aerosol properties determined via closure analyses in TARFOX and ACE 2. We exclude African dust, primarily by restricting latitudes to 25-60 N. The analyses use in situ aerosol composition measurements and air- and ship-borne sun-photometer measurements of AOD spectra. The aerosol model yields computed flux sensitivities (dFlux/dAOD) that agree with measurements by airborne flux radiometers in TARFOX. Its midvisible single-scattering albedo is 0.9. which is in the range obtained from in situ measurements of scattering and absorption in both TARFOX and ACE 2. Combining satellite-derived AOD maps with the aerosol model yields maps of 24-hour average net radiative flux changes. For simultaneous AVHRR, radiance measurements exceeded the sunphotometer AODs by about 0.04. However. shipboard sunphotometer and AVHRR AODs agreed Within 0.02 for data acquired during satellite overflights on two other days. We discuss attempts to demonstrate column closure within the MBL by comparing shipboard sunphotometer AODs and values calculated from simultaneous shipboard in-situ aerosol size distribution measurements. These comparisons were mostly unsuccessful, but they illustrate the difficulties inherent in this type of closure analysis. Specifically, AODs derived from near-surface in-situ size distribution measurements are extremely sensitive to the assumed hygroscopic growth model that itself requires an assumption of particle composition as a function of height and size, to the radiosonde-measured relative humidity, and to the vertical profile of particle number. We investigate further the effects of hygroscopic particle growth within the MBL by using shipboard lidar aerosol backscatter profiles together with the sunphotometer AOD.
The aerosol-ice formation closure pilot study (AEROICESTUDY) was a field campaign conducted at the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility at the Southern Great Plains (SGP) observatory from October 7 to 28, 2019. The purpose of this campaign was to test a field observational approach for conducting an aerosol-ice formation closure study. In other words, the goal was to predict the number of ice nucleating particles (INPs) from the measurements of physical and chemical aerosol properties and compare those to measured INP number concentrations. This project was motivated by the fact that the aerosol community has widely conducted aerosol radiative closure, aerosol-cloud condensation nuclei (CCN), and CCN-droplet closure studies to test the physical models and parameterizations that cloud-resolving and climate models rely on to perform reliable simulations of the Earth system and energy budget. However, very few closure studies related to INPs have been conducted, and to our knowledge, none using robust, size-resolved, ambient aerosol composition measurements as model inputs.
Aircraft engine combustors generally involve turbulent swirling flows in order to enhance fuel-air mixing and flame stabilization. It has long been recognized that eddy viscosity turbulence models are unable to appropriately model swirling flows. Therefore, it has been suggested that, for the modeling of these flows, a second order closure scheme should be considered because of its ability in the modeling of rotational and curvature effects. However, this scheme will require solution of many complicated second moment transport equations (six Reynolds stresses plus other scalar fluxes and variances), which is a difficult task for any CFD implementations. Also, this scheme will require a large amount of computer resources for a general combustor swirling flow. This report is devoted to the development of a cubic Reynolds stress-strain model for turbulent swirling flows, and was inspired by the work of Launder's group at UMIST. Using this type of model, one only needs to solve two turbulence equations, one for the turbulent kinetic energy k and the other for the dissipation rate epsilon. The cubic model developed in this report is based on a general Reynolds stress-strain relationship. Two flows have been chosen for model evaluation. One is a fully developed rotating pipe flow, and the other is a more complex flow with swirl and recirculation.
Closed form analytical expressions are obtained for predicting the limited rate of reactant conversion in a binary reaction of the type F + rO yields (1 + r) Product in unpremixed homogeneous turbulence. These relations are obtained by means of a single point Probability Density Function (PDF) method based on the Amplitude Mapping Closure. It is demonstrated that with this model, the maximum rate of the reactants' decay can be conveniently expressed in terms of definite integrals of the Parabolic Cylinder Functions. For the cases with complete initial segregation, it is shown that the results agree very closely with those predicted by employing a Beta density of the first kind for an appropriately defined Shvab-Zeldovich scalar variable. With this assumption, the final results can also be expressed in terms of closed form analytical expressions which are based on the Incomplete Beta Functions. With both models, the dependence of the results on the stoichiometric coefficient and the equivalence ratio can be expressed in an explicit manner. For a stoichiometric mixture, the analytical results simplify significantly. In the mapping closure, these results are expressed in terms of simple trigonometric functions. For the Beta density model, they are in the form of Gamma Functions. In all the cases considered, the results are shown to agree well with data generated by Direct Numerical Simulations (DNS). Due to the simplicity of these expressions and because of nice mathematical features of the Parabolic Cylinder and the Incomplete Beta Functions, these models are recommended for estimating the limiting rate of reactant conversion in homogeneous reacting flows. These results also provide useful insights in assessing the extent of validity of turbulence closures in the modeling of unpremixed reacting flows. Some discussions are provided on the extension of the model for treating more complicated reacting systems including realistic kinetics schemes and multi-scalar mixing with finite rate chemical reactions in more complex configurations.
This viewgraph presentation discusses project description, turbulence models, and computational engine and results for second-order closures for compressible turbulence.
The objective of the paper is to investigate the nature and values of closure parameters appearing in the proposed Large-Eddy Interaction Model for prediction of turbulent flow field. Effects of two closure parameters on predicted Reynolds stresses and other turbulence structural quantities are examined for channel flows at two Reynolds numbers.
The paper presents a multiple-time-scale turbulence model of a single point closure and a simplified split-spectrum method. Consideration is given to a class of turbulent boundary layer flows and of separated and/or swirling elliptic turbulent flows. For the separated and/or swirling turbulent flows, the present turbulence model yielded significantly improved computational results over those obtained with the standard k-epsilon turbulence model.
The implementation of the k-kL turbulence model using multiple computational uid dy- namics (CFD) codes is reported herein. The k-kL model is a two-equation turbulence model based on Abdol-Hamid's closure and Menter's modi cation to Rotta's two-equation model. Rotta shows that a reliable transport equation can be formed from the turbulent length scale L, and the turbulent kinetic energy k. Rotta's equation is well suited for term-by-term mod- eling and displays useful features compared to other two-equation models. An important di erence is that this formulation leads to the inclusion of higher-order velocity derivatives in the source terms of the scale equations. This can enhance the ability of the Reynolds- averaged Navier-Stokes (RANS) solvers to simulate unsteady ows. The present report documents the formulation of the model as implemented in the CFD codes Fun3D and CFL3D. Methodology, veri cation and validation examples are shown. Attached and sepa- rated ow cases are documented and compared with experimental data. The results show generally very good comparisons with canonical and experimental data, as well as matching results code-to-code. The results from this formulation are similar or better than results using the SST turbulence model.
Simulated data from the UCLA cumulus ensemble model are used to investigate the quasi-universal validity of closure assumptions used in existing cumulus parameterizations. A closure assumption is quasi-universally valid if it is sensitive neither to convective cloud regimes nor to horizontal resolutions of large-scale/mesoscale models. The dependency of three types of closure assumptions, as classified by Arakawa and Chen, on the horizontal resolution is addressed in this study. Type I is the constraint on the coupling of the time tendencies of large-scale temperature and water vapor mixing ratio. Type II is the constraint on the coupling of cumulus heating and cumulus drying. Type III is a direct constraint on the intensity of a cumulus ensemble. The macroscopic behavior of simulated cumulus convection is first compared with the observed behavior in view of Type I and Type II closure assumptions using 'quick-look' and canonical correlation analyses. It is found that they are statistically similar to each other. The three types of closure assumptions are further examined with simulated data averaged over selected subdomain sizes ranging from 64 to 512 km. It is found that the dependency of Type I and Type II closure assumptions on the horizontal resolution is very weak and that Type III closure assumption is somewhat dependent upon the horizontal resolution. The influences of convective and mesoscale processes on the closure assumptions are also addressed by comparing the structures of canonical components with the corresponding vertical profiles in the convective and stratiform regions of cumulus ensembles analyzed directly from simulated data. The implication of these results for cumulus parameterization is discussed.