Search NASA⌕ Search

SEARCH · Search NASA

Results for “engine CFD”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

A Scaling Study for Incompressible Multispecies Solver in Vertex-CFD

Multispecies incompressible flows occur widely in engineering and environmental applications, such as chemical reactors, fuel cells, ocean mixing, and biomedical systems. However, accurately resolving the complex transport and mixing phenomena associated with multiple interacting species remains computationally challenging, especially for large-scale problems. In this study, we present a robust, high-performance computing--enabled multispecies incompressible Navier–Stokes solver integrated within the Vertex-CFD framework. Our solver employs a fully coupled, implicit, finite element--based formulation that accurately captures the advection, diffusion, and interaction of multiple species in incompressible flows by leveraging the Kokkos library for parallel computing to achieve high computational efficiency. For pressure coupling, the entropically damped artificial compressibility method is utilized. We validated the solver against canonical test cases, including multispecies advection, diffusion, and Bateman systems; the results demonstrate second- and third-order spatial accuracy and consistent convergence. Additionally, we demonstrated the strong and weak scaling study results obtained on the leadership-class high-performance computing system, Frontier at Oak Ridge National Laboratory.

Oz, Furkan [ORNL] (ORCID:0000000265831724)↗
NASA NTRS2015

Wind Tunnel Interference Effects on Tilt Rotor Testing Using Computational Fluid Dynamics

Experimental techniques to measure rotorcraft aerodynamic performance are widely used. However, most of them are either unable to capture interference effects from bodies, or require an extremely large computational budget. The objective of the present research is to develop an XV-15 Tilt Rotor Research Aircraft rotor model for investigation of wind tunnel wall interference using a novel Computational Fluid Dynamics (CFD) solver for rotorcraft, RotCFD. In RotCFD, a mid-fidelity URANS solver is used with an incompressible flow model and a realizable k-ε turbulence model. The rotor is, however, not modeled using a computationally expensive, unsteady viscous body-fitted grid, but is instead modeled using a blade element model with a momentum source approach. Various flight modes of the XV-15 isolated rotor, including hover, tilt and airplane mode, have been simulated and correlated to existing experimental and theoretical data. The rotor model is subsequently used for wind tunnel wall interference simulations in the National Full-Scale Aerodynamics Complex (NFAC) at NASA Ames Research Center in California. The results from the validation of the isolated rotor performance showed good correlation with experimental and theoretical data. The results were on par with known theoretical analyses. In RotCFD the setup, grid generation and running of cases is faster than many CFD codes, which makes it a useful engineering tool. Performance predictions need not be as accurate as high-fidelity CFD codes, as long as wall effects can be properly simulated. For both test sections of the NFAC wall interference was examined by simulating the XV-15 rotor in the test section of the wind tunnel and with an identical grid but extended boundaries in free field. Both cases were also examined with an isolated rotor or with the rotor mounted on the modeled geometry of the Tiltrotor Test Rig (TTR). A 'quasi linear trim' was used to trim the thrust for the rotor to compare the power as a unique variable. Power differences between free field and wind tunnel cases were found from -7 % to 0 % in the 80- by 120-Foot Wind Tunnel test section and -1.6 % to 4.8 % in the 40- by 80-Foot Wind Tunnel, depending on the TTR orientation, tunnel velocity and blade setting. The TTR will be used in 2016 to test the Bell 609 rotor in a similar fashion to the research in this report.

CFD↗
NASA NTRS2016

Wind Tunnel Interference Effects on Tilt Rotor Testing Using Computational Fluid Dynamics

Experimental techniques to measure rotorcraft aerodynamic performance are widely used. However, most of them are either unable to capture interference effects from bodies, or require an extremely large computational budget. The objective of the present research is to develop an XV-15 Tiltrotor Research Aircraft rotor model for investigation of wind tunnel wall interference using a novel Computational Fluid Dynamics (CFD) solver for rotorcraft, RotCFD. In RotCFD, a mid-fidelity Unsteady Reynolds Averaged Navier-Stokes (URANS) solver is used with an incompressible flow model and a realizable k-ε turbulence model. The rotor is, however, not modeled using a computationally expensive, unsteady viscous body-fitted grid, but is instead modeled using a blade-element model (BEM) with a momentum source approach. Various flight modes of the XV-15 isolated rotor, including hover, tilt, and airplane mode, have been simulated and correlated to existing experimental and theoretical data. The rotor model is subsequently used for wind tunnel wall interference simulations in the National Full-Scale Aerodynamics Complex (NFAC) at Ames Research Center in California. The results from the validation of the isolated rotor performance showed good correlation with experimental and theoretical data. The results were on par with known theoretical analyses. In RotCFD the setup, grid generation, and running of cases is faster than many CFD codes, which makes it a useful engineering tool. Performance predictions need not be as accurate as high-fidelity CFD codes, as long as wall effects can be properly simulated. For both test sections of the NFAC wall, interference was examined by simulating the XV-15 rotor in the test section of the wind tunnel and with an identical grid but extended boundaries in free field. Both cases were also examined with an isolated rotor or with the rotor mounted on the modeled geometry of the Tiltrotor Test Rig (TTR). A "quasi linear trim" was used to trim the thrust for the rotor to compare the power as a unique variable. Power differences between free field and wind tunnel cases were found from -7 to 0 percent in the 80- by 120-Foot Wind Tunnel and -1.6 to 4.8 percent in the 40- by 80-Foot Wind Tunnel, depending on the TTR orientation, tunnel velocity, and blade setting. The TTR will be used in 2016 to test the Bell 609 rotor in a similar fashion to the research in this report.

CFD↗
NASA NTRS2016

Using RotCFD to Predict Isolated XV-15 Rotor Performance

Experimental techniques to measure rotorcraft aerodynamic performance are widely used. However, the need exists to understand the limitations of ground based testing by augmenting the analysis of experimental test results with Computational Fluid Dynamics (CFD) modeling. The immediate objective of the present research is to develop an XV-15 Tilt Rotor Research Aircraft rotor model for investigation of wind tunnel wall interference. The predicted performance of the XV-15 during various flight modes is compared to theoretical and experimental data. This research is performed to support wind tunnel tests scheduled for 2016. A mid-fidelity RANS solver, RotCFD, is used with an unsteady, incompressible flow model and a realizable k- turbulence model. The rotor is modeled using an actuator disk model or blade element model with a momentum source approach. In RotCFD the setup, grid generation and running of cases is faster than many CFD codes which makes it a useful engineering tool. Performance predictions need not be as accurate as high-fidelity CFD codes, as long as wall effects can be properly simulated. Being able to accurately predict unsteady rotorcraft performance on desktop-class computers provides a quicker analysis of highly complex flows during the initial design phase.

CFD↗
NASA NTRS1991

CFD Applications in Propulsion

An overview of various applications of Computational Fluid Dynamics (CFD) algorithms to propulsion problems is given. Topics of interest include incompressible, low speed compressible, transonic, and supersonic problems. A common family of algorithms is used for all applications and emphasis is placed on maintaining accuracy and convergence efficiency for all problems. Specific problems include pump hydrodynamics, simultaneous combustion and mixing in rocket engines, viscous nozzle flow, and CFD applications to combustion stability.

Merkle, Charles L.↗
NASA NTRS2003

Enhancement and Extension of Porosity Model in the FDNS-500 Code to Provide Enhanced Simulations of Rocket Engine Components

In the past, the design of rocket engines has primarily relied on the cold flow/hot fire test, and the empirical correlations developed based on the database from previous designs. However, it is very costly to fabricate and test various hardware designs during the design cycle, whereas the empirical model becomes unreliable in designing the advanced rocket engine where its operating conditions exceed the range of the database. The main goal of the 2nd Generation Reusable Launching Vehicle (GEN-II RLV) is to reduce the cost per payload and to extend the life of the hardware, which poses a great challenge to the rocket engine design. Hence, understanding the flow characteristics in each engine components is thus critical to the engine design. In the last few decades, the methodology of computational fluid dynamics (CFD) has been advanced to be a mature tool of analyzing various engine components. Therefore, it is important for the CFD design tool to be able to properly simulate the hot flow environment near the liquid injector, and thus to accurately predict the heat load to the injector faceplate. However, to date it is still not feasible to conduct CFD simulations of the detailed flowfield with very complicated geometries such as fluid flow and heat transfer in an injector assembly and through a porous plate, which requires gigantic computer memories and power to resolve the detailed geometry. The rigimesh (a sintered metal material), utilized to reduce the heat load to the faceplate, is one of the design concepts for the injector faceplate of the GEN-II RLV. In addition, the injector assembly is designed to distribute propellants into the combustion chamber of the liquid rocket engine. A porosity mode thus becomes a necessity for the CFD code in order to efficiently simulate the flow and heat transfer in these porous media, and maintain good accuracy in describing the flow fields. Currently, the FDNS (Finite Difference Navier-Stakes) code is one of the CFD codes which are most widely used by research engineers at NASA Marshall Space Flight Center (MSFC) to simulate various flow problems related to rocket engines. The objective of this research work during the 10-week summer faculty fellowship program was to 1) debug the framework of the porosity model in the current FDNS code, and 2) validate the porosity model by simulating flows through various porous media such as tube banks and porous plate.

Cheng, Gary↗
DOE OSTI2022

Effects of stratification and charge cooling on combustion in a gasoline direct-injection compression ignition (GDCI) engine

With the development of low temperature engine combustion strategies, performance of gasoline-type fuels under compression ignition conditions has attracted extensive research interest. Meanwhile, for the sake of co-optimization of engines and fuels for future ground transportation, identification and evaluation of general fuel properties should be a core research priority instead of endless testing of specific fuels. In this study, the roles of fuel octane sensitivity in characterizing the ignition performance of gasoline surrogates have been systematically investigated under typical gasoline direct ignition compression ignition (GDCI) engine conditions using 3D combustion CFD simulation, especially considering the subsequent in-cylinder charge stratification and charge cooling. Two different operating conditions, high boost pressure low boost temperature (beyond-RON) case and low boost pressure high boost temperature (beyond-MON) case, were considered. By comparing with our previous zero-dimensional chemical kinetic study of gasoline surrogates in advanced compression ignition (ACI) engines, the effects of stratification and charge cooling on the combustion processes are investigated. It is found that different fuel octane sensitivities lead to slight difference in equivalence ratio stratification and charge cooling due to differences in volatility. However, fuel reactivity is still the more dominant factor than the stratification and charge cooling effects in determining combustion phasing. The present results help to justify the P-T domain framework for engine autoignition analysis of overlapping pressure-temperature trajectory with ignition delay iso-contour. Here, the results also provide useful guidance to the understanding of GCI combustion process, and to the evaluation of controlling fuel properties and the selection of alternative fuels in GCI engines.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

A Study of Propane Combustion in a Spark-Ignited Cooperative Fuel Research (CFR) Engine

Liquefied petroleum gas (LPG), whose primary composition is propane, is a promising candidate for heavy-duty vehicle applications as a diesel fuel alternative due to its CO 2 reduction potential and high knock resistance. To realize diesel-like efficiencies, spark-ignited LPG engines are proposed to operate near knock-limit over a wide range of operating conditions, which necessitates an investigation of fuel-engine interactions that leads to end-gas autoignition with propane combustion. This work presents both experimental and numerical studies of stoichiometric propane combustion in a sparkignited (SI) cooperative fuel research (CFR) engine. Engine experiments are initially conducted at different compression ratio (CR) values, and the effects of CR on engine combustion are characterized. A three-pressure analysis (TPA) model based on the two-zone combustion concept is developed in GT-Power and validated using test results to estimate in-cylinder wall temperatures, residual gas fraction, etc. This model is further utilized to examine end-gas chemistry by enabling the SI turbulent flame combustion and unburned gas chemical kinetics modules. Finally, a three-dimensional (3D) computational fluid dynamic (CFD) model of the CFR engine is developed in CONVERGE, where the G-equation and SAGE detailed chemical kinetics models are implemented for combustion modeling. Here, a 153 species reduced chemical kinetics mechanism derived from the detailed NUIGMech1.1 mechanism based on the ignition delay and laminar flame speed (LFS) studies is used to generate an LFS lookup table and to describe end-gas autoignition chemistry. Multi-cycle Reynolds-averaged Navier-Stokes (RANS) simulations are then performed for the tested CRs, and the numerical model is shown to be capable of predicting the propane combustion characteristics, particularly the end-gas autoignition behavior.

20 FOSSIL-FUELED POWER PLANTS↗
NASA NTRS2011

Computational Analysis of a Low-Boom Supersonic Inlet

A low-boom supersonic inlet was designed for use on a conceptual small supersonic aircraft that would cruise with an over-wing Mach number of 1.7. The inlet was designed to minimize external overpressures, and used a novel bypass duct to divert the highest shock losses around the engine. The Wind-US CFD code was used to predict the effects of capture ratio, struts, bypass design, and angles of attack on inlet performance. The inlet was tested in the 8-ft by 6-ft Supersonic Wind Tunnel at NASA Glenn Research Center. Test results showed that the inlet had excellent performance, with capture ratios near one, a peak core total pressure recovery of 96 percent, and a stable operating range much larger than that of an engine. Predictions generally compared very well with the experimental data, and were used to help interpret some of the experimental results.

Chima, Rodrick V.↗
DOE OSTI2023

Modeling the impact of the fuel injection strategy on the combustion and performance characteristics of a heavy-duty GCI engine

Gasoline compression ignition (GCI) is a promising strategy to achieve high thermal efficiency and low emissions with limited modifications to the conventional diesel engine hardware. It is a partially premixed concept, which derives its superiority from higher volatility and longer ignition delay of gasoline-like fuels combined with higher compression ratio typical of diesel engines. Here, the present study investigates the combustion process in a GCI engine operating with different injection strategies using computational fluid dynamics (CFD). Simulations are carried out on a single cylinder of a multi cylinder heavy-duty compression ignition engine, which operates at a compression ratio of 17:1 and an engine speed of 1038 rev/min. Two different injection strategies viz., late injection (LI), and early pilot injection (EP) are investigated to understand their impact on combustion and performance of the engine. Renormalized group (RNG) k-ε model is used to describe in-cylinder turbulence and KH RT model is used to simulate the fuel spray breakup. The developed CFD methodology is validated against relevant experimental data under a wide range of operating conditions for each injection strategy. The developed CFD methodology was found to capture the engine combustion behavior quite well. Based on the validated CFD model, the differences in the progress of combustion event for the two injection strategies is highlighted. It was found that a larger pilot fuel mass fraction results in a steeper rise in the initial heat release rate which in turn influences the transition to mixing controlled combustion. In line with the experimental data, the study showed that the late pilot injection strategy with three injection pulses, results in higher performance compared to the other conditions.

33 ADVANCED PROPULSION SYSTEMS↗
NASA NTRS2004

HSCT Sector Combustor Evaluations for Demonstration Engine

In LET Task 10, critical development issues of the HSCT lean-burn low emissions combustor were addressed with a range of engineering tools. Laser diagnostics and CFD analysis were applied to develop a clearer understanding of the fuel-air premixing process and premixed combustion. Subcomponent tests evaluated the emissions and operability performance of the fuel-air premixers. Sector combustor tests evaluated the performance of the integrated combustor system. A 3-cup sector was designed and procured for laser diagnostics studies at NASA Glenn. The results of these efforts supported the earlier selection of the Cyclone Swirler as the pilot stage premixer and the IMFH (Integrated Mixer Flame Holder) tube as the main stage premixer of the LPP combustor. In the combustor system preliminary design subtask, initial efforts to transform the sector combustor design into a practical subscale engine combustor met with significant challenges. Concerns about the durability of a stepped combustor dome and the need for a removable fuel injection system resulted in the invention and refinement of the MRA (Multistage Radial Axial) combustor system in 1994. The MRA combustor was selected for the HSR Phase II LPP subscale combustor testing in the CPC Program.

Greenfield, Stuart↗
DOE OSTI2021

Fuel Property Effects on Knock Propensity and Thermal Efficiency in a Direct-Injection Spark-Ignition Engine

Engine knock remains one of the major barriers to further improvement in thermal efficiency of Direct-Injection Spark-Ignition (DISI) engines. While Research Octane Number and Motor Octane Number are often used as standard rating methods for knock resistance of fuels, the impacts of other fuel properties on knock propensity in modern engines such as heat of vaporization (HoV) and laminar flame speed (LFS) require better understanding in order to co-optimize fuels and engine designs to achieve higher thermal efficiency and lower CO2 emission. In the present study, computational fluid dynamics (CFD) is used to model a boosted DISI engine with a focus on knock prediction and fuel property effects. A level-set G-equation model is employed to capture turbulent premixed combustion, and is coupled with a transported Livengood-Wu (L-W) integral approach to predict autoignition in the unburnt region. A criterion associated with the L-W integral is developed to accurately predict knock onset and knock-limited spark-advance. This model is then applied to a sensitivity analysis of HoV and LFS on knock tendency and thermal efficiency. The pressure-temperature trajectory framework is applied and extended to study the fuel effects on auto-ignition process in the engine. An existing efficiency-based merit function, which is derived from experiments for boosted SI engines, is evaluated and improved based on the current CFD results.

DISI↗
NASA NTRSDate not supplied

Assessment of Coupling One-Dimension and Three-Dimension CFD in Conjugate Heat Transfer Modeling and Simulations for Cooling of Turbine Engine

The CHT (conjugate heat transfer) modeling and simulation work fits under a bigger area in the TTT (Transformational Tools and Technologies) project called “Cooling Tools and Technologies for Combustor and Turbine.” The existing validation data sets such as TTT GRC (Glenn Research Center) THX (Turbulent Heat FLUX) IV experiments will be used to assess the solid/fluid heat transfer development and validation. The experimental configuration was a subsonic convergent nozzle with a square exit blowing a heated exhaust over an effusion-cooled deck. The simulation solutions were obtained and compared to experimental measurements of the profiles of flow velocities and temperatures, as well as deck surface temperatures. Cooling air blowing ratios of 0, 1, and 2 were considered. The three-dimensional cooling tubes were modelled by the medial axis of the cylindrical geometries due to large length over width ratios. It represented a typical coupling of one-dimension and three-dimension simulations of conjugate heat transfer approaches. The surface temperatures on the solid surfaces and the profiles of velocities and temperatures near the surfaces compared to the experimental data were presented.

nonlinear sub-filter stresses↗
NASA NTRS1997

Bleed Hole Flow Phenomena Studied

Boundary-layer bleed is an invaluable tool for controlling the airflow in supersonic aircraft engine inlets. Incoming air is decelerated to subsonic speeds prior to entering the compressor via a series of oblique shocks. The low momentum flow in the boundary layer interacts with these shocks, growing in thickness and, under some conditions, leading to flow separation. To remedy this, bleed holes are strategically located to remove mass from the boundary layer, reducing its thickness and helping to maintain uniform flow to the compressor. The bleed requirements for any inlet design are unique and must be validated by extensive wind tunnel testing to optimize performance and efficiency. To accelerate this process and reduce cost, researchers at the NASA Lewis Research Center initiated an experimental program to study the flow phenomena associated with bleed holes. Knowledge of these flow properties will be incorporated into computational fluid dynamics (CFD) models that will aid engine inlet designers in optimizing bleed configurations before any hardware is fabricated. This ongoing investigation is currently examining two hole geometries, 90 and 20 (both with 5-mm diameters), and various flow features.

Source record↗
DOE OSTI2022

Numerical and experimental investigation of the flame kernel growth in a methane/air mixture near the lean flammability limit

Lean combustion has the potential to improve the thermal efficiency of spark-ignition engines, but it faces the significant challenge of increased cycle-to-cycle variation due to low mixture reactivity and unstable flame dynamics. Computational fluid dynamics (CFD) employing predictive models can guide engine design and optimize operating strategies for lean combustion. However, ignition and combustion models have rarely been validated at fuel-lean conditions, and a fundamental understanding of the early flame kernel growth process is also lacking for a successful sub-model development. Here, the present study develops a numerical simulation framework used to investigate early flame kernel growth in methane/air mixtures. A nanosecond-pulsed discharge (NPD) approach is employed to effectively decouple the flame kernel growth from the electrical discharge due to their difference in timescales, and equivalence ratios near the experimentally measured lean flammability limit (LFL) are selected to focus on challenging mixture conditions. Three numerical investigations, such as the choice of turbulence modeling, grid size, and grid control strategies, are examined to match both LFL and flame kernel structure measured from experiments. It is demonstrated that a quasi-direct numerical simulation (QDNS) with a fixed grid embedding of 10 μm can predict the LFL as φ CFD =0.61 and match the displacement speed of the kernel’s boundary marked in schlieren images. To predict the LFL and flame kernel shape, a fine grid (Δ≤12.5 μm) is needed to capture the consumption of formaldehyde (CH 2 O) in kernel’s reaction branches attached to the anode, and adaptive mesh refinement is replaced with the fixed embedding due to loss of simulation accuracy. Also, it is found that a large-eddy simulation (LES) using the Dynamic Structure model is not suitable for the NPD-induced flame kernel simulation because artificial sub-grid turbulent kinetic energy induced by shock dynamics alters the flow velocity calculation, resulting in divergence of LES from QDNS. Lastly, the simulation well matches the experimental data for the flame kernel evolution in three mixture conditions (φ = 0.7, 0.61, 0.55), showing toroidal flame kernel expansion and flame kernel growth/extinction.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

OC6 Phase IV: Validation of CFD Models for Stiesdal TetraSpar Floating Offshore Wind Platform

ABSTRACT With only a few floating offshore wind turbine (FOWT) farms deployed anywhere in the world, FOWT technology is still in its infancy, building on a modicum of real‐world experience to advance the nascent industry. To support further development, engineers rely heavily on modeling tools to accurately portray the behavior of these complex systems under realistic environmental conditions. This reliance creates a need for verification and validation of such tools to improve reliability of load and dynamic response prediction and analysis capabilities of FOWT systems. The Offshore Code Comparison Collaboration, Continued with Correlation and unCertainty (OC6) project was created under the framework of the International Energy Agency to address this need and considers a three‐sided verification and validation between engineering level models, computational fluid dynamics (CFD), and experimental results. In this paper, a novel floating offshore wind platform, the Stiesdal TetraSpar, is simulated using CFD under the load conditions defined by Phase IV of the OC6 project. The comparison of these CFD results against the experimental results demonstrated the ability to predict the platform response to waves when imposing the measured wave signals as input. Although validation versus experiment was largely successful, the damping behavior was impacted by uncertainties likely originating from the mooring system and sensor umbilical cable. This extensive comparison effort with multiple CFD practitioners offers insight into best practices to achieve reliable results.

17 WIND ENERGY↗
NASA NTRS1992

Application of computational fluid dynamics to the design of the film cooled STME subscale nozzle for the National Launch System

The status of computational fluid dynamics (CFD) calculations for the Space Transportation Main Engine (STME) film/dump cooled nozzle design is presented, with an emphasis on the timely impact of CFD on the design of the sub-scale nozzle coolant system. The following aspects of the sub-scale coolant delivery system were analyzed with CFD: 1) a design trade study of a mechanical flow splitting device for uniform distribution of the subsonic cavity flow, 2) a design trade study of the subsonic cavity lip to achieve film integrity, and 3) an analysis of the primary flow interaction with the core/secondary coolant streams. All design calculations were performed with the Generalized Aerodynamic Simulation Program (GASP), a 3-D, multi-block, generalized Navier-Stokes code capable of solving with frozen, finite-rate or equilibrium chemical kinetics. The initial design of the subsonic cavity flow used square posts to distribute the sonic orifice jets into a uniform flow. Calculations for this design indicated that an unacceptable mal-distribution of film occurred. Design modifications involving curved and slotted posts were computed in an effort to uniformly distribute the secondary coolant flow. Analysis of these configurations showed that although the flowfield improved in uniformity, it was still unacceptable, especially at higher feed pressures. Results from these studies were then incorporated into a design that resulted in the insertion of a porous metal ring into the subsonic cavity. Subsequent water flow model studies showed that this concept was successful in uniformly distributing flow exiting the cavity. In addition to the design of the subsonic cavity, CFD was also used to analyze the secondary coolant lip and the primary flow interaction with the core/secondary coolant streams. A series of calculations were first performed to modify the subsonic cavity lip contour. The flow over the modified lip was then computed simultaneously with the primary injectors to determine the impact of the subsonic coolant stream on the primary slot jets. Pressure, temperature, velocity, and coolant mass fraction contours are presented for these configurations. Information is given in viewgraph form.

Garrett, Joseph L.↗