Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical Simulation”

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,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Two-dimensional numerical simulation of a Stirling engine heat exchanger

The first phase of an effort to develop multidimensional models of Stirling engine components is described; the ultimate goal is to model an entire engine working space. More specifically, parallel plate and tubular heat exchanger models with emphasis on the central part of the channel (i.e., ignoring hydrodynamic and thermal end effects) are described. The model assumes: laminar, incompressible flow with constant thermophysical properties. In addition, a constant axial temperature gradient is imposed. The governing equations, describing the model, were solved using Crank-Nicloson finite-difference scheme. Model predictions were compared with analytical solutions for oscillating/reversing flow and heat transfer in order to check numerical accuracy. Excellent agreement was obtained for the model predictions with analytical solutions available for both flow in circular tubes and between parallel plates. Also the heat transfer computational results are in good agreement with the heat transfer analytical results for parallel plates.

Ibrahim, Mounir B.↗
NASA NTRS1989

Two-dimensional numerical simulation of a Stirling engine heat exchanger

The first phase of an effort to develop multidimensional models of Stirling engine components is described. The ultimate goal is to model an entire engine working space. Parallel plate and tubular heat exchanger models are described, with emphasis on the central part of the channel (i.e., ignoring hydrodynamic and thermal end effects). The model assumes laminar, incompressible flow with constant thermophysical properties. In addition, a constant axial temperature gradient is imposed. The governing equations describing the model have been solved using the Crack-Nicloson finite-difference scheme. Model predictions are compared with analytical solutions for oscillating/reversing flow and heat transfer in order to check numerical accuracy. Excellent agreement is obtained for flow both in circular tubes and between parallel plates. The computational heat transfer results are in good agreement with the analytical heat transfer results for parallel plates.

Ibrahim, Mounir↗
NASA NTRS1990

Numerical simulation of high-incidence flow over the F-18 aircraft

Navier-Stokes solutions have been obtained using the Chimera overset grid scheme for flow over the wing, fuselage, and wing leading-edge extension (LEX) of the F-18 aircraft at high incidence. Solutions are also presented for flow over the fuselage forebody at high angles of attack. The solutions are for turbulent flows at high-Reynolds number flight-test conditions, and are compared with available qualitative and quantitative experimental data. Comparisons of predicted surface flow patterns, off-surface flow visualizations, and surface-pressure distributions are in good agreement with flight-test data. The ability of the numerical method to predict the bursting of the LEX vortex as it encounters the adverse pressure gradient field of the wing is demonstrated.

Cummings, Russell M.↗
NASA NTRS1991

The role of the sea-surface temperature distribution on numerically simulated cyclogenesis during ERICA

The goal was to quantify the extent to which a sea surface temperature (SST) front can influence cyclogenesis. The approach was to use the Drexel Limited-Area Mesoscale Prediction System (LAMPS) dynamical model to simulate cyclogenesis over various SST fields. Research during the past year focused on the development and testing of a four dimensional data assimilation (FDDA) technique within LAMPS. The technique is a continuous dynamical assimilation where forcing terms are added to the governing model equations to gradually nudge the model solution toward a gridded analysis. Here, the nudging is used as a dynamic initialization tool during a 12 hour preforecast to generate model balanced initial conditions for a subsequent 24 hour numerical prediction. Tests were performed to determine which variables to nudge and how to specify the four dimensional weighting function used to scale the nudging terms. To date, optimal results were obtained by nudging the u and v components of the wind along with the potential temperature. The weighting function ranged from 0 to 1 and varies in time as a quadratic polynomial. It was initialized at 0, reached its maximum at 9 hours into the preforecast, and fell back at 0 to 12 hours. The nudging terms are included in the model equations for all grid points except those within the model predicted oceanic boundary layer. This design attempts to confine changes imposed by the specified SST field to the oceanic boundary layer during the preforecast period.

Lapenta, William M.↗
NASA NTRS1995

Numerical simulation of heat and mass transport during space crystal growth with MEPHISTO

The MEPHISTO space experiments are collaborative United States and French investigations aimed at understanding the fundamentals of crystal growth. Microgravity experiments were conducted aboard the USMP-1 and -2 missions on STS-52 and 62 in October 1992 and March 1994 respectively. MEPHISTO is a French designed and built Bridgman type furnace which uses the Seebeck technique to monitor the solid/liquid interface temperature and Peltier pulsing to mark the location and shape of the solid/liquid interface. In this paper the Bridgman growth of Sn-Bi and Bi-Sn under terrestrial and microgravity conditions is modeled using the finite element code, FIDAP*. The numerical model considers fully coupled heat and mass transport, fluid motion and solid/liquid phase changes in the crystal growth process. The primary goals of this work are: to provide a quantitative study of the thermal buoyancy-induced convection in the melt for the two flight experiments; to compare the vertical and horizontal growth configurations and systematically evaluate the effects of various gravity levels on the solute segregation. Numerical results of the vertical and horizontal Bridgman growth configurations are presented.

Yao, Minwu↗
NASA NTRS1991

Shear-induced instability and arch filament eruption - A magnetohydrodynamic (MHD) numerical simulation

A situation wherein a bipolar magnetic field embedded in a stratified solar atmosphere undergoes symmetrical shear motion at the footpoints is investigated via a 2D (nonplanar) MHD simulation. It was found that the vertical plasma flow velocities grow exponentially, leading to a new type of global MHD instability. The growth rate increases almost linearly until it reaches the same order of magnitude as the Alfven speed. Then a nonlinear MHD instability occurs beyond this point. It was found that the central loops are pinched by opposing Lorentz forces, and the outer closed loops stretch upward with the vertically-rising mass flow. The nonlinear dynamical shearing instability is illustrated by a numerical example that is given for three different values of the plasma beta that span several orders of magnitude.

Wu, S. T.↗
NASA NTRS1998

A Numerical Simulator for Three-Dimensional Flows Through Vibrating Blade Rows

The three-dimensional, multi-stage, unsteady, turbomachinery analysis, TURBO, has been extended to predict the aeroelastic and aeroacoustic response behaviors of a single blade row operating within a cylindrical annular duct. In particular, a blade vibration capability has been incorporated so that the TURBO analysis can be applied over a solution domain that deforms with a vibratory blade motion. Also, unsteady far-field conditions have been implemented to render the computational boundaries at inlet and exit transparent to outgoing unsteady disturbances. The modified TURBO analysis is applied herein to predict unsteady subsonic and transonic flows. The intent is to partially validate this nonlinear analysis for blade flutter applications, via numerical results for benchmark unsteady flows, and to demonstrate the analysis for a realistic fan rotor. For these purposes, we have considered unsteady subsonic flows through a 3D version of the 10th Standard Cascade, and unsteady transonic flows through the first stage rotor of the NASA Lewis, Rotor 67, two-stage fan.

Chuang, H. Andrew↗
DOE OSTI2014

On the equivalence of dynamically orthogonal and bi-orthogonal methods: Theory and numerical simulations

The Karhunen–Lòeve (KL) decomposition provides a low-dimensional representation for random fields as it is optimal in the mean square sense. Although for many stochastic systems of practical interest, described by stochastic partial differential equations (SPDEs), solutions possess this low-dimensional character, they also have a strongly time-dependent form and to this end a fixed-in-time basis may not describe the solution in an efficient way. Motivated by this limitation of standard KL expansion, Sapsis and Lermusiaux (2009) [26] developed the dynamically orthogonal (DO) field equations which allow for the simultaneous evolution of both the spatial basis where uncertainty ‘lives’ but also the stochastic characteristics of uncertainty. Recently, Cheng et al. (2013) [28] introduced an alternative approach, the bi-orthogonal (BO) method, which performs the exact same tasks, i.e. it evolves the spatial basis and the stochastic characteristics of uncertainty. In the current work we examine the relation of the two approaches and we prove theoretically and illustrate numerically their equivalence, in the sense that one method is an exact reformulation of the other. We show this by deriving a linear and invertible transformation matrix described by a matrix differential equation that connects the BO and the DO solutions. We also examine a pathology of the BO equations that occurs when two eigenvalues of the solution cross, resulting in an instantaneous, infinite-speed, internal rotation of the computed spatial basis. We demonstrate that despite the instantaneous duration of the singularity this has important implications on the numerical performance of the BO approach. On the other hand, it is observed that the BO is more stable in nonlinear problems involving a relatively large number of modes. Several examples, linear and nonlinear, are presented to illustrate the DO and BO methods as well as their equivalence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1996

Numerical simulation of iodine speciation in relation to water disinfection aboard manned spacecraft I. Equilibria

Elemental iodine (I2) is currently used as the drinking water disinfectant aboard the Shuttle Orbiter and will also be incorporated into the water recovery and distribution system for the International Space Station Alpha. Controlled release of I2 is achieved using the Microbial Check Valve (MCV), a flow-through device containing an iodinated polymer which imparts a bacteriostatic residual concentration of approximately 2mg/L to the aqueous stream. During regeneration of MCV canisters, I2 concentrations of approximately 300 mg/L are used. Dissolved iodine undergoes a series of hydrolytic disproportionation and related reactions which result in the formation of an array of inorganic species including: I-, I3-, HOI, OI-, IO3-, HIO3, I2OH-, I2O(-2), and H2OI+. Numerical estimation of the steady-state distribution of inorganic iodine containing species in pure water at 25 degrees C has been achieved by simultaneous solution of the multiple equilibrium expressions as a function of pH. The results are reported herein.

NASA Discipline Environmental Health↗
NASA NTRS1994

Meso-beta scale numerical simulation studies of terrain-induced jet streak mass/momentum perturbations

An in-depth analysis of observed gravity waves and their relationship to precipitation bands over the Montana mesonetwork during the 1981 CCOPE case study indicates that there were two episodes of coherent internal gravity waves. One of the fundamental unanswered questions from this research, however, concerns the dynamical processes which generated the observed waves, all of which originated from the region encompassing the borders of Montana, Idaho, and Wyoming. While geostrophic adjustment, shearing instability, and terrain where all implicated separately or in concert as possible wave generation mechanisms, the lack of upper-air data within the wave genesis region made it difficult to rigorously define the genesis processes from observations alone. In this report we employ a mesoscale numerical model to help diagnose the intricate early wave generation mechanisms during the first observed wave episode.

Lin, Yuh-Lang↗
NASA NTRS1988

Numerical simulation of the integrated space shuttle vehicle in ascent

A simulation of the flow about the integrated space shuttle vehicle in ascent mode has been undertaken for various flight conditions using the Chimera composite grid discretization approach. Overset body-conforming grids were used to represent each geometric component, and an implicit approximately factored finite-difference procedure was used to solve the three-dimensional thin-layer Navier-Stokes equations. The computational results have been compared with both wind tunnel and flight test data. Although relatively good agreement is obtained with the experimental data, further refinement and evaluation of numerical error is under way.

Buning, P. G.↗
NASA NTRS1986

Numerical simulation of electrophoresis separation processes

A new Petrov-Galerkin finite element formulation has been proposed for transient convection-diffusion problems. Most Petrov-Galerkin formulations take into account the spatial discretization, and the weighting functions so developed give satisfactory solutions for steady state problems. Though these schemes can be used for transient problems, there is scope for improvement. The schemes proposed here, which consider temporal as well as spatial discretization, provide improved solutions. Electrophoresis, which involves the motion of charged entities under the influence of an applied electric field, is governed by equations similiar to those encountered in fluid flow problems, i.e., transient convection-diffusion equations. Test problems are solved in electrophoresis and fluid flow. The results obtained are satisfactory. It is also expected that these schemes, suitably adapted, will improve the numerical solutions of the compressible Euler and the Navier-Stokes equations.

Ganjoo, D. K.↗
DOE OSTI2020

Aeroelastic oscillations of a pitching flexible wing with structural geometric nonlinearities: Theory and numerical simulation

In this paper, we focus on the derivation of an analytical model of the aeroelastic dynamics of an elastically mounted flexible wing. This equations of motion obtained serve to help understand the behaviour of the aeroelastic wind tunnel setup in question, which consists of a rectangular wing with a uniform NACA 0012 airfoil profile, whose base is free to rotate rigidly about a longitudinal axis. Of particular interest are the structural geometric nonlinearities primarily introduced by the coupling between the rigid body pitch degree-of-freedom and the continuous system. A coupled system of partial differential equations (PDEs) coupled with an ordinary differential equation (ODE) describing axial-bending-bending-torsion-pitch motion is derived using Hamilton's principle. A finite dimensional approximation of the system of coupled differential equations is obtained using the Galerkin method, leading to a system of coupled nonlinear ODEs. Subsequently, these nonlinear ODEs are solved numerically using Houbolt's method. The results that are obtained are verified by comparison with the results obtained by direct integration of the equations of motion using a finite difference scheme. Adopting a linear unsteady aerodynamic model, it is observed that the system undergoes coalescence flutter due to coupling between the rigid body pitch rotation dominated mode and the first flapwise bending dominated mode. Finally, the behaviour of the limit cycle oscillations is primarily influenced by the structural geometric nonlinear terms in the coupled system of PDEs and ODE.

42 ENGINEERING↗
DOE OSTI2019

Characterization and numerical simulation of liquid refrigerant R-134a flow emerging from a flooded evaporator tube bundle

The distribution of liquid droplets emerging from an evaporator tube bundle is characterized for refrigerant R-134a with a triangular tube arrangement with a pitch of 1.167. The purpose of this research was to improve understanding of the droplet ejection process to aid in design of evaporators typically used in larger chiller systems. A laser and camera system captured images of the evaporator headspace at varying conditions. Conventional shadowgraphy techniques were applied to recognize and match droplets for velocity calculations. The evaporator conditions varied with bundle mass fluxes of 20.3 and 40.7 kg s -1 m -2 , top-rows heat fluxes of 15.8 and 31.5 kWm -2 , and outlet saturation temperatures of 4.4 and 12.8 °C. Conditions ranged from flooded to dryout of the top rows. Droplet number, size distribution, velocity, and liquid volume fraction are presented in the headspace above the bundle. A method to numerically duplicate the droplet loading in the headspace using CFD with a Lagrangian discrete-phase model is also presented and verified, providing a powerful design tool. Liquid distribution in the headspace is found to be a strong function of all varied properties, particularly mass flux, liquid level, and saturation temperature

42 ENGINEERING↗
NASA NTRS1994

Numerical Simulation of Shock-Induced Combustion Past Blunt Bodies Using Shock-Fitting Technique

Two-dimensional axisymmetric, reacting viscous flow over blunt projectiles is computed to study shock induced combustion at Mach 5.11 and Mach 6.46 in hydrogen-air mixture. A finite-difference, shock-fitting method is used to solve the complete set of Navier Stokes and species conservation equations. In this approach, the bow shock represents a boundary of the computational domain and is treated as a discontinuity across which Rankine-Hugoniot conditions are applied. All interior details of the flow such as compression waves, reaction front, and the wall boundary layer are captured automatically in the solution. Since shock-fitting approach reduces the amount of artificial dissipation, all the intricate details of the flow are captured much more clearly than has been possible with the shock-capturing approach. This has allowed an improved understanding of the physics of shock-induced combustion over blunt projectiles and the numerical results can now be explained more readily with one dimensional wave-interaction model than before.

Ahuja, J. K.↗
NASA NTRS1994

Numerical simulation of shock-induced combustion past blunt bodies using shock-fitting technique

Two-dimensional axisymmetric, reacting viscous flow over blunt projectiles is computed to study shock-induced combustion at Mach 5.11 and Mach 6.46 in hydrogen-air mixture. A finite-difference, shock-fitting method is used to solve the complete set of Navier-Stokes and species conservation equations. In this approach, the bow shock represents a boundary of the computational domain and is treated as a discontinuity across which Rankine-Hugoniot conditions are applied. All interior details of the flow such as compression waves, reaction front, and the wall boundary layer are captured automatically in the solution. Since shock-fitting approach reduces the amount of artificial dissipation, all the intricate details of the flow are captured much more clearly than has been possible with the shock-capturing approach. This has allowed an improved understanding of the physics of shock-induced combustion over blunt projectiles and the numerical results can now be explained more readily with one-dimensional wave-interaction model than before.

Ahuja, J. K.↗
NASA NTRS1997

Preconditioning for Numerical Simulation of Low Mach Number Three-Dimensional Viscous Turbomachinery Flows

A preconditioning scheme has been implemented into a three-dimensional viscous computational fluid dynamics code for turbomachine blade rows. The preconditioning allows the code, originally developed for simulating compressible flow fields, to be applied to nearly-incompressible, low Mach number flows. A brief description is given of the compressible Navier-Stokes equations for a rotating coordinate system, along with the preconditioning method employed. Details about the conservative formulation of artificial dissipation are provided, and different artificial dissipation schemes are discussed and compared. The preconditioned code was applied to a well-documented case involving the NASA large low-speed centrifugal compressor for which detailed experimental data are available for comparison. Performance and flow field data are compared for the near-design operating point of the compressor, with generally good agreement between computation and experiment. Further, significant differences between computational results for the different numerical implementations, revealing different levels of solution accuracy, are discussed.

Tweedt, Daniel L.↗
NASA NTRS1994

Brush seal numerical simulation: Concepts and advances

The development of the brush seal is considered to be most promising among the advanced type seals that are presently in use in the high speed turbomachinery. The brush is usually mounted on the stationary portions of the engine and has direct contact with the rotating element, in the process of limiting the 'unwanted' leakage flows between stages, or various engine cavities. This type of sealing technology is providing high (in comparison with conventional seals) pressure drops due mainly to the high packing density (around 100 bristles/sq mm), and brush compliance with the rotor motions. In the design of modern aerospace turbomachinery leakage flows between the stages must be minimal, thus contributing to the higher efficiency of the engine. Use of the brush seal instead of the labyrinth seal reduces the leakage flow by one order of magnitude. Brush seals also have been found to enhance dynamic performance, cost less, and are lighter than labyrinth seals. Even though industrial brush seals have been successfully developed through extensive experimentation, there is no comprehensive numerical methodology for the design or prediction of their performance. The existing analytical/numerical approaches are based on bulk flow models and do not allow the investigation of the effects of brush morphology (bristle arrangement), or brushes arrangement (number of brushes, spacing between them), on the pressure drops and flow leakage. An increase in the brush seal efficiency is clearly a complex problem that is closely related to the brush geometry and arrangement, and can be solved most likely only by means of a numerically distributed model.

Braun, M. J.↗