Non-isothermal Multiphase 3D Reactive Transport Modeling of an Engineered Barrier System in Clay Rock.
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.
Abstract not provided.
Explore the source record for details and available documents.
We study the three-dimensional $SU(2)$ principal chiral model (PCM) using different tensor renormalization group methods based on the triad and anisotropic decomposition of the tensor. The tensor network representation is formulated based on the character expansion of the Boltzmann weight. We compare the average action obtained using these two tensor network algorithms and confirm that the resulting critical coupling and exponent are comparable with the recent estimations from the Monte Carlo methods.
Explore the source record for details and available documents.
This dataset is a model archive of the paper A bespoke model of Arctic river basins based on hillslope delineation (in prep), which introduces a watershed decomposition and parameterization method for large scale permafrost hydrology simulation. With this dataset, this study aims to address the research question: whether a computationally efficient hillslope-based modeling framework can reliably simulate discharge at Arctic river-basin scales. This dataset contains model input and output data for five modeling scenarios at a study site located in the Sagavanirktok River basin. The five modeling scenarios include three modeling cases under temperate conditions using full 3D, decomposed 3D, and decomposed 2D modeling strategies; and two modeling cases under actual Arctic conditions with permafrost using full 3D and decomposed 2D modeling strategies. Simulations were performed using the Advanced Terrestrial Simulator (ATS, v1.6 for three temperate scenarios and v1.5 for two Arctic scenarios), a physics-rich integrated surface–subsurface hydrologic model with cryo-hydrology features. For the three temperate models, simulations were conducted for the period of 10/01/1993 - 09/30/2002; and for the two Arctic models, simulations were conducted for the period of 01/01/1994 - 12/31/2002. To facilitate reproducibility of simulations, all datasets are organized hierarchically. The dataset contains: (1) Mesh files (.exo) for full 3D model, decomposed 3D models, and decomposed 2D models, located in huc/190604020802_gauge15906000/mesh/. Mesh files can be visualized through Paraview or read by Python. (2) Climate forcings (.h5) for full 3D model and decomposed 3D/2D models are located in huc/190604020802_gauge15906000/daymet_onePiece/, and huc/190604020802_gauge15906000/vp_pr_revised_daymet_1980_2006_with_wind/ separately. Accessible by Python. (3) Raw measured gage discharge (.csv) from USGS, located in huc/190604020802_gauge15906000/gaged_basin15906000_discharge_usgs/. Accessible by Python. (4) Delineated subdomain raster (.tif) and shape files (.shp), and the final parameterized results (.npy) for decomposed models, located in huc/190604020802_gauge15906000/data_preprocessed-meshing. Accessible by Python. (5) Temperate models are located in nonpermaf_huc190604020802_gauge15906000/, which includes three cases: decomposed 2D models (inside model_0*-hillslope_*), decomposed 3D models (inside model_1*-subcatchment_*), and full 3D model (inside model_2*-onepiece_*). Two step spin-up results (checkpoint_final.h5) are located in model_*1-*_spinup_steadystate and model_*2-*_spinup_cycle, separately, which are used to initialize real transient models. The input files (.xml) and output results (.dat) of the real transient models are located in model_*3-*_transient/. Especially, for two example hillslope models (ID=-11 and 11), additional h5py files are included in model_03-hillslope_transient/hillslope-11/, model_03-hillslope_transient/hillslope11, model_13-subcatchment_transient/subcatchment-11/, model_13-subcatchment_transient/subcatchment/11, respectively, which are used to plot the saturation figure (Figure 5) in the manuscript. Accessible by Python. (6) Arctic models are located in huc190604020802_gauge15906000/, which includes two cases: decomposed 2D models (inside model_04-hillslope_transient), and full 3D model (inside model_05-onepiece_transient_mannp1_ra). Three step spin-up results (checkpoint_final.h5) are located in model_01-column_freezeup/, model_02-column_spinup/, model_03-hillslope_spinup/, respectively, which are used to initialize real 2D transient hillslope models. The input files (.xml) and output results (.dat) of transient 2D hillslope models are located in model_04-hillslope_transient/. The input files (.xml) and output results (.dat) of the full 3D transient model is located in model_05-onepiece_transient_mannp1_ra/. The full 3D transient model is initialized by model_02-column_spinup/. Accessible by Python. (7) The MOSART routed discharge results (.csv) under Arctic conditions is located in huc190604020802_gauge15906000/MOSART/. Accessible by Python. (8) All Python codes (.py) used to parameterize full 3D model to decomposed 2D models are located in script/. These codes fit with watershed workflow (a watershed delineation tool) v1.4 under the branch gaob/v1.4 from https://github.com/gaobhub/watershed-workflow.git.
In an effort to miniaturize a 3D imaging system, we created two viewpoints in a single objective lens camera. This was accomplished by placing a pair of Complementary Multi-band Bandpass Filters (CMBFs) in the aperture area. Two key characteristics about the CMBFs are that the passbands are staggered so only one viewpoint is opened at a time when a light band matched to that passband is illuminated, and the passbands are positioned throughout the visible spectrum, so each viewpoint can render color by taking RGB spectral images. Each viewpoint takes a different spectral image from the other viewpoint hence yielding a different color image relative to the other. This color mismatch in the two viewpoints could lead to color rivalry, where the human vision system fails to resolve two different colors. The difference will be closer if the number of passbands in a CMBF increases. (However, the number of passbands is constrained by cost and fabrication technique.) In this paper, simulation predicting the color mismatch is reported.
A computer program that generates three-dimensional (3D) finite element models for cracked 3D solids was written. This computer program, gensurf, uses minimal input data to generate 3D finite element models for isotropic solids with elliptic or part-elliptic cracks. These models can be used with a 3D finite element program called surf3d. This report documents this mesh generator. In this manual the capabilities, limitations, and organization of gensurf are described. The procedures used to develop 3D finite element models and the input for and the output of gensurf are explained. Several examples are included to illustrate the use of this program. Several input data files are included with this manual so that the users can edit these files to conform to their crack configuration and use them with gensurf.
Additive manufacturing (AM), while commonly used for rapid prototyping and creating components with complex geometries, has not been widely adopted for critical applications across the aerospace, automotive, defense, energy, and medical industries. This is, in part, due to the challenges of controlling flaws and uncertainty in the mechanical behavior of additively manufactured components. In recent years, there has been an increase in research aimed at predicting the final mechanical properties of additively manufactured components during the printing process. To address these issues, a 3D-CNN model was trained using low-cost in situ visible-light camera data, anomaly classifications, and the chosen process parameters to predict the ultimate tensile strength (UTS), yield strength (YS), total elongation (TE), and uniform elongation (UE). The 3D-CNN layers of the model employed attention mechanisms to prioritize features in the data, thereby improving prediction accuracy. Furthermore, the effect of each process parameter and anomaly class is investigated using attention-based dynamic sigmoid weighted gates to interpret the influence each class has on the final prediction. Different combinations of the in situ data were fed into the 3D-CNN, with varying amounts of image layers, to determine the ideal combination for predicting mechanical properties in situ. Here, the 3D-CNN model achieved mean absolute percentage errors (MAPE) below 5% for both UTS and YS while using only a single camera input and under half of the available image layers.
In 2009, the Gulf of Mexico (GOM) Gas Hydrates Joint-Industry-Project (JIP) Leg II drilling program confirmed that gas hydrate occurs at high saturations within reservoir-quality sands in the GOM. A comprehensive logging-while-drilling dataset was collected from seven wells at three sites, including two wells at the Walker Ridge 313 site. By constraining the saturations and thicknesses of hydrate-bearing sands using logging-while-drilling data, two-dimensional (2D), cylindrical, r-z and three-dimensional (3D) reservoir models were simulated. The gas hydrate occurrences inferred from seismic analysis are used to delineate the areal extent of the 3D reservoir models. Numerical simulations of gas production from the Walker Ridge reservoirs were conducted using the depressurization method at a constant bottomhole pressure. Results of these simulations indicate that these hydrate deposits are readily produced, owing to high intrinsic reservoir-quality and their proximity to the base of hydrate stability. The elevated in situ reservoir temperatures contribute to high (5–40 MMscf/day) predicted production rates. The production rates obtained from the 2D and 3D models are in close agreement. To evaluate the effect of spatial dimensions, the 2D reservoir domains were simulated at two outer radii. The results showed increased potential for formation of secondary hydrate and appearance of lag time for production rates as reservoir size increases. Similar phenomena were observed in the 3D reservoir models. The results also suggest that interbedded gas hydrate accumulations might be preferable targets for gas production in comparison with massive deposits. Hydrate in such accumulations can be readily dissociated due to heat supply from surrounding hydrate-free zones. Special cases were considered to evaluate the effect of overburden and underburden permeability on production. The obtained data show that production can be significantly degraded in comparison with a case using impermeable boundaries. The main reason for the reduced productivity is water influx from the surrounding strata; a secondary cause is gas escape into the overburden. The results dictate that in order to reliably estimate production potential, permeability of the surroundings has to be included in a model.
Abstract. Electrical resistivity tomography (ERT) is a broadly accepted geophysical method for subsurface investigations. Interpretation of field ERT data usually requires the application of computationally intensive forward modeling and inversion algorithms. For large-scale ERT data, the efficiency of these algorithms depends on the robustness, accuracy, and scalability on high-performance computing resources. In this regard, we present a robust and highly scalable implementation of forward modeling and inversion algorithms for ERT data. The implementation is publicly available and developed within the framework of PFLOTRAN, an open-source, state-of-the-art massively parallel subsurface flow and transport simulation code. The forward modeling is based on a finite-volume discretization of the governing differential equations, and the inversion uses a Gauss–Newton optimization scheme. To evaluate the accuracy of the forward modeling, two examples are first presented by considering layered (1D) and 3D earth conductivity models. The computed numerical results show good agreement with the analytical solutions for the layered earth model and results from a well-established code for the 3D model. Inversion of ERT data, simulated for a 3D model, is then performed to demonstrate the inversion capability by recovering the conductivity of the model. To demonstrate the parallel performance of PFLOTRAN's ERT process model and inversion capabilities, large-scale scalability tests are performed by using up to 131 072 processes on a leadership class supercomputer. These tests are performed for the two most computationally intensive steps of the ERT inversion: forward modeling and Jacobian computation. For the forward modeling, we consider models with up to 122 ×106 degrees of freedom (DOFs) in the resulting system of linear equations and demonstrate that the code exhibits almost linear scalability on up to 10 000 DOFs per process. On the other hand, the code shows superlinear scalability for the Jacobian computation, mainly because all computations are fairly evenly distributed over each process with no parallel communication.
Based on the reduced diffusion mechanism for producing Forbush decreases (Fds) in the heliosphere, we constructed a three-dimensional (3D) diffusion barrier, and by incorporating it into a stochastic differential equation (SDE) based time-dependent, cosmic-ray transport model, a 3D numerical model for simulating Fds is built and applied to a period of relatively quiet solar activity. This SDE model generally corroborates previous Fd simulations concerning the effects of the solar magnetic polarity, the tilt angle of the heliospheric current sheet (HCS), and cosmic-ray particle energy. Because the modulation processes in this 3D model are multi-directional, the barrier’s geometrical features affect the intensity profiles of Fds differently. We find that both the latitudinal and longitudinal extent of the barrier have relatively fewer effects on these profiles than its radial extent and the level of decreased diffusion inside the disturbance. We find, with the 3D approach, that the HCS rotational motion causes the relative location from the observation point to the HCS to vary, so that a periodic pattern appears in the cosmic-ray intensity at the observing location. Correspondingly, the magnitude and recovery time of an Fd change, and the recovering intensity profile contains oscillation as well. Investigating the Fd magnitude variation with heliocentric radial distance, we find that the magnitude decreases overall and, additionally, that the Fd magnitude exhibits an oscillating pattern as the radial distance increases, which coincides well with the wavy profile of the HCS under quiet solar modulation conditions.
Two-dimensional models assuming axisymmetry are an economical way to explore the long-term evolution of black hole accretion disks, but they are only realistic if the feedback of the nonaxisymmetric turbulence on the mean momentum and magnetic fields is incorporated. Dynamo terms added to the 2D induction equation should be calibrated to 3D MHD simulations. For generality, the dynamo tensors should be calibrated as functions of local variables rather than explicit functions of spatial coordinates in a particular basis. In this paper, we study the feedback of non-axisymmetric features on the 2D mean fields using a global 3D, relativistic, Cartesian simulation from the Illinois GRMHD code. We introduce new methods for estimating overall dynamo alpha and turbulent diffusivity effects as well as measures of the dominance of non-axisymmetric components of energies and fluxes within the disk interior. We attempt closure models of the dynamo EMF using least squares fitting, considering both models where coefficient tensors are functions of space and more global, covariant models. None of these models are judged satisfactory, but we are able to draw conclusions on what sorts of generalizations are and are not promising.
A novel ground heat exchanger design, the Underground Thermal Battery (UTB), was proposed as a low-cost alternative to the conventional vertical bore ground heat exchanger (VBGHE) for the application of ground source heat pumps (GSHPs). The UTB is designed to be installed in a shallow borehole (less than 6 m deep) and thus could cost much less than the conventional VBGHE, which usually is installed in vertical bores 60 m deep. By utilizing natural convection of water and phase change materials, the UTB tempers its temperature change in response to thermal loads, which helps improve the efficiency of GSHPs. A one-dimensional (1D) model of the UTB has been developed and validated with the measured performance data of a small-scale UTB, as well as the simulation results of a more detailed three dimensional (3D) numerical model, which accounts for both heat transfer and fluid dynamics in the UTB. Further, this 1D model is computationally much more efficient than the 3D model and thus can simulate the annual performance of the UTB with acceptable computation time. The 1D model has been used to evaluate the performance of the UTB so that it can be compared with the conventional VBGHE.
A 3D radiative transfer model is used to investigate the relationship between spectral indices and fraction of absorbed photosynthetically active radiation (PAR) in horizontally heterogeneous vegetation canopies. Canopy reflection at optical wavelengths and PAR absorption are simulated. Data obtained indicate that the leaf area index of a canopy is less of an instructive parameter than the ground cover and clump leaf area index for these canopies. It is found that the relationship between the normalized difference vegetation index and fraction of absorbed PAR is almost linear and independent of spatial heterogeneity.
Shade obstructions can significantly impact the performance of photovoltaic (PV) systems. Although there are many models for partially shaded PV arrays, there is a lack of information available regarding their accuracy and uncertainty when compared with actual field performance. This work assesses the recorded performance of 46 residential PV systems, equipped with either string-level or module-level inverters, under a variety of shading conditions. We compare their energy production data to annual PV performance predictions, with a focus on the practical models developed here for National Renewable Energy Laboratory's system advisor model software. This includes assessment of shade extent on each PV system by using traditional onsite surveys and newer 3D obstruction modelling. The electrical impact of shade is modelled by either a nonlinear performance model or assumption of linear impact with shade extent, depending on the inverter type. When applied to the fleet of residential PV systems, performance is predicted with median annual bias errors of 2.5% or less, for systems with up to 20% estimated shading loss. The partial shade models are not found to add appreciable uncertainty to annual predictions of energy production for this fleet of systems but do introduce a monthly root-mean-square error of approximately 4%-9% due to seasonal effects. Here the use of a detailed 3D model results in similar or improved accuracy over site survey methods, indicating that, with proper description of shade obstructions, modelling of partially shaded PV arrays can be done completely remotely, potentially saving time and cost.
The brain-on-a-chip project aims to study brain function by developing a 3D cell culture model that mimics the cellular complexity and architecture of the brain. It will be important to identify the spatial placement of cells within the 3D cultures since cell-cell contact and interaction can influence cellular function. To characterize cellular composition in 3D cultures, immunocytochemistry (ICC) was used to identify neurons and glial cells (e.g. astrocytes, oligodendrocytes, and microglia). The optimization of ICC was conducted to improve labelling of specific cell types and their maturation state.