Search NASA⌕ Search

SEARCH · Search NASA

Results for “Statistical homogeneity”

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 397 records · Page 22

The Impact of FU Orionis Outbursts and the Solar Nebula

Protostellar systems are variable on many timescales. One of the most dramatic forms of variability known to occur in low mass stellar systems is the FU Orionis outburst (Herbig 1977). Throughout a typical outburst lasting several decades, system luminosities may be a hundred times what is typical of the quiesent state. FU Orionis outburst events are thought to have significant impact on the thermal structure of the protosolar nebula; their existence has been used to explain features in the meteoritic record from thermally induced homogenization to the formation of chondrules. Until recently, the magnitude of the likely effect from such outbursts has been largely speculative due to the lack of a detailed understanding of the outburst mechanism. Recent numerical models (Bell\& Lin 1994) have demonstrated the viability of the observational hypothesis (Hartmann\& Kenyon 1985) that the radiation observed during outburst is emitted by a luminous circumstellar disk transporting mass at a thousand times the quiesent rate. Light curves and color and line width evolution observed in FU Orionis systems are naturally explained by time dependent outbursting model disks (Bell et al. 1995). The radial temperature structure and shape of the disk during outburst derived from these models may be used to calculate the outburst's expected impact on primitive material at various radii throughout the disk. In this review, we will begin by discussing what is known about the FU Orionis outburst phenomenon from recent observations and theory including statistically deduced outburst timescales and observed peak temperatures. Unless covered by another author, we will discuss the evidence which suggests that outburst radiation is emitted by a circumstellar disk rather than by the star and will briefly review the thermal instability as a mechanism for outburst. We will then report on recent work which investigates the likely heating of solar nebula material due to FU Orionis outbursts including the following effects: (1) heating of the planet forming region by direct radiation from the hot inner nebula; (2) heating by the diffuse radiation field of a coccooning envelope; and (3) time-dependent penetration of the increased luminosity from the above sources into the optically thick nebula. Some of this work is currently in progress. The potential effects on condensation and migration in the nebula and the thermal processing of solids will be evaluated.

Bell, Robbins↗

Physics-informed machine learning of the Lagrangian dynamics of velocity gradient tensor

Reduced models describing the Lagrangian dynamics of the velocity gradient tensor (VGT) in homogeneous isotropic turbulence (HIT) are developed under the physics-informed machine learning (PIML) framework. We consider the VGT at both Kolmogorov scale and coarse-grained scale within the inertial range of HIT. Building reduced models requires resolving the pressure Hessian and subfilter contributions, which is accomplished by constructing them using the integrity bases and invariants of the VGT. The developed models can be expressed using the extended tensor basis neural network (TBNN) introduced by Ling et al. [J. Fluid Mech. 807, 155 (2016)]. Physical constraints, such as Galilean invariance, rotational invariance, and incompressibility condition, are thus embedded in the models explicitly. Our PIML models are trained on the Lagrangian data from a high-Reynolds number direct numerical simulation (DNS). To validate the results, we perform a comprehensive out-of-sample test. We observe that the PIML model provides an improved representation for the magnitude and orientation of the small-scale pressure Hessian contributions. Statistics of the flow, as indicated by the joint PDF of second and third invariants of the VGT, show good agreement with the “ground-truth” DNS data. A number of other important features describing the structure of HIT are reproduced by the model successfully. We have also identified challenges in modeling inertial range dynamics, which indicates that a richer modeling strategy is required. This helps us identify important directions for future research, in particular towards including inertial range geometry into the TBNN.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A co-kurtosis PCA based dimensionality reduction with nonlinear reconstruction using neural networks

For turbulent reacting flow systems, identification of low-dimensional representations of the thermo-chemical state space is vitally important, primarily to significantly reduce the computational cost of device-scale simulations. Principal component analysis (PCA), and its variants, are a widely employed class of methods. Recently, an alternative technique that focuses on higher-order statistical interactions, co-kurtosis PCA (CoK-PCA), has been shown to effectively provide a low-dimensional representation by capturing the stiff chemical dynamics associated with spatiotemporally localized reaction zones. While its effectiveness has only been demonstrated based on a priori analyses with linear reconstruction, in this work, we employ nonlinear techniques to reconstruct the full thermo-chemical state and evaluate the efficacy of CoK-PCA compared to PCA. Specifically, we combine a CoK-PCA-/PCA-based dimensionality reduction (encoding) with an artificial neural network (ANN) based reconstruction (decoding) and examine, a priori, the reconstruction errors of the thermo-chemical state. In addition, we evaluate the errors in species production rates and heat release rates, which are nonlinear functions of the reconstructed state, as a measure of the overall accuracy of the dimensionality reduction technique. We employ four datasets to assess CoK-PCA/PCA coupled with ANN-based reconstruction: zero-dimensional (homogeneous) reactor for autoignition of an ethylene/air mixture that has conventional single-stage ignition kinetics, a dimethyl ether (DME)/air mixture which has two-stage (low and high temperature) ignition kinetics, a one-dimensional freely propagating premixed ethylene/air laminar flame, and a two-dimensional dataset representing turbulent autoignition of ethanol in a homogeneous charge compression ignition (HCCI) engine. Finally, results from the analyses demonstrate the robustness of the CoK-PCA based low-dimensional manifold with ANN reconstruction in accurately capturing the data, specifically from the reaction zones.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An engineering model of woven composites based on micromechanics

Composites with three-dimensional woven architectures exhibit large strains to failure when compared to composites made up of the same materials but not with three-dimensional interlocking tows. The fracture mechanics of such three-dimensional architectures is a subject requiring substantial investigation and experimental testing. Classical fracture mechanics concepts (for instance, an isolated defect in a homogeneous body) will not be applicable to the woven fracture test specimen. The use of an isolated singularity to characterize an entire specimen is inadequate when the density of defects is considerable and the material is heterogeneous. Modelling of such a complex system requires a great deal of insight and consideration as well as prudent choices of model sizes to make numerical schemes feasible. The purpose of this manuscript is to review our recently acquired knowledge of damage accumulation in woven composites and to describe a practicable model of the macroscopic behavior in these and other complex composite architectures based on such knowledge. In this manuscript, discussion will be limited to uniaxial compressive loading; considerations of general loading (monotonic and cyclic) will appear in a subsequent manuscript. Our modelling efforts may be briefly described as follows: the composite is subdivided into microstructural elements (microelements) in which the micromechanical modelling is either understood rigorously or can be represented adequately by statistical parameters. There can be microstructural elements for many different types of composite components, such as the various types of warp and weft and matrix for three-dimensional woven composites. The physical dimensions of microelements are made as large as possible while the response within the element can still be represented by a single micromechatlical calculation. The various elements are linked together(sometimes by associating distinct corners and edges, sometimes by superposition) in a pattern which resembles a particular weave architecture. The model can then be loaded in any manner and the linear and nonlinear elastic responses of representative weaves can be calculated. After the elastic regime, the fracture response is determined by monitoring the damage accumulation.

Carter, W. C.↗

Scattering Properties and Brightness Temperatures Associated with Solid Precipitation

In the past few years, early solid precipitation detection and retrieval algorithms have been developed and shown to be applicable for snowing clouds and blizzards. NOAA has an operational snow versus rain classifier based on AMSU-B observations. Solid precipitation retrieval algorithms reported in the literature over the past two years include those that rely on neural nets, statistics, or physical relationships. All of the algorithms require the use of millimeter-wave radiometer observations. The millimeter-wave frequencies are especially sensitive to the scattering and emission properties of frozen particles due to the ice particle refractive index. Passive radiometric channels respond to both the integrated particle mass throughout the volume and field of view, and to the amount, location, and size distribution of the frozen (and liquid) particles with the sensitivity varying for different frequencies and hydrometeor types. This investigation probes the sensitivity of scattering and absorption coefficients, and hence computed brightness temperatures, resulting from variations in solid precipitation cloud profiles. The first study compares the single scattering, absorption, and asymmetry parameters associated with snow particles in clouds. Several methodologies are used to convert the physical characteristics (e.g., shape, size distributions, ice-air-water ratios) of ice particles to electromagnetic properties (e.g., absorption, scattering, and asymmetry factors). These methodologies include: conversion to solid ice particles, homogeneous dielectric mixing, or discrete dipole approximation. Changes in the conversion methodology can produce computed brightness temperature differences greater than 50 Kelvin.

Skofronick-Jackson, Gail M.↗

Variable-density buoyancy-driven turbulence with asymmetric initial density distribution

The effects of different initial density distributions on the evolution of buoyancy-driven homogeneous variable-density turbulence (HVDT) at low (0.05) and high (0.75) Atwood numbers are studied by using high-resolution direct numerical simulations. HVDT aims to mimic the acceleration-driven Rayleigh–Taylor and shock-driven Richtmyer–Meshkov instabilities and reveals new physics that arise from variable-density effects on the turbulent mixing. Here, the initial amounts of pure light and pure heavy flows are altered primarily to mimic the variable-density turbulence at the different locations of the Rayleigh–Taylor and Richtmyer–Meshkov instabilities’ mixing layers where the amounts of the mixing fluids are not equal. It is found that for the low Atwood number cases, the asymmetric initial density distribution has limited effects on both global and local flow evolution for HVDT. However, at high Atwood number, both global flow evolution and the local flow structures are strongly affected by the initial composition ratio. The flow composed of more light fluid reaches higher turbulent levels and the local statistics reach their fully-developed behavior earlier in the time evolution. In contrast, during the late time decay, where most of the flow is well-mixed, all parameters become independent of the initial composition ratio for both low and high Atwood number cases.

42 ENGINEERING↗

The Flux of Large Meteoroids Observed with Lunar Impact Monitoring

The flux of large meteoroids is not well determined due to relatively low number statistics, due mainly to the lack of collecting area available to meteor camera systems (10(2)-10(5) km2). Larger collecting areas are needed to provide reasonable statistics for flux calculations. The Moon, with millions of square kilometers of lunar surface, can be used as a detector for observing the population of large meteoroids in the tens of grams to kilogram mass range. This is accomplished by observing the flash of light produced when a meteoroid impacts the lunar surface, converting a portion of its kinetic energy to visible light detectable from Earth. A routine monitoring program at NASA's Marshall Space Flight Center has recorded over 300 impact flashes since early 2006. The program utilizes multiple 0.35 m (14 inch) Schmidt-Cassegrain telescopes, outfitted with video cameras using the 1/2 inch Sony EXview HAD CCDTM chip, to perform simultaneous observations of the earthshine hemisphere of the Moon when the lunar phase is between 0.1 and 0.5. This optical arrangement permits monitoring of approximately 3.8x10(6) km2 of lunar surface. A selection of 126 flashes recorded in 266.88 hours of photometric skies was analyzed, creating the largest and most homogeneous dataset of lunar impact flashes to date. Standard CCD photometric techniques outlined in [1] were applied to the video to determine the luminous energy, kinetic energy, and mass for each impactor, considering a range of luminous efficiencies. The flux to a limiting energy of 2.5x10(-6) kT TNT or 1.05×10(7) J is 1.03×10(-7) km(-2) hr(-1) and the flux to a limiting mass of 30 g is 6.14×10(-10) m(-2) yr(-1). Comparisons made with measurements and models of the meteoroid population indicate that the flux of objects in this size range is slightly lower (but within the error bars) than the power law distribution determined for the near Earth object population by [2].

Cooke, W. J.↗

Validating Mixtures of 233 U, 235 U, and 239 Pu for the Sum-of-Fractions Method

The Sum-of-Fractions method is a technique used to assure that homogeneous mixtures of fissile and fissionable isotopes are below a minimum margin of k eff or reactivity. Current work by Pacific Northwest National Laboratory examines different mixtures of 233 U, 235 U, and 239 Pu to determine critical mass limits for mixtures of transuranic actinides lacking a validation basis. To provide a validation basis for these limits, the work presented here describes the results of a sensitivity and uncertainty analysis of various mixtures of these isotopes in various concentrations moderated and reflected by light water and polyethylene. The TSUNAMI-1D sequence in the SCALE code system was used to generate sensitivity coefficients for three different concentrations of mixtures of 233 U, 235 U, and 239 Pu. The TSUNAMI-IP sequence was then used for similarity assessment (c k ) with critical benchmark experiment sensitivity data files (SDFs) from the Oak Ridge National Laboratory Verified, Archived Library of Inputs and Data and the Nuclear Energy Agency SDF database. The VADER sequence in SCALE was used for statistical testing and to generate upper subcritical limits from the data to develop a basis for validating critical mass limits.

07 ISOTOPE AND RADIATION SOURCES↗

Survey of Surveys: I. The largest compilation of radial velocities for the Galaxy

In the present-day panorama of large spectroscopic surveys, the amount, diversity, and complexity of the available data continuously increase. The overarching goal of studying the formation and evolution of our Galaxy is hampered by the heterogeneity of instruments, selection functions, analysis methods, and measured quantities. Here, we present a comprehensive catalogue, the Survey of Surveys (SoS), built by homogeneously merging the radial velocity (RV) determinations of the largest ground-based spectroscopic surveys to date, such as APOGEE, GALAH, Gaia-ESO, RAVE, and LAMOST, using Gaia as a reference. This pilot study serves to prove the concept and to test the methodology that we plan to apply in the future to the stellar parameters and abundance ratios as well. We have devised a multi-staged procedure that includes: (i) the cross match between Gaia and the spectroscopic surveys using the official Gaia cross-match algorithm, (ii) the normalisation of uncertainties using repeated measurements or the three-cornered hat method, (iii) the cross calibration of the RVs as a function of the main parameters on which depend (magnitude, effective temperature, surface gravity, metallicity, and signal-to-noise ratio) to remove trends and zero point offsets, and (iv) the comparison with external high-resolution samples, such as the Gaia RV standards and the Geneva-Copenhagen survey, to validate the homogenisation procedure and to calibrate the RV zero-point of the SoS catalogue. We provide the largest homogenised RV catalogue to date, containing almost 11 million stars, of which about half come exclusively from Gaia and half in combination with the ground-based surveys. We estimate the accuracy of the RV zero-point to be about 0.16–0.31 km s –1 and the RV precision to be in the range 0.05–1.50 km s –1 depending on the type of star and on its survey provenance. We validate the SoS RVs with open clusters from a high resolution homogeneous samples and provide the systemic velocity of 55 individual open clusters. Additionally, we provide median RVs for 532 clusters recently discovered by Gaia data. The SoS is publicly available and ready to be applied to various research projects, such as the study of star clusters, Galactic archaeology, stellar streams, or the characterisation of planet-hosting stars, to name a few. We also plan to include survey updates and more data sources in future versions of the SoS.

79 ASTRONOMY AND ASTROPHYSICS↗

Multiwavelength Observations of the AGN 1ES 0414+009 with Veritas, Fermi-LAT, Swift-XRT, and MDM

We present observations of the BL Lac object 1ES 0414+009 in the greater than 200 GeV gamma-ray band by the VERITAS array of Cherenkov telescopes. 1ES 0414+009 was observed by VERITAS between January 2008 and February 2011, resulting in 56.2 hours of good quality pointed observations. These observations resulted in a detection of 822 events from the source corresponding to a statistical significance of 6.4 standard deviations (6.4 sigma) above the background. The source flux, showing no evidence for variability, is measured as (5.2 plus or minus 1.1(sub stat) plus or minus 2.6(sub sys) x 10(exp 12) photons per square centimeter per second above 200 GeV, equivalent to approximately 2% of the Crab Nebula flux above this energy. The differential photon spectrum from 230 GeV to 850 GeV is well fit by a power law with an photon index of TAU = 3.4 plus or minus 0.5(sub stat) plus or minus 0.3(sub sys) and a flux normalization of (1.6 plus or minus 0.3(sub stat) plus or minus 0.8(sub sys) x 10(exp -11) photons per square centimeter per second at 300 GeV. We also present multiwavelength results taken in the optical (MDM), X-ray (Swift-XRT), and GeV (Fermi-LAT) bands and use these results to construct a broadband spectral energy distribution (SED). Modeling of this SED indicates that homogenous one-zone leptonic scenarios are not adequate to describe emission from the system, with a lepto-hadronic model providing a better fit to the data.

Aliu, E.↗

Development of a Time-Height Profile Assimilation Technique for Large-Eddy Simulation

Mesoscale-to-microscale coupling (MMC) aims to address the limited scope of traditional large-eddy simulations by driving the microscale flow with information concerning large-scale weather patterns provided by mesoscale models. We present a new offline MMC technique for horizontally homogeneous microscale flow conditions, in which internal forcing terms are computed based on mesoscale time–height profiles of mean-flow quantities. The advantage of such an approach is that it can be used to drive a microscale simulation with either mesoscale or observational data, and that it does not rely on specific terms in the mesoscale budget equations, which are typically not part of the default output of a mesoscale solver. The performance of the proposed profile assimilation technique is assessed based on the simulation of a typical diurnal cycle over the Scaled Wind Farm Technology site in west Texas. Results indicate that simple data assimilation techniques lead to unphysically high levels of shear and turbulence caused by the algorithm’s inability to cope with inaccuracies in the mesoscale time–height profiles. Modifying the algorithm to account for vertical coherence in the mesoscale source terms gives the microscale solver a greater ability to correct the provided mesoscale time–height profiles, leading to improved predictions of shear and turbulence statistics. The resulting turbulence statistics are in good agreement with meteorological tower observations and simulation results obtained with state-of-the-art coupling techniques using mesoscale budget components.

data assimilation↗

Experimental philosophy leading to a small scale digital data base of the conterminous United States for designing experiments with remotely sensed data

Research using satellite remotely sensed data, even within any single scientific discipline, often lacked a unifying principle or strategy with which to plan or integrate studies conducted over an area so large that exhaustive examination is infeasible, e.g., the U.S.A. However, such a series of studies would seem to be at the heart of what makes satellite remote sensing unique, that is the ability to select for study from among remotely sensed data sets distributed widely over the U.S., over time, where the resources do not exist to examine all of them. Using this philosophical underpinning and the concept of a unifying principle, an operational procedure for developing a sampling strategy and formal testable hypotheses was constructed. The procedure is applicable across disciplines, when the investigator restates the research question in symbolic form, i.e., quantifies it. The procedure is set within the statistical framework of general linear models. The dependent variable is any arbitrary function of remotely sensed data and the independent variables are values or levels of factors which represent regional climatic conditions and/or properties of the Earth's surface. These factors are operationally defined as maps from the U.S. National Atlas (U.S.G.S., 1970). Eighty-five maps from the National Atlas, representing climatic and surface attributes, were automated by point counting at an effective resolution of one observation every 17.6 km (11 miles) yielding 22,505 observations per map. The maps were registered to one another in a two step procedure producing a coarse, then fine scale registration. After registration, the maps were iteratively checked for errors using manual and automated procedures. The error free maps were annotated with identification and legend information and then stored as card images, one map to a file. A sampling design will be accomplished through a regionalization analysis of the National Atlas data base (presently being conducted). From this analysis a map of homogeneous regions of the U.S.A. will be created and samples (LANDSAT scenes) assigned by region.

Labovitz, M. L.↗

Dumb-bell galaxies in southern clusters: Catalog and preliminary statistical results

The dominant galaxy of a rich cluster is often an object whose formation and evolution is closely connected to the dynamics of the cluster itself. Hoessel (1980) and Schneider et al. (1983) estimate that 50 percent of the dominant galaxies are either of the dumb-bell type or have companions at projected distances less than 20 kpc, which is far in excess of the number expected from chance projection (see also Rood and Leir 1979). Presently there is no complete sample of these objects, with the exception of the listing of dumb-bell galaxies in BM type I and I-II clusters in the Abell statistical sample of Rood and Leir (1979). Recent dynamical studies of dumb-bell galaxies in clusters (Valentijn and Casertano, 1988) still suffer from inhomogeneity of the sample. The fact that it is a mixture of optically and radio selected objects may have introduced an unknown biases, for instance if the probability of radio emission is enhanced by the presence of close companions (Stocke, 1978, Heckman et al. 1985, Vettolani and Gregorini 1988) a bias could be present in their velocity distribution. However, this situation is bound to improve: a new sample of Abell clusters in the Southern Hemisphere has been constructed (Abell et al., 1988 hereafter ACO), which has several advantages over the original northern catalog. The plate material (IIIaJ plates) is of better quality and reaches fainter magnitudes. This makes it possible to classify the cluster types with a higher degree of accuracy, as well as to fainter magnitudes. The authors therefore decided to reconsider the whole problem constructing a new sample of dumb-bell galaxies homogeneously selected from the ACO survey. Details of the classification criteria are given.

Vettolani, G.↗

Calculation of self-diffusion coefficients in supercritical carbon dioxide using mean force kinetic theory

This paper presents an application of mean force kinetic theory (MFT) to the calculation of the self-diffusivity of CO 2 in the supercritical fluid regime. Two modifications to the typical application of MFT are employed to allow its application to a system of molecular species. Therefore, the first is the assumption that the inter-particle potential of mean force can be obtained from the molecule center-of-mass pair correlation function, which in the case of CO 2 is the C–C pair correlation function. The second is a new definition of the Enskog factor that describes the effect of correlations at the surface of the collision volume. The new definition retains the physical picture that this quantity represents a local density increase, resulting from particle correlations, relative to that in the zero density homogeneous fluid limit. These calculations are facilitated by the calculation of pair correlation functions from molecular dynamics (MD) simulations using the FEPM2 molecular CO2 model. The self-diffusivity calculated from theory is in good agreement with that from MD simulations up to and slightly beyond the density at the location of the Frenkel line. The calculation is compared with and is found to perform similarly well to other commonly used models but has a greater potential for application to systems of mixed species and to systems of particles with long range interatomic potentials due to electrostatic interactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DNS and LES of a Shear-Free Mixing Layer

The purpose of this work is twofold. First, given the computational resources available today, it is possible to reach, using DNS, higher Reynolds numbers than in Briggs et al.. In the present study, the microscale Reynolds numbers reached in the low- and high-energy homogeneous regions are, respectively, 32 and 69. The results reported earlier can thus be complemented and their robustness in the presence of increased turbulence studied. The second aim of this work is to perform a detailed and documented LES of the shear-free mixing layer. In that respect, the creation of a DNS database at higher Reynolds number is necessary in order to make meaningful LES assessments. From the point of view of LES, the shear-free mixing-layer is interesting since it allows one to test how traditional LES models perform in the presence of an inhomogeneity without having to deal with difficult numerical issues. Indeed, as argued in Briggs et al., it is possible to use a spectral code to study the shear-free mixing layer and one can thus focus on the accuracy of the modelling while avoiding contamination of the results by commutation errors etc. This paper is organized as follows. First we detail the initialization procedure used in the simulation. Since the flow is not statistically stationary, this initialization procedure has a fairly strong influence on the evolution. Although we will focus here on the shear-free mixing layer, the method proposed in the present work can easily be used for other flows with one inhomogeneous direction. The next section of the article is devoted to the description of the DNS. All the relevant parameters are listed and comparison with the Veeravalli & Warhaft experiment is performed. The section on the LES of the shear-free mixing layer follows. A detailed comparison between the filtered DNS data and the LES predictions is presented. It is shown that simple eddy viscosity models perform very well for the present test case, most probably because the flow seems to be almost isotropic in the small-scale range that is not resolved by the LES.

Knaepen, B.↗

Sustainable Development in Algeria’s Urban Areas: Population Growth and Land Consumption

We analyzed the urban development sustainability in five major urban areas of Algeria by the standard of the UN Sustainability Development Goal indicator SDG 11.3.1, which focuses on the ratio of land consumption rate to population growth rate. We utilized the annual global artificial impervious area (GAIA) dataset to characterize land-use and population data from the two censuses carried out by the National Office of Statistics (ONS) for 2008 and 2018. We discuss the prevailing relationship between urban land consumption rate and population growth rate at the smallest territorial and population census unit scale. We confirm that the indicator SDG 11.3.1 is nonlinear and that while, for example, the wilaya of Tlemcen as a whole appears to be on a sustainable path, twenty-one of its communes are not. We found that overall, and for most of its communes, the wilaya of Oran seems to have an urban land use commensurable to its population growth, but in the wilaya of Algiers, out of fifty-seven communes, only fourteen have a tendency towards sustain-able development. However, the latter wilaya hosts the country’s capital and includes government buildings that are uninhabited but are accounted for as land consumed, and as such, the relationship between urban land consumption and population growth is biased. The wilaya of Annaba showed large discrepancies in terms of land use and population growth rates, and the evolution of these quantities is not homogenous across communes and not sustainable. In the Saharan wilaya of Ghardaia, the development is not homogeneous in all communes, with smaller communes undergoing buildup increases of more than 150% over the decade. Finally, in all communes where population growth exceeded urban land growth, there will be overcrowding, an aspect neither the SDG 11.3.1 nor the impervious surface per capita indicator captures. This result, in addition to other limitations, makes SDG 11.3.1 incomplete for the determination of the sustainable development in urban areas.

Algeria↗

Parameters Impacting Columnated Granular Soil Pneumatic Seal Performance

The ability of a column of loose granular soil to form a pneumatic seal was investigated by varying the diameter of the soil column, the effective column height, and the level of compaction in the soil. The soil column diameter was tested at three levels using pipes with inner diameters measuring 5.08, 10.16, and 15.24 centimeters (2, 4, and 6 inches). Soil was filled in each pipe to form 15, 30, and 45 centimeter (6, 12, and 18 inch) tall soil columns. Each diameter/height configuration was also tested at three levels of soil compaction, compared by calculating the bulk density of the soil with mass and volume measurements. GRC-1a simulant was used, with approximate low/medium/high bulk densities of 1.6, 1.75, and 1.9 g/cc achieved with a combination of vibration and tamping. The order of tests for a given column diameter was randomized and repeated three times. With the top of the soil open to atmosphere at room temperature, compressed air was injected through a small diffuser at the column base with several small downward-facing holes. The number and size of these holes was scaled such that a constant total inlet orifice area to column cross sectional area ratio was maintained for each column diameter. Inlet air pressure was slowly increased via a precision regulator to preserve quasi-static equilibrium in the soil column to minimize the impact of dynamics. Air pressure was increased all the way through the static and bubbling regimes until slugging or turbulent behavior was observed in the soil to ensure that the entire static regime had been captured during data collection. A three-factor, three-level analysis of variance (ANOVA) statistical analysis was performed on the resulting data to determine the extent to which each physical parameter impacted the soil column seal performance. It was concluded that there is statistically significant evidence that column height, and the interaction between column height and diameter impact soil seal performance. In all other cases there was insufficient data to identify a statistically significant causal relationship. Additionally, plots were generated comparing experimental data to the predictive formula developed by Ogino et al. for fluidized beds in 1993. Because this model was developed for industrial spouted fluidized beds, the accuracy of its output prior to fluidization in the static seal ‘edge case’ is unknown, especially considering in this application the working gas was diffused across the column base rather than being injected through a spout. Further, the fidelity of the model had not yet been tested with lunar soil simulants. Plotting the Ogino et al. model alongside test data allows for a more intuitive sense of the impact of test parameters on soil seal performance, as well as providing a quick means to further tune this predictive model for more accurate use with static seals across granular lunar soil simulants. The model provided by Ogino et al. was further tuned using test data to determine the degree to which a spouted bed model could be applied to a slightly modified set of testing conditions: lunar soil simulants and a more diffuse, homogeneous application of pneumatic pressure. A Matlab script was created to test different values for the leading coefficient and exponents in Ogino’s formula. The script swept preset ranges, then iterated with higher resolutions across narrower ranges to converge on the optimal value for each parameter. The resulting adjusted model was compared to the original, as well as test data with noticeable improvements across the entire test domain.

Jack Stewart↗

Suboptimal schemes for atmospheric data assimilation based on the Kalman filter

This work is directed toward approximating the evolution of forecast error covariances for data assimilation. The performance of different algorithms based on simplification of the standard Kalman filter (KF) is studied. These are suboptimal schemes (SOSs) when compared to the KF, which is optimal for linear problems with known statistics. The SOSs considered here are several versions of optimal interpolation (OI), a scheme for height error variance advection, and a simplified KF in which the full height error covariance is advected. To employ a methodology for exact comparison among these schemes, a linear environment is maintained, in which a beta-plane shallow-water model linearized about a constant zonal flow is chosen for the test-bed dynamics. The results show that constructing dynamically balanced forecast error covariances rather than using conventional geostrophically balanced ones is essential for successful performance of any SOS. A posteriori initialization of SOSs to compensate for model - data imbalance sometimes results in poor performance. Instead, properly constructed dynamically balanced forecast error covariances eliminate the need for initialization. When the SOSs studied here make use of dynamically balanced forecast error covariances, the difference among their performances progresses naturally from conventional OI to the KF. In fact, the results suggest that even modest enhancements of OI, such as including an approximate dynamical equation for height error variances while leaving height error correlation structure homogeneous, go a long way toward achieving the performance of the KF, provided that dynamically balanced cross-covariances are constructed and that model errors are accounted for properly. The results indicate that such enhancements are necessary if unconventional data are to have a positive impact.

Todling, Ricardo↗