Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “cavity flow”

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 1,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Simulations of Cavitating Cryogenic Inducers

Simulations of cavitating turbopump inducers at their design flow rate are presented. Results over a broad range of Nss, numbers extending from single-phase flow conditions through the critical head break down point are discussed. The flow characteristics and performance of a subscale geometry designed for water testing are compared with the fullscale configuration that employs LOX. In particular, thermal depression effects arising from cavitation in cryogenic fluids are identified and their impact on the suction performance of the inducer quantified. The simulations have been performed using the CRUNCH CFD[R] code that has a generalized multi-element unstructured framework suitable for turbomachinery applications. An advanced multi-phase formulation for cryogenic fluids that models temperature depression and real fluid property variations is employed. The formulation has been extensively validated for both liquid nitrogen and liquid hydrogen by simulating the experiments of Hord on hydrofoils; excellent estimates of the leading edge temperature and pressure depression were obtained while the comparisons in the cavity closure region were reasonable.

Dorney, Dan↗
NASA NTRS2010

Design and Analysis of Outer Mold Line Close-outs for the Max Launch Abort System (MLAS) Flight Experiment

In 2007, the NASA Exploration Systems Mission Directorate (ESMD) chartered the NASA Engineering Safety Center (NESC) to demonstrate an alternate launch abort concept as risk mitigation for the Orion project's baseline "tower" design. On July 8, 2009, a full scale, passive aerodynamically stabilized Max Launch Abort System (MLAS) pad abort demonstrator was successfully launched from NASA Goddard Space Flight Center's Wallops Flight Facility. Aerodynamic close-outs were required to cover openings on the MLAS fairing to prevent aerodynamic flow-through and to maintain the MLAS OML surface shape. Two-ply duct tape covers were designed to meet these needs. The duct tape used was a high strength fiber reinforced duct tape with a rubberized adhesive that demonstrated 4.6 lb/in adhesion strength to the unpainted fiberglass fairing. Adhesion strength was observed to increase as a function of time. The covers were analyzed and experimentally tested to demonstrate their ability to maintain integrity under anticipated vehicle ascent pressure loads and to not impede firing of the drogue chute mortars. Testing included vacuum testing and a mortar fire test. Tape covers were layed-up on thin Teflon sheets to facilitate installation on the vehicle. Custom cut foam insulation board was used to fill mortar hole and separation joint cavities and provide support to the applied tape covers. Flight test results showed that the tape covers remained adhered during flight.

Woods-Vedeler, Jessica A.↗
DOE OSTI2023

Embedding ceramic components in metal structures with hybrid directed energy deposition

The combined benefit of both additive and subtractive manufacturing within the same gantry system enables hybrid directed energy deposition to create complex geometries with smooth surface finish and superior dimensional accuracy. Moreover, with layer-by-layer access to the structure during both the addition and subtraction of material, the insertion of components is now possible, assuming the components can survive the high temperatures associated with the subsequent metal deposition. Ceramic inserts are of interest for a variety of reasons including (1) to create complex interwoven ductile/brittle composites for ballistics or high-temperature applications or (2) to integrate high-temperature strain or temperature sensors protected within ceramic substrate subsumed into a larger metal structure. Here, in this work, stainless steel substrates were machined to create an internal cavity for the insertion of a ceramic component. During the investigation of several different over-the-ceramic deposition strategies, components were inserted, and different process sequences were allowed to continue to envelop the inserted ceramic with varying success. Unmelted powder was used to serve both as a thermal buffer and to provide a flush surface upon which the laser cladding could continue. Subsequent depositions were attempted with both dry and wet powder (addition of machining coolant to wet). The wet powder has previously been demonstrated to not significantly impact the mechanical properties of a final structure and provided a thermal barrier to protect the ceramic piece from the extreme temperatures of the final metal deposition. The wetting of the powder provided stability and minimized displacement caused by the powder flow from the laser cladding head. Finally, the use of an oblique angle for laser cladding allowed for the redirection of some fraction of the introduced thermal energy away from the ceramic component and, consequently, improved the survival of the ceramic inserts. With this combination of techniques, ceramic inserts survived full embedding within a 3D-printed stainless steel structure.

36 MATERIALS SCIENCE↗
NASA NTRS1988

On the origin of hot diamagnetic cavities near the earth's bow shock

The origin of hot diamagnetic cavities (HDCs) observed occasionally upstream from the earth's bow shock is investigated by examining the results of November 16, 1977, observation, when four of these events occurred on a single day, as well as plasma and field data from that day. The results suggest that HDCs may form as a result of an unusually strong interaction between shock-reflected ions and the incoming solar wind. It is proposed that this interaction stems from a temporary and localized reflection of a larger-than-normal fraction of the incident ions, which is stimulated by sudden changes in the upstream field orientation; the consequences of such a temporary overreflection are found to be consistent with many of the observed features of HDCs, including the strong slowing, deflection, and heating of the flow, as well as the localization, internal recoveries, and occasional formation upstream from the shock itself.

Thomsen, M. F.↗
NASA NTRS1999

Glacier and Ice Shelves Studies Using Satellite SAR Interferometry

Satellite radar interferometry is a powerful technique to measure the surface velocity and topography of glacier ice. On ice shelves, a quadruple difference technique separates tidal motion from the steady creep flow deformation of ice. The results provide a wealth of information about glacier grounding lines , mass fluxes, stability, elastic properties of ice, and tidal regime. The grounding line, which is where the glacier detaches from its bed and becomes afloat, is detected with a precision of a few tens of meters. Combining this information with satellite radar altimetry makes it possible to measure glacier discharge into the ocean and state of mass balance with greater precision than ever before, and in turn provide a significant revision of past estimates of mass balance of the Greenland and Antarctic Ice Sheets. Analysis of creep rates on floating ice permits an estimation of basal melting at the ice shelf underside. The results reveal that the action of ocean water in sub-ice-shelf cavities has been largely underestimated by oceanographic models and is the dominant mode of mass release to the ocean from an ice shelf. Precise mapping of grounding line positions also permits the detection of grounding line migration, which is a fine indicator of glacier change, independent of our knowledge of snow accumulation and ice melting. This technique has been successfully used to detect the rapid retreat of Pine Island Glacier, the largest ice stream in West Antarctica. Finally, tidal motion of ice shelves measured interferometrically provides a modern, synoptic view of the physical processes which govern the formation of tabular icebergs in the Antarctic.

Rignot, Eric↗
NASA NTRS2004

A New Look at Stellar Outflows: Spitzer Observations of the HH 46/47 System

We present the Early Release Observations of the HH 46/47 system and HH 46 IRS 1 source, taken with the three instruments aboard the Spitzer Space Telescope. The optically invisible southwest lobe, driven by the HH 47C bow shock, is revealed in full detail by the Infrared Array Camera (IRAC) images and displays a 'loop'-like morphology. Both of the mid-infrared outflow lobes are narrower than those of CO flow. We believe that the combination of emission by H2 rotational lines [S(11)-S(4)] and some atomic lines, which fall within the IRAC passbands, are responsible for the bulk of the observed emission, although contributions from the 3.3, 6.2, and 7.7 micron polycyclic aromatic hydrocarbon emission bands cannot be ruled out. Weak spectral features corresponding to these emitters are present in the Infrared Spectrograph spectrum of the HH 47A bow shock. The spectrum of HH 46 IRS 1 shows remarkable similarities to those of high-mass protostars, which include the presence of H2O, CO2, CH4, and possibly NH3, CH3OH, and ices. The high ice abundances and the lack of signs of thermal processing indicate that these ices in the envelope are well shielded from the powerful outflow and its cavity. Emission from the Bok globule at 24 micron is detected and displays a similar structure to that observed at 8 micron.

Herbig-Haro objects↗
NASA NTRS1988

Heating requirements and nonadiabatic surface effects for a model in the NTF cryogenic wind tunnel

A theoretical study has been made of the severity of nonadiabatic surface conditions arising from internal heat sources within a model in a cryogenic wind tunnel. Local surface heating is recognized as having an effect on the development of the boundary layer, which can introduce changes in the flow about the model and affect the wind tunnel data. The geometry was based on the NTF Pathfinder I wind tunnel model. A finite element heat transfer computer code was developed and used to compute the steady state temperature distribution within the body of the model, from which the surface temperature distribution was extracted. Particular three dimensional characteristics of the model were represented with various axisymmetric approximations of the geometry. This analysis identified regions on the surface of the model susceptible to surface heating and the magnitude of the respective surface temperatures. It was found that severe surface heating may occur in particular instances, but could be alleviated with adequate insulating material. The heat flux through the surface of the model was integrated to determine the net heat required to maintain the instrumentation cavity at the prescribed temperature. The influence of the nonadiabatic condition on boundary layer properties and on the validity of the wind tunnel simulation was also investigated.

Macha, J. M.↗
NASA NTRS2015

Development of a Premixed Combustion Capability for Scramjet Combustion Experiments

Hypersonic air-breathing engines rely on scramjet combustion processes, which involve high speed, compressible, and highly turbulent flows. The combustion environment and the turbulent flames at the heart of these engines are difficult to simulate and study in the laboratory under well controlled conditions. Typically, wind-tunnel testing is performed that more closely approximates engine testing rather than a careful investigation of the underlying physics that drives the combustion process. The experiments described in this paper, along with companion data sets being developed separately, aim to isolate the chemical kinetic effects from the fuel-air mixing process in a dual-mode scramjet combustion environment. A unique fuel injection approach is taken that produces a nearly uniform fuel-air mixture at the entrance to the combustor. This approach relies on the precombustion shock train upstream of the dual-mode scramjet combustor. A stable ethylene flame anchored on a cavity flameholder with a uniformly mixed combustor inflow has been achieved in these experiments allowing numerous companion studies involving coherent anti-Stokes Raman scattering (CARS), particle image velocimetry (PIV), and planar laser induced fluorescence (PLIF) to be performed.

Rockwell, Robert D.↗
NASA NTRS2009

In-Situ Measurement of Hall Thruster Erosion Using a Fiber Optic Regression Probe

One potential life-limiting mechanism in a Hall thruster is the erosion of the ceramic material comprising the discharge channel. This is especially true for missions that require long thrusting periods and can be problematic for lifetime qualification, especially when attempting to qualify a thruster by analysis rather than a test lasting the full duration of the mission. In addition to lifetime, several analytical and numerical models include electrode erosion as a mechanism contributing to enhanced transport properties. However, there is still a great deal of dispute over the importance of erosion to transport in Hall thrusters. The capability to perform an in-situ measurement of discharge channel erosion is useful in addressing both the lifetime and transport concerns. An in-situ measurement would allow for real-time data regarding the erosion rates at different operating points, providing a quick method for empirically anchoring any analysis geared towards lifetime qualification. Erosion rate data over a thruster s operating envelope would also be useful in the modeling of the detailed physics inside the discharge chamber. There are many different sensors and techniques that have been employed to quantify discharge channel erosion in Hall thrusters. Snapshots of the wear pattern can be obtained at regular shutdown intervals using laser profilometry. Many non-intrusive techniques of varying complexity and sensitivity have been employed to detect the time-varying presence of erosion products in the thruster plume. These include the use quartz crystal microbalances, emission spectroscopy, laser induced flourescence, and cavity ring-down spectroscopy. While these techniques can provide a very accurate picture of the level of eroded material in the thruster plume, it is more difficult to use them to determine the location from which the material was eroded. Furthermore, none of the methods cited provide a true in-situ measure of erosion at the channel surface while the thruster is in operation (i.e. none yield a continuous channel erosion measurement). A recent fundamental sensor development effort has led to a novel regression, erosion, and ablation sensor technology (REAST). The REAST sensor allows for measurement of real-time surface erosion rates at a discrete surface location. The sensor was tested using a linear Hall thruster geometry (see Fig. 1), which served as a means of producing plasma erosion of a ceramic discharge chamber. The mass flow rate, discharge voltage, and applied magnetic field strength could be varied, allowing for erosion measurements over a broad thruster operating envelope. Results are presented demonstrating the ability of the REAST sensor to capture not only the insulator erosion rates but also changes in these rates as a function of the discharge parameters.

Polzin, Kurt↗
DOE OSTI2023

Code Benchmark of the HTTF Pressurized Conduction Cooldown Test Using SAM

The High Temperature Test Facility (HTTF) at Oregon State University is an integral system test facility to simulate postulated reactor transients of prismatic high-temperature gas-cooled reactors(HTGRs). A series of test campaigns was launched, providing abundant test data that could be used to benchmark reactor system analysis codes like the System Analysis Module (SAM). In this study, a SAM model of the facility is developed based on the two-dimensional (2D) ring model approach. All components including the ceramic matrix, graphite heaters, helium coolant channels, core barrel, upcomer, pressure vessel, and reactor cavity cooling system are modeled as concentric cylindrical rings. The model is used to simulate one of the benchmark problems-Pressurized Conduction Cooldown (PCC)-within the scope of the Organisation for Economic Co-operation and Development Nuclear Energy Agency International HTTF Benchmark. The simulations consist of two parts. In the first part, operating and boundary conditions as well as thermophysical properties of materials are specified for the benchmark problem. In this work, results from the first part will be used in code-to-code comparison. In the second part, the SAM model is used to simulate Test PG-27, which is the first PCC test carried out in the HTTF, with only two of the ten heater banks activated. The results in the second part are used for code-to-data comparison. Because the helium coolant flow rate is not measured in this facility, it is estimated using the input power and inlet/outlet coolant temperatures. Additionally, radial heat flow in the ceramic blocks is complicated by hundreds of cylindrical coolant channels and heater rods embedded in them. As such, it is necessary to deduce an effective thermal conductivity for the ceramic to analyze the core thermal behavior. SAM predictions of the helium coolant and ceramic temperatures are compared with test data measured in three equivalent sectors. Overall, the SAM results agree reasonably well with test data within the variation of data among the three sectors, which demonstrates SAM's capability in capturing transient effects in HTGR using the simplified 2D ring model.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS2020

Advanced High-Lift Leading-Edge (AHLLE) Technologies for Laminar Wings

This is the Northrop Grumman final report for the Advanced High-Lift Leading-Edge Technologies for Laminar Wings conceptual design and demonstration program performed for the National Aeronautics and Space Administration Langley Research Center. Northrop Grumman conducted conceptual aerodynamic CFD and QFD analysis to downselect two leading-edge concepts for high lift. The effort included aerodynamic analysis to examine instrumentation placement and propulsion analysis to design internal cavity geometry of the high-pressure blown slot configurations. A wind tunnel model design created a 7-foot span model of the complex configuration. A wind tunnel test plan, which included baseline, deflected trailing-edge, blown deflected trailing-edge, blown leading-edge, morphed leading-edge, and blown/morphed leading-edge airfoil configurations was developed. As a result of this program approach, a successful demonstration of high-lift leading-edge technology was reached, and the target goals for improved lift were exceeded by 30% with a maximum section lift coefficient (Cl) of 5.2. Incremental lift coefficient (DCl) of 3.5 and 3.1 was achieved on the two laminar compatible airfoil configurations (at cruise). Our high-lift design yielded a massive 94% lift improvement over conventional high-lift Krueger flap configurations while providing laminar flow capability on the cruise configuration.

High-lift↗
DOE OSTI2023

Brewster angle-cavity ringdown spectroscopy for low temperature plasma measurements in multiphases

Here we report on the development of a Brewster angle-cavity ringdown spectroscopy (BA-CRDS) system for low temperature plasma diagnostics. The system can measure gas species in solutions, with a detection limit (minimum detectable absorbance) of 9.1 × 10 -5 , which is equivalent to a detection limit of 0.04 parts per billion for measuring OH radicals in water at 308 nm. With higher reflectivity ringdown mirrors and improved design of a Brewster angle cell, the detection limit can potentially be up to 10 -6 or lower. In this exploratory study, the absorption cross sections of HgBr 2 and H 2 O 2 in the aqueous phase at 256 nm are measured to be (1.8 ± 0.1) × 10 -18 cm 2 and (5.2 ± 0.5) × 10 -20 cm 2 , respectively. Furthermore, temporal profiles of absorbance from distilled water, HgBr 2 , and H 2 O 2 solutions when interacting with a helium atmospheric plasma jet are individually characterized at different plasma powers, gas flow rates, and/or solute concentrations. The observed linear temporal profiles of absorbance from the plasma-interacted water suggest formation of H 2 O 2 from plasma-generated OH radicals, while the nonlinear temporal profiles from the plasma-treated HgBr 2 solutions reveal possible removal of HgBr 2 by OH radicals. Our results demonstrate that the new BA-CRDS system is a powerful tool for quantification of reactive plasma species in multiphases or other complex settings.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

An Overview of Non-LWR Vessel Cooling Systems for Passive Decay Heat Removal (Technical Letter Final Report)

The following report provides a technical review of various reactor vessel cooling system (VCS) concepts under consideration for decay heat removal (DHR) in non-light water advanced reactor designs. This review focuses on ex-vessel designs, including the Reactor Cavity Cooling System (RCCS), the Reactor Vessel Auxiliary Cooling System (RVACS), and hybrid iterations, using both air and water cooling to achieve their heat removal function. Based on a literature review of publicly available sources published between 1979 and 2021, a technical summary is presented detailing existing and planned VCS design options, their applicability to specific reactor type, and review of authored research and development studies. Following an assessment of the availability of data and modeling tools, an evaluation was performed analyzing their likely performance during normal, degraded, and accident conditions, including reliability, stability, and longevity. With renewed consideration of air-based DHR systems by some US vendors, there is a need to fully understand the complexities inherent to natural circulation systems, and more importantly, how they may influence safety-related heat removal functions and system performance. For DHR systems that rely on air-based mode of cooling, this entails quantification of the impacts of low flow conditions during start-up, relative elevations of inlet and outlet ducts for below-grade installations, effects of the use of multiple parallel chimneys for redundancy, and impact of weather conditions at the plant site. For systems relying on water-based mode cooling, reliance on limited coolant inventory in an event of pipe break or extended duration accident scenarios, sensitivities related to stability of boiling flow and heat transfer conditions, and quantification of full-scale structural vibrations on balance of plant structures are among the complicating factors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1995

Advanced Multi-phase Flow CFD Model Development for Solid Rocket Motor Flowfield Analysis

A Navier-Stokes code, finite difference Navier-Stokes (FDNS), is used to analyze the complicated internal flowfield of the SRM (solid rocket motor) to explore the impacts due to the effects of chemical reaction, particle dynamics, and slag accumulation on the solid rocket motor (SRM). The particulate multi-phase flowfield with chemical reaction, particle evaporation, combustion, breakup, and agglomeration models are included in present study to obtain a better understanding of the SRM design. Finite rate chemistry model is applied to simulate the chemical reaction effects. Hermsen correlation model is used for the combustion simulation. The evaporation model introduced by Spalding is utilized to include the heat transfer from the particulate phase to the gase phase due to the evaporation of the particles. A correlation of the minimum particle size for breakup expressed in terms of the Al/Al2O3 surface tension and shear force was employed to simulate the breakup of particles. It is assumed that the breakup occurs when the Weber number exceeds 6. A simple L agglomeration model is used to investigate the particle agglomeration. However, due to the large computer memory requirements for the agglomeration model, only 2D cases are tested with the agglomeration model. The VOF (Volume of Fluid) method is employed to simulate the slag buildup in the aft-end cavity of the redesigned solid rocket motor (RSRM). Monte Carlo method is employed to calculate the turbulent dispersion effect of the particles. The flowfield analysis obtained using the FDNS code in the present research with finite rate chemical reaction, particle evaporation, combustion, breakup, agglomeration, and VOG models will provide a design guide for the potential improvement of the SRM including the use of materials and the shape of nozzle geometry such that a better performance of the SRM can be achieved. The simulation of the slag buildup in the aft-end cavity can assist the designer to improve the design of the RSRM geometry.

Liaw, Paul↗
NASA NTRS1987

A review of research and development on the microwave-plasma electrothermal rocket

The microwave-plasma electrothermal rocket (MWPETR) shows promise for spacecraft propulsion and maneuvering, without some of the drawbacks of competitive electric propulsion systems. In the MWPETR, the electric power is first converted to microwave-frequency radiation. In a specially-designed microwave cavity system, the electromagnetic energy of the radiation is transferred to the electrons in a plasma sustained in the working fluid. The resulting high-energy electrons transfer their energy to the atoms and molecules of the working fluid by collisions. The working fluid, thus heated, expands through a nozzle to generate thrust. In the MWPETR, no electrodes are in contact with the working fluid, the energy is transferred into the working fluid by nonthermal mechanisms, and the main requirement for the materials of construction is that the walls of the plasma chamber be insulating and transparent to microwave radiation at operating conditions. In this survey of work on the MWPETR, several experimental configurations are described and compared. Diagnostic methods used in the study are described and compared, including titration, spectroscopy, calorimetry, electric field measurements, gas-dynamic methods, and thrust measurements. Measured and estimated performance efficiencies are reported. Results of computer modeling of the plasma and of the gas flowing from the plasma are summarized.

Hawley, Martin C.↗
NASA NTRS2002

Damping of Quasi-stationary Waves Between Two Miscible Liquids

Two viscous miscible liquids with an initially sharp interface oriented vertically inside a cavity become unstable against oscillatory external forcing due to Kelvin-Helmholtz instability. The instability causes growth of quasi-stationary (q-s) waves at the interface between the two liquids. We examine computationally the dynamics of a four-mode q-s wave, for a fixed energy input, when one of the components of the external forcing is suddenly ceased. The external forcing consists of a steady and oscillatory component as realizable in a microgravity environment. Results show that when there is a jump discontinuity in the oscillatory excitation that produced the four-mode q-s wave, the interface does not return to its equilibrium position, the structure of the q-s wave remains imbedded between the two fluids over a long time scale. The damping characteristics of the q-s wave from the time history of the velocity field show overdamped and critically damped response; there is no underdamped oscillation as the flow field approaches steady state. Viscous effects serve as a dissipative mechanism to effectively damp the system. The stability of the four-mode q-s wave is dependent on both a geometric length scale as well as the level of background steady acceleration.

Duval, Walter M. B.↗
DOE OSTI2022

Challenges to Reliable Production Nitrogen Doping of Nb for SRF Accelerating Cavities

Over the last several years, alloying of the surface layer of niobium SRF cavities has been demonstrated to beneficially lower the superconducting RF surface resistance. Nitrogen, titanium, and oxygen have all been demonstrated as effective alloying agents, occupying interstitial sites in the niobium lattice within the RF penetration depth and even deeper, when allowed to thermally diffuse into the surface at appropriate temperatures. The use of nitrogen for this function has been often termed ’nitrogen doping’ and is being applied in the LCLS-II and LCLS-II HE projects. We report characterization studies of the distribution of nitrogen into the exposed niobium surface and how such distribution is affected by the quality of the vacuum furnace environment in which the doping takes place, and the complexity of nitride crystal growth on different grain orientations of surface niobium. Using state-of-the-art quantification methods by dynamic secondary ion mass spectrometry (SIMS) depth profiling in niobium, we identify several phenomena involving furnace-sourced contamination. We also highlight a potential issue with N₂ flow constraints from the flange ’caps’ used during heat treatments.

Reece, Charles↗
NASA NTRS1993

FAS multigrid calculations of three dimensional flow using non-staggered grids

Grid staggering is a well known remedy for the problem of velocity/pressure coupling in incompressible flow calculations. Numerous inconveniences occur, however, when staggered grids are implemented, particularly when a general-purpose code, capable of handling irregular three-dimensional domains, is sought. In several non-staggered grid numerical procedures proposed in the literature, the velocity/pressure coupling is achieved by either pressure or velocity (momentum) averaging. This approach is not convenient for simultaneous (block) solvers that are preferred when using multigrid methods. A new method is introduced in this paper that is based upon non-staggered grid formulation with a set of virtual cell face velocities used for pressure/velocity coupling. Instead of pressure or velocity averaging, a momentum balance at the cell face is used as a link between the momentum and mass balance constraints. The numerical stencil is limited to 9 nodes (in 2D) or 27 nodes (in 3D), both during the smoothing and inter-grid transfer, which is a convenient feature when a block point solver is applied. The results for a lid-driven cavity and a cube in a lid-driven cavity are presented and compared to staggered grid calculations using the same multigrid algorithm. The method is shown to be stable and produce a smooth (wiggle-free) pressure field.

Matovic, D.↗