Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow analysis”

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 595 records · Page 33

Tank 49H solids disturbance analysis

Tank 49H serves as the feed tank for the Salt Waste Processing Facility (SWPF). Transfers into Tank 49H may disturb any solids that have settled to the bottom of the tank, resulting in feed that may exceed the insoluble solids content limit of 1,200 mg/L of the SWPF Waste Acceptance Criteria (WAC). During a transfer into Tank 49H, material that free falls into Tank 49H through a downcomer could potentially disturb any solids on the bottom of the tank and scour or suspend solids from a settled solids layer or turbid region. A previous analysis and report evaluated the potential to disturb solids when transferring into Tank 49H through the B4 riser and recommended a minimum tank level of 120 inches to prevent disturbing any solids in the bottom of the tank. The scope of this task is to perform additional fluid flow analysis to determine whether accounting for disturbed particle settling and particle mixing and dispersion during transfer could allow the minimum liquid level to be reduced below 120 inches. The analysis utilized models from the technical literature to calculate the size and shape of the “plunging jet” as a function of input parameters such as initial velocity, initial jet diameter, elevation of the initial jet, liquid level in the tank, and solids depth. The analysis relied on the M-Star® simulations performed for the previous analysis to provide bounding estimates of the amount of solid particles disturbed and used the MStar® software to calculate the dispersion and mixing of the disturbed solids with other liquid in the tank as the solids are transported to the transfer pump. The analysis showed that with a solid particle size of 5 micron or less, a liquid level of 120 inches should be maintained to prevent significant disturbance of the solid layer at the bottom of Tank 49H.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Toward a CFD nose-to-tail capability - Hypersonic unsteady Navier-Stokes code validation

Computational fluid dynamics (CFD) research for hypersonic flows presents new problems in code validation because of the added complexity of the physical models. This paper surveys code validation procedures applicable to hypersonic flow models that include real gas effects. The current status of hypersonic CFD flow analysis is assessed with the Compressible Navier-Stokes (CNS) code as a case study. The methods of code validation discussed to beyond comparison with experimental data to include comparisons with other codes and formulations, component analyses, and estimation of numerical errors. Current results indicate that predicting hypersonic flows of perfect gases and equilibrium air are well in hand. Pressure, shock location, and integrated quantities are relatively easy to predict accurately, while surface quantities such as heat transfer are more sensitive to the solution procedure. Modeling transition to turbulence needs refinement, though preliminary results are promising.

Edwards, Thomas A.↗

Methods for Quantitative Thermal Analysis of Lithium Solid-State and Beyond Battery Safety

The use of differential scanning calorimetry (DSC) to measure the thermal behavior of individual components and electrolyte/electrode combinations is common. However, here we focus on DSC tests on an anode, cathode, and electrolyte (ACE) component combination over a temperature range that includes many of the phase transitions and key reactions (i.e., to 500 °C) that contribute to thermal runaway. This method can help quantify the complex reaction network in a full cell, thereby informing potential safety issues. Here, we used DSC heat flow data from a solid-state Li 0.43 CoO 2 +C+PVDF | LLZO | Li metal ACE sample and its components to quantify key factors affecting results. We focused on three areas: (1) ACE sample preparation and assembly in DSC pans, (2) DSC measurement parameters, and (3) heat flow analysis. Key points include the choice of component ratios (e.g., commercially relevant N:P capacity ratio), the importance of conductive carbon and binder, type of pan used, DSC ramp rate, and integration method used when dealing with broad and overlapping exothermic peaks. This work deepens the scientific basis and best practices for obtaining heat flow data from ACE samples for early-stage evaluation of solid-state and beyond battery safety.

25 ENERGY STORAGE↗

Analytical study of vortex flaps on highly swept delta wings

This paper highlights some current results from ongoing analytical studies of vortex flaps on highly swept delta wings. A brief discussion of the vortex flow analysis tools is given along with comparisons of the theories to vortex flap force and pressure data. Theoretical trends in surface pressure distribution for both angle-of-attack variation and flap deflection are correctly predicted by Free Vortex Sheet theory. Also shown are some interesting calculations for attached-flow and vortex-flow flap hinge moments that indicate flaps utilizing vortex flow may generate less hinge moment than attached flow flaps. Finally, trailing-edge flap effects on leading-edge flap thrust potential are investigated and theory-experiment comparisons made.

Frink, N. T.↗

Application of a finite element algorithm for high speed viscous flows using structured and unstructured meshes

A higher-order streamline upwinding Petrov-Galerkin finite element method is employed for high speed viscous flow analysis using structured and unstructured meshes. For a Mach 8.03 shock interference problem, successive mesh adaptation was performed using an adaptive remeshing method. Results from the finite element algorithm compare well with both experimental data and results from an upwind cell-centered method. Finite element results for a Mach 14.1 flow over a 24 degree compression corner compare well with experimental data and two other numerical algorithms for both structured and unstructured meshes.

Vemaganti, Gururaja R.↗

A title-gap flow model for use in aerodynamic loads assessment of space shuttle thermal protection system: Parallel gap faces

The problem of predicting aerodynamic loads on the insulating tiles of the space shuttle thermal protection system (TPS) is discussed and seen to require a method for predicting pressure and mass flux in the gaps between tiles. A mathematical model of the tile-gap flow is developed, based upon a slow viscous (Stokes) flow analysis, and is verified against experimental data. The tile-gap pressure field is derived from a solution of the two-dimensional Laplace equation; the mass-flux vector is then calculated from the pressure gradient. The means for incorporating this model into a lumped-parameter network analogy for porous-media flow is given. The means for incorporating this model into a lumped-parameter network analogy for porous-media flow is given. The flow model shows tile-gap mass flux to be very sensitive to the gap width indicating a need for coupling the TPS flow and tile displacement calculation. Analytical and experimental work to improve TPS flow predictions and a possible shuttle TPS hardware modification are recommended.

Dwoyer, D. L.↗

Flow through the tile gaps in the Space Shuttle Thermal Protection System

The problem of predicting aerodynamic loads on the insulating tiles of the Space Shuttle Thermal Protection System (TPS) is discussed and seen to require a method for predicting pressure and mass flux in the gaps between tiles. A mathematical model of the tile-gap flow is developed based upon a slow viscous (Stokes) flow analysis and is verified against available experimental data. This model derives the tile-gap pressure field from a solution of the two-dimensional Laplace equation; the mass flux vector is then calculated from the pressure gradient. The means for incorporating this model into a lumped-parameter network analogy for porous-media flow is also given. The flow model shows tile-gap mass flux to be very sensitive to the gap width indicating a need for coupling the TPS flow and tile displacement calculations. Finally recommendations are made concerning additional analytical and experimental work to improve TPS flow predictions.

Dwoyer, D. L.↗

Development of the METMAN Thermal Desktop Human Model for System Integration

Human thermal analysis and model development inform NASA’s space suit development, vehicle/habitat design, and survivability studies. Historically, human thermal models such as the 41-node metabolic man (METMAN) and the Wissler model have been standalone software tools developed with FORTRAN, a programming language known for its high performance in computationally intensive applications. Though efficient, these standalone programs pose challenges to coupled human-system analysis with detailed life support and thermal control subsystem models in other commercial software. This paper describes the conversion of the METMAN human thermal model from a standalone FORTRAN program to a model in Thermal Desktop (Ansys), a commonly used CAD-based simulation software for engineering that specializes in heat transfer, thermal radiation, and fluid flow analysis. This format was chosen to best facilitate model sharing and compatibility, enabling the direct integration of METMAN human thermal analysis with subsystem models across NASA programs and commercial partners.

Elizabeth Thurston↗

Least Squares Shadowing Sensitivity Analysis of Chaotic Flow Around a Two-Dimensional Airfoil

Gradient-based sensitivity analysis has proven to be an enabling technology for many applications, including design of aerospace vehicles. However, conventional sensitivity analysis methods break down when applied to long-time averages of chaotic systems. This breakdown is a serious limitation because many aerospace applications involve physical phenomena that exhibit chaotic dynamics, most notably high-resolution large-eddy and direct numerical simulations of turbulent aerodynamic flows. A recently proposed methodology, Least Squares Shadowing (LSS), avoids this breakdown and advances the state of the art in sensitivity analysis for chaotic flows. The first application of LSS to a chaotic flow simulated with a large-scale computational fluid dynamics solver is presented. The LSS sensitivity computed for this chaotic flow is verified and shown to be accurate, but the computational cost of the current LSS implementation is high.

Blonigan, Patrick J.↗

Conceptual Design of a Two Spool Compressor for the NASA Large Civil Tilt Rotor Engine

This paper focuses on the conceptual design of a two spool compressor for the NASA Large Civil Tilt Rotor engine, which has a design-point pressure ratio goal of 30:1 and an inlet weight flow of 30.0 lbm/sec. The compressor notional design requirements of pressure ratio and low-pressure compressor (LPC) and high pressure ratio compressor (HPC) work split were based on a previous engine system study to meet the mission requirements of the NASA Subsonic Rotary Wing Projects Large Civil Tilt Rotor vehicle concept. Three mean line compressor design and flow analysis codes were utilized for the conceptual design of a two-spool compressor configuration. This study assesses the technical challenges of design for various compressor configuration options to meet the given engine cycle results. In the process of sizing, the technical challenges of the compressor became apparent as the aerodynamics were taken into consideration. Mechanical constraints were considered in the study such as maximum rotor tip speeds and conceptual sizing of rotor disks and shafts. The rotor clearance-to-span ratio in the last stage of the LPC is 1.5% and in the last stage of the HPC is 2.8%. Four different configurations to meet the HPC requirements were studied, ranging from a single stage centrifugal, two axi-centrifugals, and all axial stages. Challenges of the HPC design include the high temperature (1,560deg R) at the exit which could limit the maximum allowable peripheral tip speed for centrifugals, and is dependent on material selection. The mean line design also resulted in the definition of the flow path geometry of the axial and centrifugal compressor stages, rotor and stator vane angles, velocity components, and flow conditions at the leading and trailing edges of each blade row at the hub, mean and tip. A mean line compressor analysis code was used to estimate the compressor performance maps at off-design speeds and to determine the required variable geometry reset schedules of the inlet guide vane and variable stators that would result in the transonic stages being aerodynamically matched with high efficiency and acceptable stall margins based on user specified maximum levels of rotor diffusion factor and relative velocity ratio.

Veres, Joseph P.↗

Boundary Representation Tolerance Impacts on Mesh Generation and Adaptation

The control of discretization error is critical to obtaining reliable simulation results. Re-cent progress has matured anisotropic mesh adaptation, which automates discretization error control for complex geometries. However, the meshing process can fail when geometric model Boundary REPresentation (BREP) tolerances are larger than boundary layer surface nor-mal and tangential spacing requirements. Geometry sources are often created in Mechanical Computer-Aided Design (MCAD) systems, which are designed to produce models for manufacturing and not the stricter requirements of viscous flow analysis. BREP tolerances are strained by three factors: the inherent complexity of the model (e.g., large range of scales and complex topology), the numerical difficulties arising from surface/surface intersections, and the omission of crucial data during export and translation. Manual preparation of geometry is commonly employed to enable expert-guided mesh generation, which severely inhibits work-flow automation. Accommodation of loose BREP tolerances is required for automated mesh processes, because barrier issues may not be detected until well into the solution process. To raise awareness of this class of geometry issues and their impact on simulation, a survey of these barrier issues is presented with example mitigation techniques. This awareness may also impact geometry creation workflows to prevent the introduction of these artifacts.

CAD↗

Navier-Stokes Analysis of the Flowfield Characteristics of an Ice Contaminated Aircraft Wing

An analytical study was performed as part of the NASA Lewis support of a National Transportation Safety Board (NTSB) aircraft accident investigation. The study was focused on the performance degradation associated with ice contamination on the wing of a commercial turbo-prop-powered aircraft. Based upon the results of an earlier numerical study conducted by the authors, a prominent ridged-ice formation on the subject aircraft wing was selected for detailed flow analysis using 2-dimensional (2-D), as well as, 3-dimensional (3-D) Navier-Stokes computations. This configuration was selected because it caused the largest lift decrease and drag increase among all the ice shapes investigated in the earlier study. A grid sensitivity test was performed to find out the influence of grid spacing on the lift, drag, and associated angle-of-attack for the maximum lift (C(sub lmax)). This study showed that grid resolution is important and a sensitivity analysis is an essential element of the process in order to assure that the final solution is independent of the grid. The 2-D results suggested that a severe stability and control difficulty could have occurred at a slightly higher angle-of-attack (AOA) than the one recorded by the Flight Data Recorder (FDR). This stability and control problem was thought to have resulted from a decreased differential lift on the wings with respect to the normal loading for the configuration. The analysis also indicated that this stability and control problem could have occurred whether or not natural ice shedding took place. Numerical results using an assumed 3-D ice shape showed an increase of the angle at which this phenomena occurred of about 4 degrees. As it occurred with the 2-D case, the trailing edge separation was observed but started only when the AOA was very close to the angle at which the maximum lift occurred.

Chung, J.↗

Image registration for accurate electrode deformation analysis in operando microscopy of battery materials

Operando imaging techniques have become increasingly valuable in both battery research and manufacturing. However, the reliability of these methods can be compromised by instabilities in the imaging setup and operando cells, particularly when utilizing high-resolution imaging systems. The acquired imaging data often include features arising from both undesirable system vibrations and drift, as well as the scientifically relevant deformations occurring in the battery sample during cell operation. For meaningful analysis, it is crucial to distinguish and separately evaluate these two factors. To address these challenges, we employ a suite of advanced image-processing techniques. These include fast Fourier transform analysis in the frequency domain, power spectrum-based assessments for image quality, as well as rigid and non-rigid image-registration methods. These techniques allow us to identify and exclude blurred images, correct for displacements caused by motor vibrations and sample holder drift and, thus, prevent unwanted image artifacts from affecting subsequent analyses and interpretations. Additionally, we apply optical flow analysis to track the dynamic deformation of battery electrode materials during electrochemical cycling. This enables us to observe and quantify the evolving mechanical responses of the electrodes, offering deeper insights into battery degradation. Together, these methods ensure more accurate image analysis and enhance our understanding of the chemomechanical interplay in battery performance and longevity.

Sun, Tianxiao↗

Three-dimensional boundary layer analysis program Blay and its application

The boundary layer calculation program (BLAY) is a program code which accurately analyzes the three-dimensional boundary layer of a wing with an undefined plane. In comparison with other preexisting programs, the BLAY is characterized by the following: (1) the time required for computation is shorter than any other; (2) the program is adaptable to a parallel processing computer; and (3) the program is associated with a secondary accuracy in the z-direction. As a boundary layer modification to transonic nonviscous flow analysis programs, it is used to adjust viscous and nonviscous interference problems repeatedly. Its efficiency is an important factor in cost reduction in aircraft designing.

Matsuno, K. I.↗

Quantitative Small-Signal Stability Analysis of Grid-Forming Inverters in Large-Scale Power Systems

In this paper, the challenges induced by the integration of inverter-based resources (IBRs) into existing power systems are discussed. Although the grid-forming (GFM) control has been used to enhance system performance in multiple sections of electric power systems (e.g., transmission and distribution systems, microgrids, among others), the coupling mechanism between GFM inverter and various types of other components in modern power systems should be clearly explored. Existing analyses mainly focus on qualitative investigations that verify the effectiveness of deploying GFM inverters for stability enhancement. A holistic system-level quantitative analysis is still missing to verify the detailed benefits of GFM inverters and their interactions with other generation units (e.g., synchronous generators [SGs], conventional grid-following [GFL] inverters, synchronous condensers, among others). To fill in the technical gaps, a holistic small-signal model of the entire system with both rotating generation units and IBRs is established along with the customized power flow analysis for varying steady-state operating points. Further, a quantitative method is proposed to identify the impacts of GFM inverters on system performance with an emphasis on a quantifiable stability margin and a complete modal analysis at the system level. Finally, a realistic test system is used to verify the proposed approach.

grid-following↗

The structure of dilute combusting sprays

An experimental and theoretical study of drop processes in a turbulent flame is described. The experiments involved a monodisperse (105 and 180 micro m initial diameter) stream of methanol drops injected at the base of a turbulent methane-fueled diffusion flame burning in still air. The following measurements were made: mean and fluctuating phase velocities, mean drop number flux, drop-size distributions and mean gas-phase temperatures. Measurements were compared with predictions of two separated flow models: (1) deterministic separated flow, where drop-turbulence interactions are ignored; and (2) stochastic separated flow, where drop-turbulence interactions are considered using random-walk computations. The stochastic separated flow analysis yielded best agreement with measurements, since it provides for turbulent dispersion of drops which was important for present test conditions (and probably for most combusting sprays as well). Distinguishing the presence or absence of envelope flames around the drops, however, was relatively unimportant for present test conditions, since the drops spent most of their lifetime in fuel-rich regions of the flow where this distinction is irrelevant.

Shuen, J. S.↗

User's manual for airfoil flow field computer code SRAIR

A two dimensional unsteady Navier-Stokes calculation procedure with specific application to the isolated airfoil problem is presented. The procedure solves the full, ensemble averaged Navier-Stokes equations with turbulence represented by a mixing length model. The equations are solved in a general nonorthogonal coordinate system which is obtained via an external source. Specific Cartesian locations of grid points are required as input for this code. The method of solution is based upon the Briley-McDonald LBI procedure. The manual discusses the analysis, flow of the program, control steam, input and output.

Shamroth, S. J.↗

Development of free-piston Stirling engine performance and optimization codes based on Martini simulation technique

A FORTRAN computer code is described that could be used to design and optimize a free-displacer, free-piston Stirling engine similar to the RE-1000 engine made by Sunpower. The code contains options for specifying displacer and power piston motion or for allowing these motions to be calculated by a force balance. The engine load may be a dashpot, inertial compressor, hydraulic pump or linear alternator. Cycle analysis may be done by isothermal analysis or adiabatic analysis. Adiabatic analysis may be done using the Martini moving gas node analysis or the Rios second-order Runge-Kutta analysis. Flow loss and heat loss equations are included. Graphical display of engine motions and pressures and temperatures are included. Programming for optimizing up to 15 independent dimensions is included. Sample performance results are shown for both specified and unconstrained piston motions; these results are shown as generated by each of the two Martini analyses. Two sample optimization searches are shown using specified piston motion isothermal analysis. One is for three adjustable input and one is for four. Also, two optimization searches for calculated piston motion are presented for three and for four adjustable inputs. The effect of leakage is evaluated. Suggestions for further work are given.

Martini, William R.↗