Moving from 2D to 3D Flow Modeling for the Waste Isolation Pilot Plant Performance Assessment.
Abstract not provided.
SEARCH · Search NASA
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.
Abstract not provided.
The physical phenomena associated with the modes of heat transfer within a packed pebble bed can be modelled by explicitly solving the associated partial differential equations using a finite volume method approach. Convection from the spheres to coolant, conduction inside and between spheres as well as radiation between various components of a packed bed reactor are the main modes of heat transfer, with several lesser, albeit contributing, phenomena also present. Advances both in software and computing power allow for a larger portion of a packed bed to be modeled with greater accuracy by integrating all applicable physics models in one simulation. A representative unstructured packed bed of uniform spheres was generated with a discrete element method (DEM) and the flow domain modelled using the CFD software package Star-CCM+, taking into account all primary modes of heat transfer as well as a user-specified approach to take secondary effects into account. Effects in the bulk, near-wall and wall-contact regions were investigated. By using such a fully integrated explicit model, a deeper understanding can be can be gained into the role of the underlying phenomena and their respective contributions to the overall heat transfer process. Fuel spheres were modeled explicitly with the conduction between spheres, bed and reflectors being prescribed by contact heat transfer theory as a function of temperature and compressive forces within the bed. The effective thermal conductivity (k{sub eff}), flow pressure drop, as well as temperature distributions in the spheres and coolant were calculated and compared with experimental results from published literature. Good agreement was found with experimental results and the significant contribution of radiation as well as inter-pebble conduction to the effective thermal conductivity could be observed from the results. (authors)
A computationally efficient simulation platform was developed that can provide representative synthetic ground motions from crustal earthquakes in the Stable Continental Regions of Central and Eastern US (CEUS), using 3D modeling and high-performance computing. The main objective was to use synthetic ground motion to provide constrains to refinements of exiting ergodic Ground Motion Models (GMMs), for large magnitude earthquakes and near-fault distances, for which these models are less reliable. Physics-based broadband (0-5Hz) ground motion simulations were used to estimate the near-fault ground motion amplitudes and within event and between-event variabilities associated with fault rupture characteristics. As part of a strategy for selecting a reginal velocity model and validation of developed rupture modeling technique, ground motions from the moment magnitude Mw5.0 November 7, 2016, Cushing Oklahoma, and Mw5.8 September 3, 2016, Pawnee Oklahoma earthquakes were simulated. In our simulations we used a 3D regional velocity model that was based on Saikia’s 1D velocity model. Saikia’s model demonstrated better performance in modelling high frequency regional wave propagation for CEUS region. The proposed 3D model includes lateral variations added to the 1D background model using the stochastic scheme of Pitarka and Mellors. Comparisons of the simulations with recordings of both earthquakes demonstrated the reliability of our deterministic simulation approach while emphasizing the importance of including small-scale variability in the regional velocity model needed to reproduce the observed high-frequency wave scattering effects. As part of validation analysis, comparisons with different GMMs for a Mw6.5 earthquake in the CESUS region resulted in a very good match between the simulated and empirical ground motion models. Initial investigations of within-event and between-event ground motion variabilities for Mw6.5 scenario earthquakes on a strike-slip fault, suggest that they are strongly related to spatial slip and slip rate variations, average rupture velocity, rupture area and rupture initiation location. For certain scenarios we found that the ground motion variability observed at near-fault distances (< 5 km) also persists at longer distances. Regardless of the rupture scenario, the simulated ground motion tends to fully saturate at short distances and for all periods. The near-fault saturation has to do with the attenuation of waves propagating along the fault and local rupture radiation pattern that also contribute to stronger ground motion variation at such distances. Analysis of effects of rupture initiation location suggest that the peak ground motion (PGV) and spectral acceleration (SA) can be quite variable due to rupture directivity effects. Such effects are stronger at periods longer than 1s. The effect of the 1D velocity models and surface topography on simulated ground motion were investigated by comparing three component synthetic seismograms computed at selected sites. Effect of surface topography was considered using the ratio between spectral accelerations simulated for two 1D models with flat surface topography and realistic model with surface topography. Overall, the topography slightly amplifies (by ~30%) the ground motion amplitude in the frequency range 1-3Hz. The effect of topography is more visible in the surface and coda waves portion of the seismograms.
A spectral 3D model of the stratosphere has been used to study the sensitivity of polar ozone with respect to a carbon dioxide increase. The lower stratospheric cooling associated with an imposed CO2 doubling may increase the probability of polar stratospheric cloud (PSC) formation and this affect ozone. The ozone perturbation obtained with the inclusion of a simple parameterization for heterogeneous chemistry on PSCs is compared to that relative to a pure homogeneous chemistry. In both cases the temperature perturbation is determined by a CO2 doubling, while the total chlorine content is kept at the present level. It is shown that the lower temperature may increase the depth and the extension of the ozone hole by extending the area amenable to PSC formation. It may be argued that this effect, coupled with an increasing amount of chlorine, may produce a positive feedback on the ozone destruction.
Simulations of scenario earthquake ground motions play an increasingly important role in improving seismic hazard assessment. In the San Francisco Bay Area (SFBA) the USGS has developed a 3D geologic/seismic model. Simulations of large magnitude (M) 7 Hayward Fault scenario earthquakes show dramatic differences in peak ground motions across the fault, arising from wavespeed differences, which highlights the need to evaluate the model. Our comparison of observed and simulated waveforms from M 4 SFBA earthquakes reveals bias in the travel times of direct S-waves for paths crossing the East Bay Hills, east of Hayward Fault, suggests that the shear wavespeeds in the upper crust of the USGS 3D model are too low, resulting in large delay times of S-wave arrivals at the East Bay Hills. We also retrieved noise-derived Green's functions (NGFs) and examined the path-specific bias in the USGS 3D model. We systematically compared synthetic waveforms with the 3D velocity model and observed NGFs. Our result shows that the time lags between the synthetic and observed data ranging from 2-5 seconds for the ray paths propagating the East Bay Hills. This result also suggests that the SFBA 3D model generates arrivals that are too slow. Additionally, we evaluate the Rayleigh-wave ellipticity to evaluate the subsurface velocity beneath seismic stations. A lateral variation in Rayleigh-wave ellipticity across the Hayward fault is observed, which will contribute to refine the SFBA 3D velocity model.
Not provided.
The Caerbannog Point Clouds provide point-sampled 3D models occluded in clouds of points. We synthesized the 3D point clouds from polygonal models, point-sampling the models and surrounding them in a point cloud such that the shape of the model is occluded in any 2D projection. We obscure our model-of-interest by repeatedly surrounding it with an amorphous cloud of points, giving the overall point cloud a structure of organic nature, like that of shrubbery. In a user study, participants were significantly better at identifying the models when visualized as 3D scatterplots under rotation than in axis-aligned 2D scatterplots. We provide three point-clouds, both occluded and unoccluded: the Stanford bunny, Utah Teapot, and OSG Cow.
Abstract not provided
3D models of stellar convection are necessary for accurate description of state of model near collapse. Convection is occurring at many scales with the large dominant mode near collapse. 3D instabilities can affect flow properties and mass entrainment. Mach number profiles show favorable conditions for explosion. 3D rotating progenitor models are ALSO necessary. The redistribution of AM diverges from the MESA model, with implications for remnant. There is turbulent transport of AM in convective shell regions. Multi-D models can provide input for successful CCSN models. Larger non-radial kinetic energy exists when using multi-D progenitor input. The 3D CCSN model showed prompt convection and asymmetric shock runaway. Explosion properties suggest robust impact on multi-messenger signals.
Conference presentation at International Meeting for Applied Geoscience & Energy (IMAGE), Houston, Texas, August 28 – September 1, 2023. Using field data from hydrocarbon and CO 2 applications, we illustrate the importance of a workflow and adaption to the target on hand. Verifying the geophysical acquisition and processing steps with 3D modeling and checking them against a 3D anisotropic log-derived model maintains confidence in the workflow and minimizes the influence on the data. This allows us to predict data validity and to certify the data with respect to the borehole logs.
The subsurface geology beneath the Pajarito Plateau is critical to understanding the seismic hazard of the Pajarito Fault System, yet our understanding of this geology is relatively poor. While previous 3D geologic framework models of the area have been created for the purposes of understanding hydrogeologic flow, they are inadequate for the purposes of understanding the Pajarito Fault System. The specific challenges of using oil and gas software for this purpose include: (1) the geologic complexities resulting from volcanism and tectonism; (2) a need for a high level of stratigraphic detail over a large area; (3) a near complete lack of seismic data; and (4) sparse wellbore data. Presented here is a workflow that handles these challenges of adapting commercially available software used by the oil and gas industries to this seismic hazard problem.
Using self-consistent three-dimensional (3D) N-body simulations, we investigate the physical properties of nonaxisymmetric features in a disk galaxy created by a tidal interaction with its companion. The primary galaxy consists of a stellar disk, a bulge, and a live halo, corresponding to Milky-Way-type galaxies, while the companion is represented by a halo alone. We vary the companion mass and the pericenter distance to explore situations with differing tidal strength parameterized by either the relative tidal force P or the relative imparted momentum S. We find that the formation of a tidal tail in the outer parts requires P≳0.05 or S≳0.07. A stronger interaction results in a stronger, less wound tail that forms earlier. Similarly, a stronger tidal forcing produces stronger, more loosely wound spiral arms in the inner parts. The arms are approximately logarithmic in shape, with both amplitude and pitch angle decaying with time. The derived pattern speed decreases with radius and is close to the Ω−κ/2 curve at late time, with Ω and κ denoting the angular and epicycle frequencies, respectively. This suggests that the tidally induced spiral arms are most likely kinematic density waves weakly modified by self-gravity. Compared to the razor-thin counterparts, arms in the 3D models are weaker, have a smaller pitch angle, and wind and decay more rapidly. The 3D density structure of the arms is well described by the concentrated and sinusoidal models when the arms are in the nonlinear and linear regimes, respectively. We demonstrate that dynamical friction between interacting galaxies transfers the orbital angular momentum of one galaxy to the spin angular momentum of the companion halo.
The software Fire Risk Investigation in 3D (FRI3D) has been developed over the last 2 years to integrate 3D spatial modeling with existing fire probabilistic risk assessment (PRA) models and fire simulation codes. The goal of this research and development is to automate many of the fire analysis manual tasks to reduce industry efforts in the initial fire modeling and operational costs for the model maintenance and evaluations required during normal plant operations. The tasks for Fiscal Year (FY) 2021 include first testing the FRI3D modeling capabilities by importing an industry fire model into FRI3D and making a 3D model of a complex/high-risk significant area. (For this work, the switchgear room was chosen.) Then, the second task of FY 2021 is to develop a dynamic fire PRA process that can help optimize traditional fire PRA models. The switchgear room model will be used for the dynamic fire PRA work. This report describes the work and insights learned when using FRI3D software to model both a Nuclear Regulatory Report (NUREG) example models and a full industry switchgear room.
Outline - Process-Structure-Performance Framework - Process-Specific Defects in selective laser melting (SLM) additive manufacturing (AM) process - Idealized Keyhole/Entrapped gas pores - Idealized Lack of Fusion (LoF) pores - Influence of Pore Geometry on Strain Localization - 3D Model of LoF pores - Partial validation of the 3D LoF 2 model
Many important engineering and scientific applications such as cement slurries, foams, crude oil, and granular avalanches involve the concept of yield stress. Therefore, modeling yield stress fluids in different flow configurations, including the accurate prediction of the yield surface, is important. In this paper, we present a computational model based on the finite element method to study the flow of yield stress fluids in a thin mold and compare the results with data from flow visualization experiments. Here, we use the level set method to describe the interface between the filling fluid and air. We use polypropylene glycol as a model Newtonian fluid and Carbopol for the model yield stress fluid, as the Carbopol solution demonstrates yielding without thixotropy. To describe the yielding and shear-thinning behavior, we use a generalized Newtonian constitutive equation with a Bingham–Carreau–Yasuda form. We compare the results obtained from the mold filling experiments with the results from the three-dimensional (3D) model and from a reduced-order Hele-Shaw (HS) model that is two-dimensional, including the effect of shear-thinning along the thin direction only approximately. We show that both the 3D and the HS model can capture the experimental meniscus shape reasonably well for all the fluids considered at three different flow rates. This indicates that the shape evolution is insensitive to the dimensionality of the model. However, the viscosity and yield surfaces predicted by the 3D and HS models are different. The HS model underestimates the high viscosity and unyielded regions compared to the estimation by the 3D model.
Poster presented by NETL at ICGH, 06/17 - Abstract:<p><br></p><p>We are currently constructing an Earth model of the Terrebonne Basin for comprehensive 3D modeling. Well data and a high-quality 3D seismic survey underpin the geologic conditions for the model. Additional model constraints include basal heat flow, seafloor temperature, and paleowater depth through time. Our initial investigation comprises the development of a “pseudo” 2D model centered on the deep WR313 #001 well. We are using this test case to establish the workflow for integrating BPSM with gas hydrate analyses; this will eventually be applied to the entire mini-basin. </p><p><br></p><p>A key aspect of investigating the gas hydrate system with BPSM is integrating it with the other active petroleum system(s) in the region, because the methane in hydrate deposits can derive from both biogenic and thermogenic sources. Some of the questions we may consider are: What is the history of these deposits with respect to the gas hydrate stability zone, recognizing that this area is not one characterized by simple continuous sedimentation and subsidence, but instead by periodic uplift, tilting, plate rotation, sea-level change, salt tectonics, and evolution of temperature gradients. </p><p><br></p><p>In conclusion, BPSM has been called the ‘great integrator’ in petroleum exploration. Development of a BPSM model of gas hydrates in the Terrebonne Basin of the northern GoM will provide a vehicle within which to integrate other early exploration and assessment research being conducted on gas hydrates, such as in frontier areas like Japan, as a potential resource likely to provide many decades of energy if proven to be commercially producible</p>
To fabricate porous implants with improved biocompatibility and mechanical properties that are matched to their application using selective laser melting (SLM), flow within the mold and compressive properties and performance of the porous structures must be comprehensively studied. Parametric modeling was used to build 3D models of octahedron and hexahedron structures. Finite element analysis was used to evaluate the mold flow and compressive properties of the parametric porous structures. A DiMetal-100 SLM molding apparatus was used to manufacture the porous structures and the results evaluated by light microscopy. The results showed that parametric modeling can produce robust models. Square structures caused higher blood cell adhesion than cylindrical structures. “Vortex” flow in square structures resulted in chaotic distribution of blood elements, whereas they were mostly distributed around the connecting parts in the cylindrical structures. No significant difference in elastic moduli or compressive strength was observed in square and cylindrical porous structures of identical characteristics. Hexahedron, square and cylindrical porous structures had the same stress–strain properties. For octahedron porous structures, cylindrical structures had higher stress–strain properties. Using these modeling and molding results, an important basis for designing and the direct manufacture of fixed biological implants is provided.
We present a novel three-dimensional (3D) model of coronal mass ejections (CMEs) that unifies all key evolutionary aspects of CMEs and encapsulates their 3D magnetic field configuration. This fully analytic model is capable of reproducing the global geometrical shape of a CME with all major deformations taken into account, i.e., deflection, rotation, expansion, “pancaking,” front flattening, and rotational skew. Encapsulation of 3D magnetic structure allows the model to reproduce in-situ measurements of magnetic field for trajectories of spacecraft-CME encounters of any degree of complexity. As such, the model can be used single-handedly for the consistent analysis of both remote and in-situ observations of CMEs at any heliocentric distance. We demonstrate the latter by successfully applying the model for the analysis of two CMEs.