Search NASA⌕ Search

SEARCH · Search NASA

Results for “diameter”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 361 records · Page 20

Numerical analysis of coalescence-induced bubble departure for enhanced boiling heat transfer

Boiling heat transfer plays a crucial role in a wide range of applications, such as power generation, refrigeration, electronics cooling, and pharmaceutics. Among the various factors that influence boiling heat transfer, the dynamics of vapor bubble nucleation, growth, and departure from the heated surface stand out as particularly important. An emerging phenomenon that can promote the departure of bubbles smaller than the Fritz diameter is coalescence-induced departure. If the dynamics of this process are fully understood, then surfaces can be engineered to promote faster bubble departure and substantially increase the performance of boiling heat transfer. Further, this work expands on published results by presenting a detailed numerical analysis of bubble coalescence and departure for a range of initial bubble diameters and size ratios between coalescing bubbles. Analysis of the results is focused on explaining how the release of surface energy and bubble surface dynamics lead to bubble departure, as well as fundamentally distinguishing capillary–inertial jumping and buoyant–inertial departure mechanisms across different bubble sizes and size ratios. The results show that both the initial sizes of the coalescing bubbles and the ratio between their sizes can determine whether the merged bubble will leave the surface through capillary–inertial jumping or buoyant departure. Below a certain bubble size, the release of surface energy by the merger is not sufficient to propel the merged bubble from the surface.

42 ENGINEERING↗

Thermohydraulic experiments on a supercritical carbon dioxide–air microtube heat exchanger

Heat exchangers are critical components in supercritical CO 2 Brayton cycles and provide necessary heat input, recovery, and dissipation. Tubular heat exchangers with unconventionally small tube sizes (tube diameters less than 5 mm) are promising components for supercritical CO 2 cycles and potentially provide excellent structural stability with wide scope of application. This paper provides essential design and fabrication guidelines for a compact shell-and-tube heat exchanger with microtubes (with an inner diameter equal to 1.75 mm). A heat exchanger test rig is used to evaluate the thermohydraulic performance of this heat exchanger with supercritical CO 2 and air as working fluids. Thermohydraulic data are reported for more than forty sets of experiments with varying Reynolds numbers for shell and tube flows. Critical performance metrics are calculated from the data and compared with predictions from a previously described numerical model. The average deviations between the experimental and model results fall within 10% for all critical metrics. This excellent agreement validates the numerical model for supercritical CO 2 heat exchanger optimization and scale-up.

42 ENGINEERING↗

Detailed mixing measurements from single-hole converging ECN injectors using Rayleigh scattering

Time resolved liquid and vapor fields of dodecane and oxymethylene ethers are measured from Spray A-3 and Spray D using high speed Rayleigh scattering and diffuse back illumination at the Engine Combustion Network (ECN) Spray A condition of 900 K and 22.8 kg/m 3 . Global quantities including mixture fraction, vapor and liquid penetration, as well as spreading angle are measured. The mixture fraction fields and vapor penetration profiles are well predicted by the 1-D Musculus-Kattke model. The mixture fraction field and vapor penetration from Spray A-3 are similar to those measured from Spray A in previous works. Spray D exhibits higher mixture fraction fields and vapor penetration due to its larger nozzle diameter. The quasi-steady mixture fraction fields from these injectors scale well when distance from the injector is normalized by the nozzle diameter. The turbulent dissipation structures were also analyzed based on the orientation, thickness, and magnitude of the mixing layers. Furthermore, the orientation and thickness are similar to other measurements in atmospheric gas jets, while the magnitude is lower. The thickness and magnitude are subject to uncertainties due to limitations in the imaging resolution of the system but still provide an order of magnitude as a useful reference for comparison against computational fluid dynamic simulations.

33 ADVANCED PROPULSION SYSTEMS↗

The development of natural gas and hydrogen pipeline capital cost estimating equations

Natural gas pipeline cost data collected by the Oil and Gas Journal (O&GJ) for interstate pipelines constructed from 1980 through 2017 were used to develop capital cost estimating equations that are a function of pipeline diameter, length, and U.S. region. Equations were developed for material, labor, miscellaneous, and right-of-way costs, the four cost components in the O&GJ data, for six different regions of the United States (U.S.). Each equation is a function of pipeline diameter and length. Adjustment mechanisms were then developed for converting the natural gas pipeline equations into equations for estimating the costs of hydrogen pipelines. These adjustments were based in part on an analysis completed by the National Institute for Standards and Technology (NIST). Here the results of this work were used to update cost models in the Hydrogen Delivery Scenario Analysis Model (HDSAM), developed by Argonne National Laboratory for the U.S. Department of Energy’s Hydrogen Program. Our analysis shows a wide range of pipeline cost across different U.S. regions, especially with respect to labor and right-of-way costs. The developed cost formulas for hydrogen pipelines are both important and timely as hydrogen is being considered as a zero-carbon energy carrier with the potential to decarbonize all energy sectors, and the cost of hydrogen transportation is essential for techno-economic analysis of its potential use in these sectors.

03 NATURAL GAS↗

Thermodynamic parametric analysis of refueling heavy-duty hydrogen fuel-cell electric vehicles

Reliable design and safe operation of heavy-duty hydrogen refueling stations are essential for the successful deployment of heavy-duty fuel cell electric vehicles (FCEVs). Fueling heavy-duty FCEVs is different from light-duty vehicles in terms of the dispensed hydrogen quantities and fueling rates, requiring tailored fueling station design for each vehicle class. In particular, the selection and design of the onboard hydrogen storage tank system and the fueling performance requirements influence the safe design of hydrogen fueling stations. A thermodynamic modeling and analysis are performed to evaluate the impact of various fueling parameters and boundary conditions on the fueling performance of heavy-duty FCEVs. Here, we studied the effect of dispenser pressure ramp rate and precooling temperature, initial tank temperature and pressure, ambient temperature, and onboard storage design parameters, such as onboard storage pipe diameter and length, on the fueling rate and final vehicle state-of-charge, while observing prescribed tank pressure and temperature safety limits. An important finding was the sensitivity of the temporal fueling rate profile and the final tank state of charge to the design factors impacting pressure drop between the dispenser and vehicle tank, including onboard storage pipe diameter selection, and flow coefficients of nozzle, valves, and fittings. The fueling rate profile impacts the design and cost of the hydrogen precooling unit upstream of the dispenser.

08 HYDROGEN↗

Droplet-Gas Interactions in Nanoelectrospray Multiphase Flow

Electrospray (ES) is production of a charged droplet plume that transits from an electrically biased supply capillary to a counter electrode under the influence of an electric field. Drag interactions with the surrounding gas decelerate the droplets while transferring the momentum from droplets to accelerate the surrounding gas, resulting in an induced gas flow. This work seeks to characterize the structure of the resulting gas jet and identify the key mechanisms defining the flow structure. The phenomenon of gas jets produced by momentum transfer from nano-electrospray (nES) plumes is explored with schlieren visualization, thermal anemometry, and numerical simulations. Schlieren visualization experiments provide information on the flow structure in support of simulation predictions, and the hot thermistor anemometry measurements of gas velocities outside the spray demonstrate quantitatively validated simulation results. The study reveals the formation of a moderately high velocity coaxial gas jet within the nES plume and provides insight into the overall flow structure of the induced flow. The multiphase electrohydrodynamic simulations enable numerical experimentation to explore the fundamental physics of coupled droplet-gas transport and the resulting flow structure. The simulations, confirmed by the experiments, reveal gas jetting induced by nES with a narrow (sub hundred micrometers in diameter) core originating from the nES liquid-jet breakup region and a surrounding larger-in-extent zone (several hundred micrometers in diameter) of lower velocity gas flow within the electrospray plume. This behavior is due to nES ejecting a stream of droplets from the tip of a narrow liquid cone-jet, where the combined effect of many small droplets transferring momentum to a confined region of gas yields a narrow but high-velocity gas stream. Furthermore, the micro-jets produced by nES have practical utility for mass spectrometry, 3D printing and fabrication, and thermal management.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

A parametric finite element study for determining burst strength of thin and thick-walled pressure vessels

To accurately predict the burst strength of both thin and thick-walled pressure vessels (PVs), a parametric study of PV burst strength was performed for a wide range of vessel geometries and materials using elastic-plastic finite element analysis (FEA). A valid FEA model was established through a detailed study of 2D versus 3D FEA models, the critical stress failure criterion versus the limit load criteria, and the thick-wall effect on the FEA simulations. Here, the results show that the stresses and strains at the mean diameter, rather than outside diameter, determines a more accurate burst strength for both thin and thick-walled PVs. On this basis, a parametrized FEA script using the ABAQUS Python application programming interface (API) was used to create a large database of PV burst strengths for a variety of vessel geometries and materials, demonstrating that Python scripting is a powerful technique for performing parametric studies or generating large databases. From the FEA results, using the regression method, a new burst pressure model was developed as a function of the vessel geometry (D/t ratio) and material properties (UTS and n). As validated by a large number of full-scale burst test data, the proposed burst model can very accurately predict the burst strength for both thin and thick-walled PVs.

42 ENGINEERING↗

Effect of micro-fin tube geometry parameters on evaporator material savings

A Performance Evaluation Criteria (PEC) analysis is presented that explores the influence of the micro-fin tube geometry on the required material volume. The analysis was done using a local flow boiling model, which can quantify the effect of fin parameters and the tube diameter. The PEC was done while requiring that the compared surfaces satisfy the same heat duty and pressure drop. The following tube parameters and ranges were explored: tube diameter (3 mm ≤ D o ≤ 6 mm), fin base thickness (0.12 mm ≤ t bf ≤ 0.29 mm), helix angle (6° ≤ α ≤ 30°), fin apex angle (11° ≤ β ≤ 22°), total number of fins (43 ≤ n f ≤ 70), and the fin height (0.11 mm ≤ e ≤ 0.3 mm). The tube that provided the smallest material requirements within the examined parameter ranges had: D o = 3 mm, t bf = 0.12 mm, α = 6°, β = 11°, n f = 43, and e = 0.3 mm. In general, small D o and t bf and large n f and e provided the best performance and work to minimize the required material volume to provide a given heat duty and pressure drop. The helix and the fin-tip apex angles were found to have very little influence on the performance for the examined range of parameters and angles associated with the lowest production cost should be chosen.

Evaporation↗

Evaluation of methods for characterizing the fine particulate matter emissions from aircraft and other diffusion flame combustion aerosol sources

The U. S. Environmental Protection Agency in collaboration with the U. S. Air Force Arnold Engineering Development Complex conducted the VAriable Response In Aircraft nvPM Testing (VARIAnT) 3 and 4 test campaigns to compare nvPM emissions measurements from a variety of diffusion flame combustion aerosol sources (DFCASs), including a Cummins diesel engine, a diesel powered generator, two gas turbine start carts, a J85-GE-5 turbojet engine burning multiple fuels, and a Mini-CAST soot generator. The VARIAnT research program has conducted four test campaigns to date with the VARIAnT 3 and 4 test campaigns focusing on BC mass instrument performance, use of gas turbine engine sources for BC mass instrument calibration, and the characterization of the aerosols produced from the combustion sources. Results from both campaigns revealed agreement of about 20% between the AVL Micro-Soot Sensor, the Cavity Attenuated Phase Shift (CAPS PMSSA) monitor and the thermal-optical reference method for elemental carbon (EC), independent of the calibration source used. For the Artium LII-300, the measured mass concentrations in VARIAnT3 fall within 18% and in VARIAnT4 fall within 32% of the reference EC mass concentration when calibrated on a combustor rig in VARIAnT3 and on a LGT-60 start cart in VARIAnT4, respectively. It was also found that the three mass instrument types (MSS, CAPS PMSSA, and LII-300) can exhibit different BC to reference EC ratios depending on the emission source that appear to correlate to particle geometric mean mobility diameter, morphology, or some other parameter associated with particle geometric mean diameter (GMD) with the LII-300 showing a slightly stronger apparent trend with GMD. Systematic differences in LII-300 measured mass concentrations have been reduced by calibrating with a turbine combustion source (combustor or turbine engine). With respect to the particle size measurements made, the sizing instruments (TSI SMPS, TSI EEPS, and Cambustion DMS 500) were found to be in general agreement in terms of size distributions and concentrations with some exceptions. Teflon filter measurements of the total aerosol mass produced by the various DFCASs differed from the reference EC, BC, and integrated particle sized distribution measured aerosol masses. The Teflon filter measured mass was 30 to 40% greater than the reference EC and 50 to 60% greater than the BC mass for the J85 in both test campaigns. The measurements of particle size distributions and single particle analysis by miniSPLAT indicated the presence of large particles (>100 nm) having more compact morphologies, higher effective density, and a composition dominated by OC and ash. This increased large particle fraction is also associated with higher values of single scattering albedo measured by the CAPS PMSSA instrument and higher OC/EC measurements. These measurements indicate gas turbine engine emissions can be a more heterogeneous mix of particle types beyond the original E31 assumption that it is mainly composed of black carbon.

33 ADVANCED PROPULSION SYSTEMS↗

Growth of CuFeO 2 single crystals by the optical floating-zone technique

CuFeO 2 single crystals up to 50 mm in length and up to 10 mm in diameter were grown by the optical floating-zone method. Stoichiometric polycrystalline rods with a diameter of 6–12 mm were used as feed materials to produce crystals of sufficient size to be used as substrates for the growth of thin films of delafossites. For stable growth along the c-axis, low growth rates of 0.4 mm/h are necessary. Due to the incongruent melting behavior of CuFeO 2 , a stable melt zone requires adjustment of the lamp power during growth. The melting of CuFeO 2 is not simply incongruent because the thermodynamic equilibrium includes more than two solid phases and the melt; the gas phase is also involved. Finaly, the crystals were characterized by X-ray diffraction and X-ray fluorescence measurements.

36 MATERIALS SCIENCE↗

Precious metal oxygen-evolving anodes for electrolytic reduction of metal oxides in molten LiCl-Li 2 O electrolyte

Understanding the electrochemical stability of oxygen-evolving anode materials in molten salt electrolytes is essential to enable decarbonized electrolytic reduction of metal oxides (e.g., used nuclear oxide fuels). Here, this work investigated three precious metals (Ir, Ru, and Pt) as oxygen-evolving anodes in molten LiCl-Li 2 O (99.0-1.0 wt%) at 650 °C. For consistent measurements, this work employed a three-electrode cell comprised of a two-phase Li-Bi (65-35 at%) reference electrode and a NiO counter electrode. Anodic polarization behavior of each anode was investigated via cyclic voltammetry (CV) and chronoamperometry. The onset potential for oxygen evolution was observed at E > 2.9 V (vs. Li/Li + ) for the Ir and Ru anodes and anodic current density was as high as 1.0 A cm -2 at 3.23 V. The dimensional stability of each anode was evaluated from long-term electrolysis experiments (10.0-32.1 h) at 3.23 V. Rapid consumption of the Pt anode was observed after the application of 26,453 C cm -2 with a reduction in diameter of about 15.8% due to the formation of a non-protective Li 2 PtO 3 compound. The formation of this compound was also observed during CV measurements as additional anodic waves at potentials more negative than that of oxygen evolution. In contrast, both the Ir and Ru anodes exhibited excellent dimensional stability with a reduction in diameter or thickness of less than 1.5% even after applying greater charge density of ~41,100 C cm -2 , demonstrating superior stability during oxygen evolution in the LiCl-Li 2 O electrolyte.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Simulation of radioactive plume transport in the atmosphere including dynamics of particle aggregation and breakup

Accurate prediction of the atmospheric transport of debris particles relies heavily on our knowledge of the size distribution of the particles within a debris cloud. Assuming a fixed particle size during simulations is not always viable since the size distribution of the debris can change during transport. Various microphysical processes, such as aggregation and breakup, influence debris particles and dictate any changes to the size distribution. To track those changes that can occur, a population balance model can be adopted and instituted within a model framework. Nonetheless, many of the models that simulate the transport of radioactive debris following a device-driven fission incident have historically neglected to consider these processes. As such, this work describes our effort to develop a modeling framework capable of simulating the transport and deposition of a radioactive plume generated from a fission incident with a dynamic population balance including particle aggregation and breakup. Herein the impact of aggregation and breakup, individually and collectively, on the particle size distribution is explored using the developed framework. When simulating aggregation, for example, six mechanisms, including Brownian coagulation, the convective enhancement to Brownian coagulation, van der Waals-viscous force correction for Brownian coagulation, gravitational collection, turbulent inertial motion, and turbulent shear, are considered. Brownian coagulation and its corrections have, as one would expect, a large impact on relatively small aggregates. Aggregates with a diameter that is less than or equal to 1.0 μm, for instance, comprise 50.6 vol % of all aggregates in the absence of aggregation and 31.2 vol % when Brownian coagulation and its corrections are considered. Gravitational collection and, to a much lesser extent, turbulent shear and turbulent inertial motion are, conversely, of great importance to relatively large aggregates (i.e., diameter greater than 3.0 μm). Additionally, the individual effects of atmospheric and particle parameters, such as wind speed and particle density, are examined. Of the parameters examined, turbulent energy dissipation and aggregate fractal dimension (i.e., aggregate shape with lower values representing more irregular particles) were of substantial importance since both terms directly impact aggregate stability and, by extension, the breakup rate. Large-scale transport and deposition simulations in a dry atmosphere are also presented and discussed as a proof of concept.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Impact of lead on an axisymmetric, single bead blown powder DED overhung geometry

Metal additive manufacturing can be utilized for the near net-shape manufacture of components in a layer wise technique. Traditionally, the buildup process is conducted in the direction of gravity or bottom to top in the vertical orientation; however, with the availability of commercial systems with additive heads or fixturing that can index from vertical, and the need to manufacture larger parts without modifying available systems, the limitation of manufacturing from bottom to top is removed enabling the buildup of larger components printed at an angle from vertical. The effect of gravity on the melt pool and as-deposited component quality when printing off vertical is unknown in literature, especially for axisymmetric components. This understanding is critical for the advancement of manufacturing components of increased size in 4 and 5-axis additive systems. This investigation utilizes blown powder directed energy deposition to evaluate the change in as-printed geometry when the start point is altered in relation to gravity while the part rotates to manufacture an axisymmetric component. The objectives of this work are to determine the impact of the deposition location on the geometric variability on axisymmetric components. This investigation tests the hypothesis that differences in layer height exist due to a change in catchment efficiency when the deposition location is moved to a tangent of the round geometry due to a change in the melt pool dynamics. It was found the ideal deposition location when printing at 26.6-degrees from vertical was at −90-degrees from the top center point of the component, along the tangent, where gravity was pushing the melt pool down, but the rotation of the part was pulling the deposited material towards the top center of the component. This work provides an understanding of layer height stability and catchment efficiency to guide print orientation strategy for high-aspect ratio components. It was found the −90-degree lead deposition had the best layer height stability at 0.6 % as compared to the top-center at 6.3 % and +90-degree deposition location at 9.2 % for the nominal programmed layer height. The change in layer height also effected the final diameter the greatest for the top center deposition location where the diameter diverged by 2.5 % during the overbuild condition and converged by 0.6 %. In conclusion, this finding will increase the manufacturing efficiency of axisymmetric components by increasing the passive stability of the printing process for the successful manufacture of parts without the need for perfectly calibrated manufacturing parameters.

36 MATERIALS SCIENCE↗

Fabrication of thin-walled tubes from alloy 602 CA using shear assisted processing and extrusion

Thin-walled tube is usually produced seamlessly via an expensive multi-step process or via welding of thin sheets but at reduced performance. To overcome the current process and performance inefficiencies, advanced manufacturing methods need to be investigated. In this investigation, shear assisted processing and extrusion (ShAPE) was used to fabricate thin-walled Inconel 602 tube in a single-step. Tubes measuring 0.83 m length with a 12 mm outer diameter and 1 mm wall thickness were fabricated with an average surface roughness, Ra of 1.6 µm and Rz of 15 µm. Tensile testing of tubes in the as-fabricated condition showed strength increases of 4-35% over current offerings while maintaining or improving elongation. Electron microscopy analysis revealed the recrystallized microstructure with refined inter- and intra-granular carbides. Preliminary results obtained in this investigation shows the feasibility of producing thin-wall high temperature tube in single-step via ShAPE. This achievement marks a significant step towards manufacturing larger diameter nickel alloy tubes for the U.S. Department of Energy’s Waste Treatment and Immobilization Plant.

36 MATERIALS SCIENCE↗

A improved equation of state for Xe gas bubbles in gamma U-Mo fuels

A monolithic fuel design based on a U-Mo alloy has been selected as the fuel type for conversion of the United States High-Performance Research Reactors (HPRRs). An issue with U-Mo monolithic fuel is the large amount of swelling that takes place during operation. The accurate prediction of fuel evolution under irradiation requires implementation of correct thermodynamic properties into mesoscale and continuum level fuel performance modeling codes. However, the thermodynamic properties of the fission gas bubbles (such as the relationship among bubble size, equilibrium Xe concentration, and bubble pressure) are not well known. This work studies Xe bubbles in ?U-Mo from a diameter of 3 nm up to 8.5 nm and from 400 K up to 700 K. The energetic relationship of Xe bubbles with regard to voids and Xe substitutional atoms is described. The transition is also determined for when a bubble becomes over-pressurized. Finally, an equation of state is fit to the pressure as a function of molar volume and temperature.

Beeler, Benjamin↗

Sintering effects on NpO 2 grain size and morphology: The role of precursor Np phase

Neptunium dioxide (NpO 2 ) is a key phase in nuclear material processing as a target material for the production of plutonium-238 ( 238 Pu) and has been historically synthesized via the calcination of a Np oxalate precursor. Alternative synthesis methods for NpO 2 are now more prevalent, necessitating their study and comparison with the more common oxalate route. The purpose of this work was to investigate the microstructural properties of NpO 2 synthesized from a nitrate-based Np precursor phase via the assessment of NpO 2 particle size and morphology as a function of calcination temperature and time. Scanning electron microscopy (SEM) was used to probe the primary grain size and morphology of NpO 2 after calcination at temperatures ranging from 700 to 1100 °C and hold times ranging from 1 to 10 h. Post-image analysis using ImageJ software enabled the quantification of mean particle diameter. This analysis indicated that particle diameter increases with both increasing calcination temperature and hold time. Primary particles were shown to be clumped in irregular patterns into the overall rough, blocky aggregates, but this macroscopic morphology was not affected by calcination time or temperature. Although trends in primary grain size of NpO 2 were consistent with available literature from other Np precursor phases, the macroscopic morphology of the NpO 2 aggregates was quite different than reported for other precursors. Through comparison with historical literature on Np oxalate, this work emphasizes the importance of Np precursor on the physical properties of NpO 2 .

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Optimizing Porous Transport Layer Porosity for Proton Exchange Membrane Water Electrolysis

An empirical model is presented that describes anode-side losses related to porous transport layer (PTL) morphology in proton exchange membrane water electrolysis (PEMWE). The model is based on an advanced voltage breakdown analysis that links various overpotentials to PTL morphology. Custom Ti PTLs, spanning uncommonly low porosities (22 - 31%), were fabricated and analyzed with X-ray CT to obtain pore and particle size distributions. Particle size distributions were consistent across samples with an average particle diameter of 12.0?..mu..m, whereas average pore diameters ranged from 6.0 to 7.0?..mu..m. The PTLs were tested in standard PEMWE cell assemblies with anode catalyst loadings of 0.1 mgIr cm-2 to obtain polarization curves, electrochemical impedance spectra, and augmented Tafel analysis. The PTL-dependent anode side losses were deconvoluted and assigned to excess utilization, concentration, ion transport resistance, and electrical contact resistance overpotentials. The data and model reveal an optimal 20 - 28% PTL porosity region where utilization and contact resistance overpotentials are minimized without triggering concentration and ion transport losses related to water deprivation. The optimal PTL porosity depends on the operating current density and is demonstrated at realistic PEMWE water flow rates to establish PTL design guidance for operation at scale.

08 HYDROGEN↗

Factors controlling pore network development of thermally mature Early Palaeozoic mudstones from the Baltic Basin (N Poland)

Understanding the formation of pore space, especially in low porosity shales (as source rocks and as unconventional resources), is critical to the oil and gas industry, since pores control the space available for hydrocarbon and participate in hydrocarbon transport. We examined 87 Ordovician and Silurian mudstone samples collected from four wells located in the Pomeranian part of the Baltic Basin (northern Poland), one of the primary Polish targets for hydrocarbon exploration. These samples represent the Pelplin, the Paslek, the Jantar, the Prabuty, and the Sasino Formations, which still requires more detailed porosity studies. Our study aimed to identify factors controlling porosity development, by applying bulk techniques (organic petrology and TOC analyses, quantitative mineralogy, and porosimetry) as well as nano- to microscale techniques (thin section petrography, electron microscopy). The studied samples are mainly argillaceous mudstones. The results of porosimetry measurements, combined with image analysis, indicate that the pores of all studied rocks are dominated by micropores (pores < 2 nm in diameter), mesopores (2-50 nm in diameter) and small macropores. The SEM images showed three main pore types: a) voids related to clay mineral aggregates, b) pores inside organic matter particles, and c) pores between other mineral grains. In the Jantar and Sasino mudstones, the organic matter content and its thermal maturity control porosity. The occurrence of solid bitumen in the rocks from these formations reduces samples’ mesoporosity because of the pore-clogging effect. In contrast, in the Paslek and Prabuty Formations, there is low organic matter content and specific surface area and the volume of mesopores increase with clay minerals content. In the Pelplin mudstones, there are no prevailing factors controlling porosity. As a result, we suggest that a combination of SEM image analysis and dual liquid porosity (DLP) measurements is a powerful method to assess porosity available for petroleum flow in mudstones.

03 NATURAL GAS↗