Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “statistical modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,657 records · Page 92

Modeling and parameterization of horizontally inhomogeneous cloud radiative properties

One of the fundamental difficulties in modeling cloud fields is the large variability of cloud optical properties (liquid water content, reflectance, emissivity). The stratocumulus and cirrus clouds, under special consideration for FIRE, exhibit spatial variability on scales of 1 km or less. While it is impractical to model individual cloud elements, the research direction is to model a statistical ensembles of cloud elements with mean-cloud properties specified. The major areas of this investigation are: (1) analysis of cloud field properties; (2) intercomparison of cloud radiative model results with satellite observations; (3) radiative parameterization of cloud fields; and (4) development of improved cloud classification algorithms.

Welch, R. M.↗

Structure of chalcogen overlayers on Au(111): Density functional theory and lattice-gas modeling

Ordering of different chalcogens, S, Se, and Te, on Au(111) exhibit broad similarities but also some distinct features, which must reflect subtle differences in relative values of the long-range pair and many-body lateral interactions between adatoms. We develop lattice-gas (LG) models within a cluster expansion framework, which includes about 50 interaction parameters. These LG models are developed based on density functional theory (DFT) analysis of the energetics of key adlayer configurations in combination with the Monte Carlo (MC) simulation of the LG models to identify statistically relevant adlayer motifs, i.e., model development is based entirely on theoretical considerations. The MC simulation guides additional DFT analysis and iterative model refinement. Given their complexity, development of optimal models is also aided by strategies from supervised machine learning. The model for S successfully captures ordering motifs over a broader range of coverage than achieved by previous models, and models for Se and Te capture the features of ordering, which are distinct from those for S. More specifically, the modeling for all three chalcogens successfully explains the linear adatom rows (also subtle differences between them) observed at low coverages of ~0.1 monolayer. The model for S also leads to a new possible explanation for the experimentally observed phase with a (5 × 5)-type low energy electron diffraction (LEED) pattern at 0.28 ML and to predictions for LEED patterns that would be observed with Se and Te at this coverage.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Detailed and simplified nonequilibrium helium ionization in the solar atmosphere

Helium ionization plays an important role in the energy balance of the upper chromosphere and transition region. Helium spectral lines are also often used as diagnostics of these regions. We carry out one-dimensional radiation-hydrodynamics simulations of the solar atmosphere and find that the helium ionization is set mostly by photoionization and direct collisional ionization, counteracted by radiative recombination cascades. By introducing an additional recombination rate mimicking the recombination cascades, we construct a simplified three-level helium model atom consisting of only the ground states. This model atom is suitable for modeling nonequilibrium helium ionization in three-dimensional numerical models. We perform a brief investigation of the formation of the He I 10830 and He II 304 spectral lines. Both lines show nonequilibrium features that are not recovered with statistical equilibrium models, and caution should therefore be exercised when such models are used as a basis for interpretating observations.

79 ASTRONOMY AND ASTROPHYSICS↗

Entanglement transitions from restricted Boltzmann machines

The search for novel entangled phases of matter has lead to the recent discovery of a new class of “en- tanglement transitions,” exemplified by random tensor networks and monitored quantum circuits. Most known examples can be understood as some classical ordering transitions in an underlying statistical mechanics model, where entanglement maps onto the free-energy cost of inserting a domain wall. In this paper we study the possibility of entanglement transitions driven by physics beyond such statistical mechanics mappings. Motivated by recent applications of neural-network-inspired variational Ansätze, we investigate under what conditions on the variational parameters these Ansätze can capture an entanglement transition. We study the entanglement scaling of short-range restricted Boltzmann machine (RBM) quantum states with random phases. For uncorrelated random phases, we analytically demonstrate the absence of an entanglement transition and reveal subtle finite-size effects in finite-size numerical simulations. Introducing phases with correlations decaying as 1/r α in real space, we observe three regions with a different scaling of entanglement entropy depending on the exponent α. We study the nature of the transition between these regions, finding numerical evidence for critical behavior. Furthermore, our work establishes the presence of long-range correlated phases in RBM-based wave functions as a required ingredient for entanglement transitions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Measurements of the time-dependent cosmic-ray Sun shadow with seven years of IceCube data: Comparison with the Solar cycle and magnetic field models

Observations of the time-dependent cosmic-ray Sun shadow have been proven as a valuable diagnostic for the assessment of solar magnetic field models. Here, seven years of IceCube data are compared to solar activity and solar magnetic field models. A quantitative comparison of solar magnetic field models with IceCube data on the event rate level is performed for the first time. Additionally, a first energy-dependent analysis is presented and compared to recent predictions. We use seven years of IceCube data for the moon and the Sun and compare them to simulations on data rate level. The simulations are performed for the geometrical shadow hypothesis for the moon and the Sun and for a cosmic-ray propagation model governed by the solar magnetic field for the case of the Sun. We find that a linearly decreasing relationship between Sun shadow strength and solar activity is preferred over a constant relationship at the 6.4σ level. We test two commonly used models of the coronal magnetic field, both combined with a Parker spiral, by modeling cosmic-ray propagation in the solar magnetic field. Both models predict a weakening of the shadow in times of high solar activity as it is also visible in the data. We find tensions with the data on the order of 3σ for both models, assuming only statistical uncertainties. The magnetic field model CSSS fits the data slightly better than the PFSS model. This is generally consistent with what is found previously by the Tibet AS-γ Experiment; a deviation of the data from the two models is, however, not significant at this point. Regarding the energy dependence of the Sun shadow, we find indications that the shadowing effect increases with energy during times of high solar activity, in agreement with theoretical predictions.

79 ASTRONOMY AND ASTROPHYSICS↗

Observations of the flank of earth's bow shock to -110 R(E) by ISEE 3/ICE

The ISEE 3/ICE spacecraft made over 100 crossings of earth's bow shock between x(SEC) = -20 and -110 R(E) most of them further away than -60 R(E) while resetting its course for comet Giacobini-Zinner. The parameters of the crossings are tabulated and some of their magnetic, plasma wave, and spectral profiles at low magnetosonic Mach numbers are displayed. Loci of crossings are compared with the extrapolation of a model of earlier statistical success, showing that the model remains useful in the far flank.

Greenstadt, E. W.↗

Numerical studies of dilute and dense spray characteristics

Several issues involving the improvement of physical submodels and the computational efficiency in modeling dilute and dense spray combustion are discussed. First, the implementations of a dispersion width approach accounting for turbulent dispersion within each computational parcel is discussed. This is essentially a statistical transport model and the testings of this model confirm the capability of accurately representing dispersion in nearly-homogeneous and inhomogeneous turbulent flows with improved efficiency over the delta function stochastic separated flow model. To account for the dense spray effects, an existing drop collision and coalescence model and a Taylor analogy breakup (TAB) model were employed. These models were incorporated into a state-of-the-art multiphase all-speed transient flow solution procedure. Several examples including nonevaporating, evaporating, and burning dense spray cases were studied. The numerical results show reasonably good comparisons with available experimental data in terms of spray penetration, drop sizes, and overall configuration of a spray flame.

Kim, Y. M.↗

Investigation of radiative interaction in laminar flows using Monte Carlo simulation

The Monte Carlo method (MCM) is employed to study the radiative interactions in fully developed laminar flow between two parallel plates. Taking advantage of the characteristics of easy mathematical treatment of the MCM, a general numerical procedure is developed for nongray radiative interaction. The nongray model is based on the statistical narrow band model with an exponential-tailed inverse intensity distribution. To validate the Monte Carlo simulation for nongray radiation problems, the results of radiative dissipation from the MCM are compared with two available solutions for a given temperature profile between two plates. After this validation, the MCM is employed to solve the present physical problem and results for the bulk temperature are compared with available solutions. In general, good agreement is noted and reasons for some discrepancies in certain ranges of parameters are explained.

Liu, Jiwen↗

Development of a Computational Chemical Vapor Deposition Model: Applications to Indium Nitride and Dicyanovinylaniline

A computational chemical vapor deposition (CVD) model is presented, that couples chemical reaction mechanisms with fluid dynamic simulations for vapor deposition experiments. The chemical properties of the systems under investigation are evaluated using quantum, molecular and statistical mechanics models. The fluid dynamic computations are performed using the CFD-ACE program, which can simulate multispecies transport, heat and mass transfer, gas phase chemistry, chemistry of adsorbed species, pulsed reactant flow and variable gravity conditions. Two experimental setups are being studied, in order to fabricate films of: (a) indium nitride (InN) from the gas or surface phase reaction of trimethylindium and ammonia; and (b) 4-(1,1)dicyanovinyl-dimethylaminoaniline (DCVA) by vapor deposition. Modeling of these setups requires knowledge of three groups of properties: thermodynamic properties (heat capacity), transport properties (diffusion, viscosity, and thermal conductivity), and kinetic properties (rate constants for all possible elementary chemical reactions). These properties are evaluated using computational methods whenever experimental data is not available for the species or for the elementary reactions. The chemical vapor deposition model is applied to InN and DCVA. Several possible InN mechanisms are proposed and analyzed. The CVD model simulations of InN show that the deposition rate of InN is more efficient when pulsing chemistry is used under conditions of high pressure and microgravity. An analysis of the chemical properties of DCVA show that DCVA dimers may form under certain conditions of physical vapor transport. CVD simulations of the DCVA system suggest that deposition of the DCVA dimer may play a small role in the film and crystal growth processes.

Cardelino, Carlos↗

Aero-Propulsive Modeling for eVTOL Aircraft Using Wind Tunnel Testing with Multisine Inputs

A novel approach for modeling the aero-propulsive characteristics of an electric vertical takeoff and landing (eVTOL) aircraft was developed and demonstrated in wind tunnel testing. The approach was applied to the NASA LA-8 tandem tilt-wing eVTOL aircraft, using an efficient, hybrid experiment design composed of a static I-optimal response surface design for slowly-varying test variables, and dynamic orthogonal phase-optimized multisine excitations for the control surfaces and electric propulsors. Both the static and dynamic experiment designs were executed simultaneously to collect informative data for model identification. Statistically-weighted multivariate orthogonal function modeling was used to combine local modeling results computed in the frequency domain using data collected with dynamic excitations operating on the control effectors to form an aggregate aero-propulsive model. The final identified model exhibited good predictive capability when compared to validation data acquired separately from the data used to develop the model. The required test time using these new techniques was reduced by at least a factor of five compared to previous static wind tunnel testing for the LA-8 aircraft, while providing more informative data, greater parameterization flexibility, and high-quality models.

system identification↗

Global Carbon Budget 2016

Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere the global carbon budget is important to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe data sets and methodology to quantify all major components of the global carbon budget, including their uncertainties, based on the combination of a range of data, algorithms, statistics, and model estimates and their interpretation by a broad scientific community. We discuss changes compared to previous estimates and consistency within and among components, alongside methodology and data limitations. CO2 emissions from fossil fuels and industry (EFF) are based on energy statistics and cement production data, respectively, while emissions from land-use change (ELUC), mainly deforestation, are based on combined evidence from land-cover change data, fire activity associated with deforestation, and models. The global atmospheric CO2 concentration is measured directly and its rate of growth (GATM) is computed from the annual changes in concentration. The mean ocean CO2 sink (SOCEAN) is based on observations from the 1990s, while the annual anomalies and trends are estimated with ocean models. The variability in SOCEAN is evaluated with data products based on surveys of ocean CO2 measurements. The global residual terrestrial CO2 sink (SLAND) is estimated by the difference of the other terms of the global carbon budget and compared to results of independent dynamic global vegetation models. We compare the mean land and ocean fluxes and their variability to estimates from three atmospheric inverse methods for three broad latitude bands. All uncertainties are reported as +/- 1(sigma), reflecting the current capacity to characterize the annual estimates of each component of the global carbon budget. For the last decade available (2006-2015), EFF was 9.3+/-0.5 GtC/yr, ELUC 1.0+/-0.5 GtC/yr,GATM 4.5+/-0.1 GtC/yr, SOCEAN 2.6+/-0.5 GtC/yr, and SLAND 3.1+/-0.9 GtC/yr. For year 2015 alone, the growth in EFF was approximately zero and emissions remained at 9.9+/-0.5 GtC/yr, showing a slowdown in growth of these emissions compared to the average growth of 1.8/yr that took place during 2006-2015.Also, for 2015, ELUC was 1.3+/-0.5 GtC/yr, GATM was 6.3+/-0.2 GtC/yr, SOCEAN was 3.0+/-0.5 GtC/yr, and SLAND was 1.9+/-0.9 GtC/yr. GATM was higher in 2015 compared to the past decade (2006-2015), reflecting a smaller SLAND for that year. The global atmospheric CO2 concentration reached 399.4+/-0.1 ppm averaged over 2015. For 2016, preliminary data indicate the continuation of low growth in EFF with +0.2% (range of -1.0 to +1.8% ) based on national emissions projections for China and USA, and projections of gross domestic product corrected for recent changes in the carbon intensity of the economy for the rest of the world. In spite of the low growth of EFF in 2016, the growth rate in atmospheric CO2 concentration is expected to be relatively high because of the persistence of the smaller residual terrestrial sink (SLAND) in response to El Nino conditions of 2015-2016. From this projection of EFF and assumed constant ELUC for 2016, cumulative emissions of CO2 will reach 565+/-55 GtC (2075+/-205 GtCO2) for 1870-2016, about 75% from EFF and 25% from ELUC. This living data update documents changes in the methods and data sets used in this new carbon budget compared with previous publications of this data set.

Quéré, Corinne Le↗

UNIVERSAL BEHAVIOR OF X-RAY FLARES FROM BLACK HOLE SYSTEMS

X-ray flares have been discovered in black hole systems such as gamma-ray bursts, the tidal disruption event Swift J1644+57, the supermassive black hole Sagittarius A* at the center of our Galaxy, and some active galactic nuclei. Occurrences of X-ray flares are always accompanied by relativistic jets. However, it is still unknown whether or not there is a physical analogy among such X-ray flares produced in black hole systems spanning nine orders of magnitude in mass. Here, we report observed data of X-ray flares and show that they have three statistical properties similar to solar flares, including power-law distributions of their energies, durations, and waiting times, which can be explained by a fractal-diffusive, self-organized criticality model. These statistical similarities, together with the fact that solar flares are triggered by a magnetic reconnection process, suggest that all of the X-ray flares are consistent with magnetic reconnection events, implying that their concomitant relativistic jets may be magnetically dominated.

79 ASTRONOMY AND ASTROPHYSICS↗

Statistical Analysis of Factors Riving Surface Ozone Variability over Continental South Africa

Statistical relationships between surface ozone (O3) concentration, precursor species and meteorological conditions in continental South Africa were examined from data obtained from measurement stations in north-eastern South Africa. Three multivariate statistical methods were applied in the investigation, i.e. multiple linear regression (MLR), principal component analysis (PCA) and –regression (PCR), and generalised additive model (GAM) analysis. The daily maximum 8-h moving average O3 concentrations were considered in these statistical models (dependent variable). MLR models indicated that meteorology and precursor species concentrations are able to explain ~50% of the variability in daily maximum O3 levels. MLR analysis revealed that atmospheric carbon monoxide (CO), temperature and relative humidity were the strongest factors affecting the daily O3 variability. In summer, daily O3 variances were mostly associated with relative humidity, while winter O3 levels were mostly linked to temperature and CO. PCA indicated that CO, temperature and relative humidity were not strongly collinear. GAM also identified CO, temperature and relative humidity as the strongest factors affecting the daily variation of O3. Partial residual plots found that temperature, radiation and nitrogen oxides most likely have a non-linear relationship with O3,while the relationship with relative humidity and CO is probably linear. An inter-comparison between O3 levels modelled with the three statistical models compared to measured O3 concentrations showed that the GAM model offered a slight improvement over the MLR model. These findings emphasise the critical role of regional-scale O3 precursors coupled with meteorological conditions in daily variances of O3 levels in continental South Africa.

multiple linear regression (MLR)↗

Separation of Radionuclides from a Rare Earth-Containing Solution by Zeolite Adsorption

The increasing industrial demand for rare earths requires new or alternative sources to be found. Within this context, there have been studies validating the technical feasibility of coal and coal byproducts as alternative sources for rare earth elements. Nonetheless, radioactive materials, such as thorium and uranium, are frequently seen in the rare earths’ mineralization, and causes environmental and health concerns. Consequently, there exists an urgent need to remove these radionuclides in order to produce high purity rare earths to diversify the supply chain, as well as maintain an environmentally-favorable extraction process for the surroundings. In this study, an experimental design was generated to examine the effect of zeolite particle size, feed solution pH, zeolite amount, and contact time of solid and aqueous phases on the removal of thorium and uranium from the solution. The best separation performance was achieved using 2.50 g of 12-µm zeolite sample at a pH value of 3 with a contact time of 2 h. Under these conditions, the adsorption recovery of rare earths, thorium, and uranium into the solid phase was found to be 20.43 wt%, 99.20 wt%, and 89.60 wt%, respectively. The Freundlich adsorption isotherm was determined to be the best-fit model, and the adsorption mechanism of rare earths and thorium was identified as multilayer physisorption. Further, the separation efficiency was assessed using the response surface methodology based on the development of a statistically significant model.

58 GEOSCIENCES↗

NMME Monthly / Seasonal Forecasts for NASA SERVIR Applications Science

This work details use of the North American Multi-Model Ensemble (NMME) experimental forecasts as drivers for Decision Support Systems (DSSs) in the NASA / USAID initiative, SERVIR (a Spanish acronym meaning "to serve"). SERVIR integrates satellite observations, ground-based data and forecast models to monitor and forecast environmental changes and to improve response to natural disasters. Through the use of DSSs whose "front ends" are physically based models, the SERVIR activity provides a natural testbed to determine the extent to which NMME monthly to seasonal projections enable scientists, educators, project managers and policy implementers in developing countries to better use probabilistic outlooks of seasonal hydrologic anomalies in assessing agricultural / food security impacts, water availability, and risk to societal infrastructure. The multi-model NMME framework provides a "best practices" approach to probabilistic forecasting. The NMME forecasts are generated at resolution more coarse than that required to support DSS models; downscaling in both space and time is necessary. The methodology adopted here applied model output statistics where we use NMME ensemble monthly projections of sea-surface temperature (SST) and precipitation from 30 years of hindcasts with observations of precipitation and temperature for target regions. Since raw model forecasts are well-known to have structural biases, a cross-validated multivariate regression methodology (CCA) is used to link the model projected states as predictors to the predictands of the target region. The target regions include a number of basins in East and South Africa as well as the Ganges / Baramaputra / Meghna basin complex. The MOS approach used address spatial downscaling. Temporal disaggregation of monthly seasonal forecasts is achieved through use of a tercile bootstrapping approach. We interpret the results of these studies, the levels of skill by several metrics, and key uncertainties.

Robertson, Franklin R.↗

Modular Spectral Inference Framework Applied to Young Stars and Brown Dwarfs

In practice, synthetic spectral models are imperfect, causing inaccurate estimates of stellar parameters. Using forward modeling and statistical inference, we derive accurate stellar parameters for a given observed spectrum by emulating a grid of precomputed spectra to track uncertainties. Spectral inference as applied to brown dwarfs re: Synthetic spectral models (Marley et al 1996 and 2014) via the newest grid spans a massive multi-dimensional grid applied to IGRINS spectra, improving atmospheric models for JWST. When applied to young stars(10Myr) with large starpots, they can be measured spectroscopically, especially in the near-IR with IGRINS.

From multivariate to functional data analysis: Fundamentals, recent developments, and emerging areas

Functional data analysis (FDA), which is a branch of statistics on modeling infinite dimensional random vectors resided in functional spaces, has become a major research area for Journal of Multivariate Analysis. We review some fundamental concepts of FDA, their origins and connections from multivariate analysis, and some of its recent developments, including multi-level functional data analysis, high-dimensional functional regression, and dependent functional data analysis. Here, we also discuss the impact of these new methodology developments on genetics, plant science, wearable device data analysis, image data analysis, and business analytics. Two real data examples are provided to motivate our discussions.

97 MATHEMATICS AND COMPUTING↗

Global Carbon Budget 2016

Accurate assessment of anthropogenic carbon dioxide (CO 2 ) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere – the “global carbon budget” – is important to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe data sets and methodology to quantify all major components of the global carbon budget, including their uncertainties, based on the combination of a range of data, algorithms, statistics, and model estimates and their interpretation by a broad scientific community. We discuss changes compared to previous estimates and consistency within and among components, alongside methodology and data limitations. CO 2 emissions from fossil fuels and industry ( E FF ) are based on energy statistics and cement production data, respectively, while emissions from land-use change ( E LUC ), mainly deforestation, are based on combined evidence from land-cover change data, fire activity associated with deforestation, and models. The global atmospheric CO 2 concentration is measured directly and its rate of growth ( G ATM ) is computed from the annual changes in concentration. The mean ocean CO 2 sink ( S OCEAN ) is based on observations from the 1990s, while the annual anomalies and trends are estimated with ocean models. The variability in S OCEAN is evaluated with data products based on surveys of ocean CO 2 measurements. The global residual terrestrial CO 2 sink ( S LAND ) is estimated by the difference of the other terms of the global carbon budget and compared to results of independent dynamic global vegetation models. We compare the mean land and ocean fluxes and their variability to estimates from three atmospheric inverse methods for three broad latitude bands. All uncertainties are reported as ±1 σ , reflecting the current capacity to characterise the annual estimates of each component of the global carbon budget. For the last decade available (2006–2015), E FF was 9.3 ± 0.5 GtC yr -1 , E LUC 1.0 ± 0.5 GtC yr -1 , G ATM 4.5 ± 0.1 GtC yr -1 , S OCEAN 2.6 ± 0.5 GtC yr -1 , and S LAND 3.1 ± 0.9 GtC yr -1 . For year 2015 alone, the growth in E FF was approximately zero and emissions remained at 9.9 ± 0.5 GtC yr -1 , showing a slowdown in growth of these emissions compared to the average growth of 1.8 % yr -1 that took place during 2006–2015. Also, for 2015, E LUC was 1.3 ± 0.5 GtC yr -1 , G ATM was 6.3 ± 0.2 GtC yr -1 , S OCEAN was 3.0 ± 0.5 GtC yr -1 , and S LAND was 1.9 ± 0.9 GtC yr -1 . G ATM was higher in 2015 compared to the past decade (2006–2015), reflecting a smaller S LAND for that year. The global atmospheric CO 2 concentration reached 399.4 ± 0.1 ppm averaged over 2015. For 2016, preliminary data indicate the continuation of low growth in E FF with +0.2 % (range of -1.0 to +1.8 %) based on national emissions projections for China and USA, and projections of gross domestic product corrected for recent changes in the carbon intensity of the economy for the rest of the world. In spite of the low growth of E FF in 2016, the growth rate in atmospheric CO 2 concentration is expected to be relatively high because of the persistence of the smaller residual terrestrial sink ( S LAND ) in response to El Niño conditions of 2015–2016. From this projection of E FF and assumed constant E LUC for 2016, cumulative emissions of CO 2 will reach 565 ± 55 GtC (2075 ± 205 GtCO 2 ) for 1870–2016, about 75 % from E FF and 25 % from E LUC . This living data update documents changes in the methods and data sets used in this new carbon budget compared with previous publications of this data set (Le Quéré et al., 2015b, a, 2014, 2013). All observations presented here can be downloaded from the Carbon Dioxide Information Analysis Center ( doi:10.3334/CDIAC/GCP_2016 ).

54 ENVIRONMENTAL SCIENCES↗