Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reflection data”

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 721 records · Page 40

Cyclops I: Modeling Error Study for Inertially Confined Detonation of Insensitive High Explosive (IHE) PBX-9502

A validation study is performed using the hydrodynamics code PAGOSA to assess the accuracy of reflected shock, burn front propagation, and Photon Doppler Velocimetry (PDV) predictions using several burn models and strength models. Experimental data collected from proton radiography (pRad) is compared to burn front propagation times and reflected shock patterns. Burn front isochrones are computed and the lag time for the burn front is quantified for a given burn model. The effect of strength models and damage are assessed by comparison of simulations to experimental PDV data.

97 MATHEMATICS AND COMPUTING↗

Analysis of SCALE Criticality and Sensitivity Calculations for Reflected HEU Cylinders [Slides]

The SCALE code package offers multiple nuclear data libraries supporting Monte Carlo transport, with sensitivity and uncertainty methods derived from MC transport solutions. Several libraries are multigroup, which introduce bias differing by system. Previous work has shown poor S/U results in several reflector materials: ICSBEP benchmark HMF-084 was selected to analyze biases and S/U method applicability to a variety of reflectors. Prior and ongoing work found inaccuracies in CSAS and TSUNAMI results, which were further investigated utilizing the HEU-MET-FAST-084 ICSBEP critical benchmark, chosen for its geometrical simplicity and variety of reflector materials. Perturbation of reflector thickness across various reflector materials allowed for an assortment of materials is to be tested swiftly for each sequence and method. Observation was an increasing bias of MG $k_{eff}$ relative to CE, in both direction and magnitude. IFP produced extremely reliable results. >85% of CLUTCH cases were found satisfactory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Analysis of SCALE Criticality and Sensitivity Calculations for Reflected HEU Cylinders [Abstract]

The SCALE code package offers several nuclear data libraries to support Monte Carlo (MC) transport, as well as MC-based derivation of $\kappa$ eff sensitivity and uncertainty (S/U) data. The CSAS sequence using the KENO MC code can utilize continuous-energy (CE) cross sections, or pre-generated multigroup (MG) cross section libraries. The use of MG libraries introduces bias into calculations in exchange for faster transport solutions. The TSUNAMI-3D sequence also utilizes KENO MC calculations. TSUNAMI-3D has two CE calculational methods: the Iterated Fission Probability (IFP) method, and the Contribution-Linked eigenvalue sensitivity/Uncertainty estimation via Tracklength importance CHaracterization (CLUTCH) method. Previous work has shown poor agreement between CLUTCH and confirmatory direct perturbation calculations in specific applications, e.g., fissionable and polyethylene reflectors

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NREL Solar Data, Analysis, & Tools: FY19 Research Dissemination Statistics

This report summarizes the key Fiscal Year 2019 (FY19) outreach metrics for the NREL Solar Data, Analysis, and Tools (DAT) portfolio. All metrics herein reflect the reporting period of Oct. 1, 2018, through Sept. 30, 2019. These projects are all supported by the U.S. Department of Energy Solar Energy Technologies Office Soft Costs program.

14 SOLAR ENERGY↗

A simple technique to classify diffraction data from dynamic proteins according to individual polymorphs

One often observes small but measurable differences in the diffraction data measured from different crystals of a single protein. These differences might reflect structural differences in the protein and may reveal the natural dynamism of the molecule in solution. Partitioning these mixed-state data into single-state clusters is a critical step that could extract information about the dynamic behavior of proteins from hundreds or thousands of single-crystal data sets. Mixed-state data can be obtained deliberately (through intentional perturbation) or inadvertently (while attempting to measure highly redundant single-crystal data). To the extent that different states adopt different molecular structures, one expects to observe differences in the crystals; each of the polystates will create a polymorph of the crystals. After mixed-state diffraction data have been measured, deliberately or inadvertently, the challenge is to sort the data into clusters that may represent relevant biological polystates. Here, this problem is addressed using a simple multi-factor clustering approach that classifies each data set using independent observables, thereby assigning each data set to the correct location in conformational space. This procedure is illustrated using two independent observables, unit-cell parameters and intensities, to cluster mixed-state data from chymotrypsinogen (ChTg) crystals. It is observed that the data populate an arc of the reaction trajectory as ChTg is converted into chymotrypsin.

36 MATERIALS SCIENCE↗

K-space algorithmic reconstruction (KAREN): a robust statistical methodology to separate Bragg and diffuse scattering

Diffuse scattering occurring in the Bragg diffraction pattern of a long-range-ordered structure represents local deviation from the governing regular lattice. However, interpreting the real-space structure from the diffraction pattern presents a significant challenge because of the dramatic difference in intensity between the Bragg and diffuse components of the total scattering function. In contrast to the sharp Bragg diffraction, the diffuse signal has generally been considered to be a weak expansive or continuous background signal. In this paper, using 1D and 2D models, it is demonstrated that diffuse scattering in fact consists of a complex array of high-frequency features that must not be averaged into a low-frequency background signal. To evaluate the actual diffuse scattering effectively, an algorithm has been developed that uses robust statistics and traditional signal processing techniques to identify Bragg peaks as signal outliers which can be removed from the overall scattering data and then replaced by statistically valid fill values. This method, described as a 'K-space algorithmic reconstruction' (KAREN), can identify Bragg reflections independent of prior knowledge of a system's unit cell. KAREN does not alter any data other than that in the immediate vicinity of the Bragg reflections, and reconstructs the diffuse component surrounding the Bragg peaks without introducing discontinuities which induce Fourier ripples or artifacts from underfilling `punched' voids. The KAREN algorithm for reconstructing diffuse scattering provides demonstrably better resolution than can be obtained from previously described punch-and-fill methods. The superior structural resolution obtained using the KAREN method is demonstrated by evaluating the complex ordered diffuse scattering observed from the neutron diffraction of a single plastic crystal of CBr 4 using pair distribution function analysis.

36 MATERIALS SCIENCE↗

The initial appraisal of buried DAS system in CO2CRC Otway Project: the comparison of buried standard fibre-optic and helically wound cables using 2D imaging

This study aims to assess the ability of shallow distributed acoustic sensing (DAS) to serve as a cost-effective seismic sensor array for permanent monitoring applications. To this end, as part of the CO2CRC seismic monitoring program, a fibre-optic DAS array was deployed alongside a permanently buried geophone array at the Otway Project site (Victoria, Australia). The DAS array consisted of a standard commercially available tactical fibre-optic cable, which was deployed in 0.8m deep trenches. A custom-designed helically wound (HW) cable was also deployed in one of the DAS trenches for comparison of the cable designs. Simultaneous acquisition of the seismic data was carried out using 3000 vibroseis source points and geophones, DAS standard and HW cables. For initial assessment of the seismic images acquired with DAS and to compare different cable designs, preliminary 2D seismic reflection processing is conducted on both DAS cables and geophone data along a single 2D line. The geophone data processing guided processing of the DAS data. Several shallow structures (100-450ms) and some important reflectors at 450-600ms are observed on the final DAS images. Comparison of the two different DAS cable types demonstrated that seismic imaging would benefit DAS technology. However, the benefit of utilising HW cable is modest compared with the standard cable. The workflows and results of this study pave the way for processing of the 3D seismic data set acquired with the DAS array, as well as further detailed analysis of the DAS cables and the system itself.

47 OTHER INSTRUMENTATION↗

The new silver borate Ag{sub 3}B{sub 5}O{sub 9}

Single crystals of Ag{sub 3}B{sub 5}O{sub 9} were obtained via high-pressure synthesis at 3 GPa and 600 °C, using a Walker-type multianvil high-pressure device. Ag{sub 3}B{sub 5}O{sub 9} crystalizes with a=674.7(2), b=943.5(2), c=1103.5(2) pm, V=0.7025(2) nm{sup 3}, and Z=4 in the noncentrosymmetric space group P2{sub 1}2{sub 1}2{sub 1} (no. 19). The orthorhombic structure was refined from 3740 independent reflections with R1=0.0496 and wR2=0.587 (all data). It is built up from infinite corner-sharing chains of BO{sub 4} tetrahedra along the a axis, which are interconnected by BO{sub 3} groups to form a network. In the structure, three crystallographically independent sites are occupied with Ag{sup +} cations exhibiting argentophillic interactions. The synthetic conditions as well as the results of the single crystal structure analysis are presented. - Graphical abstract: Noncentrosymmetric silver borate: During investigations in the system Ag–B–O, a new noncentrosymmetric silver borate Ag{sub 3}B{sub 5}O{sub 9} was discovered. The new structure type is built up from corner-sharing BO{sub 3} and BO{sub 4} groups, forming a network. Argentophillic interactions are clearly indicated by the Ag{sup +}⋯Ag{sup +} distances present in the structure. - Highlights: • A noncentrosymmetric borate Ag{sub 3}B{sub 5}O{sub 9} is accessible via high-pressure synthesis. • Ag{sub 3}B{sub 5}O{sub 9} is the second high-pressure silver borate. • Ag{sup +}⋯Ag{sup +} distances in Ag3B5O9 clearly indicate the presence of argentophillic interactions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Living on a Flare: Relativistic Reflection in V404 Cyg Observed by NuSTAR during Its Summer 2015 Outburst

We present first results from a series of NuSTAR observations of the black hole X-ray binary V404 Cyg obtained during its summer 2015 outburst, primarily focusing on observations during the height of this outburst activity. The NuSTAR data show extreme variability in both the flux and spectral properties of the source. This is partly driven by strong and variable line-of-sight absorption, similar to previous outbursts. The latter stages of this observation are dominated by strong flares, reaching luminosities close to Eddington. During these flares, the central source appears to be relatively unobscured and the data show clear evidence for a strong contribution from relativistic reflection, providing a means to probe the geometry of the innermost accretion flow. Based on the flare properties, analogies with other Galactic black hole binaries, and also the simultaneous onset of radio activity, we argue that this intense X-ray flaring is related to transient jet activity during which the ejected plasma is the primary source of illumination for the accretion disk. If this is the case, then our reflection modeling implies that these jets are launched in close proximity to the black hole (as close as a few gravitational radii), consistent with expectations for jet launching models that tap either the spin of the central black hole, or the very innermost accretion disk. Our analysis also allows us to place the first constraints on the black hole spin for this source, which we find to be a{sup ∗}>0.92 (99% statistical uncertainty, based on an idealized lamp-post geometry).

79 ASTRONOMY AND ASTROPHYSICS↗

The hidden treasure in your data: phasing with unexpected weak anomalous scatterers from routine data sets

Single-wavelength anomalous dispersion (SAD) utilizing anomalous signal from native S atoms, or other atoms with Z ≤ 20, generally requires highly redundant data collected using relatively long-wavelength X-rays. Here, the results from two proteins are presented where the anomalous signal from serendipitously acquired surface-bound Ca atoms with an anomalous data multiplicity of around 10 was utilized to drive de novo structure determination. In both cases, the Ca atoms were acquired from the crystallization solution, and the data-collection strategy was not optimized to exploit the anomalous signal from these scatterers. The X-ray data were collected at 0.98 Å wavelength in one case and at 1.74 Å in the other (the wavelength was optimized for sulfur, but the anomalous signal from calcium was exploited for structure solution). Similarly, using a test case, it is shown that data collected at ~1.0 Å wavelength, where the f '' value for sulfur is 0.28 e, are sufficient for structure determination using intrinsic S atoms from a strongly diffracting crystal. Interestingly, it was also observed that SHELXD was capable of generating a substructure solution from high-exposure data with a completeness of 70% for low-resolution reflections extending to 3.5 Å resolution with relatively low anomalous multiplicity. Considering the fact that many crystallization conditions contain anomalous scatterers such as Cl, Ca, Mn etc. , checking for the presence of fortuitous anomalous signal in data from well diffracting crystals could prove useful in either determining the structure de novo or in accurately assigning surface-bound atoms.

59 BASIC BIOLOGICAL SCIENCES↗

Design and performance of optimal detectors for guided wave structural health monitoring

Ultrasonic guided wave measurements in a long term structural health monitoring system are affected by measurement noise, environmental conditions, transducer aging and malfunction. This results in measurement variability which affects detection performance, especially in complex structures where baseline data comparison is required. This paper derives the optimal detector structure, within the framework of detection theory, where a guided wave signal at the sensor is represented by a single feature value that can be used for comparison with a threshold. Three different types of detectors are derived depending on the underlying structure’s complexity: (i) Simple structures where defect reflections can be identified without the need for baseline data; (ii) Simple structures that require baseline data due to overlap of defect scatter with scatter from structural features; (iii) Complex structure with dense structural features that require baseline data. The detectors are derived by modeling the effects of variabilities and uncertainties as random processes. Analytical solutions for the performance of detectors in terms of the probability of detection and false alarm are derived. A finite element model is used to generate guided wave signals and the performance results of a Monte-Carlo simulation are compared with the theoretical performance. initial results demonstrate that the problems of signal complexity and environmental variability can in fact be exploited to improve detection performance.

60 APPLIED LIFE SCIENCES↗

Structural and Stratigraphic Controls on Methane Hydrate occurrence and distribution: Gulf of Mexico, Walker Ridge 313 and Green Canyon 955: Final Report

The goal of this project was to determine structural and stratigraphic controls on hydrate occurrence and distribution in Green Canyon (GC) 955 and Walker Ridge (WR) 313 blocks using seismic and well data. Gas hydrate was discovered in these blocks in coarse- and fine-grained sediments during the 2009 Joint Industrial project (JIP) Leg 11 drilling expedition. Although the immediate interest of the exploration community is exclusively hydrate which is present in coarse–grained sediments, factors that control hydrate and free gas distribution in the two blocks and whether coarse and fine-grained hydrate-bearing units are related in any manner, formed the core of this research. The project spanned from 10/01/2012 to 07/31/2016. In the project, in both the leased blocks, the interval spanning the gas hydrate stability zone (GHSZ) was characterized using a joint analysis of sparse Ocean Bottom Seismic (OBS) and dense, surface–towed multichannel seismic (MCS) data. The project team had the luxury of calibrating their results with two well logs. Advance processing methods such as depth migration and full-waveform inversion (FWI) were used for seismic data analysis. Hydrate quantification was achieved through interpretation of the FWI velocity field using appropriate rock physics models at both blocks. The seismic modeling/inversion methodology (common to both GC955 and WR313 blocks) was as follows. First, the MCS data were depth migrated using a P-wave velocity (VP) model constructed using inversion of reflection arrival times of a few (four in both cases) key horizons carefully picked in the OBS data to farthest possible offsets. Then, the resolution of the traveltime VP model was improved to wavelength scale by inverting OBS gathers up to the highest frequency possible (21.75 Hz for GC955 and 17.5 for WR313) using FWI. Finally, the hydrate saturation (or the volume fraction) was estimated at the well location assuming one of the other hydrate morphology (filling the primary or the secondary porosity) was extrapolated out from the wells using the FWI VP as a guide. General outcomes were as follows. First and foremost, an imaging methodology using sparse seismic data, which is easily replicable at other sites with similar datasets, has been demonstrated. The end product of this methodology at both the leased blocks is quantitative estimates of hydrate distribution. Second, at both locations there is strong evidence that the base of the GHSZ, which does not appear as a clear Bottom Simulating Reflection (BSR), manifests in the VP perturbations created by FWI, suggesting that FWI is sensitive to subtle compositional changes in shallow sediments and establishes it as a valuable tool for investigations of hydrate-bearing basins. Third, through joint interpretation of the depth migrated image and the FWI VP model, how structure and stratigraphy jointly determine hydrate and free gas distribution in both blocks could be clearly visualized. The joint interpretation also suggests that the coarse and fine grained hydrate-bearing sediments at both leased are connected. Site specific results, in addition to general results, are as follows. At GC955 the overlying fine-grained hydrate-bearing unit could have been sourced from the underlying hydrate coarse-grained channel-levee complex through a chimney feature. The channel-levee system at GC955 is compartmentalized by faults, of which only a few may be impermeable. Although compartmentalized, the channel-levee system in the GC955 as a whole might be in communication except selected zones. At WR313 the overlying fine-grained fracture-filled hydrate unit appears to be sourced from below the GHSZ. The reason that only a particular fine-grained unit has hydrate, despite having lower porosity that the bounding units, could be the presence of secondary porosity (such as those formed from clay dewatering under compaction). In conclusion, the project was a pioneering effort in in joint analysis of OBS and MCS datasets for advancing the knowledge about a hydrate and free–gas system dynamics using advanced processing methods such as FWI and depth migration. Results obtained in this project can greatly advance the tools and techniques used for delineating specific hydrate prospects. Results obtained in this project can also be seamlessly incorporated into other DOE funded project on modeling the potential productivity and commercial viability of hydrate from sand-dominated reservoirs. The OBS and MCS data in this project were acquired in 2012 (after the JIP II drilling) by the USGS and therefore the results are a posteriori. Nonetheless, the seismic inversion workflow established through this project can be used to generate various what-if quantification scenarios even in absence of logs and serve as a valuable tool for guiding drilling operations. Results from this project can augment other DOE sponsored projects on determining the commercial viability of methane production from the Gulf of Mexico.

03 NATURAL GAS↗

Interpretation of HFTS 2 Microseismic Data Using Bedding-Plane-Slip Mechanism

Summary The objective of this study is to understand how microseismic events are generated during hydraulic fracturing, as well as the role of geomechanical conditions (i.e., stress and mechanical stratigraphy) in this process. In the industry, microseismic event clouds have been generally used as an “outer boundary” of the “stimulated reservoir volume” (SRV). However, by comparing with other surveillance data (low-frequency distributed acoustic sensing, or LF-DAS, strain) in the Hydraulic Fracturing Test Site (HFTS) 2 experiment, we show that this assumption is fundamentally flawed. The HFTS 2 data have three unique observations that have not been commonly observed in other datasets: (1) Due to the influence of offset pad depletion, microseismic data show that hydraulic fractures from the child well can propagate over 3,000 ft into the depleted low-stress zone. (2) By comparing microseismic and horizontal fiber LF-DAS strain data, we observe that the microseismic event cloud does not necessarily reflect the created hydraulic fracture volume. Particularly, the extent of microseismic event clouds near heel stages is much shorter than what is shown with LF-DAS strain data. (3) Microseismic event magnitudes are larger in the depleted regions. Through geomechanical analysis, we demonstrate that the “bedding-plane-slip” model is likely the mechanism for microseismic generation during hydraulic fracturing. This model successfully explains the above field observations from the HFTS 2 experiment. We also provide a quantitative relationship connecting the microseismic event magnitude with fracture width increment and layer mechanical property contrast.

Engineering↗

Measurements of Neutron Importance Distribution in {sup 235}U Fast Critical Reactor

The reactor physics related integral experiments are significant to criticality safety analysis and the design of the next generation nuclear system. Measurements of parameters such as prompt neutron decay constant α, prompt neutron generation time Λ, effective delayed neutron fraction β{sub eff}, neutron importance distribution, etc., can provide data for the verification of calculation methods, computer codes and cross-section data in criticality safety analysis and test the design of new reactors. Among these parameters, neutron importance distribution has its unique value. Beside verification and validation of nuclear data, computer codes, and simulation methods, the experimental results of neutron importance distribution are also useful to the correcting of the spatial effects in Rossi-α measurement, which is based on point-reactor neutron kinetics. The relative neutron importance distribution experiments were performed on the {sup 235}U fast critical reactor of China Academy of Engineering Physics (CAEP). The key experimental objective was to obtain the importance distribution data for nuclear data validation and criticality safety analysis. We demonstrate in this paper the experiment and the data analysis. {sup 252}Cf neutron importance spatial distributions in the diametral hole of a depleted uranium reflected {sup 235}U fast critical reactor were measured. The evaluated uncertainties of the experiment data were varying with the source location but within 1.5% ranges in all. The experiment results can provide useful data for applications like nuclear data validation and the determination of rectification factor R of Rossi-a measurements. However, the experimental uncertainties are not satisfying and can be further reduced by using a {sup 252}Cf source with much higher intensity, recording more precisely the absolute position of the neutron source, and carrying out the experiment in more shallow subcritical states in which the system multiplication is larger. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Assimilation of X-band Phased-Array Radar Data with EnKF for the Analysis and Warning Forecast of a Tornadic storm

The impact of assimilating China's operational X-band Phased-Array radar's (X-PAR) data on the analysis and warning forecast of the vortex structure and intensity of the June 8, 2018 Foshan, Guangdong province, tornadic storm was investigated for the first time using an Ensemble Kalman Filter (EnKF) data assimilation system. Both radar radial velocity (Vr) and reflectivity (Z) from two S-band operational radars and one X-PAR were assimilated. Deterministic forecasts were launched every 6 min from 05:42 UTC (20 min before the tornado touched down) to 06:00 UTC from the EnKF mean analysis field. Five experiments were conducted to examine the added capability of Z assimilation of the EnKF system, and to investigate the impact of assimilating X-PAR data on the analysis and prediction of the tornadic storm. Compared to the experiment without Z assimilation, the assimilation of Z reduced the analysis error and greatly reduced the forecast error of Z. Furthermore, the assimilation of X-PAR data greatly improved the vortex structure of the tornadic storm at low levels, and improved the intensity of the rear inflow of the tornadic storm, especially with a higher assimilation frequency. Compared to the experiments without X-PAR data assimilation, assimilating X-PAR data improved the predictability of tornadic storm.

54 ENVIRONMENTAL SCIENCES↗

The diurnal cycle of clouds and precipitation at the ARM SGP site: Cloud radar observations and simulations from the multiscale modeling framework

Millimeter Wavelength Cloud Radar (MMCR) data from December 1996 to December 2010, collected at the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) program Southern Great Plains (SGP) site, are used to examine the diurnal cycle of hydrometeor occurrence. These data are categorized into clouds (-40 dBZ e ≤ reflectivity < -10 dBZ e ), drizzle and light precipitation (-10 dBZ e ≤ reflectivity < 10 dBZ e ), and heavy precipitation (reflectivity ≥ 10 dBZ e ). The same criteria are implemented for the observation-equivalent reflectivity calculated by feeding outputs from a Multiscale Modeling Framework (MMF) climate model into a radar simulator. The MMF model consists of the National Center for Atmospheric Research Community Atmosphere Model with conventional cloud parameterizations replaced by a cloud-resolving model. We find that a radar simulator combined with the simple reflectivity categories can be an effective approach for evaluating diurnal variations in model hydrometeor occurrence. It is shown that the MMF only marginally captures observed increases in the occurrence of boundary layer clouds after sunrise in spring and autumn and does not capture diurnal changes in boundary layer clouds during the summer. Above the boundary layer, the MMF captures reasonably well diurnal variations in the vertical structure of clouds and light and heavy precipitation in the summer but not in the spring.

54 ENVIRONMENTAL SCIENCES↗

Differential equations for cosmological correlators

Cosmological fluctuations retain a memory of the physics that generated them in their spatial correlations. The strength of correlations varies smoothly as a function of external kinematics, which is encoded in differential equations satisfied by cosmological correlation functions. In this work, we provide a broader perspective on the origin and structure of these differential equations. As a concrete example, we study conformally coupled scalar fields in a power-law cosmology. The wavefunction coefficients in this model have integral representations, with the integrands being the product of the corresponding flat-space results and “twist factors” that depend on the cosmological evolution. Similar twisted integrals arise for loop amplitudes in dimensional regularization, and their recent study has led to the discovery of rich mathematical structures and powerful new tools for computing multi-loop Feynman integrals in quantum field theory. The integrals of interest in cosmology are also part of a finite-dimensional basis of master integrals, which satisfy a system of first-order differential equations. We develop a formalism to derive these differential equations for arbitrary tree graphs. The results can be represented in graphical form by associating the singularities of the differential equations with a set of graph tubings. Upon differentiation, these tubings grow in a local and predictive fashion. In fact, a few remarkably simple rules allow us to predict — by hand — the equations for all tree graphs. While the rules of this “kinematic flow” are defined purely in terms of data on the boundary of the spacetime, they reflect the physics of bulk time evolution. We also study the analogous structures in tr ϕ 3 theory, and see some glimpses of hidden structure in the sum over planar graphs. This suggests that there is an autonomous combinatorial or geometric construction from which cosmological correlations, and the associated spacetime, emerge.

Cosmological models↗