Search NASA⌕ Search

Engineering topics

Kim, Taehyun

Publications and source records attributed to Kim, Taehyun.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

DECOVALEX-2023: Task C Final Report

The Full-scale Emplacement (FE) heater experiment at the Mont Terri Underground Rock Laboratory (URL) was designed and conducted by Nagra to replicate an emplacement tunnel of Nagra’s reference repository design at 1:1 scale. Alongside testing the technical feasibility of constructing disposal tunnels, emplacing waste containers in the tunnels and then backfilling them, the main goals of the FE experiment are (1) to obtain a better understanding of the coupled effects of induced thermo-hydro-mechanical (THM) processes that may occur and (2) to validate existing coupled THM models (Müller et al., 2017). A key aspect of ensuring safety for repositories located in low-permeability rock involves minimizing any damage to the rock itself, thereby preserving its integrity and promoting a stable environment Amongst a number of processes that could damage the rock is the increase in pore pressure due to thermal loading caused by heat emitted from the waste. To reduce the potential damage of the rock, it is important to analyse the evolution of heat over time due to the heat load of the containers and assess possible consequences by coupled THM models. The aim of Task C of DECOVALEX-2023 was to build 3D numerical models of the FE experiment, focussing in particular on the heating induced pore pressure change in the Opalinus Clay. Data from a large number of sensors were available from the FE experiment for model comparison. These sensors measured temperature and relative humidity in the bentonite around the heaters, and temperature, pressure and displacement/strain in the surrounding Opalinus clay. Data were available from the start of excavation (April 2012) up to August 2020 for most sensors (more than 5 years from the start of heating in December 2014). To fulfil the overall aim of the task, the work was broken down into a number of steps, starting with simpler models to build confidence in each team’s approach and then moving to more complex models that better represent the FE experiment. Step 0 consisted of 2D benchmark models, gradually increasing the number of processes that are represented from thermal (T) only models in Step 0a, to coupled thermal hydraulic (TH) models in Step 0b with a representation of changing porosity, to coupled thermo-hydro-mechanical (THM) models in Step 0c, where porosity changes are calculated by the mechanical model. A detailed specification of processes, parameters, initial and boundary conditions was provided for this step, with the ambition that all teams would work towards close agreement in their model results, thus building confidence in the model implementations. vi It was not straightforward to achieve agreement between the teams, so additional steps (Step 0b2, 0b3, 0c2, 0c3) were added along with derivation of some analytical solutions against which the models could be compared. The reasons for the differences between teams were investigated and found to be caused primarily by different conceptual model assumptions (including temperature dependence of the thermal expansion of water), different model formulations (including porosity evolution) and differences in modelled domain sizes, boundary conditions and grid discretisation. This demonstrates that comparisons between multiple modelling teams and/or comparison with analytical results and experimental data are highly beneficial in providing an indication of uncertainty in model predictions. At the conclusion of Step 0, almost all teams had achieved a close agreement in model results and those that had not achieved an agreement knew the reason for this. Step 1 moved from 2D models to 3D models of the FE experiment without adding technical features like shotcrete or EDZ, and only considering the heating phase. Initially the 3D model was tightly specified to continue to build confidence in the model implementations (Step 1a). The results of Step 1a were compared to the data from the FE-experiment without the teams seeing the data. The teams were then provided with a sub-set of the data from the FE-experiment and invited to consider how best to use the large dataset for model comparison (Step 1b). Teams were then asked to use the data provided to calibrate their models, only changing material property values rather than adding features or processes to their models (Step 1c). In Step 1, teams were asked to only model the heating phase of the experiment, so pressure in the Opalinus Clay was reported as change in pressure since the initial conditions were specified rather than modelled. The change from 2D to 3D models was accompanied by an increase in the dispersion of results between the teams. Some of this was resolved during the task, but some remained and is potentially due to model discretisation. Calibration of parameters was useful in improving the fit of the models to the data but the remaining differences indicated that the models were missing features or processes. In Step 2, the teams were asked to update their models with additional features and processes as well as calibrating parameters to try and improve the fit of the models to the data. Teams were encouraged to represent ventilation of the open FE tunnel prior to backfilling with heaters and bentonite and in Step 2, the absolute pressure in the Opalinus Clay was compared between the teams. Teams took different approaches, but there was consideration of adding shotcrete and an EDZ into the model, representing stress change during excavation and different approaches to modelling ventilation of the FE tunnel. Overall, the documented results showed a very good agreement for temperature. The results for porewater pressure evolution showed a significant improvement for most teams compared to Step 1c with a good agreement to the measurements for several teams whereas some teams overpredicted the pressure increase and others overpredicted the drainage effect especially for the sensors close to the heater. Step 3 was an opportunity for teams to use the models developed in Step 1 and Step 2 to make predictions about the temperature and pressure changes that will be expected at the FE experiment over the next few years in light of the planned changes in thermal output of the heaters.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Stability of hydrides in sub-Neptune exoplanets with thick hydrogen-rich atmospheres

Many sub-Neptune exoplanets have been believed to be composed of a thick hydrogen-dominated atmosphere and a high-temperature heavier-element-dominant core. From an assumption that there is no chemical reaction between hydrogen and silicates/metals at the atmosphere–interior boundary, the cores of sub-Neptunes have been modeled with molten silicates and metals (magma) in previous studies. In large sub-Neptunes, pressure at the atmosphere–magma boundary can reach tens of gigapascals where hydrogen is a dense liquid. A recent experiment showed that hydrogen can induce the reduction of Fe 2+ in (Mg,Fe)O to Fe 0 metal at the pressure–temperature conditions relevant to the atmosphere–interior boundary. However, it is unclear whether Mg, one of the abundant heavy elements in the planetary interiors, remains oxidized or can be reduced by H. Our experiments in the laser-heated diamond-anvil cell found that heating of MgO + Fe to 3,500 to 4,900 K (close to or above their melting temperatures) in an H medium leads to the formation of Mg 2 FeH 6 and H 2 O at 8 to 13 GPa. At 26 to 29 GPa, the behavior of the system changes, and Mg–H in an H fluid and H 2 O were detected with separate FeH x . The observations indicate the dissociation of the Mg–O bond by H and subsequent production of hydride and water. Therefore, the atmosphere–magma interaction can lead to a fundamentally different mineralogy for sub-Neptune exoplanets compared with rocky planets. The change in the chemical reaction at the higher pressures can also affect the size demographics (i.e., “radius cliff”) and the atmosphere chemistry of sub-Neptune exoplanets.

58 GEOSCIENCES↗
DOE OSTI2023

MHz free electron laser x-ray diffraction and modeling of pulsed laser heated diamond anvil cell

A new diamond anvil cell experimental approach has been implemented at the European x-ray Free Electron Laser, combining pulsed laser heating with MHz x-ray diffraction. In this report we use this setup to determine liquidus temperatures under extreme conditions, based on the determination of time-resolved crystallization. The focus is on a Fe-Si-O ternary system, relevant for planetary cores. This time-resolved diagnostic is complemented by a finite-element model, reproducing temporal temperature profiles measured experimentally using streaked optical pyrometry. This model calculates the temperature and strain fields by including (i) pressure and temperature dependencies of material properties, and (ii) the heat-induced thermal stress, including feedback effect on material parameter variations. Making our model more realistic, these improvements are critical as they give 7000 K temperature differences compared to previous models. Laser intensities are determined by seeking minimal deviation between measured and modeled temperatures. Combining models and streak optical pyrometry data extends temperature determination below detection limit. The presented approach can be used to infer the liquidus temperature by the appearance of SiO 2 diffraction spots. In addition, temperatures obtained by the model agree with crystallization temperatures reported for Fe–Si alloys. Our model reproduces the planetary relevant experimental conditions, providing temperature, pressure, and volume conditions. Those predictions are then used to determine liquidus temperatures at experimental timescales where chemical migration is limited. This synergy of novel time-resolved experiments and finite-element modeling pushes further the interpretation capabilities in diamond anvil cell experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2021

X-ray free electron laser heating of water and gold at high static pressure

Probing of reactive materials such as H 2 O ices and fluids at the high pressures and temperatures of planetary interiors is limited by unwanted chemical reactions and confinement failure. Faster experiments can mitigate such issues, but the common approach of adiabatic compression limits the conditions achieved. This study demonstrates a fast experimental strategy for the creation and probing of selected extreme states using static compression coupled with ultrafast X-ray laser heating. Indirect X-ray heating of H 2 O through the use of a gold absorber is evidenced by sample melting inferred from textural changes in the H 2 O diffraction lines and inter-dispersion of gold and H 2 O melts. Coupled with numerical analysis of femtosecond energy absorption, thermal equilibration, and heat transfer, all evidence indicates that temperatures in excess of an electron volt have been reached in the H 2 O at high pressure. Even after repeated heating, samples stayed chemically unchanged from the starting material.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

X-ray Free Electron Laser-Induced Synthesis of ε -Iron Nitride at High Pressures

The ultrafast synthesis of ε-Fe 3 N 1+x in a diamond-anvil cell (DAC) from Fe and N 2 under pressure was observed using serial exposures of an X-ray free electron laser (XFEL). When the sample at 5 GPa was irradiated by a pulse train separated by 443 ns, the estimated sample temperature at the delay time was above 1400 K, confirmed by in situ transformation of α- to γ-iron. Ultimately, the Fe and N 2 reacted uniformly throughout the beam path to form Fe 3 N 1.33 , as deduced from its established equation of state (EOS). We thus demonstrate that the activation energy provided by intense X-ray exposures in an XFEL can be coupled with the source time structure to enable exploration of the time-dependence of reactions under high-pressure conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

The stability of subducted glaucophane with the Earth’s secular cooling

The blueschist to eclogite transition is one of the major geochemical–metamorphic processes typifying the subduction zone, which releases fluids triggering earthquakes and arc volcanism. Although glaucophane is an index hydrous mineral for the blueschist facies, its stability at mantle depths in diverse subduction regimes of contemporary and early Earth has not been experimentally determined. Here, we show that the maximum depth of glaucophane stability increases with decreasing thermal gradients of the subduction system. Along cold subduction geotherm, glaucophane remains stable down ca. 240 km depth, whereas it dehydrates and breaks down at as shallow as ca. 40 km depth under warm subduction geotherm or the Proterozoic tectonic setting. Our results imply that secular cooling of the Earth has extended the stability of glaucophane and consequently enabled the transportation of water into deeper interior of the Earth, suppressing arc magmatism, volcanism, and seismic activities along subduction zones.

58 GEOSCIENCES↗
DOE OSTI2020

DECOVALEX-2019 (Task B Final Report)

The DECOVALEX Project is an on-going international research collaboration, established in 1992, to advance the understanding and modeling of coupled Thermal (T), Hydrological (H), Mechanical (M) and Chemical (C) processes in geological in geological systems. DECOVALEX was initially motivated by the recognition that prediction of these coupled effects is an essential part of the performance and safety assessment of geologic disposal systems for radioactive waste and spent nuclear fuel. This document is the final report of Task B which was proposed and coordinated by the Swiss Federal Nuclear Safety Inspectorate (ENSI) and Lawrence Berkeley National Laboratory (LBNL), presenting the technical definitions of the problems studied, approaches applied, achievements made and outstanding issues for future research. Task B focused mainly on modelling of fault activation experiments (FS experiment) recently performed at the Mont Terri underground research laboratory, Switzerland. The FS experiment explores the coupling between fault reactivation in a clay host rock and the potential-enhanced fluid displacement through a previously low-permeability formation. Based on this experiment, the DECOVALEX-2019 Task B addresses how the change in permeability induced by the fault activation and the resulting fluid flow within the fault can be simulated including the dependence on its mechanical behavior. This will support the understanding of the processes during fault activation itself and it can help to determine consequences for contaminant transport due to potentially created permeable flow paths in otherwise low permeability argillaceous rocks.

58 GEOSCIENCES↗
DOE OSTI2020

A new hydrous iron oxide phase stable at mid-mantle pressures

The amount of hydrogen stored in the Earth's interior is important for a range of issues, including the volatile incorporation during the Earth formation and the co-evolution of the atmosphere, the hydrosphere, and the interior. Recent experiments found titanium bearing ε-FeOOH in a hydrous basaltic system at 12–19 GPa and 1300 K. Pyrite-type FeOOH was found to be stable at pressures higher than 80 GPa. These discoveries suggest possible hydrogen storage in the mantle transition zone and in the mantle below 1800 km depths, respectively. However, it remains uncertain how the potential deep hydrogen storage can be connected to the shallower storage. Here, we report a new hydrous iron oxide (η-Fe 12 O 18+x / 2 H x , x ≈ 2) stable at pressures between the stability fields of the ε- and the pyrite-type FeOOH. Our experiment also shows that the new η phase can exist together with the major lower mantle minerals including bridgmanite and periclase, making it an important hydrogen-bearing phase in the Earth's deep interior. Because of its limited H 2 O storage capacity, which is less than 1/6 of the storage capacity of the pyrite-type phase and the ε phase, the stability of the η phase would result in H 2 O loss during water transport in the mid mantle and therefore limit the amount of H 2 O potentially stored in the Fe–O–H system of the lower mantle. As a result, the large channel in the crystal structure of the η phase could provide potential storage sites for other volatile elements in the deep mantle.

58 GEOSCIENCES↗
DOE OSTI2020

Modeling of fluid injection-induced fault reactivation using coupled fluid flow and mechanical interface model

The present study is aimed at developing a numerical model to reproduce coupled hydro-mechanical processes associated with fault reactivation by fluid injection in low permeability rock, as part of the DECOVALEX-2019 project Task B. We proposed a modeling approach for simulating the processes using the TOUGH-FLAC simulator, and modeled a fault reactivation experiment conducted at Mont Terri Rock Laboratory in Switzerland. The first step of the study involved benchmark calculations considering a simplified fault plane and geometry. Fluid flow along a fault was modeled using elements of aperture-sized thickness on the basis of Darcy's law and the cubic law in TOUGH2, whereas the mechanical behavior of a single fault was represented by zero-thickness interface elements in FLAC3D upon which a slip and/or separation is allowed. A methodology to connect a TOUGH2 volume element to a FLAC3D interface element was developed for handling the hydro-mechanical interactions on the fault during fluid injection. Two different fault models for describing the evolutions of hydraulic aperture by elastic fracture opening and failure-induced aperture increase were considered in the benchmark calculations. In the coupling process, the changes in geometrical features and hydrological properties induced by mechanical deformation were continuously updated. The transient responses of the fault and host rock to stepwise pressurization were examined during the simulation. The hydro-mechanical behavior, including the injection flow rate, pressure distribution around the borehole, stress conditions, and displacements in normal and shear directions were monitored in the surrounding rock and along the fault. The results of benchmark calculations suggest that the developed model reasonably represents the hydro-mechanical behavior of a fault and the surrounding rock. This modeling approach was applied to the fault reactivation experiment of the Mont Terri Rock Laboratory. In this interpretive modeling, a parametric study was conducted to examine the effects of input parameters regarding in situ stress and fault properties on the hydro-mechanical responses of the fault to water injection. Then, an optimal parameter set to reproduce the field experiment results was chosen by trial-and-error. The injection flow rate and pressure response during fault reactivation closely matched those obtained at the site, which indicates the capability of the model to appropriately capture the progressive pathway evolution during fault reactivation tests at the site. Finally, the anchor displacements were overestimated by the model, but a fair agreement was obtained in terms of the order of magnitude and the variation tendency.

58 GEOSCIENCES↗
DOE OSTI2020

Intense Reactivity in Sulfur–Hydrogen Mixtures at High Pressure under X-ray Irradiation

Superconductivity near room temperature in the sulfur–hydrogen system arises from a sequence of reactions at high pressures, with X-ray diffraction experiments playing a central role in understanding these chemical–structural transformations and the corresponding S:H stoichiometry. Here we document X-ray irradiation acting as both a probe and as a driver of chemical reaction in this dense hydride system. We observe a reaction between molecular hydrogen (H 2 ) and elemental sulfur (S 8 ) under high pressure, induced directly by X-ray illumination, at photon energies of 12 keV using a free electron laser. The rapid synthesis of hydrogen sulfide (H 2 S) at 0.3 GPa was confirmed by optical observations, spectroscopic measurements, and microstructural changes detected by X-ray diffraction. These results document X-ray induced chemical synthesis of superconductor-forming dense hydrides, revealing an alternative production strategy and confirming the disruptive nature of X-ray exposure in studies on high-pressure hydrogen chalcogenides, from water to high-temperature superconductors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Globular cluster populations: Results including S{sup 4}G late-type galaxies

Using 3.6 and 4.5 μm images of 73 late-type, edge-on galaxies from the S{sup 4}G survey, we compare the richness of the globular cluster populations of these galaxies to those of early-type galaxies that we measured previously. In general, the galaxies presented here fill in the distribution for galaxies with lower stellar mass, M{sub *}, specifically log(M{sub ∗}/M{sub ⊙})<10, overlap the results for early-type galaxies of similar masses, and, by doing so, strengthen the case for a dependence of the number of globular clusters per 10{sup 9}M{sub ⊙} of galaxy stellar mass, T{sub N}, on M{sub *}. For 8.5<log(M{sub ∗}/M{sub ⊙})<10.5 we find the relationship can be satisfactorily described as T{sub N}=(M{sub ∗}/10{sup 6.7}){sup −0.56} when M{sub *} is expressed in solar masses. The functional form of the relationship is only weakly constrained, and extrapolation outside this range is not advised. Our late-type galaxies, in contrast to our early types, do not show the tendency for low-mass galaxies to split into two T{sub N} families. Using these results and a galaxy stellar mass function from the literature, we calculate that, in a volume-limited, local universe sample, clusters are most likely to be found around fairly massive galaxies (M{sub *} ∼ 10{sup 10.8}M{sub ⊙}) and present a fitting function for the volume number density of clusters as a function of parent-galaxy stellar mass. We find no correlation between T{sub N} and large-scale environment, but we do find a tendency for galaxies of fixed M{sub *} to have larger T{sub N} if they have converted a larger proportion of their baryons into stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

THE ODD OFFSET BETWEEN THE GALACTIC DISK AND ITS BAR IN NGC 3906

We use mid-infrared 3.6 and 4.5 μm imaging of NGC 3906 from the Spitzer Survey of Stellar Structure in Galaxies (S{sup 4}G) to understand the nature of an unusual offset between its stellar bar and the photometric center of an otherwise regular, circular outer stellar disk. We measure an offset of ∼910 pc between the center of the stellar bar and photometric center of the stellar disk; the bar center coincides with the kinematic center of the disk determined from previous HI observations. Although the undisturbed shape of the disk suggests that NGC 3906 has not undergone a significant merger event in its recent history, the most plausible explanation for the observed offset is an interaction. Given the relatively isolated nature of NGC 3906 this interaction could be with dark matter substructure in the galaxy's halo or from a recent interaction with a fast moving neighbor that remains to be identified. Simulations aimed at reproducing the observed offset between the stellar bar/kinematic center of the system and the photometric center of the disk are necessary to confirm this hypothesis and constrain the interaction history of the galaxy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

THE SPITZER SURVEY OF STELLAR STRUCTURE IN GALAXIES (S{sup 4}G): MULTI-COMPONENT DECOMPOSITION STRATEGIES AND DATA RELEASE

The Spitzer Survey of Stellar Structure in Galaxies (S{sup 4}G) is a deep 3.6 and 4.5 μm imaging survey of 2352 nearby (<40 Mpc) galaxies. We describe the S{sup 4}G data analysis pipeline 4, which is dedicated to two-dimensional structural surface brightness decompositions of 3.6 μm images, using GALFIT3.0. Besides automatic 1-component Sérsic fits, and 2-component Sérsic bulge + exponential disk fits, we present human-supervised multi-component decompositions, which include, when judged appropriate, a central point source, bulge, disk, and bar components. Comparison of the fitted parameters indicates that multi-component models are needed to obtain reliable estimates for the bulge Sérsic index and bulge-to-total light ratio (B/T), confirming earlier results. Here, we describe the preparations of input data done for decompositions, give examples of our decomposition strategy, and describe the data products released via IRSA and via our web page (www.oulu.fi/astronomy/S4G-PIPELINE4/MAIN). These products include all the input data and decomposition files in electronic form, making it easy to extend the decompositions to suit specific science purposes. We also provide our IDL-based visualization tools (GALFIDL) developed for displaying/running GALFIT-decompositions, as well as our mask editing procedure (MASK-EDIT) used in data preparation. A detailed analysis of the bulge, disk, and bar parameters derived from multi-component decompositions will be published separately.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

A CLASSICAL MORPHOLOGICAL ANALYSIS OF GALAXIES IN THE SPITZER SURVEY OF STELLAR STRUCTURE IN GALAXIES (S{sup 4}G)

The Spitzer Survey of Stellar Structure in Galaxies (S{sup 4}G) is the largest available database of deep, homogeneous middle-infrared (mid-IR) images of galaxies of all types. The survey, which includes 2352 nearby galaxies, reveals galaxy morphology only minimally affected by interstellar extinction. This paper presents an atlas and classifications of S{sup 4}G galaxies in the Comprehensive de Vaucouleurs revised Hubble-Sandage (CVRHS) system. The CVRHS system follows the precepts of classical de Vaucouleurs morphology, modified to include recognition of other features such as inner, outer, and nuclear lenses, nuclear rings, bars, and disks, spheroidal galaxies, X patterns and box/peanut structures, OLR subclass outer rings and pseudorings, bar ansae and barlenses, parallel sequence late-types, thick disks, and embedded disks in 3D early-type systems. We show that our CVRHS classifications are internally consistent, and that nearly half of the S{sup 4}G sample consists of extreme late-type systems (mostly bulgeless, pure disk galaxies) in the range Scd-Im. The most common family classification for mid-IR types S0/a to Sc is SA while that for types Scd to Sm is SB. The bars in these two type domains are very different in mid-IR structure and morphology. This paper examines the bar, ring, and type classification fractions in the sample, and also includes several montages of images highlighting the various kinds of “stellar structures” seen in mid-IR galaxy morphology.

79 ASTRONOMY AND ASTROPHYSICS↗