SEARCH · Engineering Papers

Results for “Linear Parameter Varying”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 379 records · Page 21

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Sensitivity study of hydrogen Soret transport in yttrium Hydride-Based nuclear fuel

Yttrium hydride is an excellent solid neutron moderator material for high temperature nuclear reactor applications due to its high hydrogen density and exceptional hydride stability at high temperatures. Despite these attractive characteristics, the details of how hydrogen behaves within yttrium hydride while temperature gradients exist are still not well understood. The evolution of the hydrogen composition profile resulting from a temperature gradient requires knowledge of hydrogen’s heat of transport, a critical parameter that has not yet been measured for this material. In this work, we perform hydride redistribution, hydrogen dissociation, and hydrogen leakage calculations while varying the Soret heat of transport of hydrogen in yttrium hydride to elucidate the sensitivity of hydride stability under temperature gradients to this parameter. This study analyzes hydride stability of a hypothetical uranium-yttrium hydride nuclear fuel design during operation of a high temperature liquid metal-cooled nuclear reactor. Assuming U-YH x could be fabricated in a physically stabilized manner, this fuel system can likely maintain hydride stability while operating at very high power densities and temperatures. We find that even though the hydrogen dissociation pressure in the gas gap does vary by several percent as the heat of transport temperature parameter is varied, the hydrogen content in the U-YH x fuel meat is relatively insensitive to this parameter over the course of a high burnup fuel cycle; this is due to yttrium hydride’s excellent hydrogen retention under the high temperature conditions considered here. Here, this suggests that hydride stability analyses are insensitive to the value of the Soret heat of transport in U-YH x under steady state liquid metal-cooled reactor conditions. However, the susceptibility to internal gas overpressurization-induced stress-rupture of the cladding during a high temperature transient is more sensitive to this parameter due to the non-linear dependence of hydrogen gas dissociation pressure vs. composition and temperature.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2026

Numerical Analysis of Offshore Wind-Farm-Induced Drag Effects on Coastal Upwelling Dynamics

Wind farms extract momentum from the atmospheric flow, generating wind-speed deficits both within the plant, and extending downstream. When located offshore, these deficits modulate air-sea coupling, potentially impacting coastal upwelling in sensitive regions. We investigate impacts of wind farms on coastal upwelling using kilometre-scale, three-way-coupled simulations with the coupled ocean-atmosphere-wave-sediment transport system for the US West Coast. Wind-farm effects are represented by a generalised turbine drag formulation, an idealised, height-dependent body force whose magnitude is systematically varied. This approach isolates the leading-order fluid-dynamical response in a realistic coastal configuration. The atmospheric adjustment exhibits an approximately linear relation between drag force and wind-speed deficit, with wakes that expand downstream and increase in magnitude as drag increases. An empirical orthogonal function analysis of sea-surface-temperature anomalies reveals the emergence of a canonical dipole pattern under strong drag forcing. Subsurface diagnostics show consistent shoaling of the mixed layer and suppressed upward velocities in areas near wind-farm region, accompanied by compensating enhancements of shoaling closer to the coast. These results identify turbine drag as a control parameter in assessing interactions between wind-farm wake and coastal upwelling and provide scaling relationships for understanding offshore wind-farm effects on the coastal circulation dynamics.

16 TIDAL AND WAVE POWER↗

DOE OSTI2023

The role of isotope mass and transport for H-mode access in tritium containing plasmas at JET with ITER-like wall

The required heating power, $P_{\mathrm{LH}}$, to access the high confinement regime (H-mode) in tritium containing plasmas is investigated in JET with ITER-like wall at a toroidal magnetic field of $B_t = 1.8$ T and a plasma current of $I_p = 1.7$ MA. $P_{\mathrm{LH}}$, also referred to as the L-H power threshold, is determined in plasmas of pure tritium as well as mixtures of hydrogen with tritium (H-T) and mixtures of deuterium with tritium (D-T), and is compared to the L-H power threshold in plasmas of pure hydrogen and pure deuterium. It is found that, for otherwise constant parameters, $P_{\mathrm{LH}}$ is not the same in plasmas with the same effective isotope mass, $A_{\mathrm{eff}}$, when they differ in their isotope composition. Thus, $A_{\mathrm{eff}}$ is not sufficient to describe the isotope effect of $P_{\mathrm{LH}}$ in a consistent manner for all considered isotopes and isotope mixtures. The electron temperature profiles measured at the L-H transition in the outer half of the radius are very similar for all isotopes and isotope mixtures, despite the fact that the L-H power threshold varies by a factor of about six. This finding, together with the observation of an offset linear relation between the L-H power threshold, $P_{\mathrm{LH}}$, and an effective heat diffusivity, $\chi_{\mathrm{eff}}$, indicates that the composition-dependent heat transport in the low confinement mode (L-mode) determines, how much power is needed to reach the necessary electron temperatures at the edge, and hence $P_{\mathrm{LH}}$.

tritium plasmas↗

NASA NTRS1996

Steering of Frequency Standards by the Use of Linear Quadratic Gaussian Control Theory

Linear quadratic Gaussian control is a technique that uses Kalman filtering to estimate a state vector used for input into a control calculation. A control correction is calculated by minimizing a quadratic cost function that is dependent on both the state vector and the control amount. Different penalties, chosen by the designer, are assessed by the controller as the state vector and control amount vary from given optimal values. With this feature controllers can be designed to force the phase and frequency differences between two standards to zero either more or less aggressively depending on the application. Data will be used to show how using different parameters in the cost function analysis affects the steering and the stability of the frequency standards.

Koppang, Paul↗

DOE OSTI2024

Linear Correlations of Gibbs Free Energy of REE Phosphates (Monazite, Xenotime, and Rhabdophane) and Internally Consistent Binary Mixing Properties

Rare Earth Elements (REE) phosphates (monazite, xenotime, and rhabdophane) are critical REE-bearing minerals typically formed in hydrothermal and magmatic ore deposits. The thermodynamic properties of those REE minerals are crucial to understanding the solubility, speciation, and transport of REE complexes. However, the standard-state Gibbs free energy of formation (∆G°f) values reported for these minerals in the literature vary by up to 25 kJ mol−1. Here, we present linear free energy relationships that allow the evaluation and estimation of the ∆G°f values at 25 °C and 1 bar for the three minerals from the ionic radius (rREE3+) and the non-solvation Gibbs free energy contribution to the REE3+ aqua ion (∆G°n, REE3+): ∆G°f,monazite − 399.71 rREE3+ = 1.0059 ∆G°n,REE3+ − 2522.51; ∆G°f,xenotime − 344.08 rREE3+ = 0.9909 ∆G°n,REE3+ − 2451.53; and ∆G°f,rhabdophane − 416.17 rREE3+ = 1.0067 ∆G°n, REE3+ − 2688.86. Moreover, based on the new dataset derived for REE end-members, we re-fitted the binary Margules parameter (W) from previous theoretical calculations into linear correlations: W + 0.00204 ∆G°’n,monazite = 39.3549 ∆V + 0.0641; W + 0.00255 ∆G°’n,xenotime = 25.4885 ∆V − 0.0062. The internally consistent thermodynamic properties of these REE phosphates are incorporated into the computer program Supcrtbl, which is available online at Zhu’s research website.

Pan, Ruiguang (ORCID:0000000324896545)↗

NASA NTRS1976

A study of the effect of subcritical crack growth on the geometry dependence on nonlinear fracture toughness parameters

The effect of subcritical crack growth on the geometry dependence of nonlinear fracture toughness parameters was studied by comparing the toughness values for different specimen geometries at the onset of subcritical crack growth and at the initiation of unstable crack propagation. Center-cracked thin sheet specimens of 2024-T3 and 7075-T6 aluminum alloys were tested by varying the specimen length L, width w, and crack length-to-width ratio c/w. When the onset of unstable crack propagation was selected as the critical point, the nonlinear energy toughness and the R curve toughness increased with increasing w and decreasing L and c/w. However, when the onset of subcritical crack growth was taken as the critical point, energy toughness and the linear toughness values were independent of these geometrical variables.

Jones, D. L.↗

DOE OSTI2023

Empirical radius formulas for canonical neutron stars from bidirectionally selecting features of equations of state in extended Bayesian analyses of observational data

Significant advancement in Bayesian inference of nuclear equation of state (EOS) from gravitational wave and x-ray observations of neutron stars (NSs) has been made by the nuclear astrophysics community especially since GW170817. By extending the traditional Bayesian analysis which normally ends at presenting the marginalized posterior probability distribution functions (PDFs) of individual EOS parameters and their correlations (or sometimes only the Pearson correlation coefficients which are only reliably useful when the variables are linearly correlated while they are actually often not), we search for a data-driven and robust empirical formula for the radius 𝑅 1.4 of canonical NSs in terms of the characteristic EOS parameters (features). We also identify the single most important but currently poorly known EOS parameter for determining the 𝑅 1.4 . Using three regression-model-building methodologies: bidirectional stepwise feature selection, least absolute shrinkage selection operator (LASSO) regression, and neural network regression on a large set of posterior EOSs and the corresponding 𝑅 1.4 values inferred from earlier comprehensive Bayesian analyses of NS observational data, we systematically and rigorously develop the most probable 𝑅 1.4 formulas with varying statistical accuracy and technical complexity. Here, the most important EOS parameters for determining 𝑅 1.4 are found consistently in each of the feature selection processes to be (in order of decreasing importance): curvature 𝐾 sym , slope 𝐿, skewness 𝐽 sym of nuclear symmetry energy, skewness 𝐽 0 , incompressibility 𝐾 0 of symmetric nuclear matter, and the magnitude 𝐸 sym ⁡(𝜌 0 ) of symmetry energy at the saturation density 𝜌 0 of nuclear matter.

Bayesian methods↗

NASA NTRS2023

Solar Activity Relations in Energetic Electron Events Measured By the MESSENGER Mission

Aims. We perform a statistical study of the relations between the properties of solar energetic electron (SEE) events measured by the MESSENGER mission from 2010 to 2015 and the parameters of the respective parent solar activity phenomena in order to identify the potential correlations between them. During the time of analysis, the MESSENGER heliocentric distance varied between 0.31 and 0.47 au. Methods. We used a published list of 61 SEE events measured by MESSENGER, which includes information on the near-relativistic electron peak intensities, the peak-intensity energy spectral indices, and the measured X-ray peak intensity of the flares related to the SEE events. Taking advantage of multi-viewpoint remote-sensing observations, we reconstructed, whenever possible, the associated coronal mass ejections (CMEs) and shock waves; and we determined the three-dimensional (3D) properties (location, speed, and width) of the CMEs and the maximum speed of the 3D CME-driven shocks in the corona. We used different methods (Spearman, Pearson, and a Bayesian approach, namely the Kelly method to linear regression) to estimate the correlation coefficients between the flare intensity, maximum speed at the apex of the CME-driven shock, CME speed at the apex, and CME width with the electron peak intensities and with the energy spectral indices. In this statistical study, we considered and addressed the limitations of the particle instrument on board MESSENGER (elevated background intensity level, anti-Sun pointing). Results. There is an asymmetry to the east in the range of connection angles (CAs) for which the SEE events present the highest peak intensities, where the CA is the longitudinal separation between the footpoint of the magnetic field connecting to the spacecraft and the flare location. Based on this asymmetry, we define a subsample of well-connected events as when −65° ≤ CA ≤ +33°. For the well-connected sample, we find moderate to strong correlations between the near-relativistic electron peak intensity and the 3D CME-driven shock maximum speed at the apex (Spearman: cc = 0.53 ± 0.05; Pearson: cc = 0.65 ± 0.04; Kelly: cc = 0.87 ± 0.20), the flare peak intensity (Spearman: cc = 0.63 ± 0.03; Pearson: cc = 0.59 ± 0.03; Kelly: cc = 0.74 ± 0.30), and the 3D CME speed at the apex (Spearman: cc = 0.50 ± 0.04; Pearson: cc = 0.46 ± 0.03; Kelly: cc = 0.60 ± 0.39). When including poorly connected events (full sample), the relations between the peak intensities and the solar-activity phenomena are blurred, showing lower correlation coefficients. Conclusions. Based on the comparison of the correlation coefficients presented in this study using near 0.4 au data, (1) both flare and shock-related processes may contribute to the acceleration of near relativistic electrons in large SEE events, in agreement with previous studies based on near 1 au data; and (2) the maximum speed of the CME-driven shock is a better parameter to investigate particle-acceleration-related mechanisms than the average CME speed, as suggested by the stronger correlation with the SEE peak intensities.

Sun↗

NASA NTRS1996

High Reynolds number analysis of flat plate and separated afterbody flow using non-linear turbulence models

The ability of the three-dimensional Navier-Stokes method, PAB3D, to simulate the effect of Reynolds number variation using non-linear explicit algebraic Reynolds stress turbulence modeling was assessed. Subsonic flat plate boundary-layer flow parameters such as normalized velocity distributions, local and average skin friction, and shape factor were compared with DNS calculations and classical theory at various local Reynolds numbers up to 180 million. Additionally, surface pressure coefficient distributions and integrated drag predictions on an axisymmetric nozzle afterbody were compared with experimental data from 10 to 130 million Reynolds number. The high Reynolds data was obtained from the NASA Langley 0.3m Transonic Cryogenic Tunnel. There was generally good agreement of surface static pressure coefficients between the CFD and measurement. The change in pressure coefficient distributions with varying Reynolds number was similar to the experimental data trends, though slightly over-predicting the effect. The computational sensitivity of viscous modeling and turbulence modeling are shown. Integrated afterbody pressure drag was typically slightly lower than the experimental data. The change in afterbody pressure drag with Reynolds number was small both experimentally and computationally, even though the shape of the distribution was somewhat modified with Reynolds number.

Carlson, John R.↗

DOE OSTI2021

Estimating and Forecasting Time-Varying Groundwater Recharge in Fractured Rock: A State-Space Formulation with Preferential and Diffuse Flow to the Water Table

Rapid infiltration following precipitation may result in groundwater contamination from surface contaminants or pathogens. In fractured rock, contaminants can migrate rapidly to points of groundwater withdrawals. In contrast to the temporal availability of groundwater quality chemical indicators, meteorological and groundwater level observations are available in real-time to estimate time-varying recharge, which can act as a surrogate to identify periods of rapid infiltration that may indicate contamination susceptibility. Estimating recharge using methods, such as base-flow recession, unsaturated infiltration models, or Water-Table Fluctuations (WTF), cannot capitalize on currently available technologies and telecommunication infrastructure to conduct real-time recharge estimation at scales relevant to characterizing rapid infiltration. We present a linear, physics-based State-Space (SS) model of one-dimensional infiltration to estimate recharge, which includes preferential and diffuse-flow to the water table. The model can take advantage of real-time data for water-table altitude, precipitation, and evapotranspiration. Model parameters are calibrated over an observation period, and the Kalman Filter (KF) is subsequently applied to continuously update the observed (water-table altitude) and unobserved (groundwater recharge) system states and predict future states as new data become available. The SS/KF algorithm is demonstrated at the Masser Groundwater Recharge Site in Pennsylvania, USA and comparisons are made with recharge estimates from WTF methods. Model results indicate that the frequency of observations (daily versus sub-daily) dictates the allocation between preferential and diffuse flow. Additionally, because infiltration processes encompass many nonlinearities, model parameters estimated from observation periods need to be updated at least seasonally to account for changing recharge conditions.

groundwater, Recharge, infiltration, recursive est↗

NASA NTRS1984

Characterization of Surface Reflectance Variation Effects on Remote Sensing

The use of Monte Carlo radiative transfer codes to simulate the effects on remote sensing in visible and infrared wavelengths of variables which affect classification is examined. These variables include detector viewing angle, atmospheric aerosol size distribution, aerosol vertical and horizontal distribution (e.g., finite clouds), the form of the bidirectional ground reflectance function, and horizontal variability of reflectance type and reflectivity (albedo). These simulations are used to characterize the sensitivity of observables (intensity and polarization) to variations in the underlying physical parameters both to improve algorithms for the removal of atmospheric effects and to identify techniques which can improve classification accuracy. It was necessary to revise and validate the simulation codes (CTRANS, ARTRAN, and the Mie scattering code) to improve efficiency and accommodate a new operational environment, and to build the basic software tools for acquisition and off-line manipulation of simulation results. Initial calculations compare cases in which increasing amounts of aerosol are shifted into the stratosphere, maintaining a constant optical depth. In the case of moderate aerosol optical depth, the effect on the spread function is to scale it linearly as would be expected from a single scattering model. Varying the viewing angle appears to provide the same qualitative effect as modifying the vertical optical depth (for Lambertian ground reflectance).

Pearce, W. A.↗

NASA NTRS2000

Advances in Fluorescence Sensing Systems for the Remote Assessment of Nitrogen Supply in Field Corn

The studies described herein were conducted to better define changes in fluorescence properties of leaves from field grown corn (Zea mays L.) as they relate to varying levels of nitrogen (N) fertilization. This research was directed toward: 1) providing a remote non-destructive sensing technique to aid in the determination of optimal rates of N fertilization in corn crops and, 2) defining parameters for further development of fluorescence instrumentation to be operated remotely at field canopy levels. Fluorescence imaging bands centered in the blue (450 nm), green (525 nm), red (680 nm), and far-red (740 nm) and ratios of these bands were compared with the following plant parameters: rates of photosynthesis, N:C ratio, pigment concentrations, and grain yields. Both the fluorescence and physiological measures exhibited similar curvilinear responses to N fertilization level while significant linear correlations were obtained among fluorescence bands and band ratios to certain physiological measures of plant productivity. The red / blue, red / green, far-red / blue, far-red /green fluorescence ratios are well suited for remote observation and provided high correlations to grain yield, LAI, N:C, and chlorophyll contents. The results from this investigation indicate that fluorescence technology could aid in the determination of N fertilization requirements for corn. This discussion will also address design concepts and preliminary field trials of a mobile field-based Laser Induced Fluorescence Imaging System (LIFIS) capable of simultaneously acquiring images of four fluorescence emission bands from areas of plant canopies equaling 1 sq m and greater without interference of ambient solar radiation.

Corp, L. A.↗

NASA NTRS2006

Reliability of Radioisotope Stirling Convertor Linear Alternator

Onboard radioisotope power systems being developed and planned for NASA s deep-space missions would require reliable design lifetimes of up to 14 years. Critical components and materials of Stirling convertors have been undergoing extensive testing and evaluation in support of a reliable performance for the specified life span. Of significant importance to the successful development of the Stirling convertor is the design of a lightweight and highly efficient linear alternator. Alternator performance could vary due to small deviations in the permanent magnet properties, operating temperature, and component geometries. Durability prediction and reliability of the alternator may be affected by these deviations from nominal design conditions. Therefore, it is important to evaluate the effect of these uncertainties in predicting the reliability of the linear alternator performance. This paper presents a study in which a reliability-based methodology is used to assess alternator performance. The response surface characterizing the induced open-circuit voltage performance is constructed using 3-D finite element magnetic analysis. Fast probability integration method is used to determine the probability of the desired performance and its sensitivity to the alternator design parameters.

Shah, Ashwin↗

NASA NTRS2008

Effect of Microstructure on Time Dependent Fatigue Crack Growth Behavior In a P/M Turbine Disk Alloy

A study was conducted to determine the processes which govern hold time crack growth behavior in the LSHR disk P/M superalloy. Nineteen different heat treatments of this alloy were evaluated by systematically controlling the cooling rate from the supersolvus solutioning step and applying various single and double step aging treatments. The resulting hold time crack growth rates varied by more than two orders of magnitude. It was shown that the associated stress relaxation behavior for these heat treatments was closely correlated with the crack growth behavior. As stress relaxation increased, the hold time crack growth resistance was also increased. The size of the tertiary gamma' in the general microstructure was found to be the key microstructural variable controlling both the hold time crack growth behavior and stress relaxation. No relationship between the presence of grain boundary M23C6 carbides and hold time crack growth was identified which further brings into question the importance of the grain boundary phases in determining hold time crack growth behavior. The linear elastic fracture mechanics parameter, Kmax, is unable to account for visco-plastic redistribution of the crack tip stress field during hold times and thus is inadequate for correlating time dependent crack growth data. A novel methodology was developed which captures the intrinsic crack driving force and was able to collapse hold time crack growth data onto a single curve.

Telesman, Ignacy J.↗

DOE OSTI2025

WEST L-mode record long pulses guided by predictions using Integrated Modeling

A new record was set on the WEST Tokamak, designed to operate long duration plasmas in a tungsten (W) environment, with an injected energy of 1.15 GJ and a plasma duration 364s. Scenario development was supported by integrated modeling using the High Fidelity Plasma Simulator (HFPS), the European IMAS-coupled version of JETTO/JINTRAC, which integrates physics-driven modules into a unified framework. In particular, a reduced model for Lower-Hybrid heating and Current-Drive (LHCD) and the quasi-linear turbulent transport model TGLF are crucial for long pulses predictions up to the Last Closed Flux Surface (LCFS). Using this workflow, a 100 s reference discharge was modeled and plasma kinetic profiles and loop voltage were quantitatively well matched. In preparation for the recent long duration experiments, non-inductive current-drive actuators (I P , n e , P LHCD ) were varied to determine the operational domain going towards fully non-inductive discharges. In particular, decreasing the plasma current is shown to ease the access to such conditions, with a careful monitoring of (n e , P LHCD ) to avoid machine limitations. In addition, post-prediction experiments conducted within the investigated parameter range validated the predicted dependencies and were shown to be in quantitative agreement. Exploratory work on the use of ECCD for MHD stability purpose is also introduced.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Transport in Stochastic Media with Random Chord Length Distributions

Thermal radiation transport computations in binary Markovian random mixtures rely almost exclusively on the Levermore-Pomraning (LP) model which is obtained by applying a heuristic closure to the ensemble averaged random medium transport equation. The validity of this model has been extensively tested by comparing numerical results over a broad parameter range (material types and mixing parameters) against benchmark solutions in planar geometry. The conditions under which the LP-model provides useful results and when it breaks down are now well established, but work to date has been largely restricted to homogeneous mixing statistics, i.e., the mean chord lengths of both materials are taken to be spatially constant. In recent work, this limitation was relaxed by allowing the mean chord lengths and, in a consistent fashion, the volume fractions in the LP-model to vary spatially and in a follow-up investigation benchmark solutions were obtained by ensemble averaging results over material realizations sampled from a nonhomogeneous Poisson process (NHPP). Numerical experiments in rod geometry with specifically linear and quadratic spatial dependence of chord lengths showed that the material averaged radiation intensities vary nonmonotonically with depth into the medium, in stark contrast to solutions obtained assuming uniform chord lengths. Moreover, depending on the local optical thickness and strength of scattering, the LP-model results showed locally more nuanced deviations from the benchmark solutions than was the case with constant chord lengths. These limited numerical investigations highlight the nontrivial qualitative and quantitative consequences of nonhomogeneous mixing statistics, in particular that closure approximations may not be uniformly valid or invalid over the problem domain.

42 ENGINEERING↗

NASA NTRS1995

Wavenumber selection and hysteresis in nonlinear baroclinic flow

Wavenumber transition and hysteresis in a highly unstable baroclinic flow are investigated using a high-resolution spectral numerical model. As the flow becomes more supercritical, the dominant wave gradually shifts from the most unstable wave predicted by the linear theory to a longer wave with a larger time-averaged amplitude, while the rectified mean flow attains a stronger shear at the center of the channel. The numerical results display a complex hysteresis behavior, which occurs not only between the states of different dominant wavenumbers, but also between the states of identical dominant wavenumber but of different dynamic characteristics. In a certain parameter range three stable states, each with different dominant wavenumber, are possible, and in another parameter range four stable states are possible, among them three stable states with an identical dominant wave. The numerical results suggest that a multiple weather regime exists even without external forcing in which the flow aperiodically varies between two distinct behaviors. The effects of stable higher harmonics are assessed and it is found that their presence contributes not only to the better approximation of the model solutions but also to the selection of the final equilibrium state, due to the chaotic nature of the initial transient period.

Chou, Shih-Hung↗

DOE OSTI2021

Verification of Coupled Hydraulic Fracturing Simulators Using Laboratory-Scale Experiments

In this work, we aim to verify the predictions of the numerical simulators, which are used for designing field-scale hydraulic stimulation experiments. Although a strong theoretical understanding of this process has been gained over the past few decades, numerical predictions of fracture propagation in low-permeability rocks still remains a challenge. Against this background, we performed controlled laboratory-scale hydraulic fracturing experiments in granite samples, which not only provides high-quality experimental data but also a well-characterized experimental set-up. Using the experimental pressure responses and the final fracture sizes as benchmark, we compared the numerical predictions of two coupled hydraulic fracturing simulators—CSMP and GEOS. Both the simulators reproduced the experimental pressure behavior by implementing the physics of Linear Elastic Fracture Mechanics (LEFM) and lubrication theory within a reasonable degree of accuracy. The simulation results indicate that even in the very low-porosity (1–2 %) and low-permeability (10 -18 m 2 - 10 -19 m 2 ) crystalline rocks, which are usually the target of EGS, fluid-loss into the matrix and unsaturated flow impacts the formation breakdown pressure and the post-breakdown pressure trends. Therefore, underestimation of such parameters in numerical modeling can lead to significant underestimation of breakdown pressure. The simulation results also indicate the importance of implementing wellbore solvers for considering the effect of system compressibility and pressure drop due to friction in the injection line. The varying injection rate as a result of decompression at the instant of fracture initiation affects the fracture size, while the entry friction at the connection between the well and the initial notch may cause an increase in the measured breakdown pressure.

15 GEOTHERMAL ENERGY↗