Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ensemble simulations”

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 379 records · Page 21

Analysis of instrumentation error effects on the identification accuracy of aircraft parameters

An analytical investigation is presented of the effect of unmodeled measurement system errors on the accuracy of aircraft stability and control derivatives identified from flight test data. Such error sources include biases, scale factor errors, instrument position errors, misalignments, and instrument dynamics. Two techniques (ensemble analysis and simulated data analysis) are formulated to determine the quantitative variations to the identified parameters resulting from the unmodeled instrumentation errors. The parameter accuracy that would result from flight tests of the F-4C aircraft with typical quality instrumentation is determined using these techniques. It is shown that unmodeled instrument errors can greatly increase the uncertainty in the value of the identified parameters. General recommendations are made of procedures to be followed to insure that the measurement system associated with identifying stability and control derivatives from flight test provides sufficient accuracy.

Sorensen, J. A.↗

A numerical study of the vertical transport of momentum in a tropical rainband

The vertical transport of horizontal momentum in a convective tropical rainband is studied using a two-dimensional cloud ensemble model. Twelve simulations are made under the same large-scale conditions. The vertical transports of v momentum (parallel to the rainband) are essentially the same in all of the simulations, even though the structure of the clouds is different in each of the runs. The magnitude of the v-momentum transport by clouds is fairly large. It takes only half of a day to smooth out the tropical low-level easterly jet parallel to the rainband if no other processes are operating. The vertical transports of u momentum (perpendicular to the rainband) are quite different in all of the simulations. This difference can be explained by the dissimilarities in the distributions of horizontal momentum associated with various cloud configurations. The simulated vertical transports of horizontal momentum are compared with those computed with the Schneider and Lindzen scheme. The results suggest that their scheme is basically correct and usable if some improvements are made.

Soong, S.-T.↗

Statistical comparison of galaxy formation models - The bispectrum

To discriminate between galaxy-formation models that give similar two-point correlation functions, higher moments, specifically the three-point function in transform space or bispectrum, are looked to for an ensemble of numerical simulations. Different behaviors of this statistic are found that can be used to distinguish between classes of models. Up to the uncertainty in identifying the galaxy distribution with the mass distribution in the models, the results seem to prefer the isothermal or the cold-particle models. No model, however, agrees in all detail with the observational bispectrum of the Shane-Wirtanen galaxy catalog.

Fry, J. N.↗

Numerical and experimental investigation of nonswirling and swirling confined jets

An investigation of the influence of large scale structures on the flow development for coaxial jets with sudden expansion (with and without swirl) is presented. Both an experimental study and numerical predictions were performed for a configuration corresponding to that considered by Johnson and Bennett and Roback and Johnson. The effects of large scale structures on the swirling and nonswirling flows are documented, in particular their influence on turbulence modeling and the numerical simulation. The ensemble-averaged, time-dependent Navier-Stokes equations are solved by an LBI procedure to predict the turbulent flow field. Effects of artificial dissipation and placement of the upstream boundary in the numerical computation are also discussed.

Brondum, D. C.↗

Viral coefficient and hidden mass in the galaxy groups

The purpose is the verification of the virial mass estimations for small galaxy groups. The dynamical evolution of triple and quintuple galaxies was studied by the numerical simulations. The dependence of the virial coefficient k(t) versus time was derived. Initial k(O) = O. The function k(t) has some strong oscillations from 0.02 to 0.99. Generally, these oscillations are quasiperiodical ones. Such a behavior of k(t) is caused by formation in a system of close isolated temporary double subsystems. A strong correlation between the virial coefficient and the least mutual distance in the system is observed. Such wide oscillations may add into the estimation of virial mass of the galaxy groups an uncertainty of more than one order. An additional uncertainty is introduced by the projection effect. This uncertainty for the individual estimations of the masses approach three orders. Thus any individual estimation of the virial mass is impossible for small galaxy groups. Some possibility of statistical estimation (median or average) of the total mass, including a hidden mass, is shown for the homogeneous samples. The authors propose a method for these estimations based on a comparison of the medians of dynamical parameters (a mean size in projection and a dispersion of relative radial velocities) for the simulated and observed ensembles of the galaxy groups. This method has been applied to a sample of 46 probably physical triplets of galaxies. The probable median of the hidden mass in a volume of the triplet is about 4 M, where M is the total mass of visible matter.

Anosova, Joanna P.↗

Theoretical investigation of wave-vector-dependent analytical and numerical formulations of the interband impact-ionization transition rate for electrons in bulk silicon and GaAs

The electron interband impact-ionization rate for both silicon and gallium arsenide is calculated using an ensemble Monte Carlo simulation with the expressed purpose of comparing different formulations of the interband ionization transition rate. Specifically, three different treatments of the transition rate are examined: the traditional Keldysh formula, a new k-dependent analytical formulation first derived by W. Quade, E. Scholl, and M. Rudan (1993), and a more exact, numerical method of Y. Wang and K. F. Brennan (1994). Although the completely numerical formulation contains no adjustable parameters and as such provides a very reliable result, it is highly computationally intensive. Alternatively, the Keldysh formular, although inherently simple and computationally efficient, fails to include the k dependence as well as the details of the energy band structure. The k-dependent analytical formulation of Quade and co-workers overcomes the limitations of both of these models but at the expense of some new parameterization. It is found that the k-dependent analytical method of Quade and co-workers produces very similar results to those obtained with the completely numerical model for some quantities. Specifically, both models predict that the effective threshold for impact ionization in GaAs and silicon is quite soft, that the majority of ionization events originate from the second conduction band in both materials, and that the transition rate is k dependent. Therefore, it is concluded that the k-dependent analytical model can qualitatively reproduce results similar to those obtained with the numerical model yet with far greater computational efficiency. Nevertheless, there exist some important drawbacks to the k-dependent analytical model of Quade and co-workers: These are that it does not accurately reproduce the quantum yield data for bulk silicon, it requires determination of a new parameter, related physically to the overlap intergrals of the Bloch state which can only be adjusted by comparison to experiment, and fails to account for any wave-vector dependence of the overlap integrals. As such the transition rate may be overestimated at those points for which 'near vertical,' small change in k, transitions occur.

Kolnik, Jan↗

Theoretical study of hole initiated impact ionization in bulk silicon and GaAs using a wave-vector-dependent numerical transition rate formulation within an ensemble Monte Carlo calculation

In this paper, calculations of the hole initiated interband impact ionization rate in bulk silicon and GaAs are presented based on an ensemble Monte Carlo simulation with the inclusion of a wave-vector-dependent numerical transition rate formulation. The ionization transition rate is determined for each of the three valence bands, heavy, light, and split-off, using Fermi's golden rule with a two-body, screened Coulomb interaction. The dielectric function used within the calculation is assumed to be wave-vector-dependent. Calculations of the field-dependent impact ionization rate as well as the quantum yield are presented. It is found from both the quantum yield results and examination of the hole distribution function that the effective threshold energy for hole initiated impact ionization is relatively soft, similar to that predicted for the corresponding electron initiated ionization events occur more frequently than either heavy or split-offf initiated ionization events in bulk silicon over the applied electric field strengths examined here, 250-500 kV/cm. Conversely,in GaAs, the vast majority of hole initated ionization events originate from holes within the split-off band.

Oguzman, Ismail H.↗

Ensemble Monte Carlo calculation of the hole initiated impact ionization rate in bulk GaAs and silicon using a k-dependent, numerical transition rate formulation

The hole initiated impact ionization rate in bulk silicon and GaAs is calculated using a numerical formulation of the impact ionization transition rate incorporated into an ensemble Monte Carlo simulation. The transition rate is calculated from Fermi's golden rule using a two-body screened Coulomb interaction including a wavevector dependent dielectric function. It is found that the effective threshold for hole initiated ionization is relatively soft in both materials, that the split-off band dominates the ionization process in GaAs. and that no clear dominance by any one band is observed in silicon, though the rate out of the light hole band is greatest.

Oguzman, Ismail H.↗

Semiclassical Study of the Wave Vector Dependence of the Interband Impact Ionization Rate in Bulk Silicon

We present calculations of the interband impact ionization rate calculated using a wave vector dependent (k-dependent) semiclassical formulation of the transition rate. The transition rate is determined using Fermi's golden rule from a two-body screened Coulomb interaction assuming energy and momentum conservation. The transition rate is calculated for the first two conduction bands of silicon by numerically integrating over the full Brillouin zone. The overlap integrals in the expression for the transition rate are determined numerically using a 15 band k-p calculation. It is found that the transition rate depends strongly on the initiating electron wave vector (k vector) and that the transition rate is greatest for electrons originating within the second conduction band than the first conduction band. An ensemble Monte Carlo simulation, which includes the numerically determined ionization transition rate as well as the full details of the first two conduction bands, is used to calculate the total impact ionization rate in bulk silicon. Good agreement with the experimentally determined electron ionization rate data is obtained.

Wang, Yang↗

Design Issues of GaAs and AlGaAs Delta-Doped p-i-n Quantum-Well APD's

We examine the basic design issues in the optimization of GaAs delta-doped and AlGAs delta-doped quantum-well avalanche photodiode (APD) structures using a theoretical analysis based on an ensemble Monte Carlo simulation. The devices are variations of the p-i-n doped quantum-well structure previously described in the literature. They have the same low-noise, high-gain and high-bandwidth features as the p-i-n doped quantum-well device. However, the use of delta doping provides far greater control or the doping concentrations within each stage possibly enhancing the extent to which the device can be depleted. As a result, it is expected that the proposed devices will operate at higher gain levels (at very low noise) than devices previously developed.

Wang, Yang↗

Theoretical Investigation of Wave-Vector-Dependent Analytical and Numerical Formulations of the Interband Impact-Ionization Transition Rate for Electron in Bulk Silicon and GaAs

The electron interband impact-ionization rate for both silicon and gallium arsenide is calculated using an ensemble Monte Carlo simulation with the expressed purpose of comparing different formulations of the interband ionization transition rate. Specifically, three different treatments of the transition rate are examined: the traditional Keldysh formula, a new k-dependent analytical formulation first derived by W. Quade, E Scholl, and M. Rudan, and a more exact, numerical method of Y. Wang and K. F. Brennan. Although the completely numerical formulation contains no adjustable parameters and as such provides a very reliable result, it is highly computationally intensive. Alternatively, the Keldysh formula, although inherently simple and computationally efficient, fails to include the k dependence as well as the details of the energy band structure. The k-dependent analytical formulation of Quade and co-workers overcomes the limitations of both of these models but at the expense of some new parameterization. It is found that the k-dependent analytical method of Quade and co-workers produces very similar results to those obtained with (he completely numerical model for some quantities. Specifically, both models predict that the effective threshold for impact ionization in GaAs and silicon is quite soft, that the majority of ionization events originate from the second conduction band in both materials, and that the transition rate is k dependent. Therefore, it is concluded that the k-dependent analytical model can qualitatively reproduce results similar to those obtained with the numerical model yet with far greater computational efficiency. Nevertheless, there exist some important drawbacks to the k-dependent analytical model of Quade and co-workers: These are that it does not accurately reproduce the quantum yield data for bulk silicon, it requires determination of a new parameter, related physically to (he overlap integrals of the Bloch state which can only be adjusted by comparison to experiment, and fails to account for any wave-vector dependence of the overlap integrals. As such [he transition rate may be overestimated at those points for which "near vertical," small change in k, transitions occur.

Kolnik, Jan↗

Influence of the Biosphere on Precipitation: July 1995 Studies with the ARM-CART Data

Ensemble sets of simulation experiments were conducted with a single column model (SCM) using the Goddard GEOS II GCM physics containing a recent version of the Cumulus Scheme (McRAS) and a biosphere based land-fluxes scheme (SSiB). The study used the 18 July to 5 August 1995 ARM-CART (Atmospheric Radiation Measurement-Cloud Atmospheric Radiation Test-bed) data, which was collected at the ARM-CART site in the mid-western United States and analyzed for single column modeling (SCM) studies. The new findings affirm the earlier findings that the vegetation, which increases the solar energy absorption at the surface together with soil and soil-moisture dependent processes, which modulate the surface, fluxes (particularly evapotranspiration) together help to increase the local rainfall. In addition, the results also show that for the particular study period roughly 50% of the increased evaporation over the ARM-CART site would be converted into rainfall with the Column, while the remainder would be advected out to the large-scale. Notwithstanding the limitations of only one-way interaction (i.e., the large-scale influencing the regional physics and not vice versa), the current SCM simulations show a very robust relationship. The evaporation-precipitation relationship turns out to be independent of the soil types, and soil moisture; however, it is weakly dependent on the vegetation cover because of its surface-albedo effect. Clearly, these inferences are prone to weaknesses of the SCM physics, the assumptions of the large-scale being unaffected by gridscale (SCM-scale) changes in moist processes, and other limitations of the evaluation procedures.

Sud, Y. C.↗

An optimal modification of a Kalman filter for time scales

The Kalman filter in question, which was implemented in the time scale algorithm TA(NIST), produces time scales with poor short-term stability. A simple modification of the error covariance matrix allows the filter to produce time scales with good stability at all averaging times, as verified by simulations of clock ensembles.

time scale Kalman filter↗

Investigation of Long-Term Impacts of Urbanization when Considering Global Warming for a Coastal Tropical Region

The overachieving goal of this project is to gain a better understanding of the climate impacts caused by the combined effects of land cover and land use (LCLU) changes and increasing global concentrations of green house gases (GHG) in tropical coastal areas, regions where global, regional and local climate phenomena converge, taking as the test case the densely populated northeast region of the Caribbean island of Puerto Rico. The research uses an integrated approach of high-resolution remote sensing information linked to a high resolution Regional Atmospheric Modeling System (RAMS), which was employed to perform ensembles of climate simulations (combining 2-LCLU and 2-GHG concentration scenarios). Reconstructed agricultural maps are used to define past LCLU, and combined with reconstructed sea surface temperatures (SST) for the same period form the PAST climate scenario (1951-1956); while the PRESENT scenario (2000-2004) was additionally supported by high resolution remote sensing data (10-m-res). The climate reconstruction approach is validated with available observed climate data from surface weather stations for both periods of time simulated. The selection of the past and present climate scenarios considers large-scale biases (i.e. ENSO/NAO) as reflected in the region of interest. Direct and cross comparison of the results is allowing quantifying single, combined, and competitive effects. Results indicate that global GHG have dominant effects on minimum temperatures (following regional tendencies), while urban sprawl dominates maximum temperatures. To further investigate impacts of land use the Bowen Ratio and the thermal response number (TRN) are analyzed. The Bowen ratio indicates that forestation of past agricultural high areas have an overwhelmingly mitigation effect on increasing temperatures observed in different LCLU scenarios, but when abandoned agricultural lands are located in plains, the resulting shrub/grass lands produce higher surface temperatures. The TRN (J/m^2/degC) is a surface property defined as the ratio of the surface net radiation to the rate of change in surface temperature, expresses how those fluxes are reacting to radiant energy inputs. Natural vegetated surfaces have a greater TRN than urban and barren surfaces because the net radiation processed by them is mostly used for latent heat and thermal storage heat rather than sensible heat (heating the air). Significant changes in TRN were observed in the metropolitan area of San Juan for the two analyzed periods reflecting a reduction of this variable in the present from the past consistent with increasing in thermal mass, or intense urbanization.

Gonalez, Jorge E.↗

Five Decades of Northern Land Carbon Uptake Revealed by the Interhemispheric CO2 Gradient

The global land and ocean carbon sinks have increased proportionally with increasing carbon dioxide emissions during the past decades. It is thought that Northern Hemisphere lands make a dominant contribution to the global land carbon sink; however, the long-term trend of the northern land sink remains uncertain. Here, using measurements of the interhemispheric gradient of atmospheric carbon dioxide from 1958 to 2016, we show that the northern land sink remained stable between the 1960s and the late 1980s, then increased by 0.5 ± 0.4 petagrams of carbon per year during the 1990s and by 0.6 ± 0.5 petagrams of carbon per year during the 2000s. The increase of the northern land sink in the 1990s accounts for 65% of the increase in the global land carbon flux during that period. The subsequent increase in the 2000s is larger than the increase in the global land carbon flux, suggesting a coincident decrease of carbon uptake in the Southern Hemisphere. Comparison of our findings with the simulations of an ensemble of terrestrial carbon models over the same period suggests that the decadal change in the northern land sink between the 1960s and the 1990s can be explained by a combination of increasing concentrations of atmospheric carbon dioxide, climate variability and changes in land cover. However, the increase during the 2000s is underestimated by all models, which suggests the need for improved consideration of changes in drivers such as nitrogen deposition, diffuse light and land-use change. Overall, our findings underscore the importance of Northern Hemispheric land as a carbon sink.

Carbon cycle↗

The MJO-QBO Relationship in a GCM with Stratospheric Nudging

Observational studies show a strong connection between the intraseasonal Madden-Julian oscillation (MJO) and the stratospheric quasi-biennial oscillation (QBO): the boreal winter MJO is stronger, more predictable, and has different teleconnections when the QBO in the lower stratosphere is easterly versus westerly. Despite the strength of the observed connection, global climate models do not produce an MJO-QBO link. Here the authors use a current-generation ocean-atmosphere coupled NASA Goddard Institute for Space Studies global climate model (Model E2.1) to examine the MJO-QBO link. To represent the QBO with minimal bias, the model zonal mean stratospheric zonal and meridional winds are relaxed to reanalysis fields from 1980-2017. The model troposphere, including the MJO, is allowed to freely evolve. The model with stratospheric nudging captures QBO signals well, including QBO temperature anomalies. However, an ensemble of nudged simulations still lacks an MJO-QBO connection.

Madden-Julian oscillation (MJO)↗

Role of forest regrowth in global carbon sink dynamics

Although the existence of a large carbon sink in terrestrial ecosystems is well-established, the drivers of this sink remain uncertain. It has been suggested that perturbations to forest demography caused by past land-use change, management, and natural disturbances may be causing a large component of current carbon uptake. Here we use a global compilation of forest age observations, combined with a terrestrial biosphere model with explicit modeling of forest regrowth, to partition the global forest carbon sink between old-growth and regrowth stands over the period 1981–2010. For 2001–2010 we find a carbon sink of 0.85 (0.66–0.96) Pg/year located in intact old-growth forest, primarily in the moist tropics and boreal Siberia, and 1.30 (1.03–1.96) Pg/year located in stands regrowing after past disturbance. Approaching half of the sink in regrowth stands would have occurred from demographic changes alone, in the absence of other environmental changes. These age-constrained results show consistency with those simulated using an ensemble of demographically-enabled terrestrial biosphere models following an independent reconstruction of historical land use and management. We estimate that forests will accumulate an additional 69 (44–131) Pg C in live biomass from changes in demography alone if natural disturbances, wood harvest, and reforestation continue at rates comparable to those during 1981–2010. Our results confirm that it is not possible to understand the current global terrestrial carbon sink without accounting for the sizeable sink due to forest demography. They also imply that a large portion of the current terrestrial carbon sink is strictly transient in nature.

forest↗

Comparison of Dust Optical Depth From Multi-Sensor Products and MONARCH (Multiscale Online Non-hydrostatic AtmospheRe CHemistry) Dust Reanalysis Over North Africa, the Middle East, and Europe

Aerosol reanalysis datasets are model-based, observationally constrained, continuous 3D aerosol fields with a relatively high temporal frequency that can be used to assess aerosol variations and trends, climate effects, and impacts on socioeconomic sectors, such as health. Here we compare and assess the recently published MONARCH (Multiscale Online Non-hydrostatic AtmospheRe CHemistry) high-resolution regional desert dust reanalysis over northern Africa, the Middle East, and Europe (NAMEE) with a combination of ground-based observations and space-based dust retrievals and products. In particular, we compare the total and coarse dust optical depth (DOD) from the new reanalysis with DOD products derived from MODIS (MODerate resolution Imaging Spectroradiometer), MISR (Multi-angle Imaging SpectroRadiometer), and IASI (Infrared Atmospheric Sounding Interferometer) spaceborne instruments. Despite the larger uncertainties, satellite-based datasets provide a better geographical coverage than ground-based observations, and the use of different retrievals and products allows at least partially overcoming some single-product weaknesses in the comparison. Nevertheless, limitations and uncertainties due to the type of sensor, its operating principle, its sensitivity, its temporal and spatial resolution, and the methodology for retrieving or further deriving dust products are factors that bias the reanalysis assessment. We, therefore, also use ground-based DOD observations provided by 238 stations of the AERONET (AErosol RObotic NETwork) located within the NAMEE region as a reference evaluation dataset. In particular, prior to the reanalysis assessment, the satellite datasets were evaluated against AERONET, showing moderate underestimations in the vicinities of dust sources and downwind regions, whereas small or significant overestimations, depending on the dataset, can be found in the remote regions. Taking these results into consideration, the MONARCH reanalysis assessment shows that total and coarse-DOD simulations are consistent with satellite- and ground-based data, qualitatively capturing the major dust sources in the area in addition to the dust transport patterns. Moreover, the MONARCH reanalysis reproduces the seasonal dust cycle, identifying the increased dust activity that occurred in the NAMEE region during spring and summer. The quantitative comparison between the MONARCH reanalysis DOD and satellite multi-sensor products shows that the reanalysis tends to slightly overestimate the desert dust that is emitted from the source regions and underestimate the transported dust over the outflow regions, implying that the model's removal of dust particles from the atmosphere, through deposition processes, is too effective. More specifically, small positive biases are found over the Sahara desert (0.04) and negative biases over the Atlantic Ocean and the Arabian Sea (−0.04), which constitute the main pathways of the long-range dust transport. Considering the DOD values recorded on average there, such discrepancies can be considered low, as the low relative bias in the Sahara desert (< 50 %) and over the adjacent maritime regions (< 100 %) certifies. Similarly, over areas with intense dust activity, the linear correlation coefficient between the MONARCH reanalysis simulations and the ensemble of the satellite products is significantly high for both total and coarse DOD, reaching 0.8 over the Middle East, the Atlantic Ocean, and the Arabian Sea and exceeding it over the African continent. Moreover, the low relative biases and high correlations are associated with regions for which large numbers of observations are available, thus allowing for robust reanalysis assessment.

Michail Mytilinaios↗