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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 919 records · Page 51

All sources NASA NTRS DOE OSTI

NASA NTRS1985

Numerical simulation of vortex breakdown

The breakdown of an isolated axisymmetric vortex embedded in an unbounded uniform flow is examined by numerical integration of the complete Navier-Stokes equations for unsteady axisymmetric flow. Results show that if the vortex strength is small, the solution approaches a steady flow and the vortex is stable. If the strength is large enough, the solution remains unsteady and a recirculating zone will appear near the axis, its form and internal structure resembling those of the axisymmetric breakdown bubbles with multi-cells observed by Faler and Leibovich (1978). For apppropriate combinations of flow parameters, the flow reveals quasi-periodicity. Parallel calculations with the quasi-cylindrical approximation indicate that so far as predicting of breakdown is concerned, its results coincide quite well with the results mentioned above. Both show that the vortex breakdown has little concern with the Reynolds number or with the critical classification of the upstream flow, at least for the lower range of Reynolds numbers.

Shi, X.↗

NASA NTRS1988

A 3-dimensional numerical simulation of the atmospheric injection of aerosols by a hypothetical basaltic fissure eruption

Researchers simulated the atmospheric response to a hypothetical basaltic fissure eruption using heating rates based on the Roza flow eruption. The simulation employs the Colorado State University Regional Atmospheric Model (RAMS) with scavenging effects. The numerical model is a three-dimensional non-hydrostatic time-split compressible cloud/mesoscale model. Explicit microphysics include prediction of cloud, rain, crystal, and hail precipitation types. Nucleation and phoretic scavenging are predicted assuming that the pollutant makes an effective cloud droplet nucleus. Smoke is carried as a passive tracer. Long and short wave radiation heating tendencies, including the effects of the smoke, are parameterized. The longwave emission by the lava surface is neglected in the parameterization and included as an explicit heating term instead. A regional scale domain of 100 x 100 km in the horizontal and 22 km high is used. The horizontal grid spacing is taken to be 2 km and the vertical spacing is taken to be 0.75 km. The initial atmospheric state is taken to be horizontally homogenous and based on the standard atmospheric sounding. The fissure is assumed to be 90 km long and oriented in a zig/zag pattern.

Tripoli, Gregory J.↗

NASA NTRS1990

Numerical simulation of spatially-evolving instability

A computational study of the spatial stability of plane Poiseuille flow is presented. The numerical scheme employs a time-splitting method to integrate the full Navier-Stokes equations using spectral collocation/finite-difference discretization on a non-staggered mesh. The eigenvalue decomposition procedure is applied for the solution of the Poisson equations using the capacitance matrix technique. The buffer domain method is incorporated for the outflow boundary conditions. The input perturbation velocities are obtained by solving the Orr-Sommerfeld equation for the nonlinear eigenvalue problem employing the companion matrix method. Computational results are compared with the linear theory for two-dimensional disturbances.

Danabasoglu, G.↗

NASA NTRS1992

Numerical simulation of a hovering rotor using embedded grids

The flow field for a rotor blade in hover was computed by numerically solving the compressible thin-layer Navier-Stokes equations on embedded grids. In this work, three embedded grids were used to discretize the flow field - one for the rotor blade and two to convect the rotor wake. The computations were performed at two hovering test conditions, for a two-bladed rectangular rotor of aspect ratio six. The results compare fairly with experiment and illustrates the use of embedded grids in solving helicopter type flow fields.

Duque, Earl-Peter N.↗

NASA NTRS2005

Aspects of Numerical Simulation of Circulation Control Airfoils

The mass-averaged compressible Navier-Stokes equations are solved for circulation control airfoils. Numerical solutions are computed with a multigrid method that uses an implicit approximate factorization smoother. The effects of flow conditions (e.g., free-stream Mach number, angle of attack, momentum coefficient) and mesh on the prediction of circulation control airfoil flows are considered. In addition, the impact of turbulence modeling, including curvature effects and modifications to reduce eddy viscosity levels in the wall jet (i.e., Coanda flow), is discussed. Computed pressure distributions are compared with available experimental data.

Swanson, R. C.↗

NASA NTRS1992

Numerical simulation of a 325 Mbit/s QPPM optical communication system

Optical digital receivers are being considered for intersatellite laser communication links. A demonstration system is being designed to operate at 325 MBit/s, using quaternary pulse position modulation (QPPM). Laboratory experiments have been conducted using a 50 MBit/s prototype system. A numerical model has been developed to simulate a QPPM optical receiver. The 50 MBit/s system was simulated to verify the validity of the model. The model was then used to simulate the projected behavior of the 325 MBit/s system. The model predicts a bit error rate of 10.6 at 38 incident photons per bit for 820 nm light.

Martino, Anthony J.↗

DOE OSTI2016

A novel Lagrangian approach for the stable numerical simulation of fault and fracture mechanics

The simulation of the mechanics of geological faults and fractures is of paramount importance in several applications, such as ensuring the safety of the underground storage of wastes and hydrocarbons or predicting the possible seismicity triggered by the production and injection of subsurface fluids. However, the stable numerical modeling of ground ruptures is still an open issue. The present work introduces a novel formulation based on the use of the Lagrange multipliers to prescribe the constraints on the contact surfaces. The variational formulation is modified in order to take into account the frictional work along the activated fault portion according to the principle of maximum plastic dissipation. The numerical model, developed in the framework of the Finite Element method, provides stable solutions with a fast convergence of the non-linear problem. The stabilizing properties of the proposed model are emphasized with the aid of a realistic numerical example dealing with the generation of ground fractures due to groundwater withdrawal in arid regions. - Highlights: • A numerical model is developed for the simulation of fault and fracture mechanics. • The model is implemented in the framework of the Finite Element method and with the aid of Lagrange multipliers. • The proposed formulation introduces a new contribution due to the frictional work on the portion of activated fault. • The resulting algorithm is highly non-linear as the portion of activated fault is itself unknown. • The numerical solution is validated against analytical results and proves to be stable also in realistic applications.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

NASA NTRS1985

Numerical simulations of a transverse indirect circulation and low-level jet in the exit region of an upper-level jet

A numerical study was performed of a severe weather event (tornado) which occurred on May 10, 1973 in the Ohio region. The situation was modeled with a primitive equation mesoscale dynamic formulation. Account was taken of precipitation, the planetary boundary layer parameters as bulk quantities, the vertical pressure gradient, and lateral boundary conditions based on radiosonde data. Two 12-hr simulations, adiabatic and nondivergent, respectively, were analyzed for relationships between upper and lower level jets. In the adiabatic formulation, a transverse circulation with a low level jet formed at the exit region of the upper level jet. The nondivergent situation led to similar, but weaker, phenomena. Both forms suggest that indirect circulation in the exit zone of an upper level jet is strongly influenced by the initial structure of the jet.

Brill, K. F.↗

NASA NTRS1983

Numerical simulations of counterstreaming plasmas and their relevance to interhemispheric flows

The collisionless expansion of counterstreaming plasmas has been studied in order to elucidate the basic physical processes that may be operating during the initial refilling of depleted flux tubes after a magnetic storm. The numerical technique applied is briefly described, and simulation results are presented. The simulation geometry consisted of two high-density H(+)-O(+) electron plasmas separated by a low-density H(+)-electron plasma. The temporal evolution of the expanding plasmas and the electostatic potential in the region between the two sources is described. The main interacting streams are found to be stable with respect to both the ion acoustic and ion cyclotron modes and only the suprathermal forerunner ions are unstable with respect to the ion cyclotron mode. The results also suggest that a localized potential hill can form at the equator and that this potential hill can play an important role in the subsequent trapping and thermalization of the ion streams.

Singh, N.↗

DOE OSTI2017

Numerical Simulation with Experimental Validation of the Draping Behavior of Woven Fabrics

Woven fabric composites are extensively used in molding complex geometrical shapes due to their high conformability compared to other fabrics. Preforming is an important step in the overall process. In this step, the two-dimensional fabric is draped to become the three-dimensional shape of the part prior to resin injection. During preforming, the orientation of the tows may change significantly compared to the initial orientations. Accurate prediction of the tow orientations after molding is important for evaluating the structural performance of the final part. This paper investigates the fiber angle changes for carbon fiber woven fabrics during draping over a truncated pyramid tool designed and fabricated at the General Motors Research Labs. This aspect of study is a subset of the broad study conducted under the purview of a Department of Energy project funded to GM in developing state of the art computational tools for integrated manufacturing and structural performance prediction of carbon fiber composites. Fabric bending, picture frame testing, and bias-extension evaluations were carried out to determine the material parameters for these fabrics. The PAM-FORM computer program was used to model the draping behavior of these fabrics. Following deformation, fiber angle changes at different locations on the truncated pyramid were measured experimentally. The predicted angles matched the experimental results well as measured along the centerline and at several different locations on the deformed fabric. Details of the test methods used as well as the numerical results with various simulation parameters will be provided.

42 ENGINEERING↗

NASA NTRS1984

Numerical simulation of the interaction of a vortex with stationary airfoil in transonic flow

A perturbation form of an implicit conservative, noniterative numerical algorithm for the two-dimensional thin layer Navier-Stokes and Euler equations is used to compute the interaction flow-field of a vortex with stationary airfoil. A Lamb-like analytical vortex having a finite core is chosen to interact with a thick (NACA 0012) and a thin (NACA 64A006) airfoil independently in transonic flow. Two different configurations of vortex interaction are studied, viz., (1) when the vortex is fixed at one location in the flowfield, and (2) when the vortex is convecting past the airfoil at freestream velocity. Parallel computations of this interacting flowfield are also done using a version of the Transonic Small Disturbance Code (ATRAN2). A special treatment of the leading edge region for thin airfoils is included in this code. With this, the three methods gave qualitatively similar results for the weaker interactions considered in this study. However, the strongest interactions considered proved to be beyond the capabilities of the small disturbance code. The results also show a far greater influence of the vortex on the airfoil flowfield when the vortex is stationary than when it is convecting with the flow.

Srinivasan, G. R.↗

NASA NTRS1988

Numerical simulation of mass injection for the formation of prominence magnetic field configurations. I - Asymmetric injection

A two-dimensional ideal MHD numerical model is used to investigate chromospheric mass injection into an overlying coronal dipole magnetic field. Such injection is shown to produce magnetic field deformations conducive to the formation of active region prominences. The results support a model in which an absorptive strand is formed by chromospheric asymmetric mass injection into the overlying coronal magnetic field. A necessary condition for the accumulation of the strands is that the mass injection forms a Kippenhahn-Schluter-type (1957) field configuration.

An, Chang-Hyuk↗

NASA NTRS1995

Numerical simulation of the flow about the F-18 HARV at high angle of attack

This research has been aimed at validating numerical methods for computing the flow about the complete F-18 HARV at alpha = 30 deg and alpha = 45 deg. At 30 deg angle of attack, the flow about the F-18 is dominated by the formation, and subsequent breakdown, of strong vortices over the wing leading-edge extensions (LEX). As the angle of attack is increased to alpha = 45 deg, the fuselage forebody of the F-18 contains significant laminar and transitional regions which are not present at alpha = 30 deg. Further, the flow over the LEX at alpha = 45 deg is dominated by an unsteady shedding in time, rather than strong coherent vortices. This complex physics, combined with the complex geometry of a full-aircraft configuration, provides a challenge for current computational fluid dynamics (CFD) techniques. The following sections present the numerical method and grid generation scheme that was used, a review of prior research done to numerically model the F-18 HARV, and a discussion of the current research. The current research is broken into three main topics; the effect of engine-inlet mass-flow rate on the F-18 vortex breakdown position, the results using a refined F-18 computational model to compute the flow at alpha = 30 deg and alpha = 45 deg, and research done using the simplified geometry of an ogive-cylinder configuration to investigate the physics of unsteady shear-layer shedding. The last section briefly summarizes the discussion.

Murman, Scott M.↗

DOE OSTI2016

Numerical simulations of a horizontal axis water turbine designed for underwater mooring platforms

In order to extend the operational life of Underwater Moored Platforms (UMPs), a horizontal axis water turbine is designed to supply energy for the UMPs. The turbine, equipped with controllable blades, can be opened to generate power and charge the UMPs in moored state. Three-dimensional Computational Fluid Dynamics (CFD) simulations are performed to study the characteristics of power, thrust and the wake of the turbine. Particularly, the effect of the installation position of the turbine is considered. Simulations are based on the Reynolds Averaged Navier-Stokes (RANS) equations and the shear stress transport k-ω turbulent model is utilized. The numerical method is validated using existing experimental data. The simulation results show that this turbine has a maximum power coefficient of 0.327 when the turbine is installed near the tail of the UMP. The flow structure near the blade and in the wake are also discussed.

13 HYDRO ENERGY↗

DOE OSTI2014

Numerical simulation of shock initiation of Ni/Al multilayered composites

The initiation of chemical reaction in cold-rolled Ni/Al multilayered composites by shock compression is investigated numerically. A simplified approach is adopted that exploits the disparity between the reaction and shock loading timescales. The impact of shock compression is modeled using CTH simulations that yield pressure, strain, and temperature distributions within the composites due to the shock propagation. The resulting temperature distribution is then used as initial condition to simulate the evolution of the subsequent shock-induced mixing and chemical reaction. To this end, a reduced reaction model is used that expresses the local atomic mixing and heat release rates in terms of an evolution equation for a dimensionless time scale reflecting the age of the mixed layer. The computations are used to assess the effect of bilayer thickness on the reaction, as well as the impact of shock velocity and orientation with respect to the layering. Computed results indicate that initiation and evolution of the reaction are substantially affected by both the shock velocity and the bilayer thickness. In particular, at low impact velocity, Ni/Al multilayered composites with thick bilayers react completely in 100 ms while at high impact velocity and thin bilayers, reaction time was less than 100 μs. Quantitative trends for the dependence of the reaction time on the shock velocity are also determined, for different bilayer thickness and shock orientation.

36 MATERIALS SCIENCE↗

NASA NTRSDate not supplied

Numerical Simulation of No-Vent Chill and Fill of a Large Liquid Hydrogen Tank

To enable the design of future in-space cryogenic propellant transfer vehicles, high accuracy models of the cryogenic propellant transfer process are desired. For in-space refueling applications, accurate simulations are desired of the propellant tank chilldown and no-vent fill (NVF) process in order to estimate the amount of propellant required, to determine the maximum sustainable flow rate to fill the tank, and to establish a timeline for the filling process. NVF is a complex process that involves numerous two-phase heat and mass transfer phenomena, such as boiling at the wall and evaporation and condensation at the liquid/vapor interface. In recent years, progress has been made in developing a numerical modeling technique using the Generalized Fluid System Simulation Program (GFSSP) that has been validated against several small-scale liquid nitrogen NVF tests. This paper presents GFSSP numerical model validation against three liquid hydrogen NVF tests on a large-scale propellant tank. Comparing the GFSSP model to the data for pressure, fill level, and wall temperature, the mean absolute percentage error is below 20% across the range of simulated conditions, indicating very good predictive agreement for this large-scale hydrogen NVF test series.

Alok Majumdar↗

NASA NTRSDate not supplied

Numerical Simulation of No-Vent Chill and Fill of a Large Liquid Hydrogen Tank

Alok Majumdar↗

NASA NTRS1992

Direct numerical simulation of a combusting droplet with convection

The evaporation and combustion of a single droplet under forced and natural convection was studied numerically from first principles using a numerical scheme that solves the time-dependent multiphase and multispecies Navier-Stokes equations and tracks the sharp gas-liquid interface cutting across an arbitrary Eulerian grid. The flow fields both inside and outside of the droplet are resolved in a unified fashion. Additional governing equations model the interphase mass, energy, and momentum exchange. Test cases involving iso-octane, n-hexane, and n-propanol droplets show reasonable comparison rate, and flame stand-off distance. The partially validated code is, thus, readied to be applied to more demanding droplet combustion situations where substantial drop deformation render classical models inadequate.

Liang, Pak-Yan↗