Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “equation of fluid dynamics”

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,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Mean Flow from Phase Averages in the 2D Boussinesq Equations

The atmosphere and ocean are described by highly oscillatory PDEs that challenge both our understanding of their dynamics and their numerical approximation. This paper presents a preliminary numerical study of one type of phase averaging applied to mean flows in the 2D Boussinesq equations that also has application to numerical methods. The phase averaging technique, well-known in dynamical systems theory, relies on a mapping using the exponential operator, and then an averaging over the phase. The exponential operator has connections to the Craya–Herring basis pioneered by Jack Herring to study the fluid dynamics of oscillatory, nonlinear fluid dynamics. In this paper, we perform numerical experiments to study the effect of this averaging technique on the time evolution of the solution. We explore its potential as a definition for mean flows. We also show that, as expected from theory, the phase-averaging method can reduce the magnitude of the time rate of change in the PDEs, making them potentially suitable for time stepping methods.

Wingate, Beth A.↗
NASA NTRS1983

The impact of CFD on development test facilities - A National Research Council projection

The results of a National Research Council study on the effect that advances in computational fluid dynamics (CFD) will have on conventional aeronautical ground testing are reported. Current CFD capabilities include the depiction of linearized inviscid flows and a boundary layer, initial use of Euler coordinates using supercomputers to automatically generate a grid, research and development on Reynolds-averaged Navier-Stokes (N-S) equations, and preliminary research on solutions to the full N-S equations. Improvements in the range of CFD usage is dependent on the development of more powerful supercomputers, exceeding even the projected abilities of the NASA Numerical Aerodynamic Simulator (1 BFLOP/sec). Full representation of the Re-averaged N-S equations will require over one million grid points, a computing level predicted to be available in 15 yr. Present capabilities allow identification of data anomalies, confirmation of data accuracy, and adequateness of model design in wind tunnel trials. Account can be taken of the wall effects and the Re in any flight regime during simulation. CFD can actually be more accurate than instrumented tests, since all points in a flow can be modeled with CFD, while they cannot all be monitored with instrumentation in a wind tunnel.

Korkegi, R. H.↗
NASA NTRS1989

Euler/Navier-Stokes calculations of transonic flow past fixed- and rotary-wing aircraft configurations

Computational fluid dynamics has an increasingly important role in the design and analysis of aircraft as computer hardware becomes faster and algorithms become more efficient. Progress is being made in two directions: more complex and realistic configurations are being treated and algorithms based on higher approximations to the complete Navier-Stokes equations are being developed. The literature indicates that linear panel methods can model detailed, realistic aircraft geometries in flow regimes where this approximation is valid. As algorithms including higher approximations to the Navier-Stokes equations are developed, computer resource requirements increase rapidly. Generation of suitable grids become more difficult and the number of grid points required to resolve flow features of interest increases. Recently, the development of large vector computers has enabled researchers to attempt more complex geometries with Euler and Navier-Stokes algorithms. The results of calculations for transonic flow about a typical transport and fighter wing-body configuration using thin layer Navier-Stokes equations are described along with flow about helicopter rotor blades using both Euler/Navier-Stokes equations.

Deese, J. E.↗
NASA NTRS1990

Time domain flutter analysis of cascades using a full-potential solver

A time domain approach is used to determine the dynamic aeroelastic stability of a cascade of blades. The structural model for each blade is a typical section with two degrees of freedom. The aerodynamic model is the unsteady, two-dimensional, full-potential flow through the cascade of airfoils. The unsteady equations of motion for the structure and the fluid are integrated simultaneously in time starting with the steady flowfield and a small initial disturbance applied to the airfoils. The motion of each blade is analyzed to determine the aeroelastic stability of the cascade. The effect of interblade phase angle is included in the analysis by allowing each blade to have an independent motion and considering a number of blade passages. Calculations are made using an airfoil section and structural parameters that are representative of a propfan. The results are compared with those from a separate frequency domain analysis. Good agreement between the results is observed. With the time domain approach, it is possible to consider nonlinear structural models and nonlinear force-displacement relations. The method allows a realistic simulation of the motion of the fluid and the cascade blades for a better physical understanding and it also has the potential for saving computational time when compared to the frequency domain approach for the flutter analysis of cascades.

Bakhle, Milind A.↗
NASA NTRS1990

CFD validation for aerodynamic flows - Challenge for the '90's

The process of computational fluid dynamics (CFD) validation is described from two perspectives, numerical error validation and physical model validation. Errors associated with time and space discretization, grid refinement, numerical dissipation, and level of solution convergence, are all possible sources of numerical error. The second aspect of CFD code validation is associated with the governing equations and the physical models, e.g., chemistry and turbulence models. This type of validation requires comparisons with carefully conducted experiments and is the primary aspect of the validation process discussed in this paper. Examples of validation are shown for selected Reynolds-averaged Navier-Stokes codes for conditions from transonic through hypersonic. Since code applications are becoming more complex, it no longer suffices to use data from surface or integral quantities alone to provide the required validation. Flow field surveys and experimental boundary condition measurements are emerging as critical data that must be obtained for CFD code validation. Progress in the overall level of validation is shown to be improving. However, future challenges remain, some of which are enumerated.

Marvin, Joseph G.↗
NASA NTRS1991

Active control of panel vibrations induced by boundary-layer flow

Some problems in active control of panel vibration excited by a boundary layer flow over a flat plate are studied. In the first phase of the study, the optimal control problem of vibrating elastic panel induced by a fluid dynamical loading was studied. For a simply supported rectangular plate, the vibration control problem can be analyzed by a modal analysis. The control objective is to minimize the total cost functional, which is the sum of a vibrational energy and the control cost. By means of the modal expansion, the dynamical equation for the plate and the cost functional are reduced to a system of ordinary differential equations and the cost functions for the modes. For the linear elastic plate, the modes become uncoupled. The control of each modal amplitude reduces to the so-called linear regulator problem in control theory. Such problems can then be solved by the method of adjoint state. The optimality system of equations was solved numerically by a shooting method. The results are summarized.

Chow, Pao-Liu↗
NASA NTRS1993

Computational reacting gas dynamics

In the study of high speed flows at high altitudes, such as that encountered by re-entry spacecrafts, the interaction of chemical reactions and other non-equilibrium processes in the flow field with the gas dynamics is crucial. Generally speaking, problems of this level of complexity must resort to numerical methods for solutions, using sophisticated computational fluid dynamics (CFD) codes. The difficulties introduced by reacting gas dynamics can be classified into three distinct headings: (1) the usually inadequate knowledge of the reaction rate coefficients in the non-equilibrium reaction system; (2) the vastly larger number of unknowns involved in the computation and the expected stiffness of the equations; and (3) the interpretation of the detailed reacting CFD numerical results. The research performed accepts the premise that reacting flows of practical interest in the future will in general be too complex or 'untractable' for traditional analytical developments. The power of modern computers must be exploited. However, instead of focusing solely on the construction of numerical solutions of full-model equations, attention is also directed to the 'derivation' of the simplified model from the given full-model. In other words, the present research aims to utilize computations to do tasks which have traditionally been done by skilled theoreticians: to reduce an originally complex full-model system into an approximate but otherwise equivalent simplified model system. The tacit assumption is that once the appropriate simplified model is derived, the interpretation of the detailed numerical reacting CFD numerical results will become much easier. The approach of the research is called computational singular perturbation (CSP).

Lam, S. H.↗
NASA NTRS2013

A Computational Fluid Dynamics Study of Swirling Flow Reduction by using Anti-vortex Baffle

An anti-vortex baffle is a liquid propellant management device placed adjacent to an outlet of the propellant tank. Its purpose is to substantially reduce or eliminate the formation of free surface dip and vortex, as well as prevent vapor ingestion into the outlet, as the liquid drains out through the flight. To design an effective anti-vortex baffle, Computational Fluid Dynamic (CFD) simulations were undertaken for the NASA Ares I vehicle LOX tank subjected to the simulated flight loads with and without the anti-vortex baffle. The Six Degree-Of-Freedom (6- DOF) dynamics experienced by the Crew Launch Vehicle (CLV) during ascent were modeled by modifying the momentum equations in a CFD code to accommodate the extra body forces from the maneuvering in a non-inertial frame. The present analysis found that due to large moments, the CLV maneuvering has significant impact on the vortical flow generation inside the tank. Roll maneuvering and side loading due to pitch and yaw are shown to induce swirling flow. The vortical flow due to roll is symmetrical with respect to the tank centerline, while those induced by pitch and yaw maneuverings showed two vortices side by side. The study found that without the anti-vortex baffle, the swirling flow caused surface dip during the late stage of drainage and hence early vapor ingestion. The flow can also be non-uniform in the drainage pipe as the secondary swirling flow velocity component can be as high as 10% of the draining velocity. An analysis of the vortex dynamics shows that the swirling flow in the drainage pipe during the Upper Stage burn is mainly the result of residual vortices inside the tank due to conservation of angular momentum. The study demonstrated that the swirling flow in the drainage pipe can be effectively suppressed by employing the anti-vortex baffle.

Yang, H. Q.↗
NASA NTRS2017

A Computational Fluid Dynamics Study of Swirling Flow Reduction by Using Anti-Vortex Baffle

An anti-vortex baffle is a liquid propellant management device placed adjacent to an outlet of the propellant tank. Its purpose is to substantially reduce or eliminate the formation of free surface dip and vortex, as well as prevent vapor ingestion into the outlet, as the liquid drains out through the flight. To design an effective anti-vortex baffle, Computational Fluid Dynamic (CFD) simulations were undertaken for the NASA Ares I vehicle LOX tank subjected to the simulated flight loads with and without the anti-vortex baffle. The Six Degree-Of-Freedom (6-DOF) dynamics experienced by the Crew Launch Vehicle (CLV) during ascent were modeled by modifying the momentum equations in a CFD code to accommodate the extra body forces from the maneuvering in a non-inertial frame. The present analysis found that due to large moments, the CLV maneuvering has a significant impact on the vortical flow generation inside the tank. Roll maneuvering and side loading due to pitch and yaw are shown to induce swirling flow. The vortical flow due to roll is symmetrical with respect to the tank centerline, while those induced by pitch and yaw maneuverings showed two vortices side by side. The study found that without the anti-vortex baffle, the swirling flow caused surface dip during the late stage of drainage and hence early vapor ingestion. The flow can also be non-uniform in the drainage pipe as the secondary swirling flow velocity component can be as high as 10% of the draining velocity. An analysis of the vortex dynamics shows that the swirling flow in the drainage pipe during the Upper Stage burn is mainly the result of residual vortices inside the tank due to the conservation of angular momentum. The study demonstrated that the swirling flow in the drainage pipe can be effectively suppressed by employing the anti-vortex baffle.

Yang, H. Q.↗
NASA NTRS2017

A Computational Fluid Dynamics Study of Swirling Flow Reduction by Using Anti-Vortex Baffle

An anti-vortex baffle is a liquid propellant management device placed adjacent to an outlet of the propellant tank. Its purpose is to substantially reduce or eliminate the formation of free surface dip and vortex, as well as prevent vapor ingestion into the outlet, as the liquid drains out through the flight. To design an effective anti-vortex baffle, Computational Fluid Dynamic (CFD) simulations were undertaken for the NASA Ares I vehicle LOX tank subjected to the simulated flight loads with and without the anti-vortex baffle. The Six Degree-Of-Freedom (6-DOF) dynamics experienced by the Crew Launch Vehicle (CLV) during ascent were modeled by modifying the momentum equations in a CFD code to accommodate the extra body forces from the maneuvering in a non-inertial frame. The present analysis found that due to large moments, the CLV maneuvering has a significant impact on the vortical flow generation inside the tank. Roll maneuvering and side loading due to pitch and yaw are shown to induce swirling flow. The vortical flow due to roll is symmetrical with respect to the tank centerline, while those induced by pitch and yaw maneuverings showed two vortices side by side. The study found that without the anti-vortex baffle, the swirling flow caused surface dip during the late stage of drainage and hence early vapor ingestion. The flow can also be non-uniform in the drainage pipe as the secondary swirling flow velocity component can be as high as 10% of the draining velocity. An analysis of the vortex dynamics shows that the swirling flow in the drainage pipe during the Upper Stage burn is mainly the result of residual vortices inside the tank due to the conservation of angular momentum. The study demonstrated that the swirling flow in the drainage pipe can be effectively suppressed by employing the anti-vortex baffle.

Yang, H. Q.↗
NASA NTRS2019

Full Facility Shock Frame Simulations of the Electric Arc Shock Tube

Radiative heating computations are performed for high speed lunar return experiments conducted in the Electric Arc Shock Tube (EAST) facility at NASA Ames Research Center. The nonequilibrium radiative transport equations are solved via NASA's in-house radiation code NEQAIR using flow field input from US3D flow solver. The post-shock flow properties for the 10 km/s Earth entry conditions are computed using the stagnation line of a blunt-body and a full facility CFD (Computational Fluid Dynamics) simulation of the EAST shock tube. The shocked gas in the blunt-body flow achieves a thermochemical equilibrium away from the shock front whereas EAST flow exhibits a nonequilibrium behavior due to strong viscous dissipation of the shock by boundary layer. The full-tube flow calculations capture the influence of the boundary layer on the shocked gas state and provide a realistic fluid dynamic input for the radiative predictions. The integrated radiance behind the shock is calculated in NEQAIR for wavelength regimes from Vacuum-UltraViolet (VUV) to InfraRed (IR), which are pertinent to the emission characteristics of high enthalpy shock waves in air. These radiance profiles are validated against corresponding EAST shots. The full-tube simulations successfully predict a sharp radiance peak at the shock front which gets smeared in the test data due to the spatial resolution in the measurements. The full facility based radiance behind the shock shows a slightly better match with the test data in the VUV and Red spectral regions, as compared to that from a blunt-body based predictions. The UV radiance is very similar for both geometries and under-predicts the test behavior. The IR test data matches better with the blunt-body based predictions where the full-tube simulations show a significant over-prediction.

Radiation↗
DOE OSTI2024

Investigation of inter-subject variation in ultrafine particle deposition across human nasal airways: A study involving children, adults, and the elderly

Occupational and environmental exposure to toxic nanoparticles, driven by the rapid expansion of nanotechnology, raises significant respiratory health concern. Numerous studies have explored airflow and particle dynamics in adult nasal airways, but understanding the impact of age-related anatomical changes in children and the elderly remains limited. This study systematically investigates age-related anatomical variations and associated influence on nasal airflow dynamics and ultrafine particle deposition characteristics. Using Computational Fluid-Particle Dynamics (CFPD) method, simulation was conducted under diverse inhalation conditions spanning a wide age range, including: two children (5 years old), two young adults (in their twenties), and two elderly (over 77 years old). Our results reveal distinctive variations across age groups in anatomical dimensions, which affect distribution of wall shear stress where the elderly and children display unique patterns distinct from the young adults. While total deposition efficiency differs significantly between children and adults, filtration efficiency in the subregion with most deposition, main respiratory, remains consistent. However, inter-subject differences are observed in the vestibular and olfactory regions,emphasizing nuanced impact of age-related anatomical variations. Overall and subregional empirical equations for deposition efficiency were developed by incorporating the combined diffusion parameter, Sc a Δ b , corroborating the use of geometrical characteristic parameters for each specific subject in predicting nasal deposition efficiency across age groups. Our findings contribute to predictive nanoparticle exposure analysis in nasal airways across different age groups, thereby enhancing respiratory healthcare for individuals across the life span.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Development of a conduction-based model for analyzing frozen startup of alkali-metal heat pipes

One key area of interest in heat pipe modeling/simulation is to analyze the startup behavior of the liquid-metal heat pipes (LMHPs) from a frozen state. This so-called ‘frozen startup’ process involves a complex set of nonlinear mass and heat transport phenomena, including phase transitions from solid to liquid and vapor, multiphase interactions, microporous wick flow, and compressible vapor dynamics. The complexity of these processes makes it challenging to simulate LMHP’s frozen startup using conventional numerical methods or commercial computational fluid dynamics (CFD) software. This paper presents a simplified conduction-based modeling approach that can provide practical insights into the entire LMHP frozen startup process, while alleviating the challenges of modeling its complex physics. The theoretical foundation and physical assumptions of the proposed model are based solely on heat-conduction equation, allowing for a more tractable simulation without sacrificing essential physical accuracy. The proposed model was implemented in a commercial CFD software, and its prediction was compared with the experimental data obtained from sodium heat-pipe startup experiments. The comparison highlights the proposed model's ability to capture the transient thermal behavior of LMHP during frozen startup. This study not only validates the conduction-based frozen startup modeling method but also shows its potential as a practical and efficient tool for understanding the startup performance of the LMHP systems.

Microreactor↗
DOE OSTI2026

Recurrent convolutional neural networks for modeling nonadiabatic dynamics of quantum-classical systems

Recurrent neural networks (RNNs) have recently been extensively applied to model the time evolution in fluid dynamics, weather predictions, and even chaotic systems due to their ability to capture temporal dependencies and sequential patterns in data. Here we present an RNN model based on convolutional neural networks for modeling the nonlinear nonadiabatic dynamics of hybrid quantum-classical systems. The dynamical evolution of the hybrid systems is governed by equations of motion for classical degrees of freedom and von Neumann equation for electrons. The Physics-Aware Recurrent Convolution (PARC) neural network structure incorporates a differentiator-integrator architecture that inductively models the spatiotemporal dynamics of generic physical systems. Here, we apply our RNN approach to learn the space-time evolution of a one-dimensional semiclassical Holstein model after an interaction quench. For shallow quenches (small changes in electron-lattice coupling), the deterministic dynamics can be accurately captured using a single-CNN-based recurrent network. In contrast, deep quenches induce chaotic evolution, making long-term trajectory prediction significantly more challenging. Nonetheless, we demonstrate that the PARC-CNN architecture can effectively learn the statistical climate of the Holstein model under deep-quench conditions.

Holstein model↗
DOE OSTI2023

Dynamic scaling of order parameter fluctuations in model B

Brookhaven National Laboratory has performed a dedicated set of experiments, called the Beam Energy Scan (BES), with the goal of location a possible critical point in the phase diagram of nuclear matter, This point is analogous to the endpoint of the liquid-gas phase transition in water, where the liquid phase corresponds to a phase of hadrons, and the gas phase to a plasma of quarks and gluons. In order to describe the data generated by these experiments we have studied numerical simulations of the stochastic diffusion equation with a conserved charge. Here, we focus on the dynamics in the vicinity of a critical point in the Ising universality class. The model we consider is expected to describe the critical dynamics near a possible QCD critical point if the coupling of the order parameter to the momentum density of the fluid can be neglected. The simulations are performed on a spatial lattice, and the time evolution is performed using a Metropolis algorithm. We determine the dynamical critical exponent z ≃ 3.972(2), which agrees with predictions of the epsilon expansion. We also study non-equilibrium sweeps of the reduced temperature and observe approximate Kibble-Zurek scaling.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Perovskite Sorbent Oxygen Separation Modeling with MFiX

This document chronicles the development and implementation of computational kinetic rate models that capture absorption and desorption characteristics of the National Energy Technology Laboratory (NETL) developed perovskite, Sr 1-x Ca x FeO 3-δ . Two paths to create accurate kinetic rates were followed: (1) an isothermal rate approach where thermogravimetric (TGA) data are recast as oxygen capacities through a pseudo-second order Lagergren equation (He et al., 2009); and (2) a more traditional Arrhenius approach where experimental data are fit with a power law model to derive associate activation energies (Bulfin et al., 2020a). For reference, the mathematics and associate experimental strategies that support these derivations are included in this report. In addition, computational fluid dynamics (CFD) models were developed to utilize both kinetic rate formulations and applied to simulate oxygen uptake and release in small scale scenarios. The program Multiphase Flow with interphase eXchanges (MFiX) was used to create: (1) discrete element method (DEM) simulations of a single tube of granular perovskite experiencing isothermal O 2 -absorption; and desorption and (2) two-fluid-model (TFM) non-isothermal simulations of perovskite O 2 -absorption and desorption tubes that share a wall. Conjugate heat transfer between steel walled tubes and the perovskite bed were managed with user-defined functions. As the project moves to simulating larger scale devices that will require more robust conjugate heat transfer methods, developed kinetic rates and associate methodologies have been recast for use in the ANSYS Fluent CFD program.

36 MATERIALS SCIENCE↗
NASA NTRS1982

Tenth NASTRAN User's Colloquium

The development of the NASTRAN computer program, a general purpose finite element computer code for structural analysis, was discussed. The application and development of NASTRAN is presented in the following topics: improvements and enhancements; developments of pre and postprocessors; interactive review system; the use of harmonic expansions in magnetic field problems; improving a dynamic model with test data using Linwood; solution of axisymmetric fluid structure interaction problems; large displacements and stability analysis of nonlinear propeller structures; prediction of bead area contact load at the tire wheel interface; elastic plastic analysis of an overloaded breech ring; finite element solution of torsion and other 2-D Poisson equations; new capability for elastic aircraft airloads; usage of substructuring analysis in the get away special program; solving symmetric structures with nonsymmetric loads; evaluation and reduction of errors induced by Guyan transformation.

Source record↗
NASA NTRS1987

Status of computational fluid dynamics in the United States

CFD-related progress in U.S. aerospace industries and research institutions is evaluated with respect to methods employed, their applications, and the computer technologies employed in their implementation. Goals for subsonic CFD are primarily aimed at greater fuel efficiency; those of supersonic CFD involve the achievement of high sustained cruise efficiency. Transatmospheric/hypersonic vehicles are noted to have recently become important concerns for CFD efforts. Attention is given to aspects of discretization, Euler and Navier-Stokes general purpose codes, zonal equation methods, internal and external flows, and the impact of supercomputers and their networks in advancing the state-of-the-art.

Kutler, Paul↗