Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Parameter Optimization”

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 1,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Numerical Optimization of Density Functional Tight Binding Models: Application to Molecules Containing Carbon, Hydrogen, Nitrogen, and Oxygen

New parametrizations for semiempirical density functional tight binding (DFTB) theory have been developed by the numerical optimization of adjustable parameters to minimize errors in the atomization energy and interatomic forces with respect to ab initio calculated data. Initial guesses for the radial dependences of the Slater- Koster bond integrals and overlap integrals were obtained from minimum basis density functional theory calculations. The radial dependences of the pair potentials and the bond and overlap integrals were represented by simple analytic functions. The adjustable parameters in these functions were optimized by simulated annealing and steepest descent algorithms to minimize the value of an objective function that quantifies the error between the DFTB model and ab initio calculated data. The accuracy and transferability of the resulting DFTB models for the C, H, N, and O system were assessed by comparing the predicted atomization energies and equilibrium molecular geometries of small molecules that were not included in the training data from DFTB to ab initio data. The DFTB models provide accurate predictions of the properties of hydrocarbons and more complex molecules containing C, H, N, and O.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Toward Automated Sampling of Polymorph Nucleation and Free Energies with the SGOOP and Metadynamics

Understanding the driving forces behind the nucleation of different polymorphs is of great importance for material sciences and the pharmaceutical industry. This includes understanding the reaction coordinate that governs the nucleation process and correctly calculating the relative free energies of different polymorphs. Here, we demonstrate, for the prototypical case of urea nucleation from the melt, how one can learn such a one-dimensional reaction coordinate as a function of prespecified order parameters and use it to perform efficient biased all-atom molecular dynamics simulations. The reaction coordinate is learnt as a function of the generic thermodynamic and structural order parameters using the “spectral gap optimization of order parameters (SGOOP)” approach [Tiwary, P. and Berne, B. J. Proc. Natl. Acad. Sci. U.S.A. (2016)] and is biased using well-tempered metadynamics simulations. The reaction coordinate gives insights into the role played by different structural and thermodynamics order parameters, and the biased simulations obtain accurate relative free energies for different polymorphs. This includes an accurate prediction of the approximate pressure at which urea undergoes a phase transition and one of the metastable polymorphs becomes the most stable conformation. Here, we believe the ideas demonstrated in this work will facilitate efficient sampling of nucleation in complex, generic systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Impact of SPECT corrections on 3D-dosimetry for liver transarterial radioembolization using the patient relative calibration methodology

Purpose: Many centers aim to plan liver transarterial radioembolization (TARE) with dosimetry, even without CT-based attenuation correction (AC), or with unoptimized scatter correction (SC) methods. This work investigates the impact of presence vs absence of such corrections, and limited spatial resolution, on 3D dosimetry for TARE. Methods: Three voxelized phantoms were derived from CT images of real patients with different body sizes. Simulations of {sup 99m}Tc-SPECT projections were performed with the SIMIND code, assuming three activity distributions in the liver: uniform, inside a “liver’s segment,” or distributing multiple uptaking nodules (“nonuniform liver”), with a tumoral liver/healthy parenchyma ratio of 5:1. Projection data were reconstructed by a commercial workstation, with OSEM protocol not specifically optimized for dosimetry (spatial resolution of 12.6 mm), with/without SC (optimized, or with parameters predefined by the manufacturer; dual energy window), and with/without AC. Activity in voxels was calculated by a relative calibration, assuming identical microspheres and {sup 99m}Tc-SPECT counts spatial distribution. 3D dose distributions were calculated by convolution with {sup 90}Y voxel S-values, assuming permanent trapping of microspheres. Cumulative dose-volume histograms in lesions and healthy parenchyma from different reconstructions were compared with those obtained from the reference biodistribution (the “gold standard,” GS), assessing differences for D95%, D70%, and D50% (i.e., minimum value of the absorbed dose to a percentage of the irradiated volume). γ tool analysis with tolerance of 3%/13 mm was used to evaluate the agreement between GS and simulated cases. The influence of deep-breathing was studied, blurring the reference biodistributions with a 3D anisotropic gaussian kernel, and performing the simulations once again. Results: Differences of the dosimetric indicators were noticeable in some cases, always negative for lesions and distributed around zero for parenchyma. Application of AC and SC reduced systematically the differences for lesions by 5%–14% for a liver segment, and by 7%–12% for a nonuniform liver. For parenchyma, the data trend was less clear, but the overall range of variability passed from −10%/40% for a liver segment, and −10%/20% for a nonuniform liver, to −13%/6% in both cases. Applying AC, SC with preset parameters gave similar results to optimized SC, as confirmed by γ tool analysis. Moreover, γ analysis confirmed that solely AC and SC are not sufficient to obtain accurate 3D dose distribution. With breathing, the accuracy worsened severely for all dosimetric indicators, above all for lesions: with AC and optimized SC, −38%/−13% in liver’s segment, −61%/−40% in the nonuniform liver. For parenchyma, D50% resulted always less sensitive to breathing and sub-optimal correction methods (difference overall range: −7%/13%). Conclusions: Reconstruction protocol optimization, AC, SC, PVE and respiratory motion corrections should be implemented to obtain the best possible dosimetric accuracy. On the other side, thanks to the relative calibration, D50% inaccuracy for the healthy parenchyma from absence of AC was less than expected, while the optimization of SC was scarcely influent. The relative calibration therefore allows to perform TARE planning, basing on D50% for the healthy parenchyma, even without AC or with suboptimal corrections, rather than rely on nondosimetric methods.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1978

The Incorporation of Truncated Fourier Series into Finite Difference Approximations of Structural Stability Equations

A new trigonometric approach to the finite difference calculus was applied to the problem of beam buckling as represented by virtual work and equilibrium equations. The trigonometric functions were varied by adjusting a wavelength parameter in the approximating Fourier series. Values of the critical force obtained from the modified approach for beams with a variety of boundary conditions were compared to results using the conventional finite difference method. The trigonometric approach produced significantly more accurate approximations for the critical force than the conventional approach for a relatively wide range in values of the wavelength parameter; and the optimizing value of the wavelength parameter corresponded to the half-wavelength of the buckled mode shape. It was found from a modal analysis that the most accurate solutions are obtained when the approximating function closely represents the actual displacement function and matches the actual boundary conditions.

Hannah, S. R.↗
NASA NTRS1982

Field experience with various slicing methods

Wafer slicing using internal diameter (ID) saw, multiblade slurry (MBS) saw and multiwire slurry (MWS) saw techniques were evaluated. Wafer parameters such as bow, taper, and roughness which may not be important factors for solar cell fabrication, were considerably better for ID saw than those of the MBS and MWS saw. Analysis of add-on slicing cost indicated that machine productivity seems to be a major limiting factor for ID saw, while expendible material costs are a major factor for both MBS and MWS saw. Slicing experience indicated that the most important factors controling final wafer cost are: (1) silicon cost (wafer thickness + kerf loss); (2) add-on slicing cost, and (3) mechanical yield. There is a very strong interaction between these parameters, suggesting a necessity of optimization of these parameters.

Yoo, H. I.↗
NASA NTRS1993

Improved frame synchronization schemes for INMARSAT-B/M SCPC and TDM channels

This paper proposes faster, more robust, frame synchronization schemes for various Inmarsat-B and Inmarsat-M communication and signalling channels. Equations are developed which permit frame synchronization strategies of the type specified by Inmarsat to be evaluated in terms of average true lock time, average false maintenance time, and average search time. Evaluation of the currently specified framing schemes shows that a significant performance improvement is obtained by optimizing the threshold parameters of the scheme. The optimization seeks a compromise between the conflicting requirements of maximizing true lock time and minimizing search time.

Pan, Si-Ming↗
NASA NTRSDate not supplied

Hourly Carbon Fluxes Estimation Using the GOES Advanced Baseline Imager (ABI) Data Over the Conterminous USA

Tremendous efforts by Fluxnet scientists over the past few decades have made thousands of site-years of carbon flux observations available for advancing our understanding of carbon cycling in terrestrial ecosystems. One of key Fluxnet measurements is net ecosystem exchange (NEE) as it is directly related to carbon budget of terrestrial ecosystems. However, carbon flux estimation studies using satellite remote sensing have focused mainly on daily Gross Primary Production (GPP). The satellite based carbon flux estimation used the polar orbiting satellite sensors (e.g., Moderate Resolution Imaging Spectroradiometer (MODIS)), which allow us to observe target regions only once during the day. Because daily NEE is close to zero value, the carbon flux models using the polar orbiting satellite data have not been well used for NEE estimation. The new generation of geostationary satellite sensors (e.g., GOES Advanced Baesline Imager (ABI) and Himawari Advanced Himawari Imager (AHI)) provide frequent observations, often less than every 10 minutes. Here, we use GOES ABI data to estimate hourly NEE over the conterminous USA. We used the Terrestrial Observation Prediction System (TOPS) model for estimating hourly NEE. TOPS is a diagnostic ecosystem process model that simulates the fluxes of carbon and water through vegetation in response to climate variability. For the climate input, we developed hourly climate data using the same algorithm with NASA Earth Exchange Gridded Daily Meteorology (NEX-GDM) datasets based on machine learning techniques. The hourly climate data includes precipitation, maximum temperature, minimum temperature, dew point temperature, and, in particular, solar radiation that is directly derived from the Geostationary observations. The spatial patterns of ecosystem parameters used in TOPS are optimized using satellite Solar Induced Fluorescence (SIF) data. The high frequency GPP estimations from geostationary satellite sensors make it comparable to the instantaneous SIF data than daily GPP. We also used Fluxnet data for optimization of model parameters and the validation of the output. The derived data addresses the diurnal dynamics of carbon cycling at large scales and should help in reducing the uncertainties in carbon budget studies.

Geostationary satellite↗
NASA NTRSDate not supplied

Development of Carbon Flux Model Using ABI Data Over the Conterminous US

The satellite-driven carbon flux estimation has been playing important role to estimate continental-scale carbon budget. One of the biggest recent advances in the satellite-driven carbon flux modeling is utilization of high-frequent geostationary satellites to estimate diurnal cycle in carbon fluxes. The satellite based carbon flux estimation used the polar orbiting satellite sensors (e.g., Moderate Resolution Imaging Spectroradiometer (MODIS)), which allow us to observe target regions only once during the day. The new generation of geostationary satellite sensors provide frequent observations, often less than every 10 minutes. Here, we use GOES Advanced Baseline Imager (ABI) data to estimate hourly NEE over the conterminous US. We used the Terrestrial Observation Prediction System (TOPS) model for estimating hourly NEE. TOPS is a diagnostic ecosystem process model that simulates the fluxes of carbon and water through vegetation in response to climate variability. For the climate input, we developed hourly climate data using the same algorithm with NASA Earth Exchange Gridded Daily Meteorology (NEX-GDM) datasets based on machine learning techniques. The hourly climate data includes precipitation, maximum temperature, minimum temperature, dew point temperature, and solar radiation were derived from the Geostationary observations. The spatial patterns of ecosystem parameters used in TOPS are optimized using satellite Solar Induced Fluorescence (SIF) data. The high frequency GPP estimations from geostationary satellite sensors make it comparable to the instantaneous SIF data than daily GPP. We also used Ameriflux data for optimization of model parameters and the validation of the output. The derived data addresses the diurnal dynamics of carbon cycling at large scales and should help in reducing the uncertainties in carbon budget studies.

geostationary satellite↗
NASA NTRS1995

Direct adaptive performance optimization of subsonic transports: A periodic perturbation technique

Aircraft performance can be optimized at the flight condition by using available redundancy among actuators. Effective use of this potential allows improved performance beyond limits imposed by design compromises. Optimization based on nominal models does not result in the best performance of the actual aircraft at the actual flight condition. An adaptive algorithm for optimizing performance parameters, such as speed or fuel flow, in flight based exclusively on flight data is proposed. The algorithm is inherently insensitive to model inaccuracies and measurement noise and biases and can optimize several decision variables at the same time. An adaptive constraint controller integrated into the algorithm regulates the optimization constraints, such as altitude or speed, without requiring and prior knowledge of the autopilot design. The algorithm has a modular structure which allows easy incorporation (or removal) of optimization constraints or decision variables to the optimization problem. An important part of the contribution is the development of analytical tools enabling convergence analysis of the algorithm and the establishment of simple design rules. The fuel-flow minimization and velocity maximization modes of the algorithm are demonstrated on the NASA Dryden B-720 nonlinear flight simulator for the single- and multi-effector optimization cases.

Espana, Martin D.↗
NASA NTRS2011

Packing Optimization of Sorbent Bed Containing Dissimilar and Irregular Shaped Media

The Fire Cartridge is a packed bed air filter with two different and separate layers of media designed to provide respiratory protection from combustion products after a fire event on the International Space Station (ISS). The first layer of media is a carbon monoxide catalyst and the second layer of media is universal carbon. During development of Fire Cartridge prototypes, the two media beds were noticed to have shifted inside the cartridge. The movement of media within the cartridge can cause mixing of the bed layers, air voids, and channeling, which could cause preferential air flow and allow contaminants to pass through without removal. An optimally packed bed mitigates these risks and ensures effective removal of contaminants from the air. In order to optimally pack each layer, vertical, horizontal, and orbital agitations were investigated and a packed bulk density was calculated for each method. Packed bulk density must be calculated for each media type to accommodate variations in particle size, shape, and density. Additionally, the optimal vibration parameters must be re-evaluated for each batch of media due to variations in particle size distribution between batches. For this application it was determined that orbital vibrations achieve an optimal pack density and the two media layers can be packed by the same method. Another finding was media with a larger size distribution of particles achieve an optimal bed pack easier than media with a smaller size distribution of particles.

Holland, Nathan↗
NASA NTRS2007

Oxygen production System Models for Lunar ISRU

In-Situ Resource Utilization (ISRU) seeks to make human space exploration feasible; by using available resources from a planet or the moon to produce consumables, parts, and structures that otherwise would be brought from Earth. Producing these in situ reduces the mass of such that must be launched and doing so allows more payload mass' for each mission. The production of oxygen from lunar regolith, for life support and propellant, is one of the tasks being studied under ISRU. NASA is currently funding three processes that have shown technical merit for the production of oxygen from regolith: Molten Salt Electrolysis, Hydrogen Reduction of Ilmenite, and Carbothermal Reduction. The ISRU program is currently developing system models of, the , abovementioned processes to: (1) help NASA in the evaluation process to select the most cost-effective and efficient process for further prototype development, (2) identify key parameters, (3) optimize the oxygen production process, (4) provide estimates on energy and power requirements, mass and volume.of the system, oxygen production rate, mass of regolith required, mass of consumables, and other important parameters, and (5) integrate into the overall end-to-end ISRU system model, which could be integrated with mission architecture models. The oxygen production system model is divided into modules that represent unit operations (e.g., reactor, water electrolyzer, heat exchanger). Each module is modeled theoretically using Excel and Visual Basic for Applications (VBA), and will be validated using experimental data from on-going laboratory work. This modularity (plug-n-play) feature of each unit operation allows the use of the same model on different oxygen production systems simulations resulting in comparable results. In this presentation, preliminary results for mass, power, volume will be presented along with brief description of the oxygen production system model.

Santiago-Maldonado, Edgardo↗
DOE OSTI2026

A data-driven approach to real-time vertical position estimation for NSTX-U vertical stability control

In this paper, a database of 77 996 plasma equilibrium reconstructions from 727 discharges during the initial operation of the NSTX-U spherical tokamak is analyzed to develop a statistically robust model of the plasma vertical position for real-time control. A variety of regression models are developed and tested, ranging in complexity from linear models to deep neural networks, and including input signals ranging from the four pairs of flux loops used historically on NSTX-U up to the full set of 389 real-time signals available to the plasma control system. A linear model based on 140 real-time magnetics signals is found to offer excellent accuracy, with a coefficient of determination R 2 = 0.906. The robustness of this model to limited training data, new operating scenarios, and signal errors is tested, and a procedure is demonstrated to tune the model parameters to optimize its robustness. A time-dependent plasma equilibrium solver, TokaMaker, is used to simulate vertical stability control in NSTX-U, demonstrating that it should be possible to iteratively tune the parameters of a linear vertical position model to stabilize both positive and negative triangularity plasmas in future experiments.

magnetic diagnostics↗
NASA NTRS1994

Design and aerodynamic performance evaluation of a high-work mixed flow turbine stage

As axial and radial turbine designs have been pushed to their aerothermodynamic and mechanical limits, the mixed-flow turbine (MFT) concept has been projected to offer performance and durability improvements, especially when ceramic materials are considered. The objective of this NASA/U.S. Army sponsored mixed-flow turbine (AMFT) program was to determine the level of performance attainable with MFT technology within the mechanical constraints of 1997 projected ceramic material properties. The MFT geometry is similar to a radial turbine, exhibiting a large radius change from inlet to exit, but differing in that the inlet flowpath is not purely radial, nor axial, but mixed; it is the inlet geometry that gives rise to the name 'mixed-flow'. The 'mixed' orientation of the turbine inlet offers several advantages over radial designs by allowing a nonzero inlet blade angle yet maintaining radial-element blades. The oblique inlet not only improves the particle-impact survivability of the design, but improves the aerodynamic performance by reducing the incidence at the blade inlet. The difficulty, however, of using mixed-flow geometry lies in the scarcity of detailed data and documented design experience. This paper reports the design of a MFT stage designed with the intent to maximize aerodynamic performance by optimizing design parameters such as stage reaction, rotor incidence, flowpath shape, blade shape, vane geometry, and airfoil counts using 2-D, 3-D inviscid, and 3-D viscous computational fluid dynamics code. The aerodynamic optimization was accomplished while maintaining mechanical integrity with respect to vibration and stress levels in the rotor. A full-scale cold-flow rig test was performed with metallic hardware fabricated to the specifications of the hot ceramic geometry to evaluate the stage performance.

Neri, Remo N.↗
NASA NTRS1976

Application of optimal input synthesis to aircraft parameter identification

The Frequency Domain Input Synthesis procedure is used in identifying the stability and control derivatives of an aircraft. By using a frequency-domain approach, one can handle criteria that are not easily handled by the time-domain approaches. Numerical results are presented for optimal elevator deflections to estimate the longitudinal stability and control derivatives subject to root-mean square constraints on the input. The applicability of the steady state optimal inputs to finite duration flight testing is investigated. The steady state approximation of frequency-domain synthesis is good for data lengths greater than two time cycles for the short period mode of the aircraft longitudinal motions. Phase relationships between different frequency components become important for shorter data lengths. The frequency domain inputs are shown to be much better than the conventional doublet inputs.

Gupta, N. K.↗
DOE OSTI2023

A Comparative Study of Electrical Double Layer Effects for CO Reduction Reaction Kinetics

Solvation models describe how the interactions between the solutes and solvents affect the reactivity and selectivity in electrochemical processes. In this study, we developed a framework for evaluating the effects of applied potential and electrical double layer on CO reduction (COR), comparing fully explicit, implicit, and hybrid solvation models at the standard hydrogen electrode (SHE) scale. We analyzed all crucial intermediates leading to the production of C 1 and C 2 products and found good agreement across these models. Some notable differences were observed in the implicit description of *C and *CO adsorption at higher overpotentials and overall trends in the adsorption energies within the hybrid model. Using this unified SHE framework, we built comparable microkinetic models for COR kinetics and rates. Despite small differences in thermodynamics, solvation-model-based microkinetic simulations showed good agreement for onset potentials against the benchmark experiment. Only qualitative difference was observed for C 1+ versus hydrogen evolution at high overpotentials for the implicit model. Finally, we constructed a generalized C 2 selectivity map in descriptor space (U SHE , ΔG electrolyte CO ), which highlights the limitations of a copper-based COR catalyst and guides the search for optimal descriptor parameters to maximize C 2 selectivity. In conclusion, these findings demonstrate the importance of considering electrical double layer effects in reduction reactions and offer a useful framework for comparing solvation models and predicting optimal electrochemical conditions for specific applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2014

Optimization of the Jastrow factor using the random-phase approximation and a similarity-transformed Hamiltonian: Application to band-structure calculation for some semiconductors and insulators

Based on the random-phase approximation and the transcorrelated (TC) method, we optimize the Jastrow factor together with one-electron orbitals in the Slater determinant in the correlated wave function with a new scheme for periodic systems. The TC method is one of the promising wave function theories for first-principles electronic structure calculation, where the many-body wave function is approximated as a product of a Slater determinant and a Jastrow factor, and the Hamiltonian is similarity-transformed by the Jastrow factor. Using this similarity-transformed Hamiltonian, we can optimize the one-electron orbitals without evaluating 3N-dimensional integrations for the N-electron system. In contrast, optimization of the Jastrow factor within the framework of the TC method is computationally much more expensive and has not been performed for solid-state calculations before. In this study, we also benefit from the similarity-transformation in optimizing the Jastrow factor. Our optimization scheme is tested in applications to some solids from narrow-gap semiconductors to wide-gap insulators, and it is verified that the band gap of a wide-gap insulator and the lattice constants of some solids are improved by this optimization with reasonable computational cost.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Improved Manufacturability and Throughput of Ultra- Transparent, Super-Insulating Aerogels

AeroShield Materials produces a novel silica aerogel material with exemplary thermal performance and unprecedented optical clarity, offering the potential for super-insulating fenestration and insulated glass to reduce thermal losses in the built environment by billions of dollars every year. One of the most significant challenges facing AeroShield today is the total amount of time that is required to produce large monolithic aerogel samples. This overall process can require as much as 144 hours total, which can significantly hinder scale-up to an economically viable product. The purpose of this Phase 1 research was to continue development of our novel aerogel manufacturing process to reduce material processing time by up to 10x, greatly improving product throughput and reducing cost. AeroShield’s manufacturing process can be divided into 4 major stages, each with their own distinct set of parameters and time requirements. Under this Phase 1 award, AeroShield was able to identify and optimize a number of these parameters, including molding materials, solvent rinse conditions, critical point drying time, and annealing conditions. AeroShield also performed thorough analyses on how these optimized parameters affected important final characteristics of the gels, including optical clarity, thermal conductivity, and dimensional stability. This campaign culminated in the production of laboratory scale aerogel samples using significantly lower process times of both 36 and 18 total hours, which represent Phase 1 Target and Stretch goals. In order to achieve widespread market adoption, monolithic sheets of the aerogel material must be made to industry-standard sizes (8’ x 12’) at a cost that provides 5-7 year or less breakeven energy savings for consumers (<$2 sq/ft). The work performed under this Phase 1 award shows that time and materials required to make aerogel samples, which make up a significant portion of their overall cost, can be greatly reduced without sacrificing quality. AeroShield plans to use these optimized processes to produce larger, product-relevant sized aerogels in order to achieve target final material costs.

Wilke, Kyle↗
DOE OSTI2019

Bayesian and LASSO Regressions for Comparative Permeability Modeling of Sandstone Reservoirs

In this research, Bayesian model averaging (BMA) and least absolute shrinkage and selection operator regression (LASSO) algorithms were adopted for the core permeability modeling as a function of well log and core measurements. More specifically, the core permeability (dependent factor) was modeled given the well log and core data (independent variables) and then predicted in non-cored intervals of a well in a sandstone formation. The BMA is a stochastic linear modeling and a Bayesian parameter selection. Among 50 linear models generated as a function of the independent variables and depicted in Occam’s window, the best model of the highest posterior probability is determined. The best model has the optimal subset parameters that most influences the response factor. The core permeability modeling was again conducted using the LASSO algorithm, which adopts a completely different way of modeling and subset selection by using the penalized least-squared equation. Both BMA and LASSO resulted in a very accurate prediction of the core permeability by achieving perfect measured and calculated permeability matching. Matches for the two algorithms were quantified by computing the adjusted R-squared and root-mean-squared prediction error. In addition, results of the BMA and LASSO were compared to the conventional multiple linear regression (LM). The LASSO algorithm led to an accurate matching similar to LM, but slightly better than BMA. Therefore, both BMA and LASSO represent integrated procedures for accurate permeability estimation similar to conventional regression analysis.

58 GEOSCIENCES↗