A simplified molecular model for studying vibration-dissociation coupling in fluid flows
Molecular relaxation model for vibration dissociation coupling effects in fluid flows
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.
Molecular relaxation model for vibration dissociation coupling effects in fluid flows
A coupled Boundary Element/Finite Volume Method temperature-forward/flux-hack algorithm is developed for conjugate heat transfer (CHT) applications. A loosely coupled strategy is adopted with each field solution providing boundary conditions for the other in an iteration seeking continuity of temperature and heat flux at the fluid-solid interface. The NASA Glenn Navier-Stokes code Glenn-HT is coupled to a 3-D BEM steady state heat conduction code developed at the University of Central Florida. Results from CHT simulation of a 3-D film-cooled blade section are presented and compared with those computed by a two-temperature approach. Also presented are current developments of an iterative domain decomposition strategy accommodating large numbers of unknowns in the BEM. The blade is artificially sub-sectioned in the span-wise direction, 3-D BEM solutions are obtained in the subdomains, and interface temperatures are averaged symmetrically when the flux is updated while the fluxes are averaged anti-symmetrically to maintain continuity of heat flux when the temperatures are updated. An initial guess for interface temperatures uses a physically-based 1-D conduction argument to provide an effective starting point and significantly reduce iteration. 2-D and 3-D results show the process converges efficiently and offers substantial computational and storage savings. Future developments include a parallel multi-grid implementation of the approach under MPI for computation on PC clusters.
Three model reference adaptive controllers (MRAC) with varying levels of complexity were evaluated on a high performance jet aircraft and compared along with a baseline nonlinear dynamic inversion controller. The handling qualities and performance of the controllers were examined during failure conditions that induce coupling between the pitch and roll axes. Results from flight tests showed with a roll to pitch input coupling failure, the handling qualities went from Level 2 with the baseline controller to Level 1 with the most complex MRAC tested. A failure scenario with the left stabilator frozen also showed improvement with the MRAC. Improvement in performance and handling qualities was generally seen as complexity was incrementally added; however, added complexity usually corresponds to increased verification and validation effort required for certification. The tradeoff between complexity and performance is thus important to a controls system designer when implementing an adaptive controller on an aircraft. This paper investigates this relation through flight testing of several controllers of vary complexity.
The University of California, Berkeley, extreme ultraviolet spectrometer aboard the U.S. Air Force STP 78-1 satellite measured emission features in the Earth's dayglow due to neutral and ionized species in the atmosphere, in the 35 to 140-nm range. The spectrometer was operating between March 1979 and March 1980, including the period of the magnetic storm on March 22, 1979. Some of these measurements are interpreted using the predictions of the three-dimensional time-dependent coupled ionosphere-thermosphere model developed at University College, London. The observations show a reduction in the atomic oxygen 130.4-nm airglow emission at high northern latitudes following the storm. Model simulations show that this reduction in 130.4-nm emission is associated with an increase in the O2/O ratio. Analysis of model results using electron transport and radiative transport codes show that the brightness of 130.4-nm emission at high latitudes due to resonantly scattered sunlight is approximately twice that due to photoelectron impact excitation. However, the observed decrease in the brightness at high northern latitudes is mainly due to a change in the photoelectron impact source, which contributes approximately 75% of the total, as well as its multiple scattering component; for the photoelectron impact source at 70 deg latitude and 200 km altitude, the reduction in multiple scattering is 1.5 times greater than the reduction in the initial excitation. The reduction in the airglow emission is visible only in the norther n hemisphere because the south pole was not sunlit over the storm period. The comparison of model results with observations suggests that 130.4-nm emission may be useful as a tracer for global changes in the concentration of atomic energy.
In many real-world pollution chemical reaction problems, the rate of reaction problems, the rate of reaction may be greatly affected by unmixedness. An approximate closure scheme for a chemical kinetic submodel which conforms to the principles of invariant modeling and which accounts for the effects of inhomogeneous mixing over a wide range of conditions has been developed. This submodel has been coupled successfully with invariant turbulence and diffusion models, permitting calculation of two-dimensional diffusion of two reacting (isothermally) chemical species. The initial calculations indicate the ozone reactions in the wake of stratospheric aircraft will be substantially affected by the rate of diffusion of ozone into the wake, and in the early wake, by unmixedness.
Results of research aimed at improving the predictability of off nominal internal ballistics performance of solid propellant rocket motors (SRMs) including thrust imbalance between two SRMs firing in parallel are reported. The potential effects of nozzle throat erosion on internal ballistic performance were studied and a propellant burning rate low postulated. The propellant burning rate model when coupled with the grain deformation model permits an excellent match between theoretical results and test data for the Titan IIIC, TU455.02, and the first Space Shuttle SRM (DM-1). Analysis of star grain deformation using an experimental model and a finite element model shows the star grain deformation effects for the Space Shuttle to be small in comparison to those of the circular perforated grain. An alternative technique was developed for predicting thrust imbalance without recourse to the Monte Carlo computer program. A scaling relationship used to relate theoretical results to test results may be applied to the alternative technique of predicting thrust imbalance or to the Monte Carlo evaluation. Extended investigation into the effect of strain rate on propellant burning rate leads to the conclusion that the thermoelastic effect is generally negligible for both steadily increasing pressure loads and oscillatory loads.
The effect of treatment of incoming solar radiation on the photochemistry of the troposphere is discussed. A one dimensional photochemical model of the troposphere containing the species of the nitrogen, oxygen, carbon, hydrogen, and sulfur families was developed. The vertical flux is simulated by use of the parameterized eddy diffusion coefficients. The photochemical model is coupled to a radiative transfer model that calculates the radiation field due to the incoming solar radiation which initiates much of the photochemistry of the troposphere. Vertical profiles of tropospheric species were compared with the Leighton approximation, radiative transfer, matrix inversion model. The radiative transfer code includes the effects of multiple scattering due to molecules and aerosols, pure absorption, and surface albedo on the transfer of incoming solar radiation. It is indicated that significant differences exist for several key photolysis frequencies and species number density profiles between the Leighton approximation and the profiles generated with, radiative transfer, matrix inversion technique. Most species show enhanced vertical profiles when the more realistic treatment of the incoming solar radiation field is included
The specular reflection process is shown to be a key aspect of radiation transfer by plant canopies. Polarization measurements are demonstrated as the tool for determining the specular and diffuse portions of the canopy radiance. The magnitude of the specular fraction of the reflectance is significant compared to the magnitude of the diffuse fraction. Therefore, it is necessary to consider specularly reflected light in developing and evaluating light-canopy interaction models for wheat canopies. Models which assume leaves are diffuse reflectors correctly predict only the diffuse fraction of the canopy reflectance factor. The specular reflectance model, when coupled with a diffuse leaf model, would predict both the specular and diffuse portions of the reflectance factor. The specular model predicts and the data analysis confirms that the single variable, angle of incidence of specularly reflected sunlight on the leaf, explains much of variation in the polarization data as a function of view-illumination directions.
The present use of the Rice convection model to investigate the electrodynamic coupling of the thermosphere to the inner magnetosphere encompasses the effects of EUV-driven and convection-driven neutral winds under quasi-equilibrium conditions. Convection-driven winds are included self-consistently and interactively; a steady-state wind parameterization is written analytically in terms of the electrostatic potential, which is in turn included in a closed-loop calculation for the electric potential itself. The simulations conducted show that, as the neutral system approaches a quasi-equilibrium state, the neutral winds play a much more significant role.
A one-dimensional, atmospheric boundary layer model is coupled to a thermodynamic ice model to estimate the surface turbulent fluxes over thick sea ice. The principal forcing parameters in this time-dependent model are the air temperature, humidity, and wind speed at a specified level (either at 2 m or at 850 mb) and the downwelling surface radiative fluxes. The free parameters. are the air temperature, humidity, and wind speed profiles below the specified level, the surface skin temperature, the ice temperature profile, and the surface turbulent fluxes. The goal is to determine how well we can estimate the turbulent surface heat and momentum fluxes using forcing parameters from atmospheric temperatures and radiative fluxes retrieved from the TIROS-N Operational Vertical Sounder (TOVS) data. Meteorological observations from the Lead Experiment (LeadEx, April 1992) ice camp are used to validate turbulent fluxes computed with the surface observations and the results are used to compare with estimates based on radio-sonde observations or with estimates based on TOVS data. We find that the TOVS-based estimates of the stress are significantly more accurate than those found with a constant geostrophic drag coefficient, with a root-mean-square error about half as large. This improvement is due to stratification effects included in the boundary layer model. The errors in the sensible heat flux estimates, however, are large compared to the small mean values observed during the field experiment.
This effort provides an initial glimpse at NASA capabilities available in predicting the scattering of fan noise from a non-conventional aircraft configuration. The Aircraft NOise Prediction Program, Fast Scattering Code, and the Rotorcraft Noise Model were coupled to provide increased fidelity models of scattering effects on engine fan noise sources. The integration of these codes led to the identification of several keys issues entailed in applying such multi-fidelity approaches. In particular, for prediction at noise certification points, the inclusion of distributed sources leads to complications with the source semi-sphere approach. Computational resource requirements limit the use of the higher fidelity scattering code to predict radiated sound pressure levels for full scale configurations at relevant frequencies. And, the ability to more accurately represent complex shielding surfaces in current lower fidelity models is necessary for general application to scattering predictions. This initial step in determining the potential benefits/costs of these new methods over the existing capabilities illustrates a number of the issues that must be addressed in the development of next generation aircraft system noise prediction tools.
The sensitivity of mesoscale weather prediction model to a vegetation roughness initialization is investigated for the south central United States. Three different roughness databases are employed: i) a control or standard lookup table roughness that is a function only of land cover type, ii) a spatially heterogeneous roughness database previously derived using a physically based procedure and MODIS imagery, and iii) a MODIS climatologic roughness database that possesses the same spatial heterogeneity as (i) but with mean land class values from (ii). The model used is the Weather Research and Forecast Model (WRF) coupled to the Community Land Model within the Land Information System (LIS). For each simulation, a statistical comparison is made between modeled results and ground observations from meteorological stations within the Oklahoma mesonet and surrounding region during IHOP20O2. A sensitivity analysis on the impact the MODIS-based roughness fields is also made through a time-series intercomparison of temperature bias, probability of detection (POD), average wind speed, boundary layer height, and turbulent kinetic energy (TKE) the results that, for the current replacement of the standard land-cover type based roughness values with the satellite-derived fields statistically improves model performance for most of the observed variables. Further, the satellite-based roughness enhances the surface wind speed, PBL height and TKE production on the order of 3 to l0 percent, with a lesser effect over grassland and cropland domains, and the greater effect over mixed land cover domains
The sensitivity of a mesoscale weather prediction model to a 1 km satellite-based vegetation roughness initialization is investigated for a domain within the south central United States. Three different roughness databases are employed: i) a control or standard lookup table roughness that is a function only of land cover type, ii) a spatially heterogeneous roughness database, specific to the domain, that was previously derived using a physically based procedure and Moderate Resolution Imaging Spectroradiometer (MODIS) imagery, and iii) a MODIS climatologic roughness database that like (i) is a function only of land cover type, but possesses domain specific mean values from (ii). The model used is the Weather Research and Forecast Model (WRF) coupled to the Community Land Model within the Land Information System (LIS). For each simulation, a statistical comparison is made between modeled results and ground observations within a domain including Oklahoma, Eastern Arkansas, and Northwest Louisiana during a 4-day period within IHOP 2002. Sensitivity analysis compares the impact the three roughness initializations on time-series temperature, precipitation probability of detection (POD), average wind speed, boundary layer height, and turbulent kinetic energy (TKE). Overall, the results indicate that, for the current investigation, replacement of the standard look-up table values with the satellite-derived values statistically improves model performance for most observed variables. Such natural roughness heterogeneity enhances the surface wind speed, PBL height and TKE production up to 10 percent, with a lesser effect over grassland, and greater effect over mixed land cover domains.
A microwave radiance assimilation (RA) scheme for the retrieval of snow physical state variables requires a snowpack physical model (SM) coupled to a radiative transfer model. In order to assimilate microwave brightness temperatures (Tbs) at horizontal polarization (h-pol), an SM capable of resolving melt-refreeze crusts is required. To date, it has not been shown whether an RA scheme is tractable with the large number of state variables present in such an SM or whether melt-refreeze crust densities can be estimated. In this paper, an RA scheme is presented using the CROCUS SM which is capable of resolving melt-refreeze crusts. We assimilated both vertical (v) and horizontal (h) Tbs at 18.7 and 36.5 GHz. We found that assimilating Tb at both h-pol and vertical polarization (v-pol) into CROCUS dramatically improved snow depth estimates, with a bias of 1.4 cm compared to-7.3 cm reported by previous studies. Assimilation of both h-pol and v-pol led to more accurate results than assimilation of v-pol alone. The snow water equivalent (SWE) bias of the RA scheme was 0.4 cm, while the bias of the SWE estimated by an empirical retrieval algorithm was -2.9 cm. Characterization of melt-refreeze crusts via an RA scheme is demonstrated here for the first time; the RA scheme correctly identified the location of melt-refreeze crusts observed in situ.
The FLuorescence EXplorer (FLEX) satellite mission, selected as ESA's 8th Earth Explorer, has been designed forthe measurement of sun-induced fluorescence (F) spectra emitted by plants. This will be accomplished through amulti-sensor approach by placing it in a common orbit in tandem with the Sentinel-3 (S3) mission, which willhave two optical sensors on board, OLCI (Ocean and Land Colour Instrument) and SLSTR (Sea and Land SurfaceTemperature Radiometer) to complement FLEX. These S3 instruments will be used in combination with theimaging spectrometers on board FLEX to provide data useful for atmospheric correction of FLEX data. However,a fully synergetic approach, i.e. by exploiting the spectral and directional information from all tandem missioninstruments together, is an attractive alternative which is explored in this paper. By employing all combined topof-atmosphere (TOA) spectral radiance data, one can (i) characterize the relevant optical properties of the atmosphere,(ii) retrieve biophysical canopy properties including the associated reflectance anisotropy, and (iii)retrieve a more accurate and consistent canopy F.Regarding retrieval methods, Fraunhofer Line Depth (FLD) and Spectral Fitting (SF) are well-known techniquesapplied to hyperspectral data. Both methods depend on a high spectral resolution and assume aLambertian (isotropic) canopy reflectance. However, most vegetation canopies are non-Lambertian. This impliesthat, in particular when ignoring the anisotropic surface reflection, substantial retrieval errors can occur due tothe interaction between atmospheric absorption bands and surface reflectance anisotropy. In this paper, a novelmethod based on spectral radiative transfer (RT) modeling is proposed, in which coupled RT models are used tosimulate TOA radiance spectra. These are then matched with ‘measured' spectra in order to retrieve surfacefluorescence, along with a suite of biophysical parameters, by model inversion through optimization. By applyingcoupled RT models of the soil-leaf-canopy and the surface-atmosphere systems, TOA radiance spectra canbe simulated for all optical sensors of this tandem mission. In this way, complex effects due to surface reflectanceanisotropy and the spectral sampling by the various instruments, which are difficult to compensate for in the endproducts, are properly taken into account by their incorporation in the forward modeling. Next, by model inversionof TOA radiance data via optimization, the most accurate F retrievals can be achieved in a consistentmanner, along with important canopy level biophysical parameters that may help interpret the F spectrum, suchas chlorophyll content and leaf area index (LAI). The potential of this approach has been explored in a numericalexperiment, and the results are presented in this paper. We find that, with the assumed well-characterized andplausible FLEX/S3 instrument performances, the simultaneous retrieval of biophysical canopy parameters and Fspectra would be possible with a remarkable accuracy, provided the correct atmospheric characterization isavailable.
The third primary production algorithm round robin (PPARR3) compares output from 24 models that estimate depth-integrated primary production from satellite measurements of ocean color, as well as seven general circulation models (GCMs) coupled with ecosystem or biogeochemical models. Here we compare the global primary production fields corresponding to eight months of 1998 and 1999 as estimated from common input fields of photosynthetically-available radiation (PAR), sea-surface temperature (SST), mixed-layer depth, and chlorophyll concentration. We also quantify the sensitivity of the ocean-color-based models to perturbations in their input variables. The pair-wise correlation between ocean-color models was used to cluster them into groups or related output, which reflect the regions and environmental conditions under which they respond differently. The groups do not follow model complexity with regards to wavelength or depth dependence, though they are related to the manner in which temperature is used to parameterize photosynthesis. Global average PP varies by a factor of two between models. The models diverged the most for the Southern Ocean, SST under 10 C, and chlorophyll concentration exceeding 1mg Chlm-3. Based on the conditions under which the model results diverge most, we conclude that current ocean-color-based models are challenged by high-nutrient low-chlorophyll conditions, and extreme temperatures or chlorophyll concentrations. The GCM-based models predict comparable primary production to those based on ocean color: they estimate higher values in the Southern Ocean, at low SST, and in the equatorial band, while they estimate lower values in eutrophic regions (probably because the area of high chlorophyll concentrations is smaller in the GCMs). Further progress in primary production modeling requires improved understanding of the effect of temperature on photosynthesis and better parameterization of the maximum photosynthetic rate.
Changing amplitude of the seasonal cycle of atmospheric CO2 (SCA) in the northern hemisphere is an emerging carbon cycle property. Mauna Loa (MLO) station (20°N, 156°W), which has the longest continuous northern hemisphere CO2 record, shows an increasing SCA before the 1980s (p < .01), followed by no significant change thereafter. We analyzed the potential driving factors of SCA slowing‐down, with an ensemble of dynamic global vegetation models (DGVMs) coupled with an atmospheric transport model. We found that slowing‐down of SCA at MLO is primarily explained by response of net biome productivity (NBP) to climate change, and by changes in atmospheric circulations. Through NBP, climate change increases SCA at MLO before the 1980s and decreases it afterwards. The effect of climate change on the slowing‐down of SCA at MLO is mainly exerted by intensified drought stress acting to offset the acceleration driven by CO2 fertilization. This challenges the view that CO2 fertilization is the dominant cause of emergent SCA trends at northern sites south of 40°N. The contribution of agricultural intensification on the deceleration of SCA at MLO was elusive according to land–atmosphere CO2 flux estimated by DGVMs and atmospheric inversions. Our results also show the necessity to adequately account for changing circulation patterns in understanding carbon cycle dynamics observed from atmospheric observations and in using these observations to benchmark DGVMs.
A dynamical systems approach is used to quantify the instantaneous and time-averaged predictability of a low-order moist general circulation model. Specifically, the effects on predictability of incorporating an active ocean circulation, implementing annual solar forcing, and asynchronously coupling the ocean and atmosphere are evaluated. The predictability and structure of the model attractors is compared using the Lyapunov exponents, the local divergence rates, and the correlation, fractal, and Lyapunov dimensions. The Lyapunov exponents measure the average rate of growth of small perturbations on an attractor, while the local divergence rates quantify phase-spatial variations of predictability. These local rates are exploited to efficiently identify and distinguish subtle differences in predictability among attractors. In addition, the predictability of monthly averaged and yearly averaged states is investigated by using attractor reconstruction techniques.