Large Area Crop Inventory Experiment (LACIE). Detailed description of the wheat acreage estimation procedure used in the Large Area Crop Inventory Experiment
There are no author-identified significant results in this report.
Austin, W. W.↗
SEARCH · Search NASA
Results for “large”
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.
Large deflection, large amplitude vibrations and random response of symmetrically laminated plates
Explore the source record for details and available documents.
Gray, C. E., Jr.↗
Improved large perturbation propulsion models for control system design (1988-1989) and large perturbation models of high velocity propulsion systems (1989-1990) and reduced order propulsion models for control system design (1990-1991)
Methods for modeling high speed propulsion systems will be discussed. Included in this category are internal flow propulsion systems without rotating machinery, such as inlets, ramjets, and scramjets. Among the modeling topics discussed are modeling of linear isentropic flow, heat exchange, gasdynamics, lumped parameter systems, and infinite dimensional systems. Furthermore, a generalized overview of modeling high speed propulsion systems is presented in this collection of papers.
Hartley, Tom T.↗
The distance to the Large Magellanic Cloud - Constraints from Cepheids in Large Magellanic Cloud star clusters
We have used recent observational data for the Cepheids in the rich, young LMC clusters NGC 1866 and NGC 2031 to constrain the cluster distances with the mass-equivalency (ME) method. The basis of this approach is to fix the cluster distances by requiring the Cepheid evolutionary and pulsational masses to be equal. Using evolutionary models incorporating a mild amount of core and envelope overshooting along with recent pulsational models, we derive distance moduli of 18.51 +/- 0.21 and 18.32 +/- 0.20 for NGC 1866 and NGC 2031, respectively. The quoted errors are dominated by the uncertainties in the heavy element abundances of the clusters (assumed to be 0.3 dex for both clusters), with a smaller contribution due to the apparently intrinsic spread in the masses of the Cepheids in each cluster. For the ME method, we find that Delta(m-M)(0)/Delta sub Z(1) = 0.69, where Z(1) = log (Z/0.016). This result implies that the cluster distances can be determined to better than +/- 5 percent if the cluster abundances can be measured to better than about +/- 0.15 dex.
Bertelli, Gianpaolo↗
An investigation of a wing-propeller configuration employing large-chord plain flaps and large-diameter propellers for low-speed flight and vertical take-off
Explore the source record for details and available documents.
Kuhn, Richard E↗
Supercooled Large Droplet Deformation and Breakup Near the Leading Edge of Large Transport Airfoils
No abstract available
Vargas, Mario↗
Occultation of a Large Star by the Large Plutino (28978) Ixion on 2020 October 13 UTC
Explore the source record for details and available documents.
Stephen E. Levine↗
Metrology of Large Parts
As discussed in the first chapter of this book, there are many different methods to measure a part using optical technology. Chapter 2 discussed the use of machine vision to measure macroscopic features such as length and position, which was extended to the use of interferometry as a linear measurement tool in chapter 3, and laser or other trackers to find the relation of key points on large parts in chapter 4. This chapter looks at measuring large parts to optical tolerances in the sub-micron range using interferometry, ranging, and optical tools discussed in the previous chapters. The purpose of this chapter is not to discuss specific metrology tools (such as interferometers or gauges), but to describe a systems engineering approach to testing large parts. Issues such as material warpage and temperature drifts that may be insignificant when measuring a part to micron levels under a microscope, as will be discussed in later chapters, can prove to be very important when making the same measurement over a larger part. In this chapter, we will define a set of guiding principles for successfully overcoming these challenges and illustrate the application of these principles with real world examples. While these examples are drawn from specific large optical testing applications, they inform the problems associated with testing any large part to optical tolerances. Manufacturing today relies on micrometer level part performance. Fields such as energy and transportation are demanding higher tolerances to provide increased efficiencies and fuel savings. By looking at how the optics industry approaches sub-micrometer metrology, one can gain a better understanding of the metrology challenges for any larger part specified to micrometer tolerances. Testing large parts, whether optical components or precision structures, to optical tolerances is just like testing small parts, only harder. Identical with what one does for small parts, a metrologist tests large parts and optics in particular to quantify their mechanical properties (such as dimensions, mass, etc); their optical prescription or design (i.e. radius of curvature, conic constant, vertex location, size); and their full part shape. And, just as with small parts, a metrologist accomplishes these tests using distance measuring instruments such as tape measures, inside micrometers, coordinate measuring machines, distance measuring interferometers; angle measuring instruments such as theodolites, autocollimators; and surface measuring instruments including interferometers, stylus profilers, interference microscopes, photogrammetric cameras, or other tools. However, while the methodology may be similar, it is more difficult to test a large object for the simple reason that most metrologists do not have the necessary intuition. The skills used to test small parts or optics in a laboratory do not extrapolate to testing large parts in an industrial setting any more than a backyard gardener might successfully operate a farm. But first, what is a large part? A simple definition might be the part's size or diameter. For optics and diffuse surface parts alike, the driving constraint is ability to illuminate the part's surface. For reflective convex mirrors, large is typically anything greater than 1 meter. But, for refractive optics, flats or convex mirrors, large is typically greater than 0.5 meter. While a size definition is simple, it may be less than universal. A more nuanced definition might be that a large part is any component which cannot be easily tested in a standard laboratory environment, on a standard vibration isolated table using standard laboratory infrastructure. A micro-switch or a precision lens might be easily measured to nanometer levels under a microscope in a lab, but a power turbine spline or a larger telescope mirror will not fit under that microscope and may not even fit on the table.
Stahl, H. Philip↗
Large Payload Ground Transportation and Test Considerations
Many spacecraft concepts under consideration by the National Aeronautics and Space Administration’s (NASA’s) Evolvable Mars Campaign take advantage of a Space Launch System payload shroud that may be 8 to 10 meters in diameter. Large payloads can theoretically save cost by reducing the number of launches needed--but only if it is possible to build, test, and transport a large payload to the launch site in the first place. Analysis performed previously for the Altair project identified several transportation and test issues with an 8.973 meters diameter payload. Although the entire Constellation Program—including Altair—has since been canceled, these issues serve as important lessons learned for spacecraft designers and program managers considering large payloads for future programs. A transportation feasibility study found that, even broken up into an Ascent and Descent Module, the Altair spacecraft would not fit inside available aircraft. Ground transportation of such large payloads over extended distances is not generally permitted, so overland transportation alone would not be an option. Limited ground transportation to the nearest waterway may be possible, but water transportation could take as long as 67 days per production unit, depending on point of origin and acceptance test facility; transportation from the western United States would require transit through the Panama Canal to access the Kennedy Space Center launch site. Large payloads also pose acceptance test and ground processing challenges. Although propulsion, mechanical vibration, and reverberant acoustic test facilities at NASA’s Plum Brook Station have been designed to accommodate large spacecraft, special handling and test work-arounds may be necessary, which could increase cost, schedule, and technical risk. Once at the launch site, there are no facilities currently capable of accommodating the combination of large payload size and hazardous processing such as hypergolic fuels, pyrotechnic devices, and high pressure gasses. Ironically, the limiting factor to a national heavy lift strategy may not be the rocket technology needed to throw a heavy payload, but rather the terrestrial infrastructure—roads, bridges, airframes, and buildings—necessary to transport, acceptance test, and process large spacecraft. Failure to carefully consider where and how large spacecraft are manufactured, tested, and launched could result in unforeseen cost to modify existing (or develop new) infrastructure, or incur additional risk due to increased handling operations or eliminating key verifications. Although this paper focuses on the canceled Altair spacecraft as a case study, the issues identified here have wide applicability to other large payloads, including concepts under consideration for NASA’s Evolvable Mars Campaign.
Rucker, Michelle A.↗
Small vs. Large Convective Cloud Objects from CERES Aqua Observations: Where are the Intraseasonal Variation Signals?
During inactive phases of Madden-Julian oscillation (MJO), there are plenty of deep but small convective systems and far fewer deep and large ones. During active phases of MJO, a manifestation of an increase in the occurrence of large and deep cloud clusters results from an amplification of large-scale motions by stronger convective heating. This study is designed to quantitatively examine the roles of small and large cloud clusters during the MJO life cycle. We analyze the cloud object data from Aqua CERES observations for tropical deep convective (DC) and cirrostratus (CS) cloud object types according to the real-time multivariate MJO index. The cloud object is a contiguous region of the earth with a single dominant cloud-system type. The size distributions, defined as the footprint numbers as a function of cloud object diameters, for particular MJO phases depart greatly from the combined (8-phase) distribution at large cloud-object diameters due to the reduced/increased numbers of cloud objects related to changes in the large-scale environments. The medium diameter corresponding to the combined distribution is determined and used to partition all cloud objects into "small" and "large" groups of a particular phase. The two groups corresponding to the combined distribution have nearly equal numbers of footprints. The medium diameters are 502 km for DC and 310 km for cirrostratus. The range of the variation between two extreme phases (typically, the most active and depressed phases) for the small group is 6-11% in terms of the numbers of cloud objects and the total footprint numbers. The corresponding range for the large group is 19-44%. In terms of the probability density functions of radiative and cloud physical properties, there are virtually no differences between the MJO phases for the small group, but there are significant differences for the large groups for both DC and CS types. These results suggest that the intreseasonal variation signals reside at the large cloud clusters while the small cloud clusters represent the background noises resulting from various types of the tropical waves with different wavenumbers and propagation directions/speeds.
Xu, Kuan-Man↗
Background Noises Versus Intraseasonal Variation Signals: Small vs. Large Convective Cloud Objects From CERES Aqua Observations
During inactive phases of Madden-Julian Oscillation (MJO), there are plenty of deep but small convective systems and far fewer deep and large ones. During active phases of MJO, a manifestation of an increase in the occurrence of large and deep cloud clusters results from an amplification of large-scale motions by stronger convective heating. This study is designed to quantitatively examine the roles of small and large cloud clusters during the MJO life cycle. We analyze the cloud object data from Aqua CERES (Clouds and the Earth's Radiant Energy System) observations between July 2006 and June 2010 for tropical deep convective (DC) and cirrostratus (CS) cloud object types according to the real-time multivariate MJO index, which assigns the tropics to one of the eight MJO phases each day. The cloud object is a contiguous region of the earth with a single dominant cloud-system type. The criteria for defining these cloud types are overcast footprints and cloud top pressures less than 400 hPa, but DC has higher cloud optical depths (=10) than those of CS (<10). The size distributions, defined as the footprint numbers as a function of cloud object diameters, for particular MJO phases depart greatly from the combined (8-phase) distribution at large cloud-object diameters due to the reduced/increased numbers of cloud objects related to changes in the large-scale environments. The medium diameter corresponding to the combined distribution is determined and used to partition all cloud objects into "small" and "large" groups of a particular phase. The two groups corresponding to the combined distribution have nearly equal numbers of footprints. The medium diameters are 502 km for DC and 310 km for cirrostratus. The range of the variation between two extreme phases (typically, the most active and depressed phases) for the small group is 6-11% in terms of the numbers of cloud objects and the total footprint numbers. The corresponding range for the large group is 19-44%. In terms of the probability density functions of radiative and cloud physical properties, there are virtually no differences between the MJO phases for the small group, but there are significant differences for the large groups for both DC and CS types. These results suggest that the intreseasonal variation signals reside at the large cloud clusters while the small cloud clusters represent the background noises resulting from various types of the tropical waves with different wavenumbers and propagation speeds/directions.
Xu, Kuan-Man↗
Effects of Double-Leakage Tip Clearance Flow on the Performance of a Compressor Stage with a Large Rotor Tip Gap
Effects of a large rotor tip gap on the performance of a one and half stage axial compressor are investigated in detail with a numerical simulation based on LES and available PIV data. The current paper studies the main flow physics, including why and how the loss generation is increased with the large rotor tip gap. The present study reveals that when the tip gap becomes large, tip clearance fluid goes over the tip clearance core vortex and enters into the next blade's tip gap, which is called double-leakage tip clearance flow. As the tip clearance flow enters into the adjacent blade's tip gap, a vortex rope with a lower pressure core is generated. This vortex rope breaks up the tip clearance core vortex of the adjacent blade, resulting in a large additional mixing. This double-leakage tip clearance flow occurs at all operating conditions, from design flow to near stall condition, with the large tip gap for the current compressor stage. The double-leakage tip clearance flow, its interaction with the tip clearance core vortex of the adjacent blade, and the resulting large mixing loss are the main flow mechanism of the large rotor tip gap in the compressor. When the tip clearance is smaller, flow near the end wall follows more closely with the main passage flow and this double-leakage tip clearance flow does not happen near the design flow condition for the current compressor stage. When the compressor with a large tip gap operates at near stall operation, a strong vortex rope is generated near the leading edge due to the double-leakage flow. Part of this vortex separates from the path of the tip clearance core vortex and travels from the suction side of the blade toward the pressure side of the blade. This vortex is generated periodically at near stall operation with a large tip gap. As the vortex travels from the suction side to the pressure side of the blade, a large fluctuation of local pressure forces blade vibration. Nonsynchronous blade vibration occurs due to this vortex as the frequency of this vortex generation is not the same as the rotor. The present investigation confirms that this vortex is a part of separated tip clearance vortex, which is caused by the double-leakage tip clearance flow.
flow↗
Quantifying Distribution System Resilience From Utility Data: Large Event Risk and Benefits of Investments
We focus on blackouts in electric distribution systems that have a large cost to customers. To quantify resilience to these events, we show how to calculate risk metrics from the historical outage data routinely collected by utilities' outage management systems. Risk is defined using a customer cost exceedance curve. The exceedance curve has a heavy tail that implies large fluctuations in large blackout costs, and this makes estimating the mean large cost in the usual way impractical. To avoid this problem, we use new resilience metrics describing the large event risk; these metrics are the probability of a large cost event, the annual log cost resilience index, and the average of the logarithm of the cost of large-cost events or the slope magnitude of the tail on a log–log exceedance curve. Resilience can be improved by planned investments to upgrade system components or speed up restoration. The benefits that these investments would have had if they had been made in the past can be quantified by “rerunning history” with the effects of the investment included, and then recalculating the large event risk to find the improvement in resilience. An example using utility data shows a 2% reduction in the probability of a large cost event due to 10% wind hardening and 6%–7% reduction due to 10% faster restoration in two different areas of a distribution utility. This new data-driven approach to quantify resilience and resilience investments is realistic and much easier to apply than complicated approaches based on modeling all the phases of resilience. Moreover, an appeal to improvements to past lived experience may well be persuasive to customers and regulators in making the case for resilience investments.
24 POWER TRANSMISSION AND DISTRIBUTION↗
Measurement Capabilities of Planar Doppler Velocimetry in Large-Scale Wind Tunnels
Over the past few years, Planar Doppler Velocimetry (PDT) has been shown by several laboratories to offer an attractive means for measuring three-dimensional velocity vectors everywhere in a light sheet placed in a flow. Unlike some other optical means of measuring flow velocities, PDT is particularly attractive for use in large wind tunnels where distances to the sample region may be several meters, because it does not require the spatial resolution and tracking of individual scattering particles or the alignment of crossed beams at large distances. To date, demonstrations of PDT (also called Doppler Global Velocimetry by some authors) have been made either in low speed flows without quantitative comparison to other measurements, or in supersonic flows where the Doppler shift is large and its measurement is relatively insensitive to instrumental errors. Moreover, most reported applications have relied on the use of continuous-wave lasers, which limit the measurement to time-averaged velocity fields. This work summarizes the results of two previous studies of PDT in which the use of pulsed lasers to obtain instantaneous velocity vector fields is evaluated. The objective has been to quantitatively define and demonstrate PDT capabilities for applications in large-scale wind tunnels that are intended primarily for the testing of rotorcraft and subsonic aircraft at speeds typically less than 100 m/s. For such applications, the adequate resolution of low-speed flow fields requires accurate measurements of small Doppler shifts that are obtained at distances of several meters from the sample region and with a field of view that is sufficient to encompass the entire region of interest. The use of pulsed lasers provides the unique capability to obtain not only time-averaged fields, but also their statistical fluctuation amplitudes and the spatial excursions of unsteady flow regions such as wakes, separations, and rotor-tip vortices. To accomplish the objectives of these studies, the PDT measurement process was first modeled and its performance evaluated computationally. The noise sources considered included those related to the optical and electronic properties of Charge-Coupled Device (CCD) arrays and to speckle effects associated with the coherent illumination of aerosols from pulsed lasers. The signal noise estimates were incorporated into the PDT signal analysis process and combined with computed scattering signals using a Mae scattering theory for a distributed range of aerosol particle sizes. The results are used to define the necessary instrument configuration, to show that the expected signal levels from a practical PDT system are sufficiently large to allow its useful application in large facilities, and to show that the expected velocity measurement uncertainties are small compared to the mean velocities of interest for most subsonic, large-scale wind tunnel testing, Experimental studies using several bench-top setups are then described that validate the physics of the PDT model and demonstrate the ability to obtain accurate PDT measurements using procedures that are compatible with large-scale wind tunnel operations. The PDT measurement capabilities found in this study lead to the conclusion that PDT offers significant advantages compared to other means of measuring velocity fields in large-scale wind tunnels, including those operating at low speeds.
McKenzie, Robert L.↗
Exact evaluation of large-charge correlation functions in nonrelativistic conformal field theory
The large-charge master field which generates all n -point correlation functions with an insertion of large charge Q in nonrelativistic conformal field theory is obtained. This field is used to compute Schrödinger-invariant n -point correlation functions of large-charge operators via a direct evaluation of the path integral. Conformal dimensions are found to agree with calculations based on the state-operator correspondence. The master field solution exhibits an emergent harmonic trap whose frequency is a function of the Euclidean time. The large-charge effective action with operator insertions describes a droplet of superfluid matter whose spatial size scales with the time separation of sources. The solution is used to compute Schrödinger symmetry breaking corrections in the large-charge effective field theory (EFT) due to a finite scattering length in the fundamental theory of fermions near unitarity. The scaling of these effects in the large-charge power counting scheme is established, and the size of the effects is quantified using input from quantum Monte Carlo simulations of the near-unitary gas, as well as from the large- N expansion at large charge. Published by the American Physical Society 2024
Astronomy & Astrophysics↗
On the large eddy simulation of turbulent flows in complex geometry
Application of the method of Large Eddy Simulation (LES) to a turbulent flow consists of three separate steps. First, a filtering operation is performed on the Navier-Stokes equations to remove the small spatial scales. The resulting equations that describe the space time evolution of the 'large eddies' contain the subgrid-scale (sgs) stress tensor that describes the effect of the unresolved small scales on the resolved scales. The second step is the replacement of the sgs stress tensor by some expression involving the large scales - this is the problem of 'subgrid-scale modeling'. The final step is the numerical simulation of the resulting 'closed' equations for the large scale fields on a grid small enough to resolve the smallest of the large eddies, but still much larger than the fine scale structures at the Kolmogorov length. In dividing a turbulent flow field into 'large' and 'small' eddies, one presumes that a cut-off length delta can be sensibly chosen such that all fluctuations on a scale larger than delta are 'large eddies' and the remainder constitute the 'small scale' fluctuations. Typically, delta would be a length scale characterizing the smallest structures of interest in the flow. In an inhomogeneous flow, the 'sensible choice' for delta may vary significantly over the flow domain. For example, in a wall bounded turbulent flow, most statistical averages of interest vary much more rapidly with position near the wall than far away from it. Further, there are dynamically important organized structures near the wall on a scale much smaller than the boundary layer thickness. Therefore, the minimum size of eddies that need to be resolved is smaller near the wall. In general, for the LES of inhomogeneous flows, the width of the filtering kernel delta must be considered to be a function of position. If a filtering operation with a nonuniform filter width is performed on the Navier-Stokes equations, one does not in general get the standard large eddy equations. The complication is caused by the fact that a filtering operation with a nonuniform filter width in general does not commute with the operation of differentiation. This is one of the issues that we have looked at in detail as it is basic to any attempt at applying LES to complex geometry flows. Our principal findings are summarized.
Ghosal, Sandip↗
Planar Doppler Velocimetry for Large-Scale Wind Tunnel Testing
Recently, Planar Doppler Velocimetry (PDV) has been shown by several laboratories to offer an attractive means for measuring three-dimensional velocity vectors everywhere in a light sheet placed in a flow. Unlike other optical means of measuring flow velocities, PDV is particularly attractive for use in large wind tunnels where distances to the sample region may be several meters, because it does not require the spatial resolution and tracking of individual scattering particles or the alignment of crossed beams at large distances. To date, demonstrations of PDV have been made either in low speed flows without quantitative comparison to other measurements, or in supersonic flows where the Doppler shift is large and its measurement is relatively insensitive to instrumental errors. Moreover, most reported applications have relied on the use of continuous-wave lasers, which limit the measurement to time-averaged velocity fields. This work summarizes the results of two previous studies of PDV in which the use of pulsed lasers to obtain instantaneous velocity vector fields is evaluated. The objective has been to quantitatively define and demonstrate PDV capabilities for applications in large-scale wind tunnels that are intended primarily for the production testing of subsonic aircraft. For such applications, the adequate resolution of low-speed flow fields requires accurate measurements of small Doppler shifts that are obtained at distances of several meters from the sample region. The use of pulsed lasers provides the unique capability to obtain not only time-averaged fields, but also their statistical fluctuation amplitudes and the spatial excursions of unsteady flow regions such as wakes and separations. To accomplish the objectives indicated, the PDV measurement process is first modeled and its performance evaluated computationally. The noise sources considered include those related to the optical and electronic properties of Charge-Coupled Device (CCD) arrays and to speckle effects associated with coherent illumination from pulsed lasers. The signal noise estimates are incorporated into the PDV signal analysis process and combined with computed scattering signals using a Mie scattering theory for polydisperse smoke particles. The relevant parameters incorporate a range of practical aerodynamic test conditions and facility sizes. The results define the optimum instrument configurations, show that the expected signal levels from a practical PDV system are sufficiently large to allow its useful application in large facilities, and show that the expected velocity measurement uncertainties are small compared to the mean velocities of interest for most subsonic, large-scale wind tunnel testing. Experimental studies using several experimental bench-top setups are then described that validate the physics of the PDV model and to calibrate its computed results. The validated model allows estimates of the uncertainties of PDV measurements and a complete definition of the PDV capabilities to be made with sufficient confidence to decide the viability of PDV for large-scale wind tunnel applications.
McKenzie, Robert L.↗
Efficient Implementation of an Optimal Interpolator for Large Spatial Data Sets
Scattered data interpolation is a problem of interest in numerous areas such as electronic imaging, smooth surface modeling, and computational geometry. Our motivation arises from applications in geology and mining, which often involve large scattered data sets and a demand for high accuracy. The method of choice is ordinary kriging. This is because it is a best unbiased estimator. Unfortunately, this interpolant is computationally very expensive to compute exactly. For n scattered data points, computing the value of a single interpolant involves solving a dense linear system of size roughly n x n. This is infeasible for large n. In practice, kriging is solved approximately by local approaches that are based on considering only a relatively small'number of points that lie close to the query point. There are many problems with this local approach, however. The first is that determining the proper neighborhood size is tricky, and is usually solved by ad hoc methods such as selecting a fixed number of nearest neighbors or all the points lying within a fixed radius. Such fixed neighborhood sizes may not work well for all query points, depending on local density of the point distribution. Local methods also suffer from the problem that the resulting interpolant is not continuous. Meyer showed that while kriging produces smooth continues surfaces, it has zero order continuity along its borders. Thus, at interface boundaries where the neighborhood changes, the interpolant behaves discontinuously. Therefore, it is important to consider and solve the global system for each interpolant. However, solving such large dense systems for each query point is impractical. Recently a more principled approach to approximating kriging has been proposed based on a technique called covariance tapering. The problems arise from the fact that the covariance functions that are used in kriging have global support. Our implementations combine, utilize, and enhance a number of different approaches that have been introduced in literature for solving large linear systems for interpolation of scattered data points. For very large systems, exact methods such as Gaussian elimination are impractical since they require 0(n(exp 3)) time and 0(n(exp 2)) storage. As Billings et al. suggested, we use an iterative approach. In particular, we use the SYMMLQ method, for solving the large but sparse ordinary kriging systems that result from tapering. The main technical issue that need to be overcome in our algorithmic solution is that the points' covariance matrix for kriging should be symmetric positive definite. The goal of tapering is to obtain a sparse approximate representation of the covariance matrix while maintaining its positive definiteness. Furrer et al. used tapering to obtain a sparse linear system of the form Ax = b, where A is the tapered symmetric positive definite covariance matrix. Thus, Cholesky factorization could be used to solve their linear systems. They implemented an efficient sparse Cholesky decomposition method. They also showed if these tapers are used for a limited class of covariance models, the solution of the system converges to the solution of the original system. Matrix A in the ordinary kriging system, while symmetric, is not positive definite. Thus, their approach is not applicable to the ordinary kriging system. Therefore, we use tapering only to obtain a sparse linear system. Then, we use SYMMLQ to solve the ordinary kriging system. We show that solving large kriging systems becomes practical via tapering and iterative methods, and results in lower estimation errors compared to traditional local approaches, and significant memory savings compared to the original global system. We also developed a more efficient variant of the sparse SYMMLQ method for large ordinary kriging systems. This approach adaptively finds the correct local neighborhood for each query point in the interpolation process.
Memarsadeghi, Nargess↗