Saturated Liquid Research Articles

Investigation of film boiling over tandem arrangement of elliptical cylinders in vertical cross-flow of saturated liquid

The present study investigates the phenomenon of film boiling heat transfer for the tandem arrangement of elliptical heaters placed in cross-flow of saturated liquid (water) at near critical pressure conditions (P/Pc = 0.99), using the Eulerian approach-based volume of the fluid framework. The effect of vertical spacing between the heaters (pd = 7.5a to 12.5a), their aspect ratio (ψ = 0.4 to 1), degree of wall superheat (Jav/Prv = 0.6 to 1.8) and Reynolds number for liquid flow (Rel = 10 to 40) is investigated on the bubble departure dynamics and heat transfer characteristics. It is observed that the bubble departure size from the bottom heater has remained almost unchanged for all the situations whereas its departure frequency has increased upon an increase of Reynolds number and degree of superheat. The bubble size has varied significantly for the top elliptical heater and its departure frequency has increased with the decrease of aspect ratio and increase of pitch, Reynolds number and degree of superheat. The time-averaged heat transfer coefficient has increased upon an increase in pitch distance, aspect ratio, Reynolds number and degree of superheat, however, the rate of increase has diminished at a high degree of superheat.

Enhanced internal condensation of R1233zd(E) on micro- and nanostructured copper and aluminum surfaces

In-tube condensation of refrigerants is an important process which affects thermal efficiency in many applications, ranging from refrigeration and air conditioning to electronics thermal management. In-tube heat transfer and pressure drop are important to heat exchanger sizing and design. In this work, micro- and nanostructured surfaces are applied to the internal wetted areas of copper and aluminum mini-channels to enhance the condensation heat transfer coefficient of hydrofluorocarbon R1233zd(E) refrigerant. To achieve scalable nanomanufacturing, surfaces were uniformly structured by relying on hydrochloric acid etching of aluminum and chemical oxidation of copper. The etched aluminum surfaces exhibited a 150% increase in heat transfer coefficient compared to smooth aluminum channels at specific qualities, with a 66% heat transfer coefficient improvement for complete phase change from saturated vapor to saturated liquid. Copper oxide structures showed no discernable difference in thermal-hydraulic performance when compared to smooth copper channels. Critical dimensionless parameters governing the heat transfer enhancement were identified by varying the tube internal diameter (2.3 mm to 4.7 mm), refrigerant mass flux (50 to 300 kg/(m2·s)), and refrigerant quality (0 to 1). The dimensionless parameters include the Bond number normalized to the condensate film thickness, and the Weber number modified by the vapor friction factor. The relatively small increase in pressure drop (< 10%) associated with these surface enhancements further supports the promise of this method. The scalable and cost-effective techniques used to create these aluminum microstructures may reduce manufacturing cost when compared with current enhancement approaches such as extrusion, drawing, and welding.

A Reference Equation of State with an Associating Term for the Thermodynamic Properties of Ammonia

Fluid ammonia is highly associated because of strong intermolecular hydrogen bonding. This causes different behavior of thermophysical properties from nonpolar fluids, as reflected, for example, in the heat capacity. In this work, a new equation of state explicit in the Helmholtz energy with independent variables of temperature and density was developed, which includes a new associating term. Its uncertainties in density, vapor pressure, saturated liquid and vapor densities, and caloric properties were estimated by comparisons with experimental data. The new equation of state is valid from the triple-point temperature (195.49 K) to 725 K at pressures up to 1000 MPa and densities up to 53.13 mol dm−3. Physically correct behavior within the region of validity and at extremely high temperatures and pressures, and at temperatures far below the triple-point temperature, was obtained by applying constraints on various properties. The unique physical behavior of ammonia shown in some thermodynamic properties is described in detail, which will provide a preliminary template for developing equations of state for other associating fluids.

Entropy, Heat Capacity, and Thermal Conductivity of trans-1,1,1,4,4,4-Hexafluoro-2-butene [R1336mzz(E)] Derived from the Gaseous Sound Speed Measurement

Thermophysical properties such as entropy, heat capacity, and thermal conductivity of trans-1,1,1,4,4,4-hexafluoro-2-butene, which is designated as R1336mzz(E) for refrigerant use, were derived from the sound speed measurement in the gaseous phase. The gaseous sound speed was measured by using an acoustic-microwave resonator with a cylindrical cavity, which enables simultaneous measurement of the gaseous relative dielectric permittivity. The ideal gas heat capacity was determined from the gaseous sound speed by analyzing the measurement data with the acoustic-virial equation, and the electric dipole moment was additionally determined from the gaseous relative permittivity measurements. In compliance with the fluid-property correlation models that involve parameters associated with the ideal gas heat capacity and the dipole moment under the corresponding states principle, entropies at saturated vapor and liquid, isobaric heat capacities at saturated liquid, and thermal conductivities at vapor and liquid phases were evaluated for R1336mzz(E).

Euler–Euler Multi-Scale Simulations of Internal Boiling Flow with Conjugated Heat Transfer

A numerical approach was implemented, to study a boiling flow in a horizontal serpentine tube. A NEPTUNE_CFD two-fluid model was used, to study the behavior of the refrigerant R141b in diabatic cases. The model was based on the Euler–Euler formalism of the Navier–Stokes equations, in which governing equations are solved for both phases of the fluid at each time step. The conjugate heat transfer—between the tube wall and the fluid—was considered via a coupling with the SYRTHES 4.3 software, which solves solid conduction in three dimensions. A mesh convergence study was carried out, which found that a resolution of 40 meshes per diameter was necessary for our case. The approach was validated by comparison with an experimental study of the literature, based on the faithful reproduction of the positions of two-phase flow regime transitions in the domain. Original post-processing was used, to unravel the flow characteristics. The mean and RMS fields of void fraction, temperatures and stream wise velocities in several sections were analyzed, when statistical convergence was reached. A thermal equilibrium was reached in the saturated liquid, but not in the vapor phase, due to the flow dynamic and possibly the presence of droplets. Finally, a thermal analysis of the configuration was proposed. It demonstrated the strong coupling between the temperature distribution in the solid, and the two-phase flow regimes at stake in the fluid domain.

THERMODYNAMIC PROPERTIES CALCULATION USING A MULTIOBJECTIVE APPROACH IN ESTIMATING THE PARAMETERS OF THE GEOS3C EQUATION OF STATE

Since van der Waals' first publication, a lot of equations of state have been proposed to represent the PVT behavior of pure compounds, as is the case with GEOS3C, which is a new form of generalized cubic equation of state that uses a temperature function dependent on three adjustable empirical parameters. In order to obtain the parameters that lead to the lowest errors of vapor pressure and saturated liquid volume predictions, it is possible to turn this problem into a multiobjective optimization problem. In this context, a modified MOPSO algorithm was employed. This method has the main advantage to provide a set of solutions that show the existing conflict in the attempt to minimize the properties of interest. Substances from different chemical families were used to evaluate the methodology proposed in this work and comparisons against experimental data were made. Results showed that the values obtained through multiobjective optimization for GEOS3C adjustable parameters provide better or equivalent predictions for all calculated thermodynamic properties. In addition, it was possible to observe a correlation between the parameters in the viable regions of the Pareto fronts: C1 assuming a constant value, whereas C2 and C3 are linearly dependent on each other

Experiments and qualitative analysis by artificial neural network approach on pool boiling of FC-72 on finned surfaces confined by an unheated horizontal wall

Heat transfer measurements were taken during nucleate pool boiling of FC-72 (a saturated liquid utilized for thermal control of electronic components) from extended finned surfaces placed in a horizontal channel. The influence of channel width on boiling behavior for various extended finned surfaces was experimentally investigated. These finned surfaces have an array of uniformly or non-uniformly spaced straight spines with a square cross-section. The fins had a length of 3 mm or 6 mm. Starting from a uniform configuration, for which the width spacing of the spines was 0.4 mm, non-uniform surfaces were obtained by regularly removing some rows of spines. The effects of channel width (0.5–20 mm) and spines distribution were analyzed. The case of an extended unconfined surface was investigated and used as a reference to compare the effect of the confinement. An Artificial Neural Network (ANN) with one hidden layer was trained with the experimental data to predict the heat flux. An optimization algorithm was implemented to calculate the geometry that maximizes the predicted heat flux operating near the maximum temperature limit under different constrictions as fixed fin height or fixed channel width.

Prediction of fine tailings settlement in pit lakes using a population growth model

A model for predicting the settlement trajectory in saturated fluid fine tailings deposits is described. The model uses a population growth function and leverages established theories in soil mechanics, clay–water surface interactions and biogeochemistry to derive compressibility and permeability functions for predicting the settlement behaviour of fine tailings in deep deposits, typical of end pit lakes. The method is particularly useful for tailings treated with flocculants and/or coagulants with continuously changing compressibility and permeability parameters during deposition. For oil sands or mineral sands fluid tailings treated with or without high doses of a coagulant and/or a flocculant, the consolidation parameters determined from the model are comparable with those measured using standardised large-strain consolidation methods, and the predicted settlement using a modified Gibson’s finite-strain equation closely describes the measured settlement and void ratio profiles in geocolumns.

Charles Cagniard de Latour

Charles Cagniard de la Tour (1777-1859) was a French physicist and chemist that carried a detailed study of the phase transformation of a saturated liquid under the influence of pressure and temperature and proved that inside a sealed glass tube the fluid would become first a saturated vapor, and then turn continuously into superheated vapor, no matter what the pressure. He discovered the existence of a critical point, without using this concept. Cagniard de la Tour used a microscope to investigate the phenomenon of fermentation, particularly by beer yeast, its causes and effects, and proved that it was caused by a living organism and was not a pure chemical reaction. He also studied the pyrolysis of woody matter in a sealed atmosphere and showed that it resulted in the formation of tarry and sticky carbonaceous material, accompanied by a gas atmosphere.@font-face {font-family:"Cambria Math"; panose-1:2 4 5 3 5 4 6 3 2 4; mso-font-charset:0; mso-generic-font-family:roman; mso-font-pitch:variable; mso-font-signature:-536870145 1107305727 0 0 415 0;}@font-face {font-family:Calibri; panose-1:2 15 5 2 2 2 4 3 2 4; mso-font-charset:0; mso-generic-font-family:swiss; mso-font-pitch:variable; mso-font-signature:-536859905 -1073732485 9 0 511 0;}p.MsoNormal, li.MsoNormal, div.MsoNormal {mso-style-unhide:no; mso-style-qformat:yes; mso-style-parent:""; margin:0cm; mso-pagination:widow-orphan; font-size:12.0pt; font-family:"Calibri",sans-serif; mso-ascii-font-family:Calibri; mso-ascii-theme-font:minor-latin; mso-fareast-font-family:Calibri; mso-fareast-theme-font:minor-latin; mso-hansi-font-family:Calibri; mso-hansi-theme-font:minor-latin; mso-bidi-font-family:Arial; mso-bidi-theme-font:minor-bidi;}.MsoChpDefault {mso-style-type:export-only; mso-default-props:yes; font-family:"Calibri",sans-serif; mso-ascii-font-family:Calibri; mso-ascii-theme-font:minor-latin; mso-fareast-font-family:Calibri; mso-fareast-theme-font:minor-latin; mso-hansi-font-family:Calibri; mso-hansi-theme-font:minor-latin; mso-bidi-font-family:Arial; mso-bidi-theme-font:minor-bidi;}div.WordSection1 {page:WordSection1;}

Prediction of the Thermal Conductivity for Liquid Hydrocarbons and Halogenated Hydrocarbons

The thermal conductivity of compressed liquids has received less attention in prediction models compared with saturated liquids. This work proposes a new thermal conductivity prediction model for alkanes, alkenes, and halogenated hydrocarbons for saturated and compressed liquids. The model based on dimensional analysis was developed for n-alkanes. The critical temperature, critical density, molecular weight, reduced temperature, reduced density, and acentric factor of each fluid are required. Introducing the reduced density effectively implements the thermal conductivity calculation of the compressed liquid. Then, the particle density was found to affect the thermal conductivity of halogenated and isomeric hydrocarbons compared with normal hydrocarbons. Hence, a model with a correction term related to the particle density was developed for halogenated and isomeric hydrocarbons. Good agreements between the calculated results and experimental data for n-alkanes and halogenated (or isomeric) hydrocarbons were obtained with average absolute relative deviation values of 1.74% and 3.03%, respectively.

Limiting Performance of the Ejector Refrigeration Cycle with Pure Working Fluids

An ejector refrigeration system is a promising heat-driven refrigeration technology for energy consumption. The ideal cycle of an ejector refrigeration cycle (ERC) is a compound cycle with an inverse Carnot cycle driven by a Carnot cycle. The coefficient of performance (COP) of this ideal cycle represents the theoretical upper bound of ERC, and it does not contain any information about the properties of working fluids, which is a key cause of the large energy efficiency gap between the actual cycle and the ideal cycle. In this paper, the limiting COP and thermodynamics perfection of subcritical ERC is derived to evaluate the ERC efficiency limit under the constraint of pure working fluids. 15 pure fluids are employed to demonstrate the effects of working fluids on limiting COP and limiting thermodynamics perfection. The limiting COP is expressed as the function of the working fluid thermophysical parameters and the operating temperatures. The thermophysical parameters are the specific entropy increase in the generating process and the slope of the saturated liquid, and the limiting COP increases with these two parameters. The result shows R152a, R141b, and R123 have the best performance, and the limiting thermodynamic perfections at the referenced state are 86.8%, 84.90%, and 83.67%, respectively.

Criticality points on the phase diagram of substances as scaling parameters in the law of corresponding states: alkali metals

The law of corresponding states (LCS), based on a generalized van der Waals equation of state (gvdWEoS) has been formulated by rescaling its thermodynamic variables with liquid- vapour critical—point parameters, the spinodal, the thermodynamic limit of superheat (TLS) and the supercritical (SC) point. This gvdWEoS differs from the known van der Waals equation of state (vdWEoS) by the modified attractive term. The scaling properties have been derived on the basis of the gvdWEoS. The expressions for the isothermal compressibility, a derivate of Gibbs free energy, have been derived. In all the cases of rescaling, the gvdWEoS has the similar functional form and the expression for isothermal compressibility has a similar functional form. This shows that the parameters of any spinodal-point on the phase diagram, including the TLS as a particular case, and the SC point, characteristic point on the quasispinodal, provide the alternate scales of parameterizing the thermodynamic properties of substances other than the liquid-vapour critical point. Moreover, an isothermal compressibility regularity is formulated on the basis of the gvdWEoS and it is found to satisfactorily trace the literature data for saturated liquid alkali metals.

Cryogenic processes of SuperKEKB superconducting magnet cryostat with continuous flow of saturated liquid helium

SuperKEKB is a particle accelerator to collide electron and positron beams and employs 55 superconducting (SC) magnets which are accommodated in two cryostats and cooled by subcooled liquid helium (LHe) of ∼ 20 g/s from refrigerators. One cryostat was set up in laboratory for exciting tests and cooled by continuous flow of saturated LHe (<3.0 g/s) from dewars of 1000 L. The differences of two cooling schemes of subcooled and saturated LHe were investigated and described in this paper. The required liquids to cool down the cryostat cold mass (>1700 kg) to 4.5 K were evaluated with a quantified utilization of vapor sensible heat, which fills the large gap between the maximum and minimum liquid requirement presented in most cryogenic literature. The saturated LHe flow was computed to ensure LHe accumulation in the cryostat and immerse SC magnets. The cryogenic processes were monitored by the measured cryogenic temperatures, pressures, flow rates and resistances of SC magnets, which demonstrate cryogenic stabilities with two-phase helium flow and illustrate the agreements with the thermal studies and analyses. This study provides a successful experience to cool SC magnet cryostats of large cold mass with continuous flow of saturated LHe from dewars instead of cryogenic refrigerators.

Brownian motion with radioactive decay to calculate the dynamic bulk modulus of gases saturating porous media according to Biot theory

We present a new stochastic simulation method for determining the long-wavelength effective dynamic bulk modulus of gases, such as ambient air, saturating porous media with relatively arbitrary microgeometries, i.e., simple enough to warrant Biot’s simplification that the fluid and solid motions are quasi-incompressible motions at the pore scale. The simulation method is based on the mathematical isomorphism between two different physical problems. One of them is the actual Fourier heat exchange problem between gas and solid in the context of Biot theory. The other is a diffusion-disintegration-controlled problem that considers Brownian motion of diffusing particles undergoing radioactive-type decay in the pore volume and instant decay at the pore walls. By appropriately choosing the decay time and the diffusion coefficient, the stochastic algorithm we develop to determine the average lifetime of the diffusing particles, directly gives the effective apparent modulus of the saturating fluid. We show how it leads to purely geometric stochastic constructions to determine a number of geometrical parameters. After validating the algorithm for cylindrical circular pores, its power is illustrated for the case of fibrous materials of the type used in noise control. The results agree well with a model of the effective modulus with three purely geometric parameters of the pore space: static thermal permeability divided by porosity, static thermal tortuosity, and thermal characteristic length.

Convective Flow of Blood through a Constricted Cylinder and the Effect of Cholesterol Growth Rate on the Motion in the Presence of a Magnetic Field

An inflammatory disease resulting in the pathological constriction of the intima and media of the arterial system, such as the aorta, causes symptoms like stroke, heart attack, or angina; This research investigates the convective flow of blood through a constricted cylinder and the effect of cholesterol growth rate on the motion in the presence of a magnetic field by transforming the problem into a system of time-dependent partial differential equations. The governing partial differential equations (PDEs) were transformed into dimensionless ordinary differential equations (ODEs) and solved using the Laplace method. After obtaining the analytical solution, Wolfram Mathematica was used to perform the numerical computation where the various physical parameters such as Soret number, radiation number, solutal Grashof number, Schmidt number, and Prandtl number were varied to their impact studied. The study discovered that increasing the Soret number causes an increase in blood velocity, whereas decreasing the Solutal Grashof number decreases blood velocity. Furthermore, an increase in radiation parameter increases the blood velocity, but a Soret number increase results in a decrease in the concentration of cholesterol in the fluid, a Schmidt number, and oscillatory frequency increase, causing the temperature to decrease. This research is useful for clinicians and mathematical modellers who are trying to understand the flow of cholesterol saturated fluid in the human vascular system and how best to proffer analytical and numerical solutions in synergy with laboratory investigation.

Comprehensive study on EEVs regulation characteristics of vapor-injection CO2 heat pump for electric vehicles

The regulation characteristics of electronic expansion valves (EEVs) are critical to improving the comprehensive energy efficiency of the heat pump system for electric vehicles. In this paper, a vapor-injection transcritical CO2 heat pump system for electric vehicles was first established to experimentally study the impacts of EEV1 and EEV2 openings on the overall system performance. The results show that the cooling capacities of the vapor-injection system are respectively increased by 3.13% / 5.75% with the EEV1’s opening of 15% / 25% and injection ratio of 0.0975 / 0.8305, compared against these of the basic system without vapor injection. As the injection ratio rises, the cooling capacity increment tends to increase while the cooling coefficient of performance (COP) declines from 1.7 to 1.51. Then, the regulation characteristics of EEV1 and EEV2 were quantitively analyzed, indicating that the evaporation pressure is mainly determined by EEV2’s openings. By adjusting EEV1’s opening, the COP of the CO2 vapor-injection heat pump can be maximized at 1.9, with EEV1’s opening at 20% and outdoor/indoor ambient temperature at 0°C / 20°C. When adjusting EEV2’s opening, the maximum injection pressure corresponds to the maximum COP of 1.84, with EEV2’s opening at 65% under 35°C/ 27°C. When the COP reaches the maximum, the refrigerant at the outlet of EEV1 is just saturated liquid under cooling conditions, while it is located in the two-phase region under heating conditions. Besides, the optimal COP regulation strategy of the CO2 vapor-injection heat pump was further formulated based on two modes, as the strategy in the cooling mode differs from that in the heating mode.

Evaluation of machine learning approaches for estimating thermodynamic properties of new generation refrigerant R513A

In this study, direct thermodynamic properties of R513A refrigerant such as temperature, pressure, enthalpy, and entropy were estimated using polynomial regression and GradientBoosting, XGBR, KNNRegressor, and Decision Tree as machine learning algorithms. The low-order polynomial regression results for the enthalpy and entropy of saturated vapor reached an R2 value of 0.9939 and below, while the 4th-degree polynomial regression reached an R2 value of 1.0 for the pressure of saturated vapor. Compared to other machine learning algorithms, the KNNRegressor reached an R2 value of 0.9947 and it gave better results than other methods for the saturated vapor. For enthalpy calculations of superheated vapor, R2 = 0.9984 result was reached in KNNRegressor, while in entropy calculations, XGBR and KNNRegressor gave 0.9960 and 0.9998 for R2 values, respectively. KNNRegressor gives better results in modeling with an R2 value of 0.9946 and above for the saturated liquid properties of R513A. According to the results obtained in this study, polynomial regression and machine learning methods achieve success in estimating properties of R513A refrigerant with an R2 value of 0.99 and above. The consequence of this study showed that machine learning methods can be used successfully in complicated processes such as predicting the thermodynamic properties of refrigerants.

Evidence of Nonlinear Seismic Effects in the Earth from Downhole Distributed Acoustic Sensors

Seismic velocities and elastic moduli of rocks are known to vary significantly with applied stress, which indicates that these materials exhibit nonlinear elasticity. Monochromatic waves in nonlinear elastic media are known to generate higher harmonics and combinational frequencies. Such effects have the potential to be used for broadening the frequency band of seismic sources, characterization of the subsurface, and safety monitoring of civil engineering infrastructure. However, knowledge on nonlinear seismic effects is still scarce, which impedes the development of their practical applications. To explore the potential of nonlinear seismology, we performed three experiments: two in the field and one in the laboratory. The first field experiment used two vibroseis sources generating signals with two different monochromatic frequencies. The second field experiment used a surface orbital vibrator with two eccentric motors working at different frequencies. In both experiments, the generated wavefield was recorded in a borehole using a fiber-optic distributed acoustic sensing cable. Both experiments showed combinational frequencies, harmonics, and other intermodulation products of the fundamental frequencies both on the surface and at depth. Laboratory experiments replicated the setup of the field test with vibroseis sources and showed similar nonlinear combinations of fundamental frequencies. Amplitudes of the nonlinear signals observed in the laboratory showed variation with the saturating fluid. These results confirm that nonlinear components of the wavefield propagate as body waves, are likely to generate within rock formations, and can be potentially used for reservoir fluid characterization.

Preparation of CMAS, LZS glass-ceramics and study of the glass ceramics anti-corrosion properties in liquid leadbismuth eutectic at 500°C

Liquid lead-bismuth eutectic (LBE) was an important candidate material for the fourth generation reactor coolant. Its high temperature and heavy metal corrosion conditions put forward higher requirements for the corrosion resistance of reactor structural materials. As a new kind of special material, glass-ceramic material, which possesses the excellent properties of glass and ceramics, has good mechanical property, wear-resisting property and corrosion resistance. The LBE corrosion behavior of glass-ceramics at typical operating temperatures of nuclear reactors was investigated. The lead bismuth corrosion test was carried out under the condition of oxygen saturation at 500°C for 500h. The lead bismuth corrosion resistance of the two glass ceramic samples in static LBE was analyzed by means of scanning electron microscope (SEM) and X-ray diffraction (XRD). After 500 h corrosion in oxygen saturated liquid LBE at 500°C, there was no obvious corrosion on the surface of CMAS glass ceramics. Its dense microstructure and stable phase composition effectively prevented the penetration of Pb and Bi elements, and had good resistance to corrosion; The surface of LZS glass ceramics was covered with a layer of loose and porous compounds, and the thickness of the corrosion zone was 10 μm, which occurred a certain degree of lead-bismuth corrosion, and its corrosion resistance was poor. The results showed that under the same corrosion conditions, no lead bismuth corrosion behavior was found on the surface of CaO-MgO-Al2O3-SiO2 (CMAS) glass-ceramic, while there was a corrosion layer with a thickness of about 10 μm on the surface of Li2O-ZnO-SiO2 (LZS) glass-ceramic, indicating that the lead bismuth corrosion resistance of CMAS glass ceramics was significantly better than that of LZS glass ceramics.

Pore morphology effect on elastic and fluid flow properties in Bakken formation using rock physics modeling

Unconventional geo-resources are critical due to their important contributions to energy production. In this energy transition and sustainability era, there is an increased focus on CO2-enhanced oil recovery (CO2-EOR) and geological CO2 storage (GCS) in unconventional hydrocarbon reservoirs, and the extraction of hot fluid for energy through enhanced geothermal systems. However, these energy solutions can only be achieved through efficient stimulation to develop a complex fracture network and pore structure in the host rocks to extract heat and hydrocarbon, or for CO2 storage. Using Bakken formation well data and rock physics models, this study aimed to identify the post-depositional effect of pore structure on seismic velocity, elastic moduli, and formation fluid; and further predict the best lithofacies interval for well landing, and the implications for fluid (gas, oil, and water) recovery in naturally- and often systematically-fractured geosystems. The KT and DEM models' predictions show distinct formation intervals exhibiting needle-like pores and having higher seismic velocities (Vp and Vs) and elastic moduli (K and µ), relative to other formation intervals that exhibit moldic pores. At the same fluid concentration, the needle-like pores (small aspect ratios) have a higher impact on elastic moduli, Vp, and Vs than on the moldic spherical pores with all other parameters held constant. Vp is affected more than Vs by the properties of the saturating fluid (gas, oil, or water) with Vp being greater in Bakken formation when it is water-saturated than when it is gas-saturated. Vs exhibit the reverse behavior, with Vs greater in the gas-saturated case than in the water-saturated case. Further, analyses suggest that the middle Bakken formation will have a higher susceptibility to fracturing and faulting, and hence will achieve greater fluid (oil and water) recovery. Our findings in this study provide insights that are relevant for fluid production and geo-storage in unconventional reservoirs.Article highlightsIntegrated well log data and rock physics models.Investigated the effect of changes in pore structure on elastic properties and fluid flow in shale.Increase in porosity causes a reduction in elastic moduli and seismic velocities.Vp is more affected by pore geometry than Vs depending on density and properties of saturating fluid.Lithofacies with needle−like pores are more susceptible to fracturing than lithofacies with intragranular pores.