Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mesoscale systems”

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.

412 records · Page 23

US East Coast synthetic aperture radar wind atlas for offshore wind energy

We present the first synthetic aperture radar (SAR) offshore wind atlas of the US East Coast from Georgia to the Canadian border. Images from RADARSAT-1, Envisat, and Sentinel-1A/B are processed to wind maps using the geophysical model function (GMF) CMOD5.N. Extensive comparisons with 6008 collocated buoy observations of the wind speed reveal that biases of the individual systems range from -0.8 to 0.6 m s -1 . Unbiased wind retrievals are crucial for producing an accurate wind atlas, and intercalibration of the SAR observations is therefore applied. Wind retrievals from the intercalibrated SAR observations show biases in the range of to -0.2 to 0.0 m s -1 , while at the same time improving the root-mean-squared error from 1.67 to 1.46 m s -1 . The intercalibrated SAR observations are, for the first time, aggregated to create a wind atlas at the height 10 m a.s.l. (above sea level). The SAR wind atlas is used as a reference to study wind resources derived from the Wind Integration National Dataset Toolkit (WTK), which is based on 7 years of modelling output from the Weather Research and Forecasting (WRF) model. Comparisons focus on the spatial variation in wind resources and show that model outputs lead to lower coastal wind speed gradients than those derived from SAR. Areas designated for offshore wind development by the Bureau of Ocean Energy Management are investigated in more detail; the wind resources in terms of the mean wind speed show spatial variations within each designated area between 0.3 and 0.5 m s -1 for SAR and less than 0.2m s -1 for the WTK. Our findings indicate that wind speed gradients and variations might be underestimated in mesoscale model outputs along the US East Coast.

17 WIND ENERGY↗

Mechanism-Resolved PFM of Ferroionic and Ferroelectric Responses in Thickness-Gradient Hf 0.5 Zr 0.5 O 2 Libraries

Resolving growth mechanisms and thickness evolution of functional properties is one of the key tasks in materials discovery and optimization involving thin-film materials, traditionally requiring significant experimental budgets. Here we introduce the combination of thickness-gradient libraries and automated scanning probe microscopy as a systematic pathway to elucidate growth modes and disentangle ferroelectric and electrochemical contributions in ferroelectric thin films. As a model system, we explore the Hf 0.5 Zr 0.5 O 2 (HZO) gradient thin films grown on La x Sr 1-x MnO 3 (LSMO) bottom electrode thin films. Automated piezoresponse force microscopy, spectroscopy, and lithography reveals that irreversible topographic deformation arises from electrochemical activity at the LSMO surface, whereas reversible phase inversion in HZO reflects ferroelectric switching. Automated topography height-map scans are used to further quantify nucleation density, particle-size evolution, and roughness correlations across the thickness-gradient, demonstrating that improved plume stabilization during growth suppresses interfacial reactions and promotes dense, fine-grained HZO conducive to ferroelectric phase formation. This combined materials-engineering and automated-SPM framework establishes a platform for high-throughput, mechanism-resolved characterization of ferroionic and ferroelectric responses in complex oxide films.

FOS: Physical sciences↗

Drizzle, Turbulence, and Density Currents Below Post Cold Frontal Open Cellular Marine Stratocumulus Clouds

Nine cases of postcold frontal marine stratocumulus clouds exhibiting open cellular mesoscale organization are analyzed to characterize the drizzle, turbulence, and density currents below them. Data collected by the vertically pointing Doppler cloud radar and multiple lidars part of the Atmospheric Radiation Measurement (ARM) Eastern North Atlantic (ENA) site are used in these analyses. A total of 251 drizzle shafts passed over the site with 76 shafts sampled for more than 15 min by the vertically pointing instruments. On average the drizzle shafts were ~18 km wide with below-cloud drizzle water path of ~40 g m −2 , cloud base rain rate of 7.78 mm day −1 , and cloud base drizzle modal diameter of 393 μm. The widths of the drizzle shafts did not exhibit any relationship with the cloud base rain rate, below-cloud drizzle water path, and the total water removed from the cloud. The downdrafts in the lowest 500 m within the drizzle cells strengthened with increasing cloud base rain rates, while the updrafts within drizzle cells did not exhibit this behavior. On average, the surface air density, pressure, water vapor mixing ratio, and wind speed along the background wind during the drizzle shaft were higher by 1.88 g m −3 , 8.60 Pa, 0.24 g kg −1 , and 1.22 m s −1 , respectively, from their background values. On the other hand, the surface air temperature was lower by 0.43 K during the drizzle shaft than its background value. The coincident measurements highlight the drizzle-turbulence-surface coupling in these cloud systems.

59 BASIC BIOLOGICAL SCIENCES↗

Microstructure, Thermal, and Mechanical Properties Relationships in U and UZr Alloys (Final Report)

Uranium-zirconium (U-Zr) alloys are candidate fuel systems for transmutation based reactors that can be used to burn long-lived minor actinides and fission products in fast spectrum reactors. Metallic fuels have also been gaining more recent attention for applications as accident tolerant fuels, as well as for use in small modular reactors. This research focused on a “science-based” approach to capture the connections between U and U-Zr alloys’ three-dimensional (3-D) microstructure, thermal properties, and mechanical properties through closely coordinated experiments and modeling efforts from the unirradiated to the irradiated fuels. Advanced characterization and modeling techniques were used to understand irradiation-induced microstructural evolution and its direct impact on the thermal and mechanical properties of U and U-Zr fuel. Closely coordinated experiments and modeling were performed to provide crucial data that does not currently exist. Overall, this research spanned multiple length and time scales within the models and experiments. The scope of the research encompassed the understanding of the irradiation effects in U and various U-Zr alloys with particular attention paid to three task areas: (1) microstructural evolution, (2) in-situ/ex-situ thermal and mechanical properties, and (3) multiscale modeling connections to microstructure, thermal, and mechanical properties. This research resulted in (1) the 3-D characterization of neutron irradiated U-Zr fuel in multiple phase regions to better understand fission gas swelling and constituent redistribution, (2) development of a microstructural model linking thermal and mechanical properties via in situ Raman and nanoindentation, (3) and mesoscale phase field modeling was coupled with the AEH method in the MOOSE framework was used to calculate the effective thermal conductivities of U-Zr fuels consisting of α-U and δ-UZr 2 heterogeneous microstructures.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Molecular To Mesoscale Targeting of Oxoanions with Multi-Tasking Hosts

Achieving a better understanding of anion interactions both in solution and crystalline state was the overarching goal of this project. Anions are everywhere throughout Nature and play important roles in biological and environmental processes. They can be beneficial or deleterious or both in different situations and concentrations. For either reason it is important to have molecules that can bind anions for key needs that benefit society. However, recognition of specific anions is challenging due to the diffuse nature of their negative charge(s) as well as their various shapes and sizes. Understanding the basic properties of anions and how they interact with other molecules and ions in surrounding environments is key to selective recognition. In this project multi-tasking molecules for selective binding of targeted anions were designed to achieve cooperativity and synergism in one rather than multiple host molecules, including (1) cation:anion pair hosts for anions with charges of -2 or greater; (2) pH and redox activated hosts for on-off binding and release; and (3) multiple anion capture in extended host networks. Our design strategy was to combine the use of simple inexpensive building blocks and high yield synthetic pathways to provide economically feasible scale-up for applications. Oxoanions representing multiple shapes and charges were chosen based on having the potential for significant impact on DOE separations needs. Amide/amine-based macrocycles and urea/amine-based chelates and macrocycles with multiple hydrogen bonding sites provided the basic anion-binding frameworks. Successful multi-tasking outcomes were forthcoming in all three tasks. In Task 1, successful ion pair binding for anions with multiple charges was achieved. Furthermore, the ion pair molecules were capable of extended interactions through supramolecular intertwining, like fishing nets for capturing pools of fish (also fitting with Task 3). In Task 2, molecules were synthesized possessing on-off switches. These included a pH sensitive sensor for on-off binding of anions in general, as well as an electrochemical sensor selective for sulfate capture. Three new classes of extended anion host networks capable of binding multiple ions was a major outcome of Task 3. These systems included: anion sensitive, fluorescent organogels; channel-forming macrocycles for studying anion-water including larger macrocyclic cluster sandwiches; and, the offshoot of Task 1, fishing net ion-pair networks for higher valent anions. These strategies can be expanded in the future to other ions and molecules for a better understanding of intermolecular and interionic interactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Disentangling High Harmonic Generation from Surface and Bulk States of a Topological Insulator

The discovery of topological phases has introduced a new dimension to materials science. Three-dimensional (3D) topological insulators (TIs) are a remarkable class of matter that is insulating in the bulk while hosting conductive topological surface states (TSSs) with unique charge and spin properties. High-order harmonic generation (HHG) has emerged as a powerful tool to probe condensed matter systems by providing insights into their electronic structure and dynamic behavior. Here, we investigate HHG in the prototype 3D-TI Bi$_2$Se$_3$. We demonstrate that the contributions of bulk and surface states to the harmonic emission can be controlled by tuning the thickness of thin film samples. An ultrathin (6 nm) film substantially enhances HHG from the surface states, while the bulk states dominate HHG in a thicker (50 nm) film. By applying a quasi-static terahertz perturbing field, we disentangle the bulk and surface responses and reveal the significant impact of the surface states' shift vector and Berry curvature on HHG. Our study provides effective methods for isolating the optical responses of TSSs from those of the bulk, which opens the door to resolving an ongoing debate regarding whether it is possible to reliably extract topological signatures in HHG.

Atomic Physics (physics.atom-ph)↗

Imaging a terahertz superfluid plasmon in a two-dimensional superconductor

The superconducting gap defines the fundamental energy scale for the emergence of dissipationless transport and collective phenomena in a superconductor. In layered high-temperature cuprate superconductors, in which the Cooper pairs are confined to weakly coupled two-dimensional (2D) copper–oxygen (CuO 2 ) planes, terahertz (THz) spectroscopy at subgap millielectronvolt (meV) energies has provided crucial insights into the collective superfluid response perpendicular to the superconducting layers. However, within the CuO 2 planes, the collective superfluid response manifests as plasmonic charge oscillations at energies far exceeding the superconducting gap, obscured by strong dissipation. Here, in this study, we present spectroscopic evidence of a below-gap, 2D superfluid plasmon in few-layer Bi 2 Sr 2 CaCu 2 O 8+x and spatially resolve its deeply subdiffractive THz electrodynamics. By placing the superconductor in the near field of a spintronic THz emitter, we reveal this distinct resonance—absent in bulk samples and observed only in the superconducting phase—and determine its plasmonic nature by mapping the geometric anisotropy and dispersion. Crucially, these measurements offer a direct view of the momentum-dependent and frequency-dependent superconducting transition in two dimensions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Improving the representation of shallow cumulus convection with the simplified-higher-order-closure–mass-flux (SHOC+MF v1.0) approach

Abstract. Parameterized boundary layer turbulence and moist convection remain some of the largest sources of uncertainty in general circulation models. High-resolution climate modeling aims to reduce that uncertainty by explicitly attempting to resolve deep moist convective motions. An example of such a model is the Simple Cloud-Resolving E3SM Atmosphere Model (SCREAM) with a target global resolution of 3.25 km, allowing for a more accurate representation of complex mesoscale deep convective dynamics. Yet, small-scale planetary boundary layer turbulence and shallow convection still need to be parameterized, which in SCREAM is accomplished through the turbulent-kinetic-energy-based (TKE-based) simplified higher-order closure (SHOC) – a simplified version of the assumed-double-Gaussian-PDF (probability density function) higher-order-closure method. In this paper, we implement a stochastic-multiplume mass-flux (MF) parameterization of dry and shallow convection in SCREAM to go beyond the limitations of double-Gaussian-PDF closures and couple it to SHOC (SHOC+MF). The new parameterization implemented in a single-column model type version of SCREAM produces results for two shallow cumulus convection cases (marine and continental shallow convection) that agree well with the reference data from large-eddy simulations, thus improving the general representation of the thermodynamic quantities and their turbulent fluxes as well as cloud macrophysics in the model. Furthermore, SHOC+MF parameterization shows weak sensitivity to the vertical grid resolution and model time step.

54 ENVIRONMENTAL SCIENCES↗

Lake Huron Enhances Snowfall Downwind of Lake Erie: A Modeling Study of the 2010 New Year’s Lake Effect Snowfall Event

In the Northeast and Great Lakes regions of the United States, the influence of multiple lakes on the overlying air can greatly affect lake-effect snowfall over downwind communities. To assess the impact of Lake Huron on snowfall downwind of Lake Erie, we simulated a lake-effect snow event which occurred from 1 to 6 January 2010 using the Regional Atmospheric Modeling System (RAMS). We found that the presence of Lake Huron enhances snowfall downwind of Lake Erie by almost 20% and leads to much heavier local snowfall totals than when Lake Huron is not present. This increase in snowfall is due to a lake-to-lake (L2L) convective band, as secondary circulations associated with lake-effect convection form over Lake Huron and persist overland between the lakes before reintensifying over Lake Erie. As these secondary circulations move over Lake Erie, lifting within the upward branch of the circulations accelerates the development of convection within the L2L band. Surface fluxes and convection over Lake Huron deepen the boundary layer, promoting deeper vertical development of this L2L band over Lake Erie. However, although this boundary layer modification strengthens the L2L band, we found that it actually reduced lake-effect snowfall produced by wind-parallel bands (WPB) in other parts of the lake. This indicates that boundary layer modification from upstream lakes may impact L2L bands differently than WPB which are not connected to upstream secondary circulations.

Lake effects↗

The Amazonian Low-Level Jet and Its Connection to Convective Cloud Propagation and Evolution

In this work, we describe the existence of an Amazonian low-level jet (ALLJ) that can affect the propagation and life cycle of convective systems from the northeast coast of South America into central Amazonia. Horizontal winds from reanalysis were analyzed during March–April–May (MAM) of the two years (2014–15) of the GoAmazon2014/5 field campaign. Convective system tracking was performed using GOES-13 infrared imagery and classified into days with high and weak convective activity. The MAM average winds show a nocturnal enhancement of low-level winds starting near the coast in the early evening and reaching 1600 km inland by late morning. Mean 3-hourly wind speeds maximize at 9–10 m s -1 near 900 hPa, but individual days can have nighttime low-level winds exceeding 12 m s -1 . Based on objective low-level wind criteria, the ALLJ is present 10%–40% of the time over the Amazon during MAM depending on the location and time of day. The evolution of the ALLJ across the Amazon impacts the frequency of occurrence of cloud clusters and the intensity of the moisture flux. In addition, the ALLJ is associated with the enhancement of northeasterly flow in the midtroposphere during active convective days, when vertical momentum transport may be occurring in the organized cloud clusters. During the weakly active convective period, the ALLJ is weaker near the coast but stronger across the central Amazon and appears to be linked more directly with the South American low-level jet.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Structural Changes during the Conversion Reaction of Tungsten Oxide Electrodes with Tailored, Mesoscale Porosity

In-plane tungsten oxide nanostructures, including hexagonally patterned cylinders and holes in a matrix, were fabricated via sequential infiltration synthesis (SIS) on self-assembled block copolymer templates. Using the tailored morphology and porosity of these model electrodes with in situ grazing incidence small angle X-ray scattering, the intrinsic structural change of nanoscale active materials during the conversion reaction of WO 3 + 6Li ↔ W + 3 Li 2 O was investigated at controlled electrochemical conditions. Reversible electrode volume expansion and contraction was observed during lithiation and delithiation cycles, respectively. The potential where electrode’s thickness expansion started was ~1.6 V, which is close to thermodynamically expected one for conversion reaction of WO 3 with lithium (1.65 V). The temporal evolution of the electrode volume at constant electrode potentials revealed high overpotential for bulk lithiation and slow conversion reaction kinetics, despite the tailored porosity of the SIS electrodes. Oxide cylinders showed smaller overall electrode thickness change, likely due to unconstrained lateral volume change, as compared to a matrix with holes. On the other hand, better connectivity and guided volume change of the latter electrode morphology provided improved cycling stability. Additionally, heterogeneity in an electrode, from internal pores and density gradients, was found to aggravate the fragmentation of the electrode during conversion reaction. Insights into oxide conversion reaction kinetics and the relationship between electrode mesostructure and cycling behavior obtained from this study can help guide the more rational design of conversion electrode for high performing batteries.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quantifying the Thermodynamic Impacts on the Atmospheric Boundary Layer due to the Sea Breeze in the Coastal Houston Region

The atmospheric boundary layer (ABL) is unique in coastal regions because of kinematic and thermodynamic influences from continental and marine environments. Sea-breeze (SB) circulations act to equilibrate the land–sea temperature gradient through advecting marine air onshore. The strength of the SB varies in terms of stability, temperature, and moisture advection and influences air quality and weather forecasts. The Tracking Aerosol Convection Interactions Experiment (TRACER) collected a wealth of data on coastal boundary layer evolution, including observations from uncrewed aerial systems (UASs). Vertical profiles of temperature, humidity, and winds were collected by the OU CopterSonde UAS from June to September in the coastal region of Houston. These profiles offer 5-m vertical resolution, on average, every 30 min through diurnal transitions, SB events, and nearby deep convection. During the campaign, CopterSonde observations were gathered through 17 SB events, six of which led to convection initiation. The UAS data can resolve the thermodynamic evolution and interactions between the SB and the preexisting convective boundary layer. Results show large variability across observed SBs and their impacts on temperature and moisture. The intensity of thermodynamic changes depends on the time of sea-breeze passage and influence from the Galveston Bay Breeze, a secondary marine circulation commonly observed in this region. In quantifying the spectrum of SB impacts, equivalent potential temperature θ e is used to contextualize its role in convection initiation and evolution. In conclusion, while all SBs tend to increase θ e from moisture advection, the rate and timing of the θ e rise can distinguish convective from nonconvective cases.

54 ENVIRONMENTAL SCIENCES↗

Assimilation of a Coordinated Fleet of Uncrewed Aircraft System Observations in Complex Terrain: Observing System Experiments

Uncrewed aircraft system (UAS) observations from the Lower Atmospheric Profiling Studies at Elevation–A Remotely-Piloted Aircraft Team Experiment (LAPSE-RATE) field campaign were assimilated into a high-resolution configuration of the Weather Research and Forecasting (WRF) Model. The impact of assimilating targeted UAS observations in addition to surface observations was compared to that obtained when assimilating surface observations alone using observing system experiments (OSEs) for a terrain-driven flow case and a convection initiation (CI) case observed within Colorado’s San Luis Valley (SLV). The assimilation of UAS observations in addition to surface observations results in a clear increase in skill for both flow regimes over that obtained when assimilating surface observations alone. For the terrain-driven flow case, the UAS observations improved the representation of thermal stratification across the northern SLV, which produced stronger upvalley flow over the eastern half of the SLV that better matched the observations. For the CI case, the UAS observations improved the representation of the pre-convective environment by reducing dry biases across the SLV and over the surrounding terrain. This led to earlier CI and more organized convection over the foothills that spilled outflows into the SLV, ultimately helping to increase low-level convergence and CI there. In addition, the importance of UAS capturing an outflow that originated over the Sangre de Cristo Mountains and triggered CI is discussed. These outflows and subsequent CI were not well captured in the simulation that assimilated surface observations alone. We report that observations obtained with a fleet of UAS are shown to notably improve high-resolution analyses and short-term predictions of two very different mesogamma-scale weather events.

54 ENVIRONMENTAL SCIENCES↗

Amplified Mesoscale and Submesoscale Variability and Increased Concentration of Precipitation under Global Warming over Western North America

Abstract Cold-season precipitation statistics in simulations from the storm-resolving WRF Model at 6-km and 1-h resolution over western North America are analyzed. Pseudo–global warming future simulations for the 2041–80 period, constrained by GCMs under the RCP8.5 scenario, are compared to the 1981–2020 historical simulation. The analysis focuses on the dynamical properties of precipitation time series at subdaily scales and on the morphology of storms. The statistical distribution of precipitation intensities in each pixel of the simulation domain is characterized through nonparametric statistical indicators: frequency of wet hours, mean wet-hour precipitation intensity, and Gini coefficient as a measure of the temporal concentration of the precipitation volume. Additionally, the temporal and spatial Fourier power spectra of precipitation time series and precipitation fields are analyzed. The half-power period (HPP) and half-power wavelength (HPW) are defined as spectral measures of the characteristic scales of precipitation’s temporal and spatial patterns. The results show statistically significant increases in the mean wet-hour precipitation intensity and in the Gini coefficient in 99% of the pixels, indicating that the seasonal precipitation volume becomes more concentrated within a smaller number of hours with higher precipitation intensity. The statistics of change in the frequency of wet hours are more contrasted across the simulation domain. The changes are also reflected in the power spectra, which show the spatial and temporal variability increasing proportionally more with finer spatial and temporal scales and the HPW and HPP decreasing. These projected changes are expected to have consequences, not only in terms of hydrologic impacts but also in terms of the predictability of precipitation patterns. Significance Statement The precipitation characteristics of winter storms over the western United States and southwestern Canada are analyzed in future climate simulations for the 2041–80 period. As compared to present-day climate, the most intense parts of the storms are projected to produce a higher rainfall volume, with increased concentration over smaller areas and shorter time intervals. The propensity of rainfall intensity to vary rapidly over time will be enhanced in the future according to the simulations. These model predictions imply an increased risk of rapid flooding in small basins. They also suggest that predicting several hours ahead the time and location at which a storm will produce maximum rainfall may become more challenging in the future.

Climate change↗

Application of an Embedded Fracture and Borehole Modeling Approach to the Understanding of EGS Collab Experiment 1

EGS Collab is a series of meso-scale experiments and associated numerical simulation activities being funded by the United States Department of Energy, Geothermal Technologies Office (GTO) to investigate enhanced geothermal system processes under in-situ stress and slightly elevated temperature conditions. This project is designed to provide scientists and engineers with immediate access to impermeable rock at scales larger than possible in the laboratory, but generally smaller than those for commercial production. Immediate access to rock is provided via the existing drifts of the former Homestake Gold Mine, now operated as the Sanford Underground Research Facility in Lead, South Dakota. The objectives of the EGS Collab project are to develop well controlled fracture networks between injector and producer boreholes using normal and shear stimulation for permeability enhancement. The first experimental site was located off the West Access drift on the 4850 Level (4850 feet below ground surface) in phyllite of the Precambrian Poorman formation, and involved the creation of a fracture network comprising a combination of hydraulic and natural fractures. A second experimental site is now being considered near the battery alcove on the 4100 Level in amphibolite and rhyolite of the Yates Unit. Data generated during these experiments will be compared against predictions of a suite of computer codes specifically designed to solve problems involving coupled thermal, hydrological, geomechanical, and geochemical processes. Comparisons between experimental and numerical simulation results will provide code developers with direction for improvements and verification of process models, build confidence in the suite of available numerical tools, and ultimately identify critical future development needs for the geothermal modeling community. Moreover, conducting thorough comparisons of models, modelling approaches, measurement approaches and measured data, via the EGS Collab project, will serve to identify techniques that are most likely to succeed at the Frontier Observatory for Research in Geothermal Energy (FORGE), the GTO’s flagship EGS research effort. Experiment 1 has comprised a series of successful tests, including a long-term chilled-water circulation test, but the testbed has two atypical EGS elements. Active ventilation in the adjacent drift over a 50-year period cooled the testbed rock mass and hydraulic fracturing intersected monitoring boreholes making them conduits for fluid flow. Numerical simulations executed in support of the design of the EGS Collab Experiment 1, computed a radial temperature distribution and stress gradient orthogonal to the drift axis, resulting in the forecast of an oblong hydraulic fracture geometry, extended in the direction of the drift from the stimulation borehole. This paper describes the numerical simulation of long-term chilled-water circulation test with an embedded fracture and borehole modeling approach. The principal objective of the simulation work is show agreement between experimental observations in terms of production fluid temperatures, monitoring borehole temperatures, injection and production pressures, and tracer recoveries using a single conceptual model for the fracture network, using characterization and monitoring data generally available for EGS. A secondary objective is to improve the numerical simulation result comparisons with additional data available from the broader monitoring equipment within the testbed versus deeper, hotter, and more remote EGS. One complicating factor for comparisons between the numerical simulation results and experimental observations, is the Joule-Thomson heating associated with large pressure drops across the fracture network.

numerical simulation, embedded fractures and boreh↗