Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical modelling”

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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

A model of concurrent flow flame spread over a thin solid fuel

A numerical model is developed to examine laminar flame spread and extinction over a thin solid fuel in lowspeed concurrent flows. The model provides a more precise fluid-mechanical description of the flame by incorporating an elliptic treatment of the upstream flame stabilization zone near the fuel burnout point. Parabolic equations are used to treat the downstream flame, which has a higher flow Reynolds number. The parabolic and elliptic regions are coupled smoothly by an appropriate matching of boundary conditions. The solid phase consists of an energy equation with surface radiative loss and a surface pyrolysis relation. Steady spread with constant flame and pyrolysis lengths is found possible for thin fuels and this facilitates the adoption of a moving coordinate system attached to the flame with the flame spread rate being an eigen value. Calculations are performed in purely forced flow in a range of velocities which are lower than those induced in a normal gravity buoyant environment. Both quenching and blowoff extinction are observed. The results show that as flow velocity or oxygen percentage is reduced, the flame spread rate, the pyrolysis length, and the flame length all decrease, as expected. The flame standoff distance from the solid and the reaction zone thickness, however, first increase with decreasing flow velocity, but eventually decrease very near the quenching extinction limit. The short, diffuse flames observed at low flow velocities and oxygen levels are consistent with available experimental data. The maximum flame temperature decreases slowly at first as flow velocity is reduced, then falls more steeply close to the quenching extinction limit. Low velocity quenching occurs as a result of heat loss. At low velocities, surface radiative loss becomes a significant fraction of the total combustion heat release. In addition, the shorter flame length causes an increase in the fraction of conduction downstream compared to conduction to the fuel. These heat losses lead to lower flame temperatures, and ultimately, extinction. This extinction mechanism differs from that of blowoff, where the flame is unable to be stabilized due to the high flow velocity.

Ferkul, Paul V.↗
NASA NTRS1993

The Stickney impact of Phobos - A dynamical model

The hypervelocity impact that excavated the Stickney crater on Phobos is numerically modeled in order to understand the effects this collision had on the interior and surface of Phobos. The model assumes homogeneity prior to impact based on Fujiwara's (1991) results. Constitutive relations for water ice and basalt are applied in tandem models. Fracture energetics is used to show that the impact itself is sufficient to cause the cracks observed on Phobos. It is further shown that the Stickney impact itself could not have created large voids; thus, the low density of Phobos must either be compositional in nature or else the result of porosity at a scale sufficiently small to avoid scattering the impact energy substantially.

Asphaug, E.↗
NASA NTRS2000

Performance Modeling of an Airborne Raman Water Vapor Lidar

A sophisticated Raman lidar numerical model had been developed. The model has been used to simulate the performance of two ground-based Raman water vapor lidar systems. After tuning the model using these ground-based measurements, the model is used to simulate the water vapor measurement capability of an airborne Raman lidar under both day-and night-time conditions for a wide range of water vapor conditions. The results indicate that, under many circumstances, the daytime measurements possess comparable resolution to an existing airborne differential absorption water vapor lidar while the nighttime measurement have higher resolution. In addition, a Raman lidar is capable of measurements not possible using a differential absorption system.

Whiteman, D. N.↗
NASA NTRS1976

Wind laws for shockless initialization

A system of diagnostic equations for the velocity field, or wind laws, was derived for each of a number of models of large-scale atmospheric flow. The derivation in each case is mathematically exact and does not involve any physical assumptions not already present in the prognostic equations, such as nondivergence or vanishing of derivatives of the divergence. Therefore, initial states computed by solving these diagnostic equations should be compatible with the type of motion described by the prognostic equations of the model and should not generate initialization shocks when inserted into the model. Numerical solutions of the diagnostic system corresponding to a barotropic model are exhibited. Some problems concerning the possibility of implementing such a system in operational numerical weather prediction are discussed.

Ghil, M.↗
NASA NTRS1986

Convective effects in float-zone and Czochralski melts

The hydrodynamics of crystal-growth melts is a relatively new research area. Numerical modeling of these processes is necessary. The work discussed herein is in two parts: numerical simulations of the flow in a Czochralski melt, and also of that in a float zone. In addition, for the float-zone case, energy stability theory will be used to determine stability bounds for the onset of oscillatory thermo-capillary flow. Convective effects in crystal-growth melts arise from a variety of mechanisms. Temperature gradients both in the direction of gravity and normal to it give rise to convection due to buoyancy effects. Rotation of the crucible and/or crystal causes a forced convection which may augment or oppose the buoyancy-driven flow. Finally, thermo-capillary forces (due to the variation of surface tension with temperature) drive surface motions which in turn generate convection in the bulk fluid. All of these mechanisms are present in either Czochralski or float-zone growth. The objective of the Czochralski modeling is to develop an accurate numerical simulation of the flow in a Czochralski silicon melt and to investigate the effects of various parameters on the flow properties. Like some earlier investigations, the intent is to simulate the effects of buoyancy, forced and thermo-capillary convection, including unsteady effects. Unlike earlier work, the aim is to include the effects of a variable free surface and freezing interface and, possibly incorporate nonaxisymmetric effects.

Neitzel, G. P.↗
NASA NTRS1983

Magnetohydrodynamic modeling of coronal structure and expansion

The presence of a magnetic field in the corona adds structure to the solar wind and almost certainly plays an important role in the energetics of the flow. Analytical and numerical modeling of gas-magnetic field interactions as used to compute steady, global flow are discussed. The approach used in, and results from a recent global model (Steinolfson, Suess and Wu, 1982) are discussed. Ideas on the most effective ways to improve the physical content and numerical efficiency of these models are outlined. Solutions of the MHD equations are discussed only in order to find steady-state flows, even though this often entails solving time-dependent equations.

Suess, S. T.↗
NASA NTRS1982

Initial results from a mesoscale atmospheric simulation system and comparisons with the AVE-SESAME I data set

The development of a comprehensive mesoscale atmospheric simulation system (MASS) is described in detail. The modeling system is designed for both research and real-time forecast applications. The 14-level numerical model, which has a 48 km grid mesh, can be run over most of North America and the adjacent oceanic regions. The model employs sixth-order accurate numerics, generalized similarity theory boundary-layer physics, a sophisticated cumulus parameterization scheme, and state of the art analysis and initialization techniques. Examples of model output on the synoptic and subsynoptic scales are presented for the AVE-SESAME I field experiment on 10-11 April 1979. The model output is subjectively compared to the observational analysis and the LFM II output on the synoptic scale. Subsynoptic model output is compared to analyses generated from the AVE-SESAME I data set.

Kaplan, M. L.↗
NASA NTRS2004

Development of a Computationally Efficient, High Fidelity, Finite Element Based Hall Thruster Model

This report documents the development of a two dimensional finite element based numerical model for efficient characterization of the Hall thruster plasma dynamics in the framework of multi-fluid model. Effect of the ionization and the recombination has been included in the present model. Based on the experimental data, a third order polynomial in electron temperature is used to calculate the ionization rate. The neutral dynamics is included only through the neutral continuity equation in the presence of a uniform neutral flow. The electrons are modeled as magnetized and hot, whereas ions are assumed magnetized and cold. The dynamics of Hall thruster is also investigated in the presence of plasma-wall interaction. The plasma-wall interaction is a function of wall potential, which in turn is determined by the secondary electron emission and sputtering yield. The effect of secondary electron emission and sputter yield has been considered simultaneously, Simulation results are interpreted in the light of experimental observations and available numerical solutions in the literature.

Jacobson, David↗
NASA NTRS1993

Antenna pattern control using impedance surfaces

During the period of this research project, a comprehensive study of pyramidal horn antennas was conducted. Full-wave analytical and numerical techniques were developed to analyze horn antennas with or without impedance surfaces. Based on these full-wave analytic techniques, research was conducted on the use of impedance surfaces on the walls of the horn antennas to control the antenna radiation patterns without a substantial loss of antenna gain. It was found that the use of impedance surfaces could modify the antenna radiation patterns. In addition to the analytical and numerical models, experimental models were also constructed and they were used to validate the predictions. Excellent agreement between theoretical predictions and the measured data was obtained for pyramidal horns with perfectly conducting surfaces. Very good comparisons between numerical and experimental models were also obtained for horns with impedance surfaces.

Balanis, Constantine A.↗
NASA NTRS2022

VULCAN-CFD Theory Manual: Ver. 7.2.0

VULCAN-CFD offers a comprehensive set of capabilities to enable the simulation of continuum flowfields from subsonic to hypersonic conditions. The governing equations that are employed include allowances for both chemical and thermal nonequilibrium processes, coupled with a wide variety of turbulence models for both Reynolds-averaged and large eddy simulations. A description of the physical and numerical models available in the software are presented in this document. However, it is emphasized that the descriptions provided are not intended to fully document every aspect of the models employed. Instead, the governing equations, and the models required to simulate them numerically, are presented with the goal of providing a sufficient level of detail to understand their strengths and limitations. The reader is encouraged to access the references provided throughout this document for a more complete explanation of the formulations presented.

VULCAN-CFD Theory Manual↗
NASA NTRS1986

Circuit transients due to negative bias arcs-II

Two new models of negative-bias arcing on a solar cell power system in Low Earth Orbit are presented. One is an extended, analytical model and the other is a non-linear, numerical model. The models are based on an earlier analytical model in which the interactions between solar cell interconnects and the space plasma as well as the parameters of the power circuit are approximated linearly. Transient voltages due to arcs struck at the negative thermal of the solar panel are calculated in the time domain. The new models treat, respectively, further linear effects within the solar panel load circuit and non-linear effects associated with the plasma interactions. Results of computer calculations with the models show common-mode voltage transients of the electrically floating solar panel struck by an arc comparable to the early model but load transients that differ substantially from the early model. In particular, load transients of the non-linear model can be more than twice as great as those of the early model and more than twenty times as great as the extended, linear model.

Metz, R. N.↗
NASA NTRS1988

Toward improved durability in advanced aircraft engine hot sections; Proceedings of the Thirty-third ASME International Gas Turbine and Aeroengine Congress and Exposition, Amsterdam, Netherlands, June 5-9, 1988

The present conference on durability improvement methods for advanced aircraft gas turbine hot-section components discusses NASA's 'HOST' project, advanced high-temperature instrumentation for hot-section research, the development and application of combustor aerothermal models, and the evaluation of a data base and numerical model for turbine heat transfer. Also discussed are structural analysis methods for gas turbine hot section components, fatigue life-prediction modeling for turbine hot section materials, and the service life modeling of thermal barrier coatings for aircraft gas turbine engines.

Sokolowski, Daniel E.↗
NASA NTRS1988

Diabatic initialization for a mesoscale model simulation experiment

A mesoscale model simulation of a 'surrogate' atmosphere designed for the evaluation of adiabatic initialization and moisture enhancement techniques is briefly described. The simulation provides mesoscale precipitation rates and associated diabatic heating fields which are known exactly and are internally consistent within the context of the model. The advantage of this approach is that a numerical model can easily serve as a source of mesoscale data in the usual absence of internally consistent mesoscale observed data sets.

Kalb, Michael W.↗
NASA NTRS1989

Nonlocal thermal transport in solar flares. II - Spectroscopic diagnostics

Physical parameters obtained for a flaring solar atmosphere in an earlier paper are used here to predict time-dependent emission-line profiles and integrated intensities as a function of position for two spectral lines commonly observed during solar flares: the X-ray resonance lines of Ca XIX and Mg XI. Considerations of ionization nonequilibrium during the rise phase of the flare are addressed, and the effects on the predicted spectral-line characteristics are discussed. It is concluded that some spectroscopic diagnostics favor the nonlocal model, but other long-standing discrepancies between the numerical models and the observations remain unresolved.

Karpen, Judith T.↗
NASA NTRS1989

Intensive MHD-structures penetration in the middle atmosphere initiated in the ionospheric cusp under quiet geomagnetic conditions

In connection with the recently detected quasiperiodical magnetic disturbances in the ionospheric cusp, the penetration of compressional surface magnetohydrodynamic (MHD) waves through the middle atmosphere is modelled numerically. For the COSPAR International Reference Atmosphere (CIRA) 72 model the respective energy density flux of the disturbances in the middle atmosphere is determined. On the basis of the developed model certain conclusions are reached about the height distribution of the structures (energy losses, currents, etc.) initiated by intensive magnetic cusp disturbances.

Mateev, L. N.↗
NASA NTRS1991

Numerical algorithm comparison for the accurate and efficient computation of high-incidence vortical flow

Computations from two Navier-Stokes codes, NSS and F3D, are presented for a tangent-ogive-cylinder body at high angle of attack. Features of this steady flow include a pair of primary vortices on the leeward side of the body as well as secondary vortices. The topological and physical plausibility of this vortical structure is discussed. The accuracy of these codes are assessed by comparison of the numerical solutions with experimental data. The effects of turbulence model, numerical dissipation, and grid refinement are presented. The overall efficiency of these codes are also assessed by examining their convergence rates, computational time per time step, and maximum allowable time step for time-accurate computations. Overall, the numerical results from both codes compared equally well with experimental data, however, the NSS code was found to be significantly more efficient than the F3D code.

Chaderjian, Neal M.↗
NASA NTRS1992

Approaches to ionospheric modelling, simulation and prediction

The ionosphere is a complex, multispecies, anisotropic medium that exhibits a significant variation with time, space, season, solar cycle, and geomagnetic activity. In recent years, a wide range of models have been developed in an effort to describe ionospheric behavior. The modeling efforts include: (1) empirical models based on extensive worldwide data sets; (2) simple analytical models for a restricted number of ionospheric parameters; (3) comprehensive, 3D, time-dependent models that require supercomputers; (4) spherical harmonic models based on fits to output obtained from comprehensive numerical models; and (5) ionospheric models driven by real-time magnetospheric inputs. In an effort to achieve simplicity, some of the models have been restricted to certain altitude or latitude domains, while others have been restricted to certain ionospheric parameters, such as the F-region peak density, the auroral conductivity, and the plasma temperatures. The current status of the modeling efforts is reviewed.

Schunk, R. W.↗
NASA NTRS1992

A modeling study of the time-averaged electric currents in the vicinity of isolated thunderstorms

A thorough examination of the results of a time-dependent computer model of a dipole thunderstorm revealed that there are numerous similarities between the time-averaged electrical properties and the steady-state properties of an active thunderstorm. Thus, the electrical behavior of the atmosphere in the vicinity of a thunderstorm can be determined with a formulation similar to what was first described by Holzer and Saxon (1952). From the Maxwell continuity equation of electric current, a simple analytical equation was derived that expresses a thunderstorm's average current contribution to the global electric circuit in terms of the generator current within the thundercloud, the intracloud lightning current, the cloud-to-ground lightning current, the altitudes of the charge centers, and the conductivity profile of the atmosphere. This equation was found to be nearly as accurate as the more computationally expensive numerical model, even when it is applied to a thunderstorm with a reduced conductivity thundercloud, a time-varying generator current, a varying flash rate, and a changing lightning mix.

Driscoll, Kevin T.↗