Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “space grid”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Impact of Tropical Cyclone Wind Forcing on the Global Climate in a Fully Coupled Climate Model

Tropical cyclones (TCs) alter upper-ocean temperature and influence ocean heat content via enhanced turbulent mixing. A better understanding of the role of TCs within the climate system requires a fully coupled modeling framework, where TC-induced ocean responses feed back to the atmosphere and subsequently to the climate mean state and variability. Here, we investigate the impacts of TC wind forcing on the global ocean and the associated feedbacks within the climate system using the fully coupled Community Earth System Model version 1.3 (CESM1.3). Using the low-resolution version of CESM1.3 (1° atmosphere and ocean grid spacing) with no intrinsic TCs, we conduct a suite of sensitivity experiments by inserting TC winds extracted from a high-resolution (0.25° atmosphere grid spacing) TC-permitting simulation into the low-resolution model. Results from the low-resolution TC experiment are compared to a low-resolution control simulation to diagnose TCs’ impact. We found that the added TC winds can increase ocean heat content by affecting ocean vertical mixing, air–sea enthalpy fluxes, and cloud amount. Here, the added TCs can influence mean SST, precipitation, ocean subsurface temperature, and ocean mixed layer depth. We found a strengthening of the wind-driven subtropical cells and a weakening of the Atlantic meridional overturning circulation due to the changes of surface buoyancy fluxes. TCs in the model cause anomalous equatorward ocean heat convergence in the deep tropics and an increase of poleward ocean heat transport out of the subtropics. Our modeling results provide new insights into the multiscale interactions between TCs and the coupled climate system.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

GPAW: An open Python package for electronic structure calculations

We review the GPAW open-source Python package for electronic structure calculations. GPAW is based on the projector-augmented wave method and can solve the self-consistent density functional theory (DFT) equations using three different wave-function representations, namely real-space grids, plane waves, and numerical atomic orbitals. The three representations are complementary and mutually independent and can be connected by transformations via the real-space grid. This multi-basis feature renders GPAW highly versatile and unique among similar codes. By virtue of its modular structure, the GPAW code constitutes an ideal platform for the implementation of new features and methodologies. Moreover, it is well integrated with the Atomic Simulation Environment (ASE), providing a flexible and dynamic user interface. In addition to ground-state DFT calculations, GPAW supports many-body GW band structures, optical excitations from the Bethe–Salpeter Equation, variational calculations of excited states in molecules and solids via direct optimization, and real-time propagation of the Kohn–Sham equations within time-dependent DFT. A range of more advanced methods to describe magnetic excitations and non-collinear magnetism in solids are also now available. In addition, GPAW can calculate non-linear optical tensors of solids, charged crystal point defects, and much more. Recently, support for graphics processing unit (GPU) acceleration has been achieved with minor modifications to the GPAW code thanks to the CuPy library. We end the review with an outlook, describing some future plans for GPAW.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1990

The molecular spiral arms of NGC 6946

From CO-12(J=1 to 0) observations at 45 seconds resolution Tacconi and Young (1989) have found evidence for enhancements in both the CO emissivity and the massive star formation efficiency (MSFE) on optical spiral arms of the bright spiral galaxy NGC 6946. In the optically luminous and well-defined spiral arm in the NE quadrant, there are enhancements in both the H2 surface density and MSFE relative to the interarm regions. In contrast, a poorly defined arm in the SW shows no arm-interarm contrast in the MSFE. To further investigate the molecular gas content of these two spiral arms, researchers have made CO-12 J=2 to 1 and 3 to 2 observations with the James Clerk Maxwell Telescope. In the J=2 to 1 line, they made observations of the NE and SW spiral arm and interarm regions in 4 x 9 10 seconds spaced grids (36 points per grid). Because of decreased sensitivity in the J=3 to 2 line, they were limited to mapping the two arm regions in 2 x 3 10 seconds spaced grids (6 points per grid). The centers of each of the grids lie 2.4 minutes to the NE and 2.3 minutes to the SW of the nucleus of NGC 6946. With the CO J=2 to 1 data researchers are able to fully resolve the two observed spiral arms in NGC 6946. In both cases the CO emission is largely confined to the optical spiral arm regions with the peak observed T asterisk sub A being up to 4 times higher on the spiral arms than in the interarm regions. Researchers are currently estimating massive star formation efficiencies on and off the spiral arms through direct comparison of the CO maps with an H alpha image. They are also comparing the CO J=2 to 1 data with an HI map made at similar resolution. Thus, they will be able to determine structure in all components of the IS on scales of less than 20 inches.

Tacconi, L. J.↗
NASA NTRS2008

Calculations of High-Temperature Jet Flow Using Hybrid Reynolds-Average Navier-Stokes Formulations

Two multiscale-type turbulence models are implemented in the PAB3D solver. The models are based on modifying the Reynolds-averaged Navier Stokes equations. The first scheme is a hybrid Reynolds-averaged- Navier Stokes/large-eddy-simulation model using the two-equation k(epsilon) model with a Reynolds-averaged-Navier Stokes/large-eddy-simulation transition function dependent on grid spacing and the computed turbulence length scale. The second scheme is a modified version of the partially averaged Navier Stokes model in which the unresolved kinetic energy parameter f(sub k) is allowed to vary as a function of grid spacing and the turbulence length scale. This parameter is estimated based on a novel two-stage procedure to efficiently estimate the level of scale resolution possible for a given flow on a given grid for partially averaged Navier Stokes. It has been found that the prescribed scale resolution can play a major role in obtaining accurate flow solutions. The parameter f(sub k) varies between zero and one and is equal to one in the viscous sublayer and when the Reynolds-averaged Navier Stokes turbulent viscosity becomes smaller than the large-eddy-simulation viscosity. The formulation, usage methodology, and validation examples are presented to demonstrate the enhancement of PAB3D's time-accurate turbulence modeling capabilities. The accurate simulations of flow and turbulent quantities will provide a valuable tool for accurate jet noise predictions. Solutions from these models are compared with Reynolds-averaged Navier Stokes results and experimental data for high-temperature jet flows. The current results show promise for the capability of hybrid Reynolds-averaged Navier Stokes and large eddy simulation and partially averaged Navier Stokes in simulating such flow phenomena.

Abdol-Hamid, Khaled S.↗
NASA NTRS2005

Numerical Study of High-Temperature Jet Flow Using RANS/LES and PANS Formulations

Two multi-scale-type turbulence models are implemented in the PAB3D solver. The models are based on modifying the Reynolds Averaged Navier-Stokes (RANS) equations. The first scheme is a hybrid RANS/LES model utilizing the two-equation (k(epsilon)) model with a RANS/LES transition function dependent on grid spacing and the computed turbulence length scale. The second scheme is a modified version of the Partially Averaged Navier-Stokes (PANS) model, where the unresolved kinetic energy parameter (f(sub k)) is allowed to vary as a function of grid spacing and the turbulence length scale. This parameter is estimated based on a novel two-stage procedure to efficiently estimate the level of scale resolution possible for a given flow on a given grid for Partial Averaged Navier-Stokes (PANS). It has been found that the prescribed scale resolution can play a major role in obtaining accurate flow solutions. The parameter f(sub k) varies between zero and one and equal to one in the viscous sub layer, and when the RANS turbulent viscosity becomes smaller than the LES viscosity. The formulation, usage methodology, and validation examples are presented to demonstrate the enhancement of PAB3D's time-accurate and turbulence modeling capabilities. The accurate simulations of flow and turbulent quantities will provide valuable tool for accurate jet noise predictions. Solutions from these models are compared to RANS results and experimental data for high-temperature jet flows. The current results show promise for the capability of hybrid RANS/LES and PANS in simulating such flow phenomena.

Abdol-Hamid, Khaled S.↗
NASA NTRS1993

Computations of spray, fuel-air mixing, and combustion in a lean-premixed-prevaporized combustor

A code was developed for computing the multidimensional flow, spray, combustion, and pollutant formation inside gas turbine combustors. The code developed is based on a Lagrangian-Eulerian formulation and utilizes an implicit finite-volume method. The focus of this paper is on the spray part of the code (both formulation and algorithm), and a number of issues related to the computation of sprays and fuel-air mixing in a lean-premixed-prevaporized combustor. The issues addressed include: (1) how grid spacings affect the diffusion of evaporated fuel, and (2) how spurious modes can arise through modelling of the spray in the Lagrangian computations. An upwind interpolation scheme is proposed to account for some effects of grid spacing on the artificial diffusion of the evaporated fuel. Also, some guidelines are presented to minimize errors associated with the spurious modes.

Dasgupta, A.↗
DOE OSTI2024

Application of the AI2 Climate Emulator to E3SMv2's Global Atmosphere Model, With a Focus on Precipitation Fidelity

Abstract Can the current successes of global machine learning‐based weather simulators be generalized beyond 2‐week forecasts to stable and accurate multiyear runs? The recently developed AI2 Climate Emulator (ACE) suggests this is feasible, based upon 10‐year simulations with a network trained on output from a physics‐based global atmosphere model using a grid spacing of approximately 110 km and forced by a repeating annual cycle of sea‐surface temperature. Here we show that ACE, without modification, can be trained to emulate another major atmospheric model, EAMv2, run at a comparable grid spacing for at least 10 years with similarly small climate biases—a prerequisite to wider applicability. With an analysis that combines multiple temporal, spatial, and frequency domain perspectives, we show that ACE faithfully represents the spatiotemporal structure of EAMv2 precipitation and related variables. Finally, we show that a pretrained ACE network is able to adapt to a new global climate model simulation data set with 10 fewer training steps than when starting from random initialization, all while still maintaining low levels of climate bias. Further analysis of these fine‐tuning experiments reveal ACE's intriguing ability to interpolate between distinct global climate models.

Duncan, James P. C.↗
DOE OSTI2020

A Baseline for Global Weather and Climate Simulations at 1 km Resolution

In an attempt to advance the understanding of the Earth's weather and climate by representing deep convection explicitly, we present a global, four-month simulation (November 2018 to February 2019) with ECMWF's hydrostatic Integrated Forecasting System (IFS) at an average grid spacing of 1.4 km. The impact of explicitly simulating deep convection on the atmospheric circulation and its variability is assessed by comparing the 1.4 km simulation to the equivalent well-tested and calibrated global simulations at 9 km grid spacing with and without parametrized deep convection. The explicit simulation of deep convection at 1.4 km results in a realistic large-scale circulation, better representation of convective storm activity, and stronger convective gravity wave activity when compared to the 9 km simulation with parametrized deep convection. Comparison of the 1.4 km simulation to the 9 km simulation without parametrized deep convection shows that switching off deep convection parametrization at a too coarse resolution (i.e., 9 km) generates too strong convective gravity waves. Based on the limited statistics available, improvements to the Madden-Julian Oscillation or tropical precipitation are not observed at 1.4 km, suggesting that other Earth system model components and/or their interaction are important for an accurate representation of these processes and may well need adjusting at deep convection resolving resolutions. Overall, the good agreement of the 1.4 km simulation with the 9 km simulation with parametrized deep convection is remarkable, despite one of the most fundamental parametrizations being turned off at 1.4 km resolution and despite no adjustments being made to the remaining parametrizations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Improving the precision of forces in real-space pseudopotential density functional theory

The high-order finite difference real-space pseudopotential density functional theory (DFT) approach is a valuable method for large-scale, massively parallel DFT calculations. A significant challenge in the approach is the oscillating “egg-box” error introduced by aliasing associated with a coarse grid spacing. To address this issue while minimizing computational cost, we developed a finite difference interpolation (FDI) scheme [Roller et al., J. Chem. Theory Comput. 19, 3889 (2023)] as a means of exploiting the high resolution of the pseudopotential to reduce egg-box effects systematically. Here, we show an implementation of this method in the PARSEC code and examine the practical utility of the combination of FDI with additional methods for improving force precision and/or reducing its computational cost, including orbital-based forces, compensating charges (namely, adding and subtracting a judiciously chosen charge density such that the total density is unaltered), and a modified spatial domain in which the real-space grid is defined. Using selected small molecules, as well as metallic Li, as test cases, we show that a combination of all four aspects leads to a significant reduction in computational cost while retaining a high level of precision that supports accurate structures and vibrational spectra, as well as stable and accurate molecular dynamics runs.

Chemistry↗
NASA NTRS2004

HiPEP Ion Optics System Evaluation Using Gridlets

Experimental measurements are presented for sub-scale ion optics systems comprised of 7 and 19 aperture pairs with geometrical features that are similar to the HiPEP ion optics system. Effects of hole diameter and grid-to-grid spacing are presented as functions of applied voltage and beamlet current. Recommendations are made for the beamlet current range where the ion optics system can be safely operated without experiencing direct impingement of high energy ions on the accelerator grid surface. Measurements are also presented of the accelerator grid voltage where beam plasma electrons backstream through the ion optics system. Results of numerical simulations obtained with the ffx code are compared to both the impingement limit and backstreaming measurements. An emphasis is placed on identifying differences between measurements and simulation predictions to highlight areas where more research is needed. Relatively large effects are observed in simulations when the discharge chamber plasma properties and ion optics geometry are varied. Parameters investigated using simulations include the applied voltages, grid spacing, hole-to-hole spacing, doubles-to-singles ratio, plasma potential, and electron temperature; and estimates are provided for the sensitivity of impingement limits on these parameters.

Willliams, John D.↗
DOE OSTI2022

Theory, observations, and simulations of kinetic entropy in a magnetotail electron diffusion region

In this work, we examine velocity-space kinetic entropy, a spatially local measure of entropy for systems out of thermal equilibrium, during an encounter of an electron diffusion region at a magnetic reconnection site in Earth’s magnetotail by the Magnetospheric Multiscale (MMS) mission. We start by generalizing the theory of kinetic entropy to the case of non-uniform velocity space grids and transforming the equations into spheri- cal energy coordinates useful to experimental plasma detectors. The theory is then applied to MMS data and compared to particle-in-cell simulations of reconnection. We demonstrate that the entropy-based non-Maxwellianity measure from the MMS data is of sufficiently high precision to reliably identify non-Maxwellian distributions and therefore the measurements when kinetic effects are most significant. By comparing two different non-Maxwellian measures, we show that total entropy density suffers from “information loss” because it lacks a dependence on the velocity space grid, and so has lost information about how well a distribution function is resolved. Local velocity-space kinetic entropy density recovers this information. We quantify information loss and argue that the considerations needed to minimize it are crucial for instruments designed to measure distribution functions in situ.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2004

Numerical Study of Flow Past a Circular Cylinder Using RANS, Hybrid RANS/LES and PANS Formulations

Two multiscale type turbulence models are implemented in the PAB3D solver. The models are based on modifying the Reynolds Averaged Navier-Stokes (RANS) equations. The first scheme is a hybrid RANS/LES model utilizing the two-equation (k(sub epsilon)) model with a RANS/LES transition function dependent on grid spacing and the computed turbulence length scale. The second scheme is a modified version of the partially averaged Navier-Stokes (PANS) model, where the unresolved kinetic energy parameter (f(sub k)) is allowed to vary as a function of grid spacing and the turbulence length scale. Solutions from these models are compared to RANS results and experimental data for a stationary and rotating cylinder. The parameter f(sub k) varies between zero and one and has the characteristic to be equal to one in the viscous sub layer, and when the RANS turbulent viscosity becomes smaller than the LES viscosity. The formulation, usage methodology, and validation example are presented to demonstrate the enhancement of PAB3D's time-accurate and turbulence modeling capabilities. The models are compared to RANS results and experimental data for turbulent separated flows (TS) and laminar separated flows (LS) around stationary and rotating cylinders. For a stationary cylinder, the TS case is accurately simulated using the general two-equation k(sub epsilon) turbulence model (eddy-viscosity model). PAB3D accurately predicts the drag coefficient (CD), lift coefficient (CL) and the Strouhal number (St). The LS case was a challenge for the RANS computation with an eddy-viscosity turbulence model. The RANS/LES and PANS performed well and showed marked improvements over the RANS solution. The modified PANS model was the most accurate. For the rotating cylinder, the spin ratio varied from zero to one, and the PANS results were in good agreement with published experimental data. RANS/LES and PANS capture both temporal and spatial fluctuations and produce large-scale structures that do not occur in the RANS simulation. The current results show promise for the capability of PANS in simulating unsteady and complex flow phenomena.

Elmiligui, Alaa A.↗
NASA NTRS2013

Downscaling Reanalysis over Continental Africa with a Regional Model: NCEP Versus ERA Interim Forcing

Five annual climate cycles (1998-2002) are simulated for continental Africa and adjacent oceans by a regional atmospheric model (RM3). RM3 horizontal grid spacing is 0.44deg at 28 vertical levels. Each of 2 simulation ensembles is driven by lateral boundary conditions from each of 2 alternative reanalysis data sets. One simulation downs cales National Center for Environmental Prediction reanalysis 2 (NCPR2) and the other the European Centre for Medium Range Weather Forecasts Interim reanalysis (ERA-I). NCPR2 data are archived at 2.5deg grid spacing, while a recent version of ERA-I provides data at 0.75deg spacing. ERA-I-forced simulations are recomrp. ended by the Coordinated Regional Downscaling Experiment (CORDEX). Comparisons of the 2 sets of simulations with each other and with observational evidence assess the relative performance of each downscaling system. A third simulation also uses ERA-I forcing, but degraded to the same horizontal resolution as NCPR2. RM3-simulated pentad and monthly mean precipitation data are compared to Tropical Rainfall Measuring Mission (TRMM) data, gridded at 0.5deg, and RM3-simulated circulation is compared to both reanalyses. Results suggest that each downscaling system provides advantages and disadvantages relative to the other. The RM3/NCPR2 achieves a more realistic northward advance of summer monsoon rains over West Africa, but RM3/ERA-I creates the more realistic monsoon circulation. Both systems recreate some features of JulySeptember 1999 minus 2002 precipitation differences. Degrading the resolution of ERA-I driving data unrealistically slows the monsoon circulation and considerably diminishes summer rainfall rates over West Africa. The high resolution of ERA-I data, therefore, contributes to the quality of the downscaling, but NCPR2laterai boundary conditions nevertheless produce better simulations of some features.

Reanalysis↗
NASA NTRS2001

Effect of Under-Resolved Grids on High Order Methods

There has been much discussion on verification and validation processes for establishing the credibility of CFD simulations. Since the early 1990s, many of the aeronautical and mechanical engineering related reference journals mandated that any accepted articles in numerical simulations (without known solutions to compared with) need to perform a minimum of one level of grid refinement and time step reduction. Due to the difficulty in analysis, the effect of under-resolved grids and the nonlinear behavior of available spatial discretizations, are scarcely discussed in the literature. Here, an under-resolved numerical simulation is one where the grid spacing being used is too coarse to resolve the smallest physically relevant scales of the chosen continuum governing equations that are of interest to the numerical modeler. With the advent of new developments in fourth-order or higher spatial schemes, it has become common to regard high order schemes as more accurate, reliable and require less grid points. The danger comes when one tries to perform computations with the coarsest grid possible while still hoping to maintain numerical results sufficiently accurate for complex flows, and especially, data-limited problems. On one hand, high order methods when applies to highly coupled multidimensional complex nonlinear problems might have different stability, convergence and reliability behavior than their well studied low order counterparts, especially for nonlinear schemes such as TVD, MUSCL with limiters, ENO, WENO and discrete Galerkin. On the other hand, high order methods involve more operation counts and systematic grid convergence study can be time consuming and prohibitively expansive. At the same time it is difficult to fully understand or categorize the different nonlinear behavior of finite discretizations, especially at the limits of under-resolution when different types of bifurcation phenomena might occur, depending on the combination of grid spacings, time steps, initial conditions and numerical treatments of boundary conditions.

Yee, H. C.↗
DOE OSTI2022

Using aircraft measurements to characterize subgrid-scale variability of aerosol properties near the Atmospheric Radiation Measurement Southern Great Plains site

Abstract. Complex distributions of aerosol properties evolve in space and time as a function of emissions, new particle formation, coagulation, condensational growth, chemical transformation, phase changes, turbulent mixing and transport, removal processes, and ambient meteorological conditions. The ability of chemical transport models to represent the multi-scale processes affecting the life cycle of aerosols depends on their spatial resolution since aerosol properties are assumed to be constant within a grid cell. Subgrid-scale-dependent processes that affect aerosol populations could have a significant impact on the formation of particles, their growth to cloud condensation nuclei (CCN) sizes, aerosol–cloud interactions, dry deposition and rainout and hence their burdens, lifetimes, and radiative forcing. To address this issue, we characterize subgrid-scale variability in terms of measured aerosol number, size, composition, hygroscopicity, and CCN concentrations made by repeated aircraft flight paths over the Atmospheric Radiation Measurement (ARM) program's Southern Great Plains (SGP) site during the Holistic Interactions of Shallow Clouds, Aerosols and Land Ecosystem (HI-SCALE) campaign. Subgrid variability is quantified in terms of both normalized frequency distributions and percentage difference percentiles using grid spacings of 3, 9, 27, and 81 km that represent those typically used by cloud-system-resolving models as well as the current and next-generation climate models. Even though the SGP site is a rural location, surprisingly large horizontal gradients in aerosol properties were frequently observed. For example, 90 % of the 3, 9, and 27 km cell mean organic matter concentrations differed from the 81 km cell around the SGP site by as much as ∼ 46 %, large spatial variability in aerosol number concentrations and size distributions were found during new particle formation events, and consequently 90 % of the 3, 9, and 27 km cell mean CCN number concentrations differed from the 81 km cell mean by as much as ∼ 38 %. The spatial variability varied seasonally for some aerosol properties, with some having larger spatial variability during the spring and others having larger variability during the late summer. While measurements at a single surface site cannot reflect the surrounding variability of aerosol properties at a given time, aircraft measurements that are averaged within an 81 km cell were found to be similar to many, but not all, aerosol properties measured at the ground SGP site. This analysis suggests that it is reasonable to directly compare most ground SGP site aerosol measurements with coarse global climate model predictions. In addition, the variability quantified by the aircraft can be used as an uncertainty range when comparing the surface point measurements with model predictions that use coarse grid spacings.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1990

Topographic distribution of clear-sky radiation over the Konza Prairie, Kansas

A method for exploring the topographic modulation of clear-sky incoming radiation using a two-stream atmospheric model together with digital elevation data is presented to obtain the temporal and spatial distribution of this energy. Utilizing this methodology for the Konza Prairie, it is observed that the variance and spatial autocorrelation of the simulated radiation data changed with sun angle and elevation grid spacing. An analytical formula that describes the behavior of the variance as a function of optical depth, sun angle, and mean terrain slope can be developed by considering direct solar radiation variability on a constant-slope terrain which has a uniform albedo and a uniform distribution of slopes. The decrease in mean slope and the increase in slope and autocorrelation is attributed to a loss of some of the shorter-wavelength terrain features at coarser grid spacings.

Dubayah, Ralph↗
NASA NTRS1994

Inducer analysis/pump model development

Current design of high performance turbopumps for rocket engines requires effective and robust analytical tools to provide design information in a productive manner. The main goal of this study was to develop a robust and effective computational fluid dynamics (CFD) pump model for general turbopump design and analysis applications. A finite difference Navier-Stokes flow solver, FDNS, which includes an extended k-epsilon turbulence model and appropriate moving zonal interface boundary conditions, was developed to analyze turbulent flows in turbomachinery devices. In the present study, three key components of the turbopump, the inducer, impeller, and diffuser, were investigated by the proposed pump model, and the numerical results were benchmarked by the experimental data provided by Rocketdyne. For the numerical calculation of inducer flows with tip clearance, the turbulence model and grid spacing are very important. Meanwhile, the development of the cross-stream secondary flow, generated by curved blade passage and the flow through tip leakage, has a strong effect on the inducer flow. Hence, the prediction of the inducer performance critically depends on whether the numerical scheme of the pump model can simulate the secondary flow pattern accurately or not. The impeller and diffuser, however, are dominated by pressure-driven flows such that the effects of turbulence model and grid spacing (except near leading and trailing edges of blades) are less sensitive. The present CFD pump model has been proved to be an efficient and robust analytical tool for pump design due to its very compact numerical structure (requiring small memory), fast turnaround computing time, and versatility for different geometries.

Cheng, Gary C.↗