Search NASA⌕ Search

SEARCH · Search NASA

Results for “surface reconstructions”

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 1,549 records · Page 86
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Intrinsic Lyman Alpha Profiles of High-Velocity G, K, and M Dwarfs

Observations of H I Lyman alpha, the brightest UV emission line of late-type stars, are critical for understanding stellar chromospheres and transition regions, modeling photochemistry in exoplanet atmospheres, and measuring the abundances of neutral hydrogen and deuterium in the interstellar medium. Yet, Lyman alpha observations are notoriously challenging due to severe attenuation from interstellar gas, hindering our understanding of this important emission line's basic morphology. We present high-resolution far- and near-UV spectroscopy of five G, K, and M dwarfs with radial velocities large enough to Doppler shift the stellar Lyman alpha emission line away from much of the interstellar attenuation, allowing the line core to be directly observed. We detect self-reversal in the Lyman alpha emission line core for all targets, and we show that the self-reversal depth decreases with increasing surface gravity. Mg II self-reversed emission line profiles provide some useful information to constrain the Lyman alpha line core, but the differences are significant enough that Mg II cannot be used directly as an intrinsic Lyman alpha template during reconstructions. We show that reconstructions that neglect self-reversal could overestimate intrinsic Lyman alpha fluxes by 60-100% for G and K dwarfs and 40-170% for M dwarfs.

Allison Ann Youngblood↗
NASA NTRS2021

Advances in Land Surface Models and Indicators for Drought Monitoring and Prediction

Millions of people across the globe are affected by droughts every year, and recent droughts have highlighted the considerable agricultural impacts and economic costs of these events. Monitoring the state of droughts depends on integrating multiple indicators that each capture particular aspects of hydrologic impact and various types and phases of drought. As the capabilities of land-surface models and remote sensing have improved, important physical processes such as dynamic, interactive vegetation phenology, groundwater, and snow pack evolution now support a range of drought indicators that better reflect coupled water, energy and carbon cycle processes. In this work, we discuss these advances, including newer classes of indicators that can be applied to improve the characterization of drought onset, severity and duration. We utilize a new model-based drought reconstruction to illustrate the role of dynamic phenology and groundwater in drought assessment. Further, through case studies on flash droughts, snow droughts, and drought recovery, we illustrate the potential advantages of advanced model physics and observational capabilities, especially from remote sensing, in characterizing droughts.

Drought↗
NASA NTRS1999

Solar Particle Event Exposures and Local Tissue Environments in Free Space and on Martian Surface

Solar particle events (SPEs) are a concern to space missions outside Earth s geomagnetic field. The September 29, 1989 SPE is the largest ground-level event since February 23, 1956. It is an iron-rich event for which the spectra are well measured. Because ten times this event matches the ground level data of the February 1956 SPE, it is suggested that an event with ten-times the scaled spectra of the September 29, 1989 SPE be used as a worst case SPE for spacecraft design. For the worst case SPE, the input spectra were reconstructed using Nymmik's (1995) model for protons, the O and Fe ion spectra of Tylka et al. (1997) to evaluate the iron enhancement ratio, and the Solar Energetic Particle Baseline (SEPB) composition of McGuire et al. (1986) for the heavy ions. The necessary transport properties of the shielding materials and the astronaut s body tissues are evaluated using the HZETRN code. Three shield configurations (assumed to be aluminum) are considered: space suit taken as 0.3 g/sq cm, helmet/pressure vessel as 1 g/sq cm, and equipment room of 5 g/sq cm. A shelter is taken as 10 g/sq cm on the Martian surface. The effect of shielding due to the Martian atmosphere is included. The astronaut geometry is taken from the computerized anatomical man (CAM) model.

Kim, M. Y.↗
DOE OSTI2019

Pt 8 cluster on alumina under a pressure of hydrogen: Support-dependent reconstruction from first-principles global optimization

Alumina supported Pt nano clusters under hydrogen environment play a crucial role for many heterogeneous catalysis applications. As such, we conducted grand canonical genetic algorithm (GCGA) simulations for supported Pt 8 clusters in hydrogen gas environment to study the intra-cluster, cluster-support, and cluster-adsorbate interactions. Two alumina surfaces, α-Al 2 O 3 (0001) and $γ$-Al 2 O 3 (100), and two conditions T=600°C, p H2 = 0.1bar and T=25°C, p H2 = 1.0 bar were considered corresponding to low and high hydrogen chemical potential $μ$ H , respectively. The low free energy ensemble of Pt 8 is decorated by a medium (2–12 H), resp. high (20–30 H), number of hydrogen atoms under equilibrium at low $μ$ H , respectively, high $μ$ H , and undergoes different morphology transformations on the two surfaces. On α-Al2O3(0001), Pt8 is mostly 3D but very fluxional in structure at low μH and converts to open one-layer 2D structures with minimal fluxionality at high μH, whereas on $γ$-Al 2 O 3 (100), the exact opposite occurs, and Pt 8 clusters present one-layer 2D shapes at low $μ$ H and switch into compact 3D shapes under high $μ$ H , during which the Pt 8 cluster preserves moderate fluxionality. Further analysis reveals similar Pt–Pt bond length increase when switching from low $μ$ H to high $μ$ H on both surfaces, though morphology transformations are different. Electronic structure analysis shows the existence of bonding interactions between Pt and Lewis acidic Al 3+ sites along with the Pt – O interaction, which implies the necessity to include Al neighbors to discuss the electronic structure of small Pt clusters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Deglaciation of the Greenland and Laurentide ice sheets interrupted by glacier advance during abrupt coolings

Understanding patterns of ice-sheet deglaciation is key for predicting the rate of future ice-sheet melt, yet the processes underlying deglaciation remain elusive. The early Holocene (11.7 ka to 8.2 ka; Greenlandian Stage) represents the most recent period when the Laurentide and Greenland ice sheets underwent large-scale recession. Moreover, this ice-sheet recession occurred under the backdrop of regional temperatures that were similar to or warmer than today, and comparable to those projected for the upcoming centuries. Reconstructing Laurentide and Greenland ice sheet behavior during the early Holocene, and elucidating the mechanisms dictating this behavior may serve as a partial analog for future Greenland ice-sheet change in a warming world. Here, we present 123 new 10 Be surface exposure ages from two sites on Baffin Island and southwestern Greenland that constrain the behavior of the Laurentide and Greenland ice sheets, and an independent alpine glacier during the early Holocene. On Baffin Island, we focus on a unique area where moraines deposited by the Laurentide Ice Sheet rest directly adjacent to moraines deposited by an independent alpine glacier. Sixty-one 10 Be ages reveal that advances and/or stillstands of the Laurentide Ice Sheet and an alpine glacier occurred in unison around 11.8 ka, 10.3 ka, and 9.2 ka. Sixty-two 10 Be ages from southwestern Greenland indicate that the Greenland Ice Sheet margin experienced re-advances or stillstands around 11.6 ka, 10.4 ka, 9.1 ka, 8.1 ka, and 7.3 ka. Our results reveal that ice sheets respond to climate perturbations on the same centennial timescale as small alpine glaciers. Further, we hypothesize that during the warming climate of the early Holocene, freshening of the North Atlantic Ocean induced by ice-sheet melt resulted in regional cooling and brief periods of ice-sheet stabilization superimposed on net glacier recession. These observations point to a negative feedback mechanism inherent to melting ice sheets in the Baffin Bay region that slows ice-sheet recession during intervals of otherwise rapid deglaciation.

58 GEOSCIENCES↗
DOE OSTI2013

IMAGING STARSPOT EVOLUTION ON KEPLER TARGET KIC 5110407 USING LIGHT-CURVE INVERSION

The Kepler target KIC 5110407, a K-type star, shows strong quasi-periodic light curve fluctuations likely arising from the formation and decay of spots on the stellar surface rotating with a period of 3.4693 days. Using an established light-curve inversion algorithm, we study the evolution of the surface features based on Kepler space telescope light curves over a period of two years (with a gap of .25 years). At virtually all epochs, we detect at least one large spot group on the surface causing a 1%-10% flux modulation in the Kepler passband. By identifying and tracking spot groups over a range of inferred latitudes, we measured the surface differential rotation to be much smaller than that found for the Sun. We also searched for a correlation between the 17 stellar flares that occurred during our observations and the orientation of the dominant surface spot at the time of each flare. No statistically significant correlation was found except perhaps for the very brightest flares, suggesting that most flares are associated with regions devoid of spots or spots too small to be clearly discerned using our reconstruction technique. While we may see hints of long-term changes in the spot characteristics and flare statistics within our current data set, a longer baseline of observation will be needed to detect the existence of a magnetic cycle in KIC 5110407.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Reducing model error using optimized galaxy selection: weak lensing cluster mass estimation

Galaxy clusters are one of the most powerful probes to study extensions of General Relativity and the Standard Cosmological Model. Upcoming surveys like the Vera Rubin Observatory’s Legacy Survey of Space and Time are expected to revolutionise the field, by enabling the analysis of cluster samples of unprecedented size and quality. To reach this era of high-precision cluster cosmology, the mitigation of sources of systematic error is crucial. A particularly important challenge is bias in cluster mass measurements induced by the mismodelling of photometric redshift estimates of source galaxies. This work proposes a method to optimise the source sample selection in cluster weak lensing analyses drawn from wide-field survey lensing catalogs to reduce the bias on reconstructed cluster masses. We use a combinatorial optimisation scheme and methods from variational inference to select galaxies in latent space to produce a probabilistic galaxy source sample catalog for highly accurate cluster mass estimation. We show that our method reduces the critical surface mass density Σ crit relative modelling bias on the 60-70% level, while maintaining up to 90% of galaxies. We highlight that our methodology has applications beyond cluster mass estimation as an approach to jointly combine galaxy selection and model inference under sources of systematics.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1990

Loop current eddy shedding estimated using Geosat altimeter data

A method for extracting periodic signals from altimeter data is applied in the Gulf of Mexico using data from the first two years of the Geosat Exact Repeat Mission. The least squares fitting of selected frequencies of sine and cosine functions to the time series of sea surface height at each point of the Geosat ground track is used to estimate the dominant frequency in the gulf for this time period. Optimal interpolation of the sine and cosine coefficients at the dominant frequency to a uniformly spaced grid allows reconstruction of the time series in the gulf and shows eddies shedding from the Loop Current with a period of approximately 10.5 months, consistent with shedding periods determined from numerical models and analysis of satellite infrared data. This technique is a novel method for investigation of periodic signals in Geosat altimetry and is not severely affected by the data outages present in the time series.

Jacobs, Gregg A.↗
NASA NTRS2006

Processing of Mars Exploration Rover Imagery for Science and Operations Planning

The twin Mars Exploration Rovers (MER) delivered an unprecedented array of image sensors to the Mars surface. These cameras were essential for operations, science, and public engagement. The Multimission Image Processing Laboratory (MIPL) at the Jet Propulsion Laboratory was responsible for the first-order processing of all of the images returned by these cameras. This processing included reconstruction of the original images, systematic and ad hoc generation of a wide variety of products derived from those images, and delivery of the data to a variety of customers, within tight time constraints. A combination of automated and manual processes was developed to meet these requirements, with significant inheritance from prior missions. This paper describes the image products generated by MIPL for MER and the processes used to produce and deliver them.

Instruments↗
NASA NTRS2002

Interference from the Robledo DSN Transmitters to Central Madrid IMT-2000/UMTS System through Terrain Diffraction at S-Band

This study evaluates the possible interference from DSN Robledo 70-m transmitter with Madrid IMT-2000/UMTS wireless users in Spain as both systems will share the same frequency band. Using the effective earth radius, the 50 km terrain profile between Robledo and Madrid is modified and reconstructed. The diffraction propagation losses due to mountain peaks are calculated for the receivers in Madrid urban area. The mountains along the path are simplified into a rounded knife-edge and a rounded obstacle. The results show that for a near surface receiver (1.5 m above the ground) in Madrid, interference signal powers received are less than -135 dBm, which is far below the -109 dBm, the IMT-2000 wireless phone threshold. When a receiver is located at about 40 m above the ground (e.g., the top of Clock Tower of Cibeles Palace), diffraction will generate interference power less than -115 dBm. We find that our calculation results are basically consistent with those from the Longley-Rice model, while the latter has smaller loss because of the low resolution terrain profile used. As a comparison, we also find that the measurements of interference powers of -121.2 dBm at the top of Clock tower is in the range of the estimation. We conclude that the interference through the diffraction mechanism will not cause any problem to IMT-2000/UMTS users at near the surface of Madrid urban area.

telecommunications↗
NASA NTRS2009

Comparison of Node-Centered and Cell-Centered Unstructured Finite-Volume Discretizations: Viscous Fluxes - Part 1

Discretization of the viscous terms in current finite-volume unstructured-grid schemes are compared using node-centered and cell-centered approaches in two dimensions. Accuracy and efficiency are studied for six nominally second-order accurate schemes: a node-centered scheme, cell-centered node-averaging schemes with and without clipping, and cell-centered schemes with unweighted, weighted, and approximately mapped least-square face gradient reconstruction. The grids considered range from structured (regular) grids to irregular grids composed of arbitrary mixtures of triangles and quadrilaterals, including random perturbations of the grid points to bring out the worst possible behavior of the solution. Two classes of tests are considered. The first class of tests involves smooth manufactured solutions on both isotropic and highly anisotropic grids with discontinuous metrics, typical of those encountered in grid adaptation. The second class concerns solutions and grids varying strongly anisotropically over a curved body, typical of those encountered in high-Reynolds number turbulent flow simulations. Results from the first class indicate the face least-square methods, the node-averaging method without clipping, and the node-centered method demonstrate second-order convergence of discretization errors with very similar accuracies per degree of freedom. The second class of tests are more discriminating. The node-centered scheme is always second order with an accuracy and complexity in linearization comparable to the best of the cell-centered schemes. In comparison, the cell-centered node-averaging schemes are less accurate, have a higher complexity in linearization, and can fail to converge to the exact solution when clipping of the node-averaged values is used. The cell-centered schemes using least-square face gradient reconstruction have more compact stencils with a complexity similar to the complexity of the node-centered scheme. For simulations on highly anisotropic curved grids, the least-square methods have to be amended either by introducing a local mapping of the surface anisotropy or modifying the scheme stencil to reflect the direction of strong coupling.

Diskin, Boris↗
NASA NTRS2012

Mars Science Laboratory Navigation Results

The Mars Science Laboratory (MSL), carrying the Curiosity rover to Mars, was launched on November 26, 2011, from Cape Canaveral, Florida. The target for MSL was selected to be Gale Crater, near the equator of Mars, with an arrival date in early August 2012. The two main interplanetary navigation tasks for the mission were to deliver the spacecraft to an entry interface point that would allow the rover to safely reach the landing area, and to tell the spacecraft where it entered the atmosphere of Mars, so it could guide itself accurately to close proximity of the landing target. MSL used entry guidance as it slowed down from the entry speed to a speed low enough to allow for a successful parachute deployment, and this guidance allowed shrinking the landing ellipse to a 99% conservative estimate of 7 by 20 kilometers. Since there is no global positioning system in Mars, achieving this accuracy was predicated on flying a trajectory that closely matched the reference trajectory used to design the guidance algorithm, and on initializing the guidance system with an accurate Mars-relative entry state that could be used as the starting point to integrate the inertial measurement unit data during entry and descent. The pre-launch entry flight path angle (EFPA) delivery requirement was +/- 0.20 deg, but after launch a smaller threshold of +/- 0.05 deg was used as the criteria for late trajectory correction maneuver (TCM) decisions. The pre-launch requirement for entry state knowledge was 2.8 kilometers in position error and 2 meters per second in velocity error, but also smaller thresholds were defined after launch to evaluate entry state update opportunities. The biggest challenge for the navigation team was to accurately predict the trajectory of the spacecraft, so the estimates of the entry conditions could be stable, and late trajectory correction maneuvers or entry parameter updates could be waved off. As a matter of fact, the prediction accuracy was such that the last TCM performed was a small burn executed eight days before landing, and the entry state that was calculated just 36 hours after that TCM, and that was uploaded to the spacecraft the same day, did not need to be updated. The final EFPA was 0.013 deg shallower than the -15.5 deg target, and the on-board entry state was just 200 meters in position and 0.11 meters per second in velocity from the post-landing reconstructed entry state. Overall the entry delivery and knowledge requirements were fulfilled with a margin of more than 90% with respect to the pre-launch thresholds. This excellent accuracy contributed to a very successful and accurate entry, descent, and landing, and surface mission.

navigation↗
NASA NTRS2016

Mars 2020 Entry, Descent and Landing Instrumentation 2 (MEDLI2)

The Mars Entry Descent and Landing Instrumentation 2 (MEDLI2) sensor suite will measure aerodynamic, aerothermodynamic, and TPS performance during the atmospheric entry, descent, and landing phases of the Mars 2020 mission. The key objectives are to reduce design margin and prediction uncertainties for the aerothermal environments and aerodynamic database. For MEDLI2, the sensors are installed on both the heatshield and backshell, and include 7 pressure transducers, 17 thermal plugs, and 3 heat flux sensors (including a radiometer). These sensors will expand the set of measurements collected by the highly successful MEDLI suite, collecting supersonic pressure measurements on the forebody, a pressure measurement on the aftbody, direct heat flux measurements on the aftbody, a radiative heating measurement on the aftbody, and multiple near-surface thermal measurements on the thermal protection system (TPS) materials on both the forebody and aftbody. To meet the science objectives, supersonic pressure transducers and heat flux sensors are currently being developed and their qualification and calibration plans are presented. Finally, the reconstruction targets for data accuracy are presented, along with the planned methodologies for achieving the targets.

Aerodynamics↗
DOE OSTI2024

A multi-model assessment of the early last deglaciation (PMIP4 LDv1): a meltwater perspective

Abstract. The last deglaciation (∼20–11 ka BP) is a period of a major, long-term climate transition from a glacial to interglacial state that features multiple centennial- to decadal-scale abrupt climate variations whose root cause is still not fully understood. To better understand this time period, the Paleoclimate Modelling Intercomparison Project (PMIP) has provided a framework for an internationally coordinated endeavour in simulating the last deglaciation whilst encompassing a broad range of models. Here, we present a multi-model intercomparison of 17 transient simulations of the early part of the last deglaciation (∼20–15 ka BP) from nine different climate models spanning a range of model complexities and uncertain boundary conditions and forcings. The numerous simulations available provide the opportunity to better understand the chain of events and mechanisms of climate changes between 20 and 15 ka BP and our collective ability to simulate them. We conclude that the amount of freshwater forcing and whether it follows the ice sheet reconstruction or induces an inferred Atlantic meridional overturning circulation (AMOC) history, heavily impacts the deglacial climate evolution for each simulation rather than differences in the model physics. The course of the deglaciation is consistent between simulations except when the freshwater forcing is above 0.1 Sv – at least 70 % of the simulations agree that there is warming by 15 ka BP in most places excluding the location of meltwater input. For simulations with freshwater forcings that exceed 0.1 Sv from 18 ka BP, warming is delayed in the North Atlantic and surface air temperature correlations with AMOC strength are much higher. However, we find that the state of the AMOC coming out of the Last Glacial Maximum (LGM) also plays a key role in the AMOC sensitivity to model forcings. In addition, we show that the response of each model to the chosen meltwater scenario depends largely on the sensitivity of the model to the freshwater forcing and other aspects of the experimental design (e.g. CO2 forcing or ice sheet reconstruction). The results provide insight into the ability of our models to simulate the first part of the deglaciation and how choices between uncertain boundary conditions and forcings, with a focus on freshwater fluxes, can impact model outputs. We can use these findings as helpful insight in the design of future simulations of this time period.

Snoll, Brooke (ORCID:0000000230914394)↗
NASA NTRS2014

Neutron Imaging for Selective Laser Melting Inconel Hardware with Internal Passages

Additive Manufacturing is showing great promise for the development of new innovative designs and large potential life cycle cost reduction for the Aerospace Industry. However, more development work is required to move this technology into space flight hardware production. With selective laser melting (SLM), hardware that once consisted of multiple, carefully machined and inspected pieces, joined together can be made in one part. However standard inspection techniques cannot be used to verify that the internal passages are within dimensional tolerances or surface finish requirements. NASA/MSFC traveled to Oak Ridge National Lab's (ORNL) Spallation Neutron Source to perform some non-destructive, proof of concept imaging measurements to assess the capabilities to understand internal dimensional tolerances and internal passages surface roughness. This presentation will describe 1) the goals of this proof of concept testing, 2) the lessons learned when designing and building these Inconel 718 test specimens to minimize beam time, 3) the neutron imaging test setup and test procedure to get the images, 4) the initial results in images, volume and a video, 4) the assessment of using this imaging technique to gather real data for designing internal flow passages in SLM manufacturing aerospace hardware, and lastly 5) how proper cleaning of the internal passages is critically important. In summary, the initial results are very promising and continued development of a technique to assist in SLM development for aerospace components is desired by both NASA and ORNL. A plan forward that benefits both ORNL and NASA will also be presented, based on the promising initial results. The initial images and volume reconstruction showed that clean, clear images of the internal passages geometry are obtainable. These clear images of the internal passages of simple geometries will be compared to the build model to determine any differences. One surprising result was that a new cleaning process was used on these simply geometric specimens that resulted in what appears to be very smooth internal surfaces, when compared to other aerospace hardware cleaning methods.

Tramel, Terri L.↗
NASA NTRS2013

Modeling Insights into Deuterium Excess as an Indicator of Water Vapor Source Conditions

Deuterium excess (d) is interpreted in conventional paleoclimate reconstructions as a tracer of oceanic source region conditions, such as temperature, where precipitation originates. Previous studies have adopted co-isotopic approaches to estimate past changes in both site and oceanic source temperatures for ice core sites using empirical relationships derived from conceptual distillation models, particularly Mixed Cloud Isotopic Models (MCIMs). However, the relationship between d and oceanic surface conditions remains unclear in past contexts. We investigate this climate-isotope relationship for sites in Greenland and Antarctica using multiple simulations of the water isotope-enabled Goddard Institute for Space Studies (GISS) ModelE-R general circulation model and apply a novel suite of model vapor source distribution (VSD) tracers to assess d as a proxy for source temperature variability under a range of climatic conditions. Simulated average source temperatures determined by the VSDs are compared to synthetic source temperature estimates calculated using MCIM equations linking d to source region conditions. We show that although deuterium excess is generally a faithful tracer of source temperatures as estimated by the MCIM approach, large discrepancies in the isotope-climate relationship occur around Greenland during the Last Glacial Maximum simulation, when precipitation seasonality and moisture source regions were notably different from present. This identified sensitivity in d as a source temperature proxy suggests that quantitative climate reconstructions from deuterium excess should be treated with caution for some sites when boundary conditions are significantly different from the present day. Also, the exclusion of the influence of humidity and other evaporative source changes in MCIM regressions may be a limitation of quantifying source temperature fluctuations from deuterium excess in some instances.

simulation↗
DOE OSTI2013

Modifications of Carbonate Fracture Hydrodynamic Properties by CO 2 -Acidified Brine Flow

Acidic reactive flow in fractures is relevant in subsurface activities such as CO{sub 2} geological storage and hydraulic fracturing. Understanding reaction-induced changes in fracture hydrodynamic properties is essential for predicting subsurface flows such as leakage, injectability, and fluid production. In this study, x-ray computed tomography scans of a fractured carbonate caprock were used to create three dimensional reconstructions of the fracture before and after reaction with CO{sub 2}-acidified brine (Ellis et al., 2011, Greenhouse Gases: Sci. Technol., 1:248-260). As expected, mechanical apertures were found to increase substantially, doubling and even tripling in some places. However, the surface geometry evolved in complex ways including ‘comb-tooth’ structures created from preferential dissolution of calcite in transverse sedimentary bands, and the creation of degraded zones, i.e. porous calcite-depleted areas on reacted fracture surfaces. These geometric alterations resulted in increased fracture roughness, as measured by surface Z{sub 2} parameters and fractal dimensions D{sub f}. Computational fluid dynamics (CFD) simulations were conducted to quantify the changes in hydraulic aperture, fracture transmissivity and permeability. The results show that the effective hydraulic apertures are smaller than the mechanical apertures, and the changes in hydraulic apertures are nonlinear. Overestimation of flow rate by a factor of two or more would be introduced if fracture hydrodynamic properties were based on mechanical apertures, or if hydraulic aperture is assumed to change proportionally with mechanical aperture. The differences can be attributed, in part, to the increase in roughness after reaction, and is likely affected by contiguous transverse sedimentary features. Hydraulic apertures estimated by the 1D statistical model and 2D local cubic law (LCL) model are consistently larger than those calculated from the CFD simulations. In addition, a novel ternary segmentation method was devised to handle the degraded zones, allowing for a bounding analysis of the effects on hydraulic properties. We found that the degraded zones account for less than 15% of the fracture volume, but cover 70% to 80% of the fracture surface. When the degraded zones are treated as part of the fracture, the fracture transmissivities are two to four times larger because the fracture surfaces after reaction are not as rough as they would be if one considers the degraded zone as part of the rock. Therefore, while degraded zones created during geochemical reactions may not significantly increase mechanical aperture, this type of feature cannot be ignored and should be treated with prudence when predicting fracture hydrodynamic properties.

roughness↗