Search NASA⌕ Search

SEARCH · Search NASA

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 487 records · Page 27

Numerical simulation of X-ray-heated winds in Seyfert galaxies. I - The case of zero angular momentum

Numerical hydrodynamics simulations are presented of X-ray heated winds believed to form the 'mirrors' which reflect light from obscured Seyfert galaxy nuclei into our line of sight. Above a likely range of parameters, it is shown that steady state flows are created when a cool torus surrounds an AGN. Some of the mass evaporated off the inner edge of the torus is captured by the central gravitational well, while the rest is driven off in a warm wind. When the central gravity is strong relative to the radiative heating rate, i.e., L/L(E) is small, nearly all the injected matter is captured. When the gravity is weaker, a given combination of heating rate and injected pressure leads to a maximum mass loss rate; when the rate at which mass is injected exceeds this amount, the excess is captured. The dividing line between these two cases comes at L/L(E) = about 0.08.

Balsara, Dinshaw S.↗

Improvements of an ice-phase microphysics parameterization for use in numerical simulations of tropical convection

It is important to properly simulate the extent and ice water content of tropical anvil clouds in numerical models that explicitly include cloud formation because of the significant effects that these clouds have on the radiation budget. For this reason, a commonly used bulk ice-phase microphysics parameterization was modified to more realistically simulate some of the microphysical processes that occur in tropical anvil clouds. Cloud ice growth by the Bergeron process and the associated formation of snow were revised. The characteristics of graupel were also modified in accord with a previous study. Numerical simulations of a tropical squall line demonstrate that the amount of cloud ice and the extent of anvil clouds are increased to more realistic values by the first two changes.

Krueger, Steven K.↗

Numerical simulation of homogeneous, isotropic, compressible turbulence

A new numerical algorithm is developed and applied to the direct simulation of compressible, homogeneous turbulent flows at low Mach numbers. A split method in time first solves a subset of the equations explicitly, followed by an implicit treatment of the pressure terms. The definition of an average sound speed over the entire field allows the implicit equations to be solved analytically, while relaxing the severe time limit imposed by the large sound speeds. Results from direct simulations on 96(3) grids provide a data base against which a new subgrid-scale model for compressible homogeneous turbulence is tested. This model reduces to the linear combination model in the absence of compressibility.

Erlebacher, G.↗

Numerical simulations of unsteady airfoil-vortex interactions

Numerical methods based on the thin layer Navier-Stokes, Euler and transonic small disturbance equations, all three of which use the perturbation form of an implicit numerical algorithm, are used to formulate the unsteady interactions of a moving vortex with a stationary airfoil. Results are presented for both transonic and subsonic flows, as well as for weak and strong interactions. In general, the results exhibit an overwhelming influence of the vortex on the flow field around the airfoil, especially when the convecting vortex is within one chord of the airfoil. It is noted that even the ostensibly weak interactions produced vortex/shock-induced boundary layer separation.

Srinivasan, G. R.↗

First-Principles Modeling of IAT-Driven Anomalous Resistivity in Hollow Cathode Discharges II: Numerical Simulations and Comparisons with Experiments

We present a model that quantifies the magnitude of the ion-acoustic turbulence (IAT) in the plume of hollow cathodes and its effect on the resistivity and ion heating. The model takes the form of a partial differential equation (PDE) that can be solved concurrently with the equations of motion for a partially ionized plasma already included in our numerical code for the simulation of the plasma discharge in hollow cathodes, OrCa2D. We also determine that self-induced magnetic fields are not negligible in hollow cathodes operating at large discharge currents and implement in our code Ampere’s law and modifications to Ohm’s law that account for this effect. Numerical simulations that employed these models show large improvements in our agreement with experimental measurements with respect to a previous model, which assumed complete saturation of the IAT and did not account for the growth stage of the waves. In particular, the model is able to accurately predict the location and magnitude of the maximum resitivity to the electron current along the cathode centerline.

Lopez Ortega, Alejandro↗

Numerical simulation of electrothermal de-icing systems

Transient simulations of de-icing of composite aircraft components by electrothermal heating have been computed for both one and two-dimensional rectangular geometries. The implicit Crank-Nicolson formulation is used to insure stability of the finite-differenced heat conduction equations and the phase change in the ice layer is simulated using the Enthalpy method. Numerical solutions illustrating de-icer performance for various composite aircraft blades and environmental conditions are presented. Comparisons are made with previous studies and with available experimental data. Initial results using a coordinate mapping technique to describe the actual blade geometry are discussed.

De Witt, K. J.↗

Numerical simulations of the bending of narrow-angle-tail radio jets by ram pressure or pressure gradients

Three-dimensional numerical hydrodynamic simulations are used to study the bending of radio jets. The simulations are compared with observations of jets in narrow-angle-tail radio sources. Two mechanisms for the observed bending are considered: direct bending of quasi-continuous jets by ram pressure from intergalactic gas and bending by pressure gradients in the interstellar gas of the host galaxy, the pressure gradients themselves being the result of ram pressure by intergalactic gas. It is shown that the pressure gradients are much less effective in bending jets, implying that the jets have roughly 30 times lower momentum fluxes if they are bent by this mechanism. Ram-pressure bending produces jets with 'kidney-shaped' cross sections; when observed from the side, these jets appear to have diffuse extensions on the downstream side. On the other hand, pressure-gradient bending causes the jets to be densest near their upstream side.

Soker, Noam↗

Numerically Simulating an Expanding Continuum Jet into a Surrounding Non-Continuum Region

A numerical technique that simulates the interaction between a continuum jet and a free-molecular, rarefied or transitional flow field is presented and applied. The technique is based on selecting a boundary between the expanding jet plume and interacting flow field defined by the Bird breakdown parameter, P. By properly choosing the boundary, the jet plume computed by computational fluid dynamics (CFD) is uncoupled from the interaction region that is simulated by direct simulation Monte Carlo (DSMC). The mechanics of uncoupling the continuum and non-continuum regions is discussed first. Also, a discussion of expanding continuum plume breakdown and the formulation of the Bird breakdown parameter is presented. To show the usefulness, examples of applying the CFD-DSMC technique are given. Many of the examples were taken from projects worked by the staff of the Aerothermodynamics Branch at the NASA Langley Research Center. These examples show a variety of applications of the technique and allow those examples not previously documented to be formally presented.

Jet interaction↗

Consistent large-eddy simulation of a temporal mixing layer laden with evaporating drops : Part 1. Direct numerical simulation formulation and a priori analysis

This paper is organized as follows. In § 2, we present the DNS formulation. This formulation consists of a set of coupled conservation equations for the gas phase and the liquid phase (i.e. the drops). The coupling describing the interaction between the two phases is provided in the gas-phase conservation equations by source terms. In § 3, we analyse the DNS results; in particular, we evaluate the adopted form of the heat and mass fluxes. The gas-phase LES equations are developed in § 4 and assumptions that may simplify these equations are examined using the DNS database. Thereafter, the contribution of the source terms is examined by considering the budget of the LES equations. The final form of the LES equations are then presented, in which only the SGS fluxes and the FSTs need be modelled. A parallel perspective on the importance of the source terms is presented in § 5, wherein the irreversible entropy production (i.e. the dissipation) expression for two-phase flow with phase change is derived and evaluated for transitional DNS and the filtered flow fields. Models for the FSTs appearing in the LES equations are investigated in § 6, while several constantcoefficient SGS models are proposed and calibrated in § 7. The present study aims to create simple models so that the modelling approach can be validated by means of an SGS models are proposed and calibrated in § 7. The present study aims to create simple models so that the modelling approach can be validated by means of an a posteriori study.

"Okong'o, Nora A.", 'Bellan, Josette↗

Numerical simulation of reconnection in an emerging magnetic flux region

Numerical solutions in two dimensions are presented for the resistive MHD equations of an initial boundary value problem, simulating reconnection between an emerging magnetic flux region and an overlying coronal magnetic field. This numerical model displays four main phases, which are interpreted in terms of (1) a slowly evolving quasi-steady phase during which most of the magnetic flux emerges, with approximate equilibrium between magnetic and pressure forces; (2) an impulsive phase in which either the reconnection or continuing emergence of the first phase comes to disrupt the equilibrium, and extensive acceleration occurs as the high pressure region within the emerging region drives the fluid upwards and outwards; (3) a second quasi-steady phase; and (4) a potential-static phase in which continuing reconnection, ohmic dissipation, and fluid transport through boundaries depletes the system of all currents and flows.

Forbes, T. G.↗

Direct numerical simulations of a spatially developing plane wake

In the present paper, direct numerical methods by which to simulate the spatially developing free shear flows in the transitional region are described and the numerical results of a spatially developing plane wake are presented. The incompressible time-dependent Navier-Stokes equations were solved using Pade finite difference approximations in the streamwise direction, a mapped pseudospectral Fourier method in the cross-stream direction, and a third-order compact Runge-Kutta scheme for time advancement. The unstable modes of the Orr-Sommerfeld equations were used to perturb the inlet of the wake. Statistical analyses were performed and some numerical results were compared with experimental measurements. When only the fundamental mode is forced, the energy spectra show amplification of the fundamental and its higher harmonics. In this case, unperturbed alternate vortices develop in the saturation region of the wake. The phase jitter around the fundamental frequency plays a critical role in generating vortices of random shape and spacing. Large- and small-scale distortions of the fundamental structure are observed. Pairing of vortices of the same sign is observed, as well as vortex coupling of vortices of the opposite sign.

Maekawa, Hiroshi↗

Numerical simulation of turbulent flows

Computational models of turbulence in incompressible Newtonian fluids governed by the Navier-Stokes equations are reviewed. The governing equations are presented, and both direct and large-eddy-simulations are examined. Resolution requirements and numerical techniques of spatial representation, definition of initial and boundary conditions, and time advancement are considered. Results of simulations of homogeneous turbulence in uniform shear, the evolution of a turbulent mixing layer, and turbulent channel flow are presented graphically and discussed.

Rogallo, R. S.↗

Numerical simulation of boundary layer excitation by surface heating/cooling

The concept of active control of growing disturbances in an unstable compressible flow by using time periodic, localized surface heating is studied numerically. The simulations are calculated by a fourth-order accurate solution of the compressible, laminar Navier-Stokes equations. Fourth-order accuracy is particularly important for this problem because the solution must be computed over many wavelengths. The numerical results demonstrate the growth of an initially small fluctuation into the nonlinear regime where a local breakdown into smaller scale disturbances can be observed. It is shown that periodic surface heating over a small strip can reduce the level of the fluctuation provided that the phase of the heating current is properly chosen.

Bayliss, A.↗

Numerical simulation of transonic flow over porous airfoils

A numerical study was made to examine the effect of a porous surface on the aerodynamic performance of a transonic airfoil. The pressure jump across the normal shock wave on the upper surface of the airfoil was reduced by making the surface below the shock porous. The weakened shock is preceded by an oblique shock at the upstream end of the porous surface where air is blown out of the cavity. The lambda shock structure shown in the numerical result qualitatively agrees with that observed in the wind tunnel. According to the present analysis, the porous airfoil has a smaller drag and a higher lift than the solid airfoil.

Chow, C.-Y.↗

Higher-Order Compact Schemes for Numerical Simulation of Incompressible Flows

A higher order accurate numerical procedure has been developed for solving incompressible Navier-Stokes equations for 2D or 3D fluid flow problems. It is based on low-storage Runge-Kutta schemes for temporal discretization and fourth and sixth order compact finite-difference schemes for spatial discretization. The particular difficulty of satisfying the divergence-free velocity field required in incompressible fluid flow is resolved by solving a Poisson equation for pressure. It is demonstrated that for consistent global accuracy, it is necessary to employ the same order of accuracy in the discretization of the Poisson equation. Special care is also required to achieve the formal temporal accuracy of the Runge-Kutta schemes. The accuracy of the present procedure is demonstrated by application to several pertinent benchmark problems.

Wilson, Robert V.↗

Implicit and Multigrid Method for Ideal Multigrid Convergence: Direct Numerical Simulation of Separated Flow Around NACA 0012 Airfoil

Numerical investigation of flow separation over a NACA 0012 airfoil at large angles of attack has been carried out. The numerical calculation is performed by solving the full Navier-Stokes equations in generalized curvilinear coordinates. The second-order LU-SGS implicit scheme is applied for time integration. This scheme requires no tridiagonal inversion and is capable of being completely vectorized, provided the corresponding Jacobian matrices are properly selected. A fourth-order centered compact scheme is used for spatial derivatives. In order to reduce numerical oscillation, a sixth-order implicit filter is employed. Non-reflecting boundary conditions are imposed at the far-field and outlet boundaries to avoid possible non-physical wave reflection. Complex flow separation and vortex shedding phenomenon have been observed and discussed.

Liu, Chao-Qun↗

Numerical simulation of boundary-layer excitation by surface heating/cooling

The concept of active control of growing disturbances in an unstable compressible flow by using time periodic, localized surface heating is studied numerically. The simulations are calculated by a fourth-order accurate solution of the compressible, laminar Navier-Stokes equations. Fourth-order accuracy is particularly important for this problem because the solution must be computed over many wavelengths. The numerical results demonstrate the growth of an initially small fluctuation into the nonlinear regime where a local breakdown into smaller scale disturbances can be observed. It is shown that periodic surface heating over a small strip can reduce the level of the fluctuation provided that the phase of the heating current is properly chosen.

Bayliss, A.↗