Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Weak form”

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 451 records · Page 25

Ultra compact Bragg grating devices with broadband selectivity

Current silicon waveguide Bragg gratings typically introduce perturbation to the optical mode in the form of modulation of the waveguide width or cladding. However, since such a perturbation approach is limited to weak perturbations to avoid intolerable scattering loss and higher-order modal coupling, it is difficult to produce ultra-wide stopbands. In this Letter, we report an ultra-compact Bragg grating device with strong perturbations by etching nanoholes in the waveguide core to enable an ultra-large stopband with apodization achieved by proper location of the nanoholes. With this approach, a 15 µm long device can generate a stopband as wide as 110 nm that covers the entire C + L band with a 40 dB extinction ratio and over a 10 dB sidelobe suppression ratio (SSR). Similar structures can be further optimized to achieve higher SSR of > <#comment/> 17 d B for a stopband of about 80 nm.

Li, Ang (ORCID:0000000232251280)↗

Flutter analysis using transversality theory

A new method of calculating flutter boundaries of undamped aeronautical structures is presented. The method is an application of the weak transversality theorem used in catastrophe theory. In the first instance, the flutter problem is cast in matrix form using a frequency domain method, leading to an eigenvalue matrix. The characteristic polynomial resulting from this matrix usually has a smooth dependence on the system's parameters. As these parameters change with operating conditions, certain critical values are reached at which flutter sets in. Our approach is to use the transversality theorem in locating such flutter boundaries using this criterion: at a flutter boundary, the characteristic polynomial does not intersect the axis of the abscissa transversally. Formulas for computing the flutter boundaries and flutter frequencies of structures with two degrees of freedom are presented, and extension to multi-degree of freedom systems is indicated. The formulas have obvious applications in, for instance, problems of panel flutter at supersonic Mach numbers.

Afolabi, D.↗

Effect of high-latitude ionospheric convection on Sun-aligned polar caps

A coupled magnetospheric-ionospheric (M-I) magnetohydrodynamic (MHD) model has been used to simulate the formation of Sun-aligned polar cap arcs for a variety of interplanetary magnetic field (IMF) dependent polar cap convection fields. The formation process involves launching an Alfven shear wave from the magnetosphere to the ionosphere where the ionospheric conductance can react self-consistently to changes in the upward currents. We assume that the initial Alfven shear wave is the result of solar wind-magnetosphere interactions. The simulations show how the E region density is affected by the changes in the electron precipitation that are associated with the upward currents. These changes in conductance lead to both a modified Alfven wave reflection at the ionosphere and the generation of secondary Alfven waves in the ionosphere. The ensuing bouncing of the Alfven waves between the ionosphere and magnetosphere is followed until an asymptotic solution is obtained. At the magnetosphere the Alfven waves reflect at a fixed boundary. The coupled M-I Sun-aligned polar cap arc model of Zhu et al.(1993a) is used to carry out the simulations. This study focuses on the dependence of the polar cap arc formation on the background (global) convection pattern. Since the polar cap arcs occur for northward and strong B(sub y) IMF conditions, a variety of background convection patterns can exist when the arcs are present. The study shows that polar cap arcs can be formed for all these convection patterns; however, the arc features are dramatically different for the different patterns. For weak sunward convection a relatively confined single pair of current sheets is associated with the imposed Alfven shear wave structure. However, when the electric field exceeds a threshold, the arc structure intensifies, and the conductance increases as does the local Joule heating rate. These increases are faster than a linear dependence on the background electric field strength. Furthermore, above the threshold, the single current sheet pair splits into multiple current sheet pairs. For the fixed initial ionospheric and magnetospheric conditions used in this study, the separation distance between the current pairs was found to be almost independent of the background electric field strength. For either three-cell or distorted two-cell background convection patterns the arc formation favored the positive B(sub y) case in the northern hemisphere.

Sojka, J. J.↗

Geometry of interplanetary magnetic clouds

Two dimensional magnetohydrodynamic simulations are presented of the distortion of a magnetic flux rope that is being accelerated through ambient solar wind plasma. The flux rope magnetic field has an axial component parallel to the solar wind field and an azimuthal component, which lies in the simulation plane. As the flux rope moves through the solar wind plasma, vortices form on its trailing edge and couple strongly to its interior. If the flux rope azimuthal field is weak, it deforms into an elongated banana-like shape a few Alfven transit times. A strong azimuthal field component tends to inhibit this distortion. If the flux rope is taken to model a magnetic cloud, it is suggested that the shape of the cloud at 1 AU is determined by its distortion in the inner solar wind. Distortion timescales beyond 1 AU are estimated as many days. It is estimated that effective drag coefficients somewhat greater than unity are appropriate for modelling flux rope propagation.

Cargill, P. J.↗

Kinetic structure of intermediate shocks: Implications for the magnetopause

A general study of the structure and stability of intermediate shocks (IS) in an isotropic plasma is presented using a hybrid as well as a resistive Hall MHD code. Special emphasis is put on the question of whether the rotational layers observed at the magnetopause can be intermediate shocks. The shocks are formed dynamically by the interaction between a flowing plasma and a stationary piston. Coplanar ISs (both strong and weak) are found to be stable in a collisionless plasma. The existence of slow shocks in a high beta plasma is also established for the first time. Noncoplanar ISs are found to be time-dependent, evolving toward a rotational discontinuity (RD) after some characteristic time tau which can be quite long (1000 Omega(exp -1), where Omega is the ion gyrofrequency). The value tau is larger the closer the rotation angle is to 180 deg. Rotations larger than 180 deg are found to be unstable, decaying into a state of minimum shear (i.e., rotation angle less than 180 deg). There are various length scales associated with an IS in the kinetic regime. The shortest scale is found to be the length scale over which rotation of the transverse component of the magnetic field takes place. This scale can have a half width as small as one ion inertial length (c/omega(sub p)) for electron sense rotations and 3c/omega(sub p) for ion sense rotations, for an upstream ion beta of unity. Both of these scales are consistent with the observed thickness at the magnetopause and identical to the corresponding RD scales. A detailed study of the mode conversion of the Alfven ion cyclotron waves (A/IC) waves across both slow and intermediate shocks and the resulting downstream wave spectrum are presented. The possibility that the large number of relfected ions observed at the magnetopause may be due to the presence of strong ISs is considered. The identification of strong ISs and their distinction from RDs should be possible in observations due to significant differences that exist between jump conditions and overall structure of the two discontinuities. The jumps in the plasma parameters across a weak IS are typically small. This together with the fact that the weak ISs and RDs have very similar thickness and other overall properties makes the distinction between weak ISs and RDs in the observations largely inconsequential. However, at large noncoplanarity angles the weak IS approaches the RD limit in a relatively short time (approximately less than 100 Omega (exp -1)). Thus, magnetopause rotations with large noncoplanarity angles are most likely either RDs or strong ISs. Finally, direct comparisons between fluid (resistive Hall MHD) and kinetic simulations show that fluid theory is not applicable to study of ISs in a collisionless plasma.

Karimabadi, H.↗

Mercury's Magnetosphere

Among the major discoveries made by the Mariner 10 mission to the inner planets was the existence of an intrinsic magnetic field at Mercury with a dipole moment of approx. 300 nT R(sup 3, sub M). This magnetic field is sufficient to stand off the solar wind at an altitude of about 1 R(sub M) (i.e. approx. 2439 km). Hence, Mercury possesses a 'magnetosphere' from which the so]ar wind plasma is largely excluded and within which the motion of charged particles is controlled by the planetary magnetic field. Despite its small size relative to the magnetospheres of the other planets, a Mercury orbiter mission is a high priority for the space physics community. The primary reason for this great interest is that Mercury unlike all the other planets visited thus far, lacks a significant atmosphere; only a vestigial exosphere is present. This results in a unique situation where the magnetosphere interacts directly with the outer layer of the planetary crust (i.e. the regolith). At all of the other planets the topmost regions of their atmospheres become ionized by solar radiation to form ionospheres. These planetary ionospheres then couple to electrodynamically to their magnetospheres or, in the case of the weakly magnetized Venus and Mars, directly to the solar wind. This magnetosphere-ionosphere coupling is mediated largely through field-aligned currents (FACs) flowing along the magnetic field lines linking the magnetosphere and the high-latitude ionosphere. Mercury is unique in that it is expected that FACS will be very short lived due to the low electrical conductivity of the regolith. Furthermore, at the earth it has been shown that the outflow of neutral atmospheric species to great altitudes is an important source of magnetospheric plasma (following ionization) whose composition may influence subsequent magnetotail dynamics. However, the dominant source of plasma for most of the terrestrial magnetosphere is the 'leakage'of solar wind across the magnetopause and more direct entry through the northern and southern cusps. Although Mariner 10 did not return plasma composition measurements, the Hermean magnetosphere should be ideal for measuring the manner and rate of solar wind plasma entry due to the lack of strong internal atmospheric sources. Finally, the solar wind conditions experienced by Mercury as it orbits the Sun at 0.31 to 0.47 AU are quite different from those typically encountered by the Earth. This may allow for new understanding of the external factors affecting the transfer of mass, momentum and energy from the solar wind to planetary magnetospheres. This article provides a brief overview of what is now known about Mercury's magnetosphere and why it is a priority target for future planetary missions.

Slavin, J. A.↗

Wave-Particle Interactions Associated with Nongyrotropic Distribution Functions: A Hybrid Simulation Study

Nongyrotropic plasma distribution functions can be formed in regions of space where guiding center motion breaks down as a result of strongly curved and weak ambient magnetic fields. Such are the conditions near the current sheet in the Earth's middle and distant magnetotail, where observations of nongyrotropic ion distributions have been made. Here a systematic parameter study of nongyrotropic proton distributions using electromagnetic hybrid simulations is made. We model the observed nongyrotropic distributions by removing a number of arc length segments from a cold ring distribution and find significant differences with the results of simulations that initially have a gyrotropic ring distribution. Model nongyrotropic distributions with initially small perpendicular thermalization produce growing fluctuations that diffuse the ions into a stable Maxwellian-like distribution within a few proton gyro periods. The growing waves produced by nongyrotropic distributions are similar to the electromagnetic proton cyclotron waves produced by a gyrotropic proton ring distribution in that they propagate parallel to the background magnetic field and occur at frequencies on the order of the proton gyrofrequency, The maximum energy of the fluctuating magnetic field increases as the initial proton distribution is made more nongyrotropic, that is, more highly bunched in perpendicular velocity space. This increase can be as much as twice the energy produced in the gyrotropic case.

Convery, P. D.↗

Cosmology from cross-correlation of ACT-DR4 CMB lensing and DES-Y3 cosmic shear

Cross-correlation between weak lensing of the Cosmic Microwave Background (CMB) and weak lensing of galaxies offers a way to place robust constraints on cosmological and astrophysical parameters with reduced sensitivity to certain systematic effects affecting individual surveys. We measure the angular cross-power spectrum between the Atacama Cosmology Telescope (ACT) DR4 CMB lensing and the galaxy weak lensing measured by the Dark Energy Survey (DES) Y3 data. Our baseline analysis uses the CMB convergence map derived from ACT-DR4 and Planck data, where most of the contamination due to the thermal Sunyaev Zel’dovich effect is removed, thus avoiding important systematics in the cross-correlation. In our modelling, we consider the nuisance parameters of the photometric uncertainty, multiplicative shear bias and intrinsic alignment of galaxies. The resulting cross-power spectrum has a signal-to-noise ratio = 7.1 and passes a set of null tests. We use it to infer the amplitude of the fluctuations in the matter distribution (S_8 ≡ σ_8(Ω_m/0.3)^0.5 = 0.782 ± 0.059) with informative but well-motivated priors on the nuisance parameters. We also investigate the validity of these priors by significantly relaxing them and checking the consistency of the resulting posteriors, finding them consistent, albeit only with relatively weak constraints. This cross-correlation measurement will improve significantly with the new ACT-DR6 lensing map and form a key component of the joint 6×2pt analysis between DES and ACT.

79 ASTRONOMY AND ASTROPHYSICS↗

Weak magnetohydrodynamic turbulence theory revisited

Two recent papers, P. H. Yoon and G. Choe, Phys. Plasmas 28, 082306 (2021) and Yoon et al., Phys. Plasmas 29, 112303 (2022), utilized in the derivation of the kinetic equation for the intensity of turbulent fluctuations the assumption that the wave spectra are isotropic, that is, the ensemble-averaged magnetic field tensorial fluctuation intensity is given by the isotropic diagonal form, ⟨δBiδBj⟩k=⟨δB2⟩kδij. However, it is more appropriate to describe the incompressible magnetohydrodynamic turbulence involving shear Alfvénic waves by modeling the turbulence spectrum as being anisotropic. That is, the tensorial fluctuation intensity should be different in diagonal elements across and along the direction of the wave vector, ⟨δBiδBj⟩k=12 ⟨δB⊥2⟩k(δij−kikj/k2)+⟨δB∥2⟩k(kikj/k2). In the present paper, we thus reformulate the weak magnetohydrodynamic turbulence theory under the assumption of anisotropy and work out the form of nonlinear wave kinetic equation.

Physics↗

The black hole weak gravity conjecture with multiple charges

We study the effect of higher-derivative corrections on asymptotically flat, fourdimensional, dyonic black holes in low-energy models of gravity coupled to N U(1) gauge fields. For large extremal black holes, the leading O (1/Q 2 ) correction to the extremality bound is calculated from the most general low-energy effective action containing operators with up to four derivatives. Motivated by the multi-charge generalization of the Weak Gravity Conjecture, we analyze the necessary kinematic conditions for an asymptotically large extremal black hole to decay into a multi-particle state of extremal black holes. In the large black hole regime, we show that the convex hull condition degenerates to the requirement that a certain quartic form constructed from the Wilson coefficients of the fourderivative effective operators, is everywhere positive. Using on-shell unitarity methods, we show that higher-derivative operators are renormalized at one-loop only if they generate local, on-shell matrix elements that are invariant tensors of the electromagnetic duality group U(N). The one-loop logarithmic running of the four-derivative Wilson coefficients is calculated and shown to imply the positivity of the extremality form at some finite value of Q 2 . This result generalizes an argument recently given by Charles [1], and shows that under the given assumptions the multi-charge Weak Gravity Conjecture is not a Swampland criterion.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Forming electrical interconnections through semiconductor wafers

An information processing system based on CMOS/SOS technology is being developed by NASA to process digital image data collected by satellites. An array of holes is laser drilled in a semiconductor wafer, and a conductor is formed in the holes to fabricate electrical interconnections through the wafers. Six techniques are used to form conductors in the silicon-on-sapphire (SOS) wafers, including capillary wetting, wedge extrusion, wire intersection, electroless plating, electroforming, double-sided sputtering and through-hole electroplating. The respective strengths and weaknesses of these techniques are discussed and compared, with double-sided sputtering and the through-hole plating method achieving best results. In addition, hollow conductors provided by the technique are available for solder refill, providing a natural way of forming an electrically connected stack of SOS wafers.

Anthony, T. R.↗

Massive star-forming galaxies have converted most of their halo gas into stars

In the local Universe, the efficiency for converting baryonic gas into stars is very low. In dark matter halos where galaxies form and evolve, the average efficiency varies with galaxy stellar mass and has a maximum of about 20% for Milky-Way-like galaxies. The low efficiency at higher mass is believed to be the result of some quenching processes, such as the feedback from active galactic nuclei. We perform an analysis of weak lensing and satellite kinematics for SDSS central galaxies. Our results reveal that the efficiency is much higher, more than 60%, for a large population of massive star-forming galaxies around 10 11 M ⊙ . This suggests that these galaxies acquired most of the gas in their halos and converted it into stars without being significantly affected by quenching processes. This population of galaxies is not reproduced in current galaxy formation models, indicating that our understanding of galaxy formation is incomplete. The implications of our results on circumgalactic media, star-formation quenching, and disk galaxy rotation curves are discussed. We also examine systematic uncertainties in halo-mass and stellar-mass measurements that might influence our results.

79 ASTRONOMY AND ASTROPHYSICS↗

Advances in heavy alkaline earth chemistry provide insight into complexation of weakly polarizing Ra 2+ , Ba 2+ , and Sr 2+ cations

Numerous technologies—with catalytic, therapeutic, and diagnostic applications—would benefit from improved chelation strategies for heavy alkaline earth elements: Ra 2+ , Ba 2+ , and Sr 2+ . Unfortunately, chelating these metals is challenging because of their large size and weak polarizing power. We found 18-crown-6-tetracarboxylic acid (H 4 COCO) bound Ra 2+ , Ba 2+ , and Sr 2+ to form M(H x COCO) x–2 . Upon isolating radioactive 223 Ra from its parent radionuclides ( 227 Ac and 227 Th), 223 Ra 2+ reacted with the fully deprotonated COCO 4- chelator to generate Ra(COCO) 2- (aq) (log K Ra(COCO)2- = 5.97 ± 0.01), a rare example of a molecular radium complex. Comparative analyses with Sr 2+ and Ba 2+ congeners informed on what attributes engendered success in heavy alkaline earth complexation. Chelators with high negative charge [-4 for Ra(COCO) 2- (aq) ] and many donor atoms [≥11 in Ra(COCO) 2- (aq) ] provided a framework for stable complex formation. These conditions achieved steric saturation and overcame the weak polarization powers associated with these large dicationic metals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Radiant Extinction of Gaseous Diffusion Flames

The absence of buoyancy-induced flows in microgravity (mu-g) and the resulting increase in the reactant residence time significantly alters the fundamentals of many combustion processes. Substantial differences between normal gravity (ng) and mu-g flames have been reported in experiments on candle flames, flame spread over solids, droplet combustion, and others. These differences are more basic than just in the visible flame shape. Longer residence times and higher concentration of combustion products in the flame zone create a thermochemical environment that changes the flame chemistry and the heat and mass transfer processes. Processes such as flame radiation, that are often ignored in ng, become very important and sometimes even controlling. Furthermore, microgravity conditions considerably enhance flame radiation by: (1) the build-up of combustion products in the high-temperature reaction zone which increases the gas radiation; and (2) longer residence times make conditions appropriate for substantial amounts of soot to form which is also responsible for radiative heat loss. Thus, it is anticipated that radiative heat loss may eventually extinguish the "weak" (low burning rate per unit flame area) mu-g diffusion flame. Yet, space shuttle experiments on candle flames show that in an infinite ambient atmosphere, the hemispherical candle flame in mu-g will burn indefinitely. This may be because of the coupling between the fuel production rate and the flame via the heat-feedback mechanism for candle flames, flames over solids and fuel droplet flames. Thus, to focus only on the gas-phase phenomena leading to radiative extinction, aerodynamically stabilized gaseous diffusion flames are examined. This enables independent control of the fuel flow rate to help identify conditions under which radiative extinction occurs. Also, spherical geometry is chosen for the mu-g experiments and modeling because: (1) It reduces the complexity by making the problem one-dimensional; (2) The spherical diffusion flame completely encloses the soot which is formed on the fuel rich side of the reaction zone. This increases the importance of flame radiation because now both soot and gaseous combustion products co-exist inside the high temperature spherical diffusion flame; (3) For small fuel injection velocities, as is usually the case for a pyrolyzing solid, the diffusion flame in mu-g around the solid naturally develops spherical symmetry. Thus, spherical diffusion flames are of interest to fires in mu-g and identifying conditions that lead to radiation-induced extinction is important for spacecraft fire safety.

Berhan, Sean↗

Unsteady Spherical Diffusion Flames in Microgravity

The absence of buoyancy-induced flows in microgravity (mu-g) and the resulting increase in the reactant residence time significantly alters the fundamentals of many combustion processes. Substantial differences between normal gravity (ng) and (mu-g) flames have been reported in experiments on candle flames, flame spread over solids, droplet combustion, and others. These differences are more basic than just in the visible flame shape. Longer residence times and higher concentration of combustion products in the flame zone create a thermochemical environment that changes the flame chemistry and the heat and mass transfer processes. Processes such as flame radiation, that are often ignored in ng, become very important and sometimes even controlling. Furthermore, microgravity conditions considerably enhance flame radiation by: (i) the build-up of combustion products in the high-temperature reaction zone which increases the gas radiation, and (ii) longer residence times make conditions appropriate for substantial amounts of soot to form which is also responsible for radiative heat loss. Thus, it is anticipated that radiative heat loss may eventually extinguish the "weak" (low burning rate per unit flame area) mu-g diffusion flame. Yet, space shuttle experiments on candle flames show that in an infinite ambient atmosphere, the hemispherical candle flame in mu-g will burn indefinitely. This may be because of the coupling between the fuel production rate and the flame via the heat-feedback mechanism for candle flames, flames over solids and fuel droplet flames. Thus, to focus only on the gas-phase phenomena leading to radiative extinction, aerodynamically stabilized gaseous diffusion flames are examined. This enables independent control of the fuel flow rate to help identify conditions under which radiative extinction occurs. Also, spherical geometry is chosen for the mu-g experiments and modeling because: (i) It reduces the complexity by making the problem one-dimensional; (ii) The spherical diffusion flame completely encloses the soot which is formed on the fuel rich side of the reaction zone. This increases the importance of flame radiation because now both soot and gaseous combustion products co-exist inside the high temperature spherical diffusion flame. (iii) For small fuel injection velocities, as is usually the case for a pyrolyzing solid, the diffusion flame in mu-g around the solid naturally develops spherical symmetry. Thus, spherical diffusion flames are of interest to fires in mu-g and identifying conditions that lead to radiation-induced extinction is important for spacecraft fire safety.

Atreya, Arvind↗

Neutron Skins: Weak Elastic Scattering and Neutron Stars

The recently completed PREX-2 campaign measured the density distribution of neutrons in the lead nucleus as a function of momentum transfer (the form factor), confirmed a relatively large extent of the neutrons beyond the protons in the nucleus (the neutron skin), and provided a precise determination of the density of protons and neutrons at the center of a heavy nucleus. In turn, the measured form factor can be related to various nuclear and neutron star properties. The NICER X-ray telescope has inferred the masses and radii of some X-ray pulsars (neutron stars), although complications arise when determining these quantities independently. Further improvements in NICER have enabled simultaneous mass–radius determinations that had not previously been possible. During the next decade, measurements in astrophysics, gravitational-wave astronomy, and nuclear physics are expected to provide a wealth of more precise data. In this review, we present an overview of the current state of neutron skin measurements and offer insights into prospects for the future.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Heterogeneous Ice Growth in Micron-Sized Water Droplets Due to Spontaneous Freezing

Understanding how ice nucleates and grows into larger crystals is of crucial importance for many research fields. The purpose of this study was to shed light on the phase and structure of ice once a nucleus is formed inside a metastable water droplet. Wide-angle X-ray scattering (WAXS) was performed on micron-sized droplets evaporatively cooled to temperatures where homogeneous nucleation occurs. We found that for our weak hits ice grows more cubic compared to the strong hits that are completely hexagonal. Due to efficient heat removal caused by evaporation, we propose that the cubicity of ice at the vicinity of the droplet’s surface is higher than for ice formed within the bulk of the droplet. Moreover, the Bragg peaks were classified based on their geometrical shapes and positions in reciprocal space, which showed that ice grows heterogeneously with a significant population of peaks indicative of truncation rods and crystal defects. Frequent occurrences of the (100) reflection with extended in-planar structure suggested that large planar ice crystals form at the droplet surface, then fracture into smaller domains to accommodate to the curvature of the droplets. Planar faulting due to misaligned domains would explain the increased cubicity close to the droplet surface.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A statistical and simulation-informed model for estimating permeability from pore size distribution in saturated geomaterials

Accurate permeability estimation is essential across subsurface engineering applications but remains challenging due to the complex pore structures of natural geomaterials. Traditional empirical methods and simplified theoretical models often inadequately capture the role of pore size distribution and connectivity. Here, this study develops a statistical and simulation-informed permeability model that collapses pore-scale complexity into a compact scaling of the form k = αϕμ d 2 , where ϕ is porosity, μ d is mean pore size, and α is a weakly varying coefficient. By combining pore network simulations with statistical analysis of unimodal and bimodal pore size distributions, we identify three key findings: (i) permeability is much more sensitive to mean pore size than to porosity; (ii) across extensive datasets, the ratio σ d /μ d (standard deviation to mean) clusters around a characteristic value ∼0.4, allowing the effects of the full pore size distribution to be represented by μ d and a narrowly varying α ≈ 0.05; and (iii) for bimodal systems, there exists a critical fraction of small pores ∼0.78 above which flow becomes small-pore dominated, enabling the definition of an effective flow-controlling pore population and facilitating simplified permeability estimation for such systems. The resulting model, which requires only porosity and a representative mean pore size as inputs, is validated against comprehensive experimental datasets (>1700 samples) spanning diverse soils and rocks and achieves good predictive accuracy. Overall, this work provides a physically grounded yet practically simple permeability estimator suitable for subsurface engineering, environmental protection, and resource management applications.

Permeability↗