Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow lines”

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 847 records · Page 47

Validation of Universal Cryogenic Flow Boiling Correlations in Thermal Desktop for Liquid Hydrogen

Developing accurate models of two-phase cryogenic flow will reduce risk and reduce margins for future NASA vehicles such as the Nuclear Thermal Propulsion (NTP) system. Currently there is a need for more accurate, direct cryogenic data-anchored models for various boiling and two-phase phenomena. The focus of the current work is on modeling steady state cryogenic flow boiling in the transfer line that connects a propellant storage tank to an engine or customer receiver tank. This paper presents Thermal Desktop model validation results against several historical liquid hydrogen (LH2) heated tube experiments that cover a wide range of inlet conditions, mass flux, and heat flux. New universal cryogenic flow boiling correlations were recently developed and anchored to the largest cryogenic flow boiling database assembled in the world to-date. The new universal correlations are patched together and implemented into Thermal Desktop to compare model performance against the built-in flow boiling correlations. Axial wall temperature predictions as a function of preponderant parameters like heat flux, mass flux, inlet pressure, inlet quality, tube diameter, etc. using the two sets of flow boiling correlations are compared against the experimental data. Overall, the Thermal Desktop model with the new universal cryogenic flow boiling correlations demonstrates an improvement in predictive performance over the Thermal Desktop model using built-in correlations for both wall temperature and location of critical heat flux compared to the data for LH2 flow boiling in heated tubes.

cryogenic propellant transfer↗

Multiple Velocity Profile Measurements in Hypersonic Flows using Sequentially-Imaged Fluorescence Tagging

Nitric-oxide planar laser-induced fluorescence (NO PLIF) was used to perform velocity measurements in hypersonic flows by generating multiple tagged lines which fluoresce as they convect downstream. For each laser pulse, a single interline, progressive scan intensified CCD camera was used to obtain separate images of the initial undelayed and delayed NO molecules that had been tagged by the laser. The CCD configuration allowed for sub-microsecond acquisition of both images, resulting in sub-microsecond temporal resolution as well as sub-mm spatial resolution (0.5-mm x 0.7-mm). Determination of axial velocity was made by application of a cross-correlation analysis of the horizontal shift of individual tagged lines. Quantification of systematic errors, the contribution of gating/exposure duration errors, and influence of collision rate on fluorescence to temporal uncertainty were made. Quantification of the spatial uncertainty depended upon the analysis technique and signal-to-noise of the acquired profiles. This investigation focused on two hypersonic flow experiments: (1) a reaction control system (RCS) jet on an Orion Crew Exploration Vehicle (CEV) wind tunnel model and (2) a 10-degree half-angle wedge containing a 2-mm tall, 4-mm wide cylindrical boundary layer trip. The experiments were performed at the NASA Langley Research Center's 31-inch Mach 10 wind tunnel.

Bathel, Brett F.↗

Automated Extraction of Secondary Flow Features

The use of Computational Fluid Dynamics (CFD) has become standard practice in the design and development of the major components used for air and space propulsion. To aid in the post-processing and analysis phase of CFD many researchers now use automated feature extraction utilities. These tools can be used to detect the existence of such features as shocks, vortex cores and separation and re-attachment lines. The existence of secondary flow is another feature of significant importance to CFD engineers. Although the concept of secondary flow is relatively understood there is no commonly accepted mathematical definition for secondary flow. This paper will present a definition for secondary flow and one approach for automatically detecting and visualizing secondary flow.

Dorney, Suzanne M.↗

An Observational Evaluation of RKW Theory over the U.S. Southern Great Plains

The theory of Rotunno et al. (“RKW” theory) addresses the behavior of squall-line cold pools in vertically sheared flows. It predicts that, within a given thermodynamic environment, a balance between baroclinic vorticity generation by the cold pool and low-level environmental vertical wind shear induces an upright updraft along the gust front that maximizes the initiation of new convective cells. Although this theory has been evaluated numerically, its applicability to observed systems remains unclear and is limited by a lack of critical measurements, including high-frequency thermodynamic and wind profiles across the gust front. Herein, observations from the Atmospheric Radiation Measurement Southern Great Plains (ARM-SGP) observatory near Lamont, Oklahoma, are used to evaluate RKW theory for 10 well-observed squall lines over a 11-yr period. For this evaluation, RKW parameters including cold-pool intensity (c), low-level ambient, line-normal vertical shear (ΔV n ), subcloud and cloud-layer updraft tilts, and multiple measures of system intensity are estimated. Furthermore, the c estimates rely on thermodynamic retrievals from the Atmosphere Emitted Radiance Interferometer (AERI), which are uncertain but verify reasonably well against independent observations. As predicted by the theory, for c/ΔV n ≥ 1, c/ΔV n correlates positively with updraft tilt and negatively with system intensity, but these results are not always statistically significant and are also sensitive to the method by which ΔV n is evaluated. Specifically, ΔV n evaluations that extend above the cold-pool top yield greater consistency with RKW predictions. Also, some measures of intensity correlate more strongly with standard moist instability metrics than with RKW parameters.

Cold pools↗

Development of a recombinant anti-VEGFR2-EPCAM bispecific antibody to improve antiangiogenic efficiency

Highlights: • Targeted therapy can be a useful therapeutic method along with traditional methods for breast cancer treatment. • Decreasing Akt/Nf-Kβ signaling pathway is crucial target in order to decreasing angiogenesis in TNBC. • The recombinant antibody can have inhibitory effects on production IL6,8, MMP2, and MMP9. Tumor progression and metastasis, especially in invasive cancers (such as triple-negative breast cancer [TNBC]), depend on angiogenesis, in which vascular epithelial growth factor (VEGF)/vascular epithelial growth factor receptor [1] has a decisive role, followed by the metastatic spread of cancer cells. Although some studies have shown that anti-VEGFR2/VEGF monoclonal antibodies demonstrated favorable results in the clinic, this approach is not efficient, and further investigations are needed to improve the quality of cancer treatment. Besides, the increased expression of epithelial cell adhesion molecule (EpCAM) in various cancers, for instance, invasive breast cancer, contributes to angiogenesis, facilitating the migration of tumor cells to other parts of the body. Thus, the main goal of our study was to target either VEGFR2 or EpCAM as pivotal players in the progression of angiogenesis in breast cancer. Regarding cancer therapy, the production of bispecific antibodies is easier and more cost-effective compared to monoclonal antibodies, targeting more than one antigen or receptor; for this reason, we produced a recombinant antibody to target cells expressing EpCAM and VEGFR2 via a bispecific antibody to decrease the proliferation and metastasis of tumor cells. Following the cloning and expression of our desired anti-VEGFR2/EPCAM sequence in E. coli, the accuracy of the expression was confirmed by Western blot analysis, and its binding activities to VEGFR2 and EPCAM on MDA-MB-231 and MCF-7 cell lines were respectively indicated by flow cytometry. Then, its anti-proliferative potential was indicated by 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) and apoptosis assay to evaluate inhibitory effects of the antibody on tumor cells. Subsequently, the data indicated that migration, invasion, and angiogenesis were inhibited in breast cancer cell lines via the bispecific antibody. Furthermore, cytokine analysis indicated that the bispecific antibody could moderate interleukin 8 (IL-8) and IL-6 as key mediators in angiogenesis progression in breast cancer. Thus, our bispecific antibody could be considered as a promising candidate tool to decrease angiogenesis in TNBC.

60 APPLIED LIFE SCIENCES↗

Millennial–Scale Vulnerability of the Antarctic Ice Sheet to Regional Ice Shelf Collapse

The Antarctic Ice Sheet (AIS) remains the largest uncertainty in projections of future sea level rise. A likely climate–driven vulnerability of the AIS is thinning of floating ice shelves resulting from surface–melt–driven hydrofracture or incursion of relatively warm water into subshelf ocean cavities. The resulting melting, weakening, and potential ice shelf collapse reduces shelf buttressing effects. Upstream ice flow accelerates, causing thinning, grounding–line retreat, and potential ice sheet collapse. While high–resolution projections have been performed for localized Antarctic regions, full–continent simulations have typically been limited to low–resolution models. Here we quantify the vulnerability of the entire present–day AIS to regional ice shelf collapse on millennial timescales treating relevant ice flow dynamics at the necessary ~1–km resolution. Collapse of any of the ice shelves dynamically connected to the West Antarctic Ice Sheet (WAIS) is sufficient to trigger ice sheet collapse in marine–grounded portions of the WAIS. As a result, vulnerability elsewhere appears limited to localized responses.

54 ENVIRONMENTAL SCIENCES↗

Damping inter-area oscillation using reinforcement learning controlled TCSC

Appropriate damping control of the inter-area oscillation would ensure the grid stability and maintain the tie-line power transfer capability. In this paper, we proposed a novel reinforcement learning (RL) based power oscillation damping (POD) controller that uses Thyristor Controlled Series Compensators (TCSC) to damp inter-area oscillations. By leveraging the unbiased gradient direction estimation of the natural evolution strategy (NES), we successfully regulated the power flows on the tie-lines and damped inter-area oscillations through dynamically modulating the inserted reactance of the TCSC. The proposed RL-based POD controller has been tested on both two-area four-machine system and North American Minni-WECC system. Extensive studies have demonstrated the excellent performance of the proposed RL-based TCSC POD controller in damping inter-area oscillations.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Transmutation Scoping Studies for a Chloride Molten Salt Reactor

Over the past few years, there has been strong renewed interest from private industry, mostly from start-up enterprises, in molten salt reactor (MSR) technologies because of the unique properties of this class of reactors. These are reactors in which the fuel is homogeneously mixed with the coolant in the form of liquid salts and is circulated continuously into and out of the active core region with on-line fuel management, salt treatment, and salt processing. In response to such wide-spread interest, Argonne National Laboratory is expanding its well-established reactor modelling and simulation expertise and infrastructure to enable detailed analysis and design of MSRs. The tools being developed are able to simulate the continuous fuel flow, the complex on-line fuel management and elemental removal processes (e.g., fission product removal) using depletion steps representative of a real MSR system. Leveraging these capabilities, a parametric study on the transmutation performance of a simplified actinide-burning MSR concept that uses a chloride-based salt was performed. This type of salt has attracted attention over the more commonly discussed fluoride-based salts since no tritium is produced as a result of irradiation and it is compatible with a fast neutron spectrum. The studies discussed in this paper examine the performance of a burner MSR design with a fixed core size and power density over a range of possible fuel salt molar ratios with NaCl-MgCl2 as the carrier salt. The intent is to quantify the impact on the required transuranics content of the make-up fuel, the actinide transmutation rates, and other performance characteristics for typical burner MSR designs.

Heidet, Florent↗

Assessing Uncertainties in Boundary Layer Transition Predictions for HIFiRE-1 at Non-zero Angles of Attack

Boundary layer stability was analyzed for the HIFiRE-1 flight vehicle geometry for ground tests conducted at the CUBRC LENS I hypersonic shock test facility and the Langley Research Center (LaRC) 20- inch Mach 6 Tunnel. Boundary layer stability results were compared to transition onset location obtained from discrete heat transfer measurements from thin film gauges during the CUBRC test and spatially continuous heat transfer measurements from thermal phosphor paint data during the LaRC test. The focus of this analysis was on conditions at non-zero angles of attack as stability analysis has already been performed at zero degrees angle of attack. Also, the transition onset data obtained during flight testing was at nonzero angles of attack, so this analysis could be expanded in the future to include the results of the flight test data. Stability analysis was performed using the 2D parabolized stability software suite STABL (Stability and Transition Analysis for Hypersonic Boundary Layers) developed at the University of Minnesota and the mean flow solutions were computed using the DPLR finite volume Navier-Stokes computational fluid dynamics (CFD) solver. A center line slice of the 3D mean flow solution was used for the stability analysis to incorporate the angle of attack effects while still taking advantage of the 2D STABL software suite. The N-factors at transition onset and the value of Re(sub theta)/M(sub e), commonly used to predict boundary layer transition onset, were compared for all conditions analyzed. Ground test data was analyzed at Mach 7.2 and Mach 6.0 and angles of attack of 1deg, 3deg and 5deg. At these conditions, the flow was found to be second mode dominant for the HIFiRE-1 slender cone geometry. On the leeward side of the vehicle, a strong trend of transition onset location with angle of attack was observed as the boundary layer on the leeward side of the vehicle developed inflection points at streamwise positions on the vehicle that correlated to angle of attack. Inflection points are a strong instability mechanism that lead to rapid breakdown and transition to turbulence. The transition onset location on the windward side of the vehicle displayed no trend with angle of attack or freestream Reynolds number and transition was observed farther down the vehicle than observed on the leeward side of the vehicle. In analysis of both windward and leeward sides of the vehicle, use of the N factor methodology to develop trends to predict boundary layer transition onset showed improvements over the Re(sub theta)/M(sub e) empirical correlation methodology. Stronger correlations and less scatter in the data were observed when using the N factor method for these cases.

Marek, Lindsay C.↗

A Modular Approach to Model Oscillating Control Surfaces Using Navier Stokes Equations

The use of active controls for rotorcraft is becoming more important for modern aerospace configurations. Efforts to reduce the vibrations of helicopter blades with use of active-controls are in progress. Modeling oscillating control surfaces using the linear aerodynamics theory is well established. However, higher-fidelity methods are needed to account for nonlinear effects, such as those that occur in transonic flow. The aeroelastic responses of a wing with an oscillating control surface, computed using the transonic small perturbation (TSP) theory, have been shown to cause important transonic flow effects such as a reversal of control surface effectiveness that occurs as the shock wave crosses the hinge line. In order to account for flow complexities such as blade-vortex interactions of rotor blades higher-fidelity methods based on the Navier-Stokes equations are used. Reference 6 presents a procedure that uses the Navier-Stokes equations with moving-sheared grids and demonstrates up to 8 degrees of control-surface amplitude, using a single grid. Later, this procedure was extended to accommodate larger amplitudes, based on sliding grid zones. The sheared grid method implemented in EulerlNavier-Stokes-based aeroelastic code ENS AERO was successfully applied to active control design by industry. Recently there are several papers that present results for oscillating control surface using Reynolds Averaged Navier-Stokes (RANS) equations. References 9 and 10 report 2-D cases by filling gaps with overset grids. Reference 9 compares integrated forces with the experiment at low oscillating frequencies whereas Ref. 10 reports parametric studies but with no validation. Reference II reports results for a 3D case by modeling the gap region with a deformed grid and compares force results with the experiment only at the mid-span of flap. In Ref. II grid is deformed to match the control surface deflections at the section where the measurements are made. However, there is no indication in Ref. II that the gaps are explicitly modeled as in Ref. 6. Computations using overset grids are reported in Ref. 12 for a case by adding moving control surface to an existing blade but with no validation either with an experiment or another computation.

Navier↗

Synthetic Absorption Lines for a Clumpy Medium: A Spectral Signature for Cloud Acceleration in AGN?

There is increasing evidence that the highly ionized multiphase components of AGN disc winds may be due to thermal instability. The ions responsible for forming the observed X-ray absorption lines may only exist in relatively cool clumps that can be identified with the so-called warm absorbers. Here we calculate synthetic absorption lines for such warm absorbers from first principles by combining 2D hydrodynamic solutions of a two-phase medium with a dense grid of photoionization models to determine the detailed ionization structure of the gas. Our calculations reveal that cloud disruption, which leads to a highly complicated velocity field (i.e. a clumpy flow), will only mildly affect line shapes and strengths when the warm gas becomes highly mixed but not depleted. Prior to complete disruption, clouds that are optically thin to the driving UV resonance lines will cause absorption at an increasingly blueshifted line-of-sight velocity as they are accelerated. This behavior will imprint an identifiable signature on the line profile if warm absorbers are enshrouded in an even broader absorption line produced by a high column of intercloud gas. Interestingly, we show that it is possible to develop a spectral diagnostic for cloud acceleration by differencing the absorption components of a doublet line, a result that can be qualitatively understood using a simple partial covering model. Our calculations also permit us to comment on the spectral differences between cloud disruption and ionization changes driven by flux variability. Notably, cloud disruption offers another possibility for explaining absorption line variability.

line: formation – galaxies: nuclei – quasars: ↗

DoE Plasma Center for Momentum Transport and Flow Self-Organization in Plasmas: Non-linear Emergent Structure Formation in magnetized Plasmas and Rotating Magnetofluids

This report covers the UW-Madison activities that took place within a larger DoE Center Administered and directed by Professor George Tynan at the University of California, San Diego. The work at Wisconsin will also be covered in the final reporting for the entire center, which will be submitted by UCSD. There were two main activities, one experimental and one that was theoretical in nature, as part of the Center activities at the University of Wisconsin, Madison. First, the Center supported an experimentally focused postdoc (Chris Cooper) to carry out fundamental studies of momentum transport in rotating and weakly magnetized plasma. His experimental work was done on the Plasma Couette Experiment, a cylindrical plasma confinement device, with a plasma flow created through electromagnetically stirring plasma at the plasma edge facilitated by arrays of permanent magnets. Cooper's work involved developing optical techniques to measure the ion temperature and plasma flow through Doppler-shifted line radiation from the plasma argon ions. This included passive emission measurements and development of a novel ring summing Fabry-Perot spectroscopy system, and the active system involved using a diode laser to induce fluorescence. On the theoretical side, CMTFO supported a postdoc (Johannes Pueschel) to carry out a gyrokinetic extension of residual zonal flow theory to the case with magnetic fluctuations, showing that magnetic stochasticity disrupts zonal flows. The work included a successful comparison with gyrokinetic simulations. This work and its connection to the broader CMTFO will be covered more thoroughly in the final CMTFO report from Professor Tynan.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multiple flow velocities in the transition region

Observations with the High Resolution Telescope and Spectrograph (HRTS) of the C IV lines at 1550 A above sunspots reveal the coexistence of several distinct high velocity gas components within the instrument resolution element of 1 x 1 arcsec. A further investigation of the HRTS material demonstrates that such multiple gas velocities are not restricted to the temperature region around 100,000 K, where the C IV lines are formed, but that the gas flows with essentially unchanged speed through all temperature layers from 250,000 K to 10,000 K. Furthermore, the phenomenon is not restricted to sunspots, but occurs also in other solar regions, often most easily distinguished in active regions. The characteristics of the velocity fields are described. The mass fluxes are analyzed, and the requirements they place on the resolving power of future space instruments are discussed.

Brekke, P.↗

A time accurate finite volume method for propulsion chamber flows

An implicit three-dimensional time-accurate method for propulsion chamber flows is proposed which uses line Gauss-Seidel relaxation and multiple axial sweeps for the convergence of each time step. The general time-integration algorithm employed includes such schemes as the Euler implicit method. The results of spatial and temporal accuracy tests reveal that Roe's (1981) flux difference splitting provides excellent tracking of acoustic wave speeds. In comparison with other methods, no low mean flow Mach number convergence limitation or Courant number stabilization restriction is observed.

Beddini, R. A.↗

Computational Aerodynamic Simulations of an 840 ft/sec Tip Speed Advanced Ducted Propulsor Fan System Model for Acoustic Methods Assessment and Development

Computational Aerodynamic simulations of an 840 ft/sec tip speed, Advanced Ducted Propulsor fan system were performed at five different operating points on the fan operating line, in order to provide detailed internal flow field information for use with fan acoustic prediction methods presently being developed, assessed and validated. The fan system is a sub-scale, lownoise research fan/nacelle model that has undergone extensive experimental testing in the 9- by 15- foot Low Speed Wind Tunnel at the NASA Glenn Research Center, resulting in quality, detailed aerodynamic and acoustic measurement data. Details of the fan geometry, the computational fluid dynamics methods, the computational grids, and various computational parameters relevant to the numerical simulations are discussed. Flow field results for three of the five operating conditions simulated are presented in order to provide a representative look at the computed solutions. Each of the five fan aerodynamic simulations involved the entire fan system, excluding a long core duct section downstream of the core inlet guide vane. As a result, only fan rotational speed and system bypass ratio, set by specifying static pressure downstream of the core inlet guide vane row, were adjusted in order to set the fan operating point, leading to operating points that lie on a fan operating line and making mass flow rate a fully dependent parameter. The resulting mass flow rates are in good agreement with measurement values. The computed blade row flow fields for all five fan operating points are, in general, aerodynamically healthy. Rotor blade and fan exit guide vane flow characteristics are good, including incidence and deviation angles, chordwise static pressure distributions, blade surface boundary layers, secondary flow structures, and blade wakes. Examination of the computed flow fields reveals no excessive boundary layer separations or related secondary-flow problems. A few spanwise comparisons between computational and measurement data in the bypass duct show that they are in good agreement, thus providing a partial validation of the computational results.

Tweedt, Daniel L.↗

Computational Aerodynamic Simulations of a 1484 ft/sec Tip Speed Quiet High-Speed Fan System Model for Acoustic Methods Assessment and Development

Computational Aerodynamic simulations of a 1484 ft/sec tip speed quiet high-speed fan system were performed at five different operating points on the fan operating line, in order to provide detailed internal flow field information for use with fan acoustic prediction methods presently being developed, assessed and validated. The fan system is a sub-scale, low-noise research fan/nacelle model that has undergone experimental testing in the 9- by 15-foot Low Speed Wind Tunnel at the NASA Glenn Research Center. Details of the fan geometry, the computational fluid dynamics methods, the computational grids, and various computational parameters relevant to the numerical simulations are discussed. Flow field results for three of the five operating points simulated are presented in order to provide a representative look at the computed solutions. Each of the five fan aerodynamic simulations involved the entire fan system, which includes a core duct and a bypass duct that merge upstream of the fan system nozzle. As a result, only fan rotational speed and the system bypass ratio, set by means of a translating nozzle plug, were adjusted in order to set the fan operating point, leading to operating points that lie on a fan operating line and making mass flow rate a fully dependent parameter. The resulting mass flow rates are in good agreement with measurement values. Computed blade row flow fields at all fan operating points are, in general, aerodynamically healthy. Rotor blade and fan exit guide vane flow characteristics are good, including incidence and deviation angles, chordwise static pressure distributions, blade surface boundary layers, secondary flow structures, and blade wakes. Examination of the computed flow fields reveals no excessive or critical boundary layer separations or related secondary-flow problems, with the exception of the hub boundary layer at the core duct entrance. At that location a significant flow separation is present. The region of local flow recirculation extends through a mixing plane, however, which for the particular mixing-plane model used is now known to exaggerate the recirculation. In any case, the flow separation has relatively little impact on the computed rotor and FEGV flow fields.

Tweedt, Daniel L.↗

A CFD-based aerodynamic design procedure for hypersonic wind-tunnel nozzles

A new procedure which unifies the best of current classical design practices, computational fluid dynamics (CFD), and optimization procedures is demonstrated for designing the aerodynamic lines of hypersonic wind-tunnel nozzles. The new procedure can be used to design hypersonic wind tunnel nozzles with thick boundary layers where the classical design procedure has been shown to break down. An efficient CFD code, which solves the parabolized Navier-Stokes (PNS) equations using an explicit upwind algorithm, is coupled to a least-squares (LS) optimization procedure. A LS problem is formulated to minimize the difference between the computed flow field and the objective function, consisting of the centerline Mach number distribution and the exit Mach number and flow angle profiles. The aerodynamic lines of the nozzle are defined using a cubic spline, the slopes of which are optimized with the design procedure. The advantages of the new procedure are that it allows full use of powerful CFD codes in the design process, solves an optimization problem to determine the new contour, can be used to design new nozzles or improve sections of existing nozzles, and automatically compensates the nozzle contour for viscous effects as part of the unified design procedure. The new procedure is demonstrated by designing two Mach 15, a Mach 12, and a Mach 18 helium nozzles. The flexibility of the procedure is demonstrated by designing the two Mach 15 nozzles using different constraints, the first nozzle for a fixed length and exit diameter and the second nozzle for a fixed length and throat diameter. The computed flow field for the Mach 15 least squares parabolized Navier-Stokes (LS/PNS) designed nozzle is compared with the classically designed nozzle and demonstrates a significant improvement in the flow expansion process and uniform core region.

Korte, John J.↗

Optimization of Time-Dependent Particle Tracing Using Tetrahedral Decomposition

An efficient algorithm is presented for computing particle paths, streak lines and time lines in time-dependent flows with moving curvilinear grids. The integration, velocity interpolation and step-size control are all performed in physical space which avoids the need to transform the velocity field into computational space. This leads to higher accuracy because there are no Jacobian matrix approximations or expensive matrix inversions. Integration accuracy is maintained using an adaptive step-size control scheme which is regulated by the path line curvature. The problem of cell-searching, point location and interpolation in physical space is simplified by decomposing hexahedral cells into tetrahedral cells. This enables the point location to be done analytically and substantially faster than with a Newton-Raphson iterative method. Results presented show this algorithm is up to six times faster than particle tracers which operate on hexahedral cells yet produces almost identical particle trajectories.

Kenwright, David↗