Search NASA⌕ Search

SEARCH · Search NASA

Results for “Exponential Matrix”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Gold Nanoparticles (AuNPs) as Reactive Matrix for Detection of Trace Levels of HCN in Air by Laser Desorption/Ionization Mass Spectrometry (LDI-MS)

Under direct laser desorption/ionization mass spectrometric conditions, the irradiation of target spots made of gold nanoparticle residues generates a series of peaks at m/z 197, 394, 591… representing Au{sub n}{sup −} ions (n = 1–3). In contrast, spectra recorded from gold nanoparticles directly mixed with an alkali cyanide exhibited an additional peak at m/z 249, indicating an abundant generation of gaseous [Au(CN){sub 2}]{sup −} ions upon irradiation. The relative intensity of the m/z 249 peak surged when the amount of cyanide in the mixture was increased. Most remarkably, a peak at m/z 249 was observed even from neat AuNPs upon irradiation, if a nearby spot, which was not irradiated, happened to bear a cyanide sample. We postulated that traces of HCN emanating from the headspace of aqueous cyanide solution during the sample-plate preparation is sufficient to convert gold to AuCN, which is subsequently detected as [Au(CN){sub 2}]{sup −}. Further experiments demonstrated that the relative intensity of the m/z 249 peak diminishes exponentially as the AuNP spot becomes more distant from the putative HCN source. Eventually, the method was developed as an efficient procedure to detect HCN or alkali cyanides. Using KCN, the detection limits were determined to be below 10 pg of CN{sup −} per spot. The method also demonstrated that, upon crushing, the seeds or roots of certain fruits and vegetables such as apple, peach, radish, and cassava, but not carrot, release HCN in amounts detectable by this method. .

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2014

Enhanced exchange bias effect in size modulated Sm{sub 0.5}Ca{sub 0.5}MnO{sub 3} phase separated manganite

The effect of grain size modulation on exchange bias effect in CE-type antiferromagnetic Sm{sub 0.5}Ca{sub 0.5}MnO{sub 3} phase separated manganite is reported here. With the reduction of particle size, ferromagnetic clusters are found to form in the charge ordered antiferromagnetic matrix and gradually become larger. The horizontal and vertical shifts of the magnetic hysteresis loops in the field cooled magnetization process clearly indicate the size dependent exchange bias effect and it can be tuned with the reduction of particle sizes. The values of exchange bias parameter, i.e., exchange bias field (H{sub E}), coercivity (H{sub C}), remanence asymmetry (M{sub E}), and magnetic coercivity (M{sub C}) are found to depend strongly on the particle size. The variations of H{sub E} follow non-monotonic dependencies with reduction in particle size and show maximum (1205 Oe) at particle size of 150 nm at T = 5 K, which can be ascribed due to the changes in uncompensated surface spins. The values of H{sub E} and M{sub E} are found to decrease exponentially with increasing temperature below the spin- or cluster-glass like freezing temperature. The spin relaxation model has been employed for analysis of large magnetic training effect. The linear relationship between H{sub E} and M{sub E} further confirms the role of uncompensated surface spins. In view of spintronics application of manganites, the present observation of large exchange bias shift in this half-doped manganite may have great technological importance.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2004

Random Matrix Approach to Quantum Adiabatic Evolution Algorithms

We analyze the power of quantum adiabatic evolution algorithms (Q-QA) for solving random NP-hard optimization problems within a theoretical framework based on the random matrix theory (RMT). We present two types of the driven RMT models. In the first model, the driving Hamiltonian is represented by Brownian motion in the matrix space. We use the Brownian motion model to obtain a description of multiple avoided crossing phenomena. We show that the failure mechanism of the QAA is due to the interaction of the ground state with the "cloud" formed by all the excited states, confirming that in the driven RMT models. the Landau-Zener mechanism of dissipation is not important. We show that the QAEA has a finite probability of success in a certain range of parameters. implying the polynomial complexity of the algorithm. The second model corresponds to the standard QAEA with the problem Hamiltonian taken from the Gaussian Unitary RMT ensemble (GUE). We show that the level dynamics in this model can be mapped onto the dynamics in the Brownian motion model. However, the driven RMT model always leads to the exponential complexity of the algorithm due to the presence of the long-range intertemporal correlations of the eigenvalues. Our results indicate that the weakness of effective transitions is the leading effect that can make the Markovian type QAEA successful.

Boulatov, Alexei↗
DOE OSTI2024

Two-fermion negativity and confinement in the Schwinger model

We consider the fermionic (logarithmic) negativity between two fermionic modes in the Schwinger model. Recent results pointed out that fermionic systems can exhibit stronger entanglement than bosonic systems, exhibiting a negativity that decays only algebraically. The Schwinger model is described by fermionic excitations at short distances, while its asymptotic spectrum is the one of a bosonic theory. We show that the two-mode negativity detects this confining, fermion-to-boson transition, shifting from an algebraic decay to an exponential decay at distances of the order of the de Broglie wavelength of the first excited state. We derive analytical expressions in the massless Schwinger model and confront them with tensor network simulations. We also perform tensor network simulations in the massive model, which is not solvable analytically, and close to the Ising quantum critical point of the Schwinger model, where we show that the negativity behaves as its bosonic counterpart. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Scalable Circuits for Preparing Ground States on Digital Quantum Computers: The Schwinger Model Vacuum on 100 Qubits

The vacuum of the lattice Schwinger model is prepared on up to 100 qubits of IBM’s Eagle-processor quantum computers. A new algorithm to prepare the ground state of a gapped translationally invariant system on a quantum computer is presented, which we call “scalable circuits ADAPT-VQE” (SC-ADAPT-VQE). This algorithm uses the exponential decay of correlations between distant regions of the ground state, together with ADAPT-VQE, to construct quantum circuits for state preparation that can be scaled to arbitrarily large systems. These scalable circuits can be determined with use of classical computers, avoiding the challenging task of optimizing parameterized circuits on a quantum computer. SC-ADAPT-VQE is applied to the Schwinger model, and is shown to be systematically improvable, with an accuracy that converges exponentially with circuit depth. Both the structure of the circuits and the deviations of prepared wave functions are found to become independent of the number of spatial sites, L . This allows a controlled extrapolation of the circuits, determined with use of small or modest-sized systems, to arbitrarily large L . The circuits for the Schwinger model are determined on lattices up to L = 14 (28 qubits) with the Qiskit classical simulator, and are subsequently scaled up to prepare the L = 50 (100 qubits) vacuum on IBM’s 127-superconducting-qubit quantum computers ibm_brisbane and ibm_cusco. After introduction of an improved error-mitigation technique, which we call “operator decoherence renormalization”, the chiral condensate and charge-charge correlators obtained from the quantum computers are found to be in good agreement with classical matrix product state simulations. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Selection and characterization of ssDNA aptamers specifically recognizing pathogenic Vibrio alginolyticus

Abstract Vibrio alginolyticus ( V. alginolyticus ) is a major opportunistic pathogen to both marine animals and humans, which has also caused heavy economic losses to mariculture. The aim of this study was to develop highly specific aptamers for V. alginolyticus . Single‐stranded DNA (ssDNA) aptamers with high binding affinity to viable V. alginolyticus were generated by Systematic Evolution of Ligands by Exponential Enrichment (SELEX) and identified by flow cytometric analysis in this study. The selected aptamers showed high specificity for V. alginolyticus and low apparent binding for other bacteria. The aptamers formed distinct stem‐loop structures, which could form the basis of aptamers’ specific binding to the target V. alginolyticus . Aptamer VA2 and VA8 showed particularly high binding affinity constant (K d ) of 14.31 ± 4.26 and 90.00 ± 13.51 nM, respectively. The aptamers produced no cytotoxic effects in vitro and in vivo. ssDNA aptamers were successfully selected against the viable bacteria pathogen V. alginolyticus by SELEX. The aptamers selected in this study could be not only applied as specific chemical molecular probes for studying V. alginolyticus pathogenesis to Trachinotus ovatus , but also developing rapid convenient diagnosis assay for V. alginolyticus infection, even when applied to the complex sample matrix, such as food and environment samples.

Yu, Qing↗
DOE OSTI2024

Microbunch rotation in an x-ray free-electron laser using a first-order achromatic bend

Electrons in an x-ray free electron laser (XFEL) develop periodic density fluctuations, known as microbunches, which enable the exponential gain of x-ray power in an XFEL. When an electron beam microbunched at a hard x-ray wavelength is kicked, microbunches are often washed out due to the dispersion and R 56 of the bend. An achromatic (dispersion-free) bend with a small R 56 , however, can preserve microbunches, which rotate to follow the new trajectory of the electron bunch. Rotated microbunches can subsequently interact in a repointed undulator to produce a new beam of off-axis x rays. In this work, we demonstrate hard x-ray multiplexing in the Linac Coherent Light Source hard x-ray undulator line using microbunch rotation through a 10 μ rad first-order-achromatic bend created by transversely offsetting quadrupole magnets in the FODO lattice. Quadrupole offsets are determined analytically from beam-matrix theory. We also discuss the application of microbunch rotation to out-coupling a cavity-based XFEL. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗
DOE OSTI2014

Enhancing phonon flow through one-dimensional interfaces by impedance matching

We extend concepts from microwave engineering to thermal interfaces and explore the principles of impedance matching in 1D. The extension is based on the generalization of acoustic impedance to nonlinear dispersions using the contact broadening matrix Γ(ω), extracted from the phonon self energy. For a single junction, we find that for coherent and incoherent phonons, the optimal thermal conductance occurs when the matching Γ(ω) equals the Geometric Mean of the contact broadenings. This criterion favors the transmission of both low and high frequency phonons by requiring that (1) the low frequency acoustic impedance of the junction matches that of the two contacts by minimizing the sum of interfacial resistances and (2) the cut-off frequency is near the minimum of the two contacts, thereby reducing the spillage of the states into the tunneling regime. For an ultimately scaled single atom/spring junction, the matching criterion transforms to the arithmetic mean for mass and the harmonic mean for spring constant. The matching can be further improved using a composite graded junction with an exponential varying broadening that functions like a broadband antireflection coating. There is, however, a trade off as the increased length of the interface brings in additional intrinsic sources of scattering.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1985

Models of the upper atmosphere

Several computer models of the thermosphere are discussed. J70MM, the Jacchia 1970 model with matrix and matrix mean output, and J703X, the Jacchia 1970, 1971, and 1977 models, are of primary interest. The subroutines in these programs were studied in detail, and several mistakes were found and corrected. It is proposed that a weighted average of the three-hour geomagnetic index be used in the models rather than a single index value. The densities (for a given date and time) generated by the J703X program (all three models), using both a single value of the index and weighted averages for three different lag times, are tabulated and discussed. Most of the equations used in the models are the results of empirical curve fits. An attempt was made to generate a theoretical prediction of the thermospheric temperature profile based on a simple physical model of atmospheric heat conduction in a spherically-symmetric shell. The exponential heating theory and the temperature dependence of k are discussed.

Davis, T. M.↗
DOE OSTI2023

Effective description of sub-maximal chaos: stringy effects for SYK scrambling

It has been proposed that the exponential decay and subsequent power law saturation of out-of-time-order correlation functions can be universally described by collective ‘scramblon’ modes. We develop this idea from a path integral perspective in several examples, thereby establishing a general formalism. After reformulating previous work on the Schwarzian theory and identity conformal blocks in two-dimensional CFTs relevant for systems in the infinite coupling limit with maximal quantum Lyapunov exponent, we focus on theories with sub-maximal chaos: we study the large-q limit of the SYK quantum dot and chain, both of which are amenable to analytical treatment at finite coupling. In both cases we identify the relevant scramblon modes, derive their effective action, and find bilocal vertex functions, thus constructing an effective description of chaos. The final results can be matched in detail to stringy corrections to the gravitational eikonal S-matrix in holographic CFTs, including a stringy Regge trajectory, bulk to boundary propagators, and multi-string effects that are unexplored holographically.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2011

Quantum state tomography with fully symmetric measurements and product measurements

We introduce random-matrix theory to study the tomographic efficiency of a wide class of measurements constructed out of weighted 2-designs, including symmetric informationally complete (SIC) probability operator measurements (POMs). In particular, we derive analytic formulas for the mean Hilbert-Schmidt distance and the mean trace distance between the estimator and the true state, which clearly show the difference between the scaling behaviors of the two error measures with the dimension of the Hilbert space. We then prove that the product SIC POMs, the multipartite analog of the SIC POMs, are optimal among all product measurements in the same sense as the SIC POMs are optimal among all joint measurements. We further show that, for bipartite systems, there is only a marginal efficiency advantage of the joint SIC POMs over the product SIC POMs. In marked contrast, for multipartite systems, the efficiency advantage of the joint SIC POMs increases exponentially with the number of parties.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRSDate not supplied

The System Complexity Metric (SCM) Explains Systems Design and is Correlated with Cost and Failure Rate

The human short term memory span and working capacity is limited to three to five items, especially if they are organized complex “chunks” of information. The impression of complexity occurs when a system is simply difficult to understand, where there is no apparent pattern to predict its behavior. Hierarchical systems design can reduce perceived complexity and increase the amount of information that can be managed. The SCM was developed to measure complexity and help compare proposed overall system architectures before detailed design information is available. The SCM is defined as the sum of the number of major nodes, N, in the system block diagram plus the number of one-way interactions, I, between the nodes. SCM = N + I. SCM’s are easily determined by direct inspection of high-level block diagrams of life support systems. Axiomatic design develops a hierarchy of subsystem requirements and designs together in a top-down, back-and-forth process. A coupling matrix is used to control the relationships between the subsystem functions and design concepts. Axiomatic design can improve system design by decoupling requirements and designs. Axiomatic design was applied to the planning of a closed life support system, similar to that used on the International Space Station. A materially open as opposed to a closed system design was created by removing the interconnections required to close the system. The open system had the same number of designed subsystems as the closed system, but it had many fewer interconnections and its SCM was lower by about half. The costs were estimated and the MTBF (Mean Time Before Failure) tabulated for open and closed space life support systems. The estimated costs were linearly proportional to SCM for the wide variations of SCM in life support, but small differences may not be significant. The flight and preflight MTBF’s both declined exponentially with increasing MTBF, faster than MTBF-2, even though the preflight estimated MTBF’s were about ten times higher than the flight MTBF’s.

System Complexity Metric (SCM)↗
DOE OSTI2011

Numerical and Experimental Study of Mixing Processes Associated with Hydrogen and High Hydrogen Content Fuels

As simulation capability improves exponentially with increasingly more cost effective CPUs and hardware, it can be used “routinely” for engineering applications. Many commercial products are available and they are marketed as increasingly powerful and easy to use. The question remains as to the overall accuracy of results obtained. To support the validation of the CFD, a hierarchical experiment was established in which the type of fuel injection (radial, axial) as well as level of swirl (non-swirling, swirling) could be systematically varied. The effort was limited to time efficient approaches (i.e., generally RANS approaches) although limited assessment of time resolved methods (i.e., unsteady RANS and LES) were considered. Careful measurements of the flowfield velocity and fuel concentration were made using both intrusive and non-intrusive methods. This database was then used as the basis for the assessment of the CFD approach. The numerical studies were carried out with a statistically based matrix. As a result, the effect of turbulence model, fuel type, axial plane, turbulent Schmidt number, and injection type could be studied using analysis of variance. The results for the non-swirling cases could be analyzed as planned, and demonstrate that turbulence model selection, turbulence Schmidt number, and the type of injection will strongly influence the agreement with measured values. Interestingly, the type of fuel used (either hydrogen or methane) has no influence on the accuracy of the simulations. For axial injection, the selection of proper turbulence Schmidt number is important, whereas for radial injection, the results are relatively insensitive to this parameter. In general, it was found that the nature of the flowfield influences the performance of the predictions. This result implies that it is difficult to establish a priori the “best” simulation approach to use. However, the insights from the relative orientation of the jet and flow do offer some guidance for which approach to take. Overall, the results underscore the importance of model “anchoring” (i.e., “tuning” the model to provide “reasonable” agreement with a well characterized geometry/flow). Finally, the results obtained have been carefully compiled into a standalone database following a standard format that is contained in an Appendix. This database is thus available for use by others for CFD modeling evaluations.

08 HYDROGEN↗
DOE OSTI2021

A theoretical approach for transient shock strengthening in high-energy-density laser compression experiments

In high-energy-density shock compression experiments, the desired state of compression is typically achieved by passing shock waves through a sequence of different materials. Here, a theoretical approach for transiently strengthening such shocks passing through interfaces of experimentally relevant impedance ratios is examined. A semi-analytical method based on characteristics analysis is used to solve the problem of a shock passing from one material to another through an intermediate region of non-uniform impedance; by appropriately designing this region, a greater shock strength can be achieved in the second material for a finite duration than in the absence of this region. When a shock passes into a material of higher impedance, the shock strengthening increases to a maximum before decreasing to an asymptotic value as the strength of the incident shock is increased. For shocks passing into materials of decreasing impedance, the shock strengthening increases monotonically as the strength of the incident shock is increased when the impedance ratio is above a critical threshold. Incorporating multiple intermediate materials can further increase the strength of the transmitted shock, with an exponential discretization of intermediate material impedances being the most effective distribution for strengthening strong shocks. The results suggest that up to 25% and 9% increases in pressure behind the leading shock can be achieved for materials of increasing and decreasing impedance, respectively. The technique is applied to the design of laser-driven dynamic compression experiments, and the results of the analysis are verified via comparison to simulations performed with the HYADES hydrodynamics code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Grain detachment and transport clogging during mineral dissolution in carbonate rocks with permeable grain boundaries

In this paper, we explore the effect of permeable grain boundaries on carbonate rock dissolution in order to improve our understanding of grain detachment and migration during reactive flow. To do so we investigated the effects of grain size and enhanced reactivity along grain boundaries on global dissolution kinetics. Variations of permeability and porosity were calculated, and exponential relationships were observed in fractured rocks. Our model employed a reactive transport framework based on the Darcy-Brinkman equation to simulate calcite dissolution in carbonate rocks composed of microporous grains and, included fluid transport along grain boundaries. The model includes fluid flow, solute transport by advection–diffusion, heterogeneous reaction between fluid and minerals and grain detachment with subsequent grain transport in macropores. The migration of solid particles due to dissolution was based on cluster analysis and local movement, and the results show that grain detachment can lead to significant decreases in permeability due to the clogging of transport pathways. Microporous media with smaller grain sizes (fine grains in this study) showed a higher average reaction rate than those with coarser grains. In addition, as the overall rates of geochemical processes are commonly affected by the presence of texture heterogeneities such as fractures, a single fracture (macropore) introduced into a microporous matrix composed of permeable grains and grain boundaries was modelled as a large channel connecting the inlet and outlet of the simulation domain. It was found that macropore clogging by grain detachment temporary decreased permeability, lowering the long-term global reaction rates. Lastly, we observed that increases in flow rate can reduce detachment and fracture clogging by reducing local dissolution along grain boundaries.

58 GEOSCIENCES↗
DOE OSTI2018

Spectroscopic and Microscopic Evidence of Biomediated HgS Species Formation from Hg(II)–Cysteine Complexes: Implications for Hg(II) Bioavailability

In this work we investigated the chemistry of Hg(II) during exposure of exponentially growing bacteria ( Escherichia coli , Bacillus subtilis , and Geobacter sulfurreducens ) to 50 nM, 500 nM, and 5 μM total Hg(II) with and without added cysteine. With X-ray absorption spectroscopy, we provide direct evidence of the formation of cell-associated HgS for all tested bacteria. The addition of cysteine (100–1000 μM) promotes HgS formation (>70% of total cell-associated Hg(II)) as a result of the biodegradation of added cysteine to sulfide. Cell-associated HgS species are also detected when cysteine is not added as a sulfide source. Two phases of HgS, cinnabar (α-HgS) and metacinnabar (β-HgS), form depending on the total concentration of Hg(II) and sulfide in the exposure medium. However, α-HgS exclusively forms in assays that contain an excess of cysteine. Scanning transmission electron microscopy images reveal that nanoparticulate HgS (s) is primarily located at the cell surface/extracellular matrix of Gram-negative E. coli and G. sulfurreducens and in the cytoplasm/cell membrane of Gram-positive B. subtilis . Intracellular Hg(II) was detected even when the predominant cell-associated species was HgS. This study shows that HgS species can form from exogenous thiol-containing ligands and endogenous sulfide in Hg(II) biouptake assays under nondissimilatory sulfate reducing conditions, providing new considerations for the interpretation of Hg(II) biouptake results.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2011

Al-Ca and Al-Fe metal-metal composite strength, conductivity, and microstructure relationships

Deformation processed metal-metal composites (DMMC’s) are composites formed by mechanical working (i.e., rolling, swaging, or wire drawing) of two-phase, ductile metal mixtures. Since both the matrix and reinforcing phase are ductile metals, the composites can be heavily deformed to reduce the thickness and spacing of the two phases. Recent studies have shown that heavily drawn DMMCs can achieve anomalously high strength and outstanding combinations of strength and conductivity. In this study, Al-Fe wire composite with 0.07, 0.1, and 0.2 volume fractions of Fe filaments and Al-Ca wire composite with 0.03, 0.06, and 0.09 volume fractions of Ca filaments were produced in situ, and their mechanical properties were measured as a function of deformation true strain. The Al-Fe composites displayed limited deformation of the Fe phase even at high true strains, resulting in little strengthening effect in those composites. Al-9vol%Ca wire was deformed to a deformation true strain of 13.76. The resulting Ca second-phase filaments were deformed to thicknesses on the order of one micrometer. The ultimate tensile strength increased exponentially with increasing deformation true strain, reaching a value of 197 MPa at a true strain of 13.76. This value is 2.5 times higher than the value predicted by the rule of mixtures. A quantitative relationship between UTS and deformation true strain was determined. X-ray diffraction data on transformation of Al + Ca microstructures to Al + various Al-Ca intermetallic compounds were obtained at the Advanced Photon Source at Argonne National Laboratory. Electrical conductivity was measured over a range of true strains and post-deformation heat treatment schedules.

36 MATERIALS SCIENCE↗
NASA NTRS2003

Flutter and Divergence Analysis using the Generalized Aeroelastic Analysis Method

The Generalized Aeroelastic Analysis Method (GAAM) is applied to the analysis of three well-studied checkcases: restrained and unrestrained airfoil models, and a wing model. An eigenvalue iteration procedure is used for converging upon roots of the complex stability matrix. For the airfoil models, exact root loci are given which clearly illustrate the nature of the flutter and divergence instabilities. The singularities involved are enumerated, including an additional pole at the origin for the unrestrained airfoil case and the emergence of an additional pole on the positive real axis at the divergence speed for the restrained airfoil case. Inconsistencies and differences among published aeroelastic root loci and the new, exact results are discussed and resolved. The generalization of a Doublet Lattice Method computer code is described and the code is applied to the calculation of root loci for the wing model for incompressible and for subsonic flow conditions. The error introduced in the reduction of the singular integral equation underlying the unsteady lifting surface theory to a linear algebraic equation is discussed. Acknowledging this inherent error, the solutions of the algebraic equation by GAAM are termed 'exact.' The singularities of the problem are discussed and exponential series approximations used in the evaluation of the kernel function shown to introduce a dense collection of poles and zeroes on the negative real axis. Again, inconsistencies and differences among published aeroelastic root loci and the new 'exact' results are discussed and resolved. In all cases, aeroelastic flutter and divergence speeds and frequencies are in good agreement with published results. The GAAM solution procedure allows complete control over Mach number, velocity, density, and complex frequency. Thus all points on the computed root loci can be matched-point, consistent solutions without recourse to complex mode tracking logic or dataset interpolation, as in the k and p-k solution methods.

Edwards, John W.↗