Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “STAGNATION PRESSURE”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
NASA NTRS1974

Rearward-facing steps in laminar supersonic flows with and without suction

An experimental investigation of heat-transfer and pressure distributions within regions of laminar separated flows produced by two-dimensional rearward-facing steps has been carried out at freestream Mach numbers of around 4 in the range of step height-to-boundary layer thickness varying from 0.1 to 2.4. With no suction from the separated area, the ratio of the maximum post-step heat transfer to the attached-flow values was less than unity. The maximum heating-rate region was located far downstream of the reattachment plate stagnation point. Mass suction from the separated area increased the local heating rates, this effect was however relatively weak for purely laminar flow conditions and the competing effect of the step height clearly predominated. At step heights comparable with boundary-layer thickness, even removing the entire approaching boundary layer was not sufficient to raise the post-step heating rates above the flat-plate values.

Jakubowski, A. K.↗
NASA NTRS1999

Numerical Simulations of Acoustically Driven, Burning Droplets

This computational study focuses on understanding and quantifying the effects of external acoustical perturbations on droplet combustion. A one-dimensional, axisymmetric representation of the essential diffusion and reaction processes occurring in the vicinity of the droplet stagnation point is used here in order to isolate the effects of the imposed acoustic disturbance. The simulation is performed using a third order accurate, essentially non-oscillatory (ENO) numerical scheme with a full methanol-air reaction mechanism. Consistent with recent microgravity and normal gravity combustion experiments, focus is placed on conditions where the droplet is situated at a velocity antinode in order for the droplet to experience the greatest effects of fluid mechanical straining of flame structures. The effects of imposed sound pressure level and frequency are explored here, and conditions leading to maximum burning rates are identified.

Kim, H.-C.↗
NASA NTRS1956

The Compressible Laminar Boundary Layer with Heat Transfer and Arbitrary Pressure Gradient

An approximate method for the calculation of the compressible laminar boundary layer with heat transfer and arbitrary pressure gradient, based on Thwaites' correlation concept, is presented. With the definition of dimensionless shear and heat-transfer parameters and an assumed correlation of these parameters in terms of a momentum parameter, a complete system of relations for calculating skin friction and heat transfer results. Knowledge of velocity or temperature profiles is not necessary in using this calculation method. When the method is applied to a convergent-divergent, axially symmetric rocket nozzle, it shows that high rates of heat transfer are obtained at the initial stagnation point and at the throat of the nozzle. Also indicated are negative displacement thicknesses in the convergent portion of the nozzle; these occur because of the high density within the lower portions of the cooled boundary layer. (author)

Cohen, Clarence B↗
NASA NTRS1989

Heat transfer to throat tubes in a square-chambered rocket engine at the NASA Lewis Research Center

A gaseous H2/O2 rocket engine was constructed at the NASA-Lewis to provide a high heat flux source representative of the heat flux to the blades in the high pressure fuel turbopump (HPFTP) during startup of the space shuttle main engines. The high heat flux source was required to evaluate the durability of thermal barrier coatings being investigated for use on these blades. The heat transfer, and specifically, the heat flux to tubes located at the throat of the test rocket engine was evaluated and compared to the heat flux to the blades in the HPFTP during engine startup. Gas temperatures, pressures and heat transfer coefficients in the test rocket engine were measured. Near surface metal temperatures below thin thermal barrier coatings were also measured at various angular orientations around the throat tube to indicate the angular dependence of the heat transfer coefficients. A finite difference model for a throat tube was developed and a thermal analysis was performed using the measured gas temperatures and the derived heat transfer coefficients to predict metal temperatures in the tube. Near surface metal temperatures of an uncoated throat tube were measured at the stagnation point and showed good agreement with temperatures predicted by the thermal model. The maximum heat flux to the throat tube was calculated and compared to that predicted for the leading edge of an HPFTP blade. It is shown that the heat flux to an uncooled throat tube is slightly greater than the heat flux to an HPFTP blade during engine startup.

Nesbitt, James A.↗
NASA NTRS2014

Experimental Study of Slat Noise from 30P30N Three-Element High-Lift Airfoil in JAXA Hard-Wall Low-Speed Wind Tunnel

Aeroacoustic measurements associated with noise radiation from the leading edge slat of the canonical, unswept 30P30N three-element high-lift airfoil configuration have been obtained in a 2 m x 2 m hard-wall wind tunnel at the Japan Aerospace Exploration Agency (JAXA). Performed as part of a collaborative effort on airframe noise between JAXA and the National Aeronautics and Space Administration (NASA), the model geometry and majority of instrumentation details are identical to a NASA model with the exception of a larger span. For an angle of attack up to 10 degrees, the mean surface Cp distributions agree well with free-air computational fluid dynamics predictions corresponding to a corrected angle of attack. After employing suitable acoustic treatment for the brackets and end-wall effects, an approximately 2D noise source map is obtained from microphone array measurements, thus supporting the feasibility of generating a measurement database that can be used for comparison with free-air numerical simulations. Both surface pressure spectra obtained via KuliteTM transducers and the acoustic spectra derived from microphone array measurements display a mixture of a broad band component and narrow-band peaks (NBPs), both of which are most intense at the lower angles of attack and become progressively weaker as the angle of attack is increased. The NBPs exhibit a substantially higher spanwise coherence in comparison to the broadband portion of the spectrum and, hence, confirm the trends observed in previous numerical simulations. Somewhat surprisingly, measurements show that the presence of trip dots between the stagnation point and slat cusp enhances the NBP levels rather than mitigating them as found in a previous experiment.

Murayama, Mitsuhiro↗
NASA NTRSDate not supplied

A Model for Ice Accretion Roughness Evolution and Spatial Variations

Over the past decade, multiple investigations of ice accretion roughness and spatial variations have been performed in the Icing Research Tunnel (IRT) at the NASA Glenn Research Center. The early investigations used models of NACA 0012 airfoils with different chord sizes and focused on temporal scaling and primary cloud scaling parameters such as stagnation point collection efficiency. Subsequent investigations included the effects of model sweep, airfoil shape, airfoil lifting condition, and freestream static temperature. To develop a predictive model for roughness evolution in generalized icing situations, the maximum roughness values for the airfoils and conditions used in the angle of attack and freestream temperature investigations were scaled to eliminate the temporal variations. LEWICE simulations were employed to identify the local collection efficiency at the locations of maximum roughness, and a two-dimensional panel-method code was used to identify the local static pressure at the location of the maximum roughness. Because of the stochastic nature of ice accretion roughness, a physics-directed approach was employed to develop a multi-dimensional correlation based on the local pressure coefficient, the local total temperature, the cloud properties, and the cloud exposure time. The resulting correlation predictions are compared to ice shapes measured in the IRT for both a 21-in. NACA 0012 airfoil model and a 60-in. HAARP-II model. The resulting predictions indicate that the local freezing fraction must be included in the roughness predictive model. Implications regarding icing heat transfer predictions using the resulting roughness model are also discussed.

Icing↗
DOE OSTI2022

Probing atomic physics at ultrahigh pressure using laser-driven implosions

Abstract Spectroscopic measurements of dense plasmas at billions of atmospheres provide tests to our fundamental understanding of how matter behaves at extreme conditions. Developing reliable atomic physics models at these conditions, benchmarked by experimental data, is crucial to an improved understanding of radiation transport in both stars and inertial fusion targets. However, detailed spectroscopic measurements at these conditions are rare, and traditional collisional-radiative equilibrium models, based on isolated-atom calculations and ad hoc continuum lowering models, have proved questionable at and beyond solid density. Here we report time-integrated and time-resolved x-ray spectroscopy measurements at several billion atmospheres using laser-driven implosions of Cu-doped targets. We use the imploding shell and its hot core at stagnation to probe the spectral changes of Cu-doped witness layer. These measurements indicate the necessity and viability of modeling dense plasmas with self-consistent methods like density-functional theory, which impact the accuracy of radiation transport simulations used to describe stellar evolution and the design of inertial fusion targets.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS2022

Design of Hypervelocity Flow Generator (HFG) and Its Diagnostics

Ground facilities for hypersonic research are a key element for successful development of aerodynamically proven hypervelocity vehicles. Design concepts and diagnostics of a hypervelocity flow generator (HFG) were made as a test platform for hypersonic/hypervelocity spacecraft models at the NASA Langley Research Center. The HFG is a hypersonic flow field generator using optically heated gas which is blown into an 80 cu m vacuum chamber. The vacuum chamber is kept at a stable vacuum pressure with a combination of three large vacuum pumps, while the HFG is in the test mode. The HFG provides a relatively small test section with approximately a 20 cm window. This facility was designed to generate 2.45 km/sec of flow speed, and potentially generate a continuous flow with the nozzle and vacuum system. The window in test section provides a direct view of the shock wave around a model in order to measure temperature, pressure, and density profile within the shock layer. One of the key test goals under this project is to understand why the emission spectra from the standing shockwave plasma predicted by Lora-Loran codes are significantly different from the measured emission spectra from the Flight Investigation Reentry Environment (FIRE II) Flight. The correct estimate of the thermal loading on the leading edge of hypersonic vehicles greatly affects the aerodynamic design, the material selection for the vehicle, and the cooling requirement and can be obtained by the precise modeling of emission spectra from shockwave plasma. However, the estimation of thermal loading is not an easy task due to complex non-equilibrium radiative process within high temperature shock layers that still falls into a category of cold plasma. Direct flight experiments are the most desirable, but not a cost-effective approach. Analysis by computational fluid dynamics (CFD) offers many test flexibilities. However, the CFD codes must be fully tested and validated with experimental data before the codes are effectively used for practical design. Large discrepancies between experiments and codes appear in hypersonic/hypervelocity flow regimes at high altitudes of 60 km ~ 90 km. This HFG facility offers some important parameters for CFD code validation, such as collision cross-sections, relaxation times, reaction rate coefficients and transportation coefficients. The HFG test facility is based on the ejection flow of high temperature gas heated over to 3500 K through a nozzle. The tungsten gas chamber of the HFG is heated up to a desired temperature by a 60 kW optical power beam source. This system consists of an optical power source, a thermal chamber, an expansion nozzle, a test section, and an 80 cu m vacuum tank. 60 kW optical input power is obtained from the 150 kW Vortek arc lamp system (by Vortek Industries, Vancouver, Canada). This optical beam is focused to heat the gas chamber within which a flow media is heated. The maximum achievable temperature of the flow medium reached approximately 3500 K or even higher but is limited by the melting point of the chamber material used. The exhaust velocity through the nozzle was determined by the stagnation temperature and the molecular weight of the working medium at the test section. To provide design parameters, a NASA Chemical Equilibrium with Application (CEA) computer program is used for the simulation of aerothermal data. This CEA program can calculate chemical equilibrium and properties of complex mixtures using shock tube parameters. For nitrogen gas at 2666 K stagnation temperature, the maximum achievable velocity at the test section is approximately 2.45 km/sec which is within the range of the thermal velocity of 8000 K shock layer. Based on the calculation through the CEA program, the design parameters of a HFG were determined and implemented for the test section that includes an expansion nozzle. The installed test section of the experimental facility can sustain a condition of a re-entry vehicle from the Space at an altitude of 60-90 km. After installation of the HFG, the system was fully tested and its operational parameters were measured. An 80 cu m vacuum chamber of HFG was set at 1 torr level to keep a stable downstream condition. The pumping time to reach the minimum vacuum pressure (~ 1 torr) at the test chamber from atmospheric pressure was approximately an hour using Kinney (MBV-14000/MB -1600/KT-300) and Stokes (Model 1772 and 412) pumps. Such a setup condition allowed a continued stable operation of the HFG experiment with flowrates through 1-, 3-, or 5-mm diameters of nozzle throat. The flow characteristics of the HFG for various operating conditions were performed using a focused Schlieren method. At 0.8 torr chamber pressure, a barrel shock was observed at the test section. With a sphere obstacle of 12 mm diameter in the flow, a bow shock (~ 2-mm thickness) was observed by a focused Schlieren visualization method.

Hypervelocity flow↗
NASA NTRS1972

Thermodynamic properties of hydrogen-helium plasmas.

Calculation of the thermodynamic properties of an atomic hydrogen-helium plasma for postulated conditions present in a stagnation shock layer of a spacecraft entering the atmosphere of Jupiter. These properties can be used to evaluate transport properties, to calculate convective heating, and to investigate nonequilibrium behavior. The calculations have been made for temperatures from 10,000 to 100,000 K, densities of 10 to the minus 7th and .00001 g cu cm, and three plasma compositions: pure hydrogen, 50% hydrogen/50% helium, and pure helium. The shock layer plasma consists of electrons, protons, atomic hydrogen, atomic helium, singly ionized helium, and doubly atomized helium. The thermodynamic properties which have been investigated are: pressure, average molecular weight, internal energy, enthalpy, entropy, specific heat, and isentropic speed of sound. A consistent model was used for the reduction of the ionization potential in the calculation of the partition functions.

Nelson, H. F.↗
NASA NTRSDate not supplied

MARS 2020 Backshell Radiative Heating Measurement and Shock Tube Verification

In February 2021, the Perseverance rover was brought to the surface of Mars by the Mars 2020 mission. A feature of the Mars 2020 capsule was instrumentation to measure its entry, descent and landing (EDL) with the so-called Mars EDL Instrumentation 2 (MEDLI2) [1]. The MEDLI2 introduced, among other things, backshell instrumentation, including a broadband radiometer. The radiometer was mounted on the leeside of the vehicle next to thermocouple plugs and a heat flux gauge. The data returned by the leeside MEDLI2 heat flux gauge is largely analogous to that measured by the COMARS gauge flown on the ExoMars Schiaparelli [2] entry in that it is measuring heat flux in an area that is entirely dominated by radiative heating. The COMARS measurement provided excellent validation of backshell radiative heating models for Mars entry, albeit at a limited number of points [3]. The MEDLI2 heat flux gauge measurement effectively confirmed the quality of the prediction, extended over the full trajectory, as will be presented in this paper. The MEDLI2 radiometer,however, was blocked by ablation products and suffered a loss of half of its signal. The second backshell heat flux gauge installed on the windside of the vehicle was also well predicted, although the heating had both radiative and convective contributions. It was desired to reproduce the conditions of the Mars 2020 entry via ground testing in the Electric Arc Shock Tube (EAST) at NASA Ames. Tests to verify stagnation line heating were previously reported in EAST, confirming the presence of shock layer radiation as the major discrepancy in heatshield temperature modelling [4]. Therefore, tests for stagnation line heating were not repeated. Instead, the shock tube informed bias method [5, 6] was used to identify test conditions that may produce similarity to streamlines that pass around the backside of the vehicle and are responsible for the radiation observed at the two heat flux gauge locations. This method was used to identify a range of velocities and densities in the shock tube that are relevant for confirming the radiative environment encountered. This paper reports the results obtained in the 10 cm diameter EAST shock tube, corresponding to later trajectory points at ambient pressures of 1.1-2.0 Torr and velocities from 1.2-3.5 km/s. The test series employed two primary diagnostics: emission spectroscopy and tunable diode laser absorption spectroscopy (TDLAS). The emission spectroscopy performed broadband measurements of the radiative emission of the 4.3 and 2.7 m bands of CO2 at flight similar conditions, obtaining both spectral and spatial data corresponding to the relaxation behind the shock front. The TDLAS measured the absorption of several lines of CO and CO2 and obtains species number densities and temperatures as a function of time behind the shock front. This paper will review highlights of this test series and analyses of the emission and absorption data. While the datasets generally show good agreement with predictions, a few discrepancies and items for additional investigation are identified and will be discussed

Brett A Cruden↗
NASA NTRS1995

Measurements of scramjet thrust in shock tunnels

By using results obtained in tests on supersonic combustion of hydrogen in air, the conditions governing model size and operating pressure levels for shock tunnel experiments on models of flight vehicles with scramjet propulsion are established. It is seen that large models are required. The development of the stress wave force balance is then described, and its use as a method of measuring thrust/drag on such models is discussed. Test results on a simple, fully integrated scramjet model, with intakes, combustion chambers, thrust surfaces and exterior surfaces, using a 13 percent silane 87 percent hydrogen fuel mixture, showed that a steady state with thrust generation could be achieved within the shock tunnel test time, and the thrust could be measured. Results are presented for a range of stagnation enthalpies, and show that the scramjet model produces net positive thrust at velocities up to 2.4 km/sec.

Stalker, R. J.↗
NASA NTRS1996

Adaptive grid methods for RLV environment assessment and nozzle analysis

Rapid access to highly accurate data about complex configurations is needed for multi-disciplinary optimization and design. In order to efficiently meet these requirements a closer coupling between the analysis algorithms and the discretization process is needed. In some cases, such as free surface, temporally varying geometries, and fluid structure interaction, the need is unavoidable. In other cases the need is to rapidly generate and modify high quality grids. Techniques such as unstructured and/or solution-adaptive methods can be used to speed the grid generation process and to automatically cluster mesh points in regions of interest. Global features of the flow can be significantly affected by isolated regions of inadequately resolved flow. These regions may not exhibit high gradients and can be difficult to detect. Thus excessive resolution in certain regions does not necessarily increase the accuracy of the overall solution. Several approaches have been employed for both structured and unstructured grid adaption. The most widely used involve grid point redistribution, local grid point enrichment/derefinement or local modification of the actual flow solver. However, the success of any one of these methods ultimately depends on the feature detection algorithm used to determine solution domain regions which require a fine mesh for their accurate representation. Typically, weight functions are constructed to mimic the local truncation error and may require substantial user input. Most problems of engineering interest involve multi-block grids and widely disparate length scales. Hence, it is desirable that the adaptive grid feature detection algorithm be developed to recognize flow structures of different type as well as differing intensity, and adequately address scaling and normalization across blocks. These weight functions can then be used to construct blending functions for algebraic redistribution, interpolation functions for unstructured grid generation, forcing functions to attract/repel points in an elliptic system, or to trigger local refinement, based upon application of an equidistribution principle. The popularity of solution-adaptive techniques is growing in tandem with unstructured methods. The difficultly of precisely controlling mesh densities and orientations with current unstructured grid generation systems has driven the use of solution-adaptive meshing. Use of derivatives of density or pressure are widely used for construction of such weight functions, and have been proven very successful for inviscid flows with shocks. However, less success has been realized for flowfields with viscous layers, vortices or shocks of disparate strength. It is difficult to maintain the appropriate mesh point spacing in the various regions which require a fine spacing for adequate resolution. Mesh points often migrate from important regions due to refinement of dominant features. An example of this is the well know tendency of adaptive methods to increase the resolution of shocks in the flowfield around airfoils, but in the incorrect location due to inadequate resolution of the stagnation region. This problem has been the motivation for this research.

Thornburg, Hugh J.↗
NASA NTRS1995

Investigation Flow Uniformity in a Supersonic Duct with High Enthalpy Flows

Flow uniformity in a high enthalpy facility is investigated. The scramjet research facility is composed of a rectangular combustor duct connected to a 100 MW electric arc air heater. The Mach 3.3 flow is accelerated through a two-dimensional contoured nozzle. Instream measurements were made with water-cooled Pitot probes and stagnation point heat flux gages at stream enthalpy levels ranging from 4 to 7 Mj/kg. Flow surveys were made on the flow centerline and off centerline in order to measure the three dimensional uniformity of the flow in the rectangular duct. Measurements indicated that although the flow in the aspect ratio 6:1 duct was relatively uniform on the centerline, three dimensional viscous effects were apparent near the corners. Flow through the nozzle and constant area duct was modeled computationally using a two dimensional, Navier-Stokes, reacting gas code. The computations predict that the flow in the test section is in vibrational equilibrium. The computed and measured Pitot pressure and heat flux profiles are in reasonable agreement with the experimental data.

Balboni, John↗
NASA NTRS1994

Turbulence modeling for separated flow

Two projects are described in this report. The first involves assessing turbulence models in separated flow. The second addresses the anomalous behavior of certain turbulence models in stagnation point flow. The primary motivation for developing turbulent transport models is to provide tools for computing non-equilibrium, or complex, turbulent flows. Simple flows can be analyzed using data correlations or algebraic eddy viscosities, but in more complicated flows such as a massively separated boundary layer, a more elaborate level of modeling is required. It is widely believed that at least a two-equation transport model is required in such cases. The transport equations determine the evolution of suitable velocity and time-scales of the turbulence. The present study included assessment of second-moment closures in several separated flows, including sharp edge separation; smooth wall, pressure driven separation; and unsteady vortex shedding. Flows with mean swirl are of interest for their role in enhancing mixing both by turbulent and mean motion. The swirl can have a stabilizing effect on the turbulence. An axi-symmetric extension to the INS-2D computer program was written adding the capability of computing swirling flow. High swirl can produce vortex breakdown on the centerline of the jet and it occurs in various combustors.

Durbin, Paul A.↗
NASA NTRS1991

Velocity and pressure characteristics of a model SSME high pressure fuel turbopump

Under the present effort an experiment rig has been constructed, an instrumentation package developed and a series of mean and rms velocity and pressure measurements made in a turbopump which modelled the first stage of the Space Shuttle Main Engine (SSME) High Pressure Fuel Turbopump. The rig was designed so as to allow initial experiments with a single configuration consisting of a bell-mouth inlet, a flight impeller, a vaneless diffuser and a volute. Allowance was made for components such as inlet guide vanes, exit guide vanes, downstream pumps, etc. to be added in future experiments. This flexibility will provide a clear baseline set of experiments and allow evaluation in later experiments of the effect of adding specific components upon the pump performance properties. The rotational speed of the impeller was varied between 4260 and 7680 rpm which covered the range of scaled SSME rotation speeds when due allowance is made for the differing stagnation temperature, model to full scale. The results at the inlet obtained with rotational speeds of 4260, 6084 and 7680 rpm showed that the axial velocity at the bell-mouth inlet remained roughly constant at 2.2 of the bulk velocity at the exit of the turbopump near the center of the inlet, but it decreased rapidly with increasing radius at all three speeds. Reverse flow occurred at a radius greater than 0.9 R for all three speeds and the maximum negative velocity reduced from 1.3 of the bulk velocity at the exit of the turbopump at 4260 rpm to 0.35 at 7680 rpm, suggesting that operating at a speed closer to the design condition of 8700 rpm improved the inlet characteristics. The reverse flow caused positive prerotation at the impeller inlet which was negligibly small near the center but reached 0.7 of the impeller speed at the outer annulus. The results in the diffuser and the volute obtained at 7680 rpm show that the hub and shroud walls of the diffuser were characterized by regions of transient reverse flow with negative revolution-averaged velocity of 8 percent of the maximum forward revolution-averaged velocity at the center of the diffuser passage near the shroud wall.

Tse, D. G-N.↗
NASA NTRS1974

Solutions of the Navier-Stokes equations for vortex breakdown

Steady solutions of the Navier-Stokes equations, in terms of velocity and pressure, for breakdown in an unconfined viscous vortex are obtained numerically using the artificial compressibility technique of Chorin combined with an ADI finite-difference scheme. Axisymmetry is assumed and boundary conditions are carefully applied at the boundaries of a large finite region in an axial plane while resolution near the axis is maintained by a coordinate transformation. The solutions, which are obtained for Reynolds numbers up to 200 based on the free-stream axial velocity and a characteristic core radius, show that breakdown results from the diffusion and convection of vorticity away from the vortex core which, because of the strong coupling between the circumferential and axial velocity fields in strongly swirling flows, can lead to stagnation and reversal of the axial flow near the axis.

Grabowski, W. J.↗
NASA NTRS1988

Numerical investigations in three-dimensional internal flows

Renewed interest in simulation of high enthalpy flows was brought about by interest in the design and testing of hypersonic propulsion components. This renewed interest has led to refurbishment of high speed ground based wind tunnel facilities. Simulation of the high enthalpy flows occurring in the NASP propulsion system is of particular interest. Studying the combustor presentes some of the more challenging problems in the ground based simulation of the fluid flow. A difficulty in simulating the required conditions at the entrance to the combustor is related to achieveing a sufficiently high static pressure such that combustion would be possible at a high Mach number. An analysis of the flow within the original nozzle was not performed previously using numerical simulations, although it was operated experimentally. A Navier-Stokes code was chosen in the present study to examine analytically the flow within representative selected nozzles for use in combustion experiments. The numerical simulation of these fluid flows is popular and may be able to shed light on the direction that redesign of the shock tube nozzles should take. The numerical codes were modified to account for the axially symmetric flow field, for the inclusion of equilibrium real gas effects and for ideal gases whose thermodynamic properties are frozen at the stagnation conditions. Numerical simulation of flows within the original and two proposed nozzles, including the energy lost to the relatively cold walls, are investigated.

Rose, William C.↗
DOE OSTI2022

Design and execution of a verification, validation, and uncertainty quantification plan for a numerical model of left ventricular flow after LVAD implantation

Left ventricular assist devices (LVADs) are implantable pumps that act as a life support therapy for patients with severe heart failure. Despite improving the survival rate, LVAD therapy can carry major complications. Particularly, the flow distortion introduced by the LVAD in the left ventricle (LV) may induce thrombus formation. While previous works have used numerical models to study the impact of multiple variables in the intra-LV stagnation regions, a comprehensive validation analysis has never been executed. The main goal of this work is to present a model of the LV-LVAD system and to design and follow a verification, validation and uncertainty quantification (VVUQ) plan based on the ASME V&V40 and V&V20 standards to ensure credible predictions. The experiment used to validate the simulation is the SDSU cardiac simulator, a bench mock-up of the cardiovascular system that allows mimicking multiple operation conditions for the heart-LVAD system. The numerical model is based on Alya, the BSC’s in-house platform for numerical modelling. Alya solves the Navier-Stokes equation with an Arbitrary Lagrangian-Eulerian (ALE) formulation in a deformable ventricle and includes pressure-driven valves, a 0D Windkessel model for the arterial output and a LVAD boundary condition modeled through a dynamic pressure-flow performance curve. The designed VVUQ plan involves: (a) a risk analysis and the associated credibility goals; (b) a verification stage to ensure correctness in the numerical solution procedure; (c) a sensitivity analysis to quantify the impact of the inputs on the four quantities of interest (QoIs) (average aortic root flow $Q^{avg}_{Ao}$, maximum aortic root flow $Q^{max}_{Ao}$, average LVAD flow $Q^{avg}_{VAD}$, and maximum LVAD flow $Q^{max}_{VAD}$; (d) an uncertainty quantification using six validation experiments that include extreme operating conditions. Numerical code verification tests ensured correctness of the solution procedure and numerical calculation verification showed a grid convergence index (GCI)95% <3.3%. The total Sobol indices obtained during the sensitivity analysis demonstrated that the ejection fraction, the heart rate, and the pump performance curve coefficients are the most impactful inputs for the analysed QoIs. The Minkowski norm is used as validation metric for the uncertainty quantification. It shows that the midpoint cases have more accurate results when compared to the extreme cases. The total computational cost of the simulations was above 100 [core-years] executed in around three weeks time span in Marenostrum IV supercomputer. This work details a novel numerical model for the LV-LVAD system, that is supported by the design and execution of a VVUQ plan created following recognised international standards. We present a methodology demonstrating that stringent VVUQ according to ASME standards is feasible but computationally expensive.

59 BASIC BIOLOGICAL SCIENCES↗