Numerical simulation of fog formation and liquid water content on polydisperse multi-composition aerosols due to combustion-related pollutants
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
The characteristics of unsteady flow motions have long been a serious concern in the study of various fluid dynamic and combustion problems. With the advancement of computer resources, numerical approaches to these problems appear to be feasible. The objective of this paper is to assess the accuracy of several numerical schemes for unsteady flow calculations. In the present study, Fourier error analysis is performed for various numerical schemes based on a two-dimensional wave equation. Four methods sieved from the error analysis are then adopted for further assessment. Model problems include unsteady quasi-one-dimensional inviscid flows, two-dimensional wave propagations, and unsteady two-dimensional inviscid flows. According to the comparison between numerical and exact solutions, although second-order upwind scheme captures the unsteady flow and wave motions quite well, it is relatively more dissipative than sixth-order central difference scheme. Among various numerical approaches tested in this paper, the best performed one is Runge-Kutta method for time integration and six-order central difference for spatial discretization.
The GCM-Reality Intercomparison Project for SPARC (GRIPS) is assessing and monitoring the performance of state-of-the-art general circulation models (GCMs). A wide variety of tasks have been initiated. These are designed to: (1) assess the ability of the GCMs to represent the current climatological structure of the troposphere and middle atmosphere,(2) to compare their response to imposed forcing anomalies, and (3) to estimate the certainty with which future climate perturbations can be predicted. This paper is concerned with assessments of the climatological states in the GCM simulations. Comparing the simulations with observational datasets reveals considerable discrepancies in the modelled fields. While it might be anticipated that certain types of biases in the model simulations might be related to the formulation of different aspects of the numerical package (dynamical schemes, cloud schemes, radiation transfer, inclusion of gravity wave drag), there is no clear relationship between these features. This paper attempts to draw a more comprehensive picture of the GCMs'performance than has previously been shown, by comparing the dominant forcing mechanisms in the models with observational estimates, and relating model deficiencies to the differences in the physical mechanisms in the GCMS.
The computation of unsteady shock waves, which contribute significantly to noise generation in supersonic jet flows, is investigated. The paper focuses on the difficulties of computing slowly moving shock waves. Numerical error is found to manifest itself principally as a spurious entropy wave. Calculations presented are performed using a third-order essentially nonoscillatory scheme. The effect of stencil biasing parameters and of two versions of numerical flux formulas on the magnitude of spurious entropy are investigated. The level of numerical error introduced in the calculation is quantified as a function of shock pressure ratio, shock speed, Courant number, and mesh density. The spurious entropy relative to the entropy jump across a static shock decreases with increasing shock strength and shock velocity relative to the grid, but is insensitive to Courant number. The structure of the spurious entropy wave is affected by the choice of flux formulas and algorithm biasing parameters. The effect of the spurious numerical waves on the calculation of sound amplification by a shock wave is investigated. For this class of problem, the acoustic pressure waves are relatively unaffected by the spurious numerical phenomena.
The computation of unsteady shock waves, which contribute significantly to noise generation in supersonic jet flows, is investigated. This paper focuses on the difficulties of computing slowly moving shock waves. Numerical error is found to manifest itself principally as a spurious entropy wave. Calculations presented are performed using a third order essentially nonoscillatory scheme. The effect of stencil biasing parameters and of two versions of numerical flux formulas on the magnitude of spurious entropy are investigated. The level of numerical error introduced in the calculation in quantified as a function of shock pressure ratio, shock speed, Courant number, and mesh density. The spurious entropy relative to the entropy jump across a static shock decreases with increasing shock strength and shock velocity relative to the grid, but is insensitive to Courant number. The structure of the spurious entropy wave is affected by the choice of flux formulas and algorithm biasing parameters. The effect of the spurious numerical waves on the calculation of sound amplification by a shock wave is investigated. For this class of problem, the acoustic pressure waves are relatively unaffected by the spurious numerical phenomena.
In several papers, the solar cycle (SC) effect in the lower atmosphere has been linked observationally to the Quasi-biennial Oscillation (QBO) of the zonal circulation, which is generated primarily by small-scale gravity waves (GW). Salby and Callaghan (2000) in particular analyzed the QBO, covering more than 40 years, and discovered that it contains a large SC signature at 20 km. With our Numerical Spectral Model (NSM), we conducted a 3D study to describe the QBO under the influence of the SC, and some results have been published (Mayr et al., GRL, 2005,2006). For a SC period of 10 years, the relative amplitude of radiative forcing is taken to vary exponentially with height, i.e., 0.2% at the surface, 2% at 50 km, 20% at 100 km and above. Applying spectral analysis to filter out and identify the SC signature, the model generates a relatively large modulation of the QBO, which reproduces the observations qualitatively. Our numerical results demonstrate that the modulation of the QBO, with constant phase relative to the SC, persist at least for 60 years. The same model run generates in the seasonal variations a hemispherically symmetric Equatorial Annual Oscillation (EAO, with 12-month period), which is confined to low latitudes like the QBO and is also modulated by the SC. Although the amplitude of the EAO is relatively small, its SC modulation is large, and it is in phase with that of the QBO. The SC modulated EAO is evidently the pathway and pacemaker for the solar influence on the QBO. To shed light on the dynamical processes involved, we present model results that show how the seasonal cycle induces the SC modulations of the EAO and QBO. Our analysis further demonstrates that the SC modulations of the QBO and EAO are amplified by the GW interaction with the flow. The GW momentum source clearly shows a SC modulation that is in phase with the corresponding modulations of the QBO and EAO. By tapping the momentum from the upward propagating GWs, the QBO and EAO apparently serve as conduits to amplify and transfer to lower altitudes the larger SC variations in the UV absorbed in the mesosphere. Our model also produces in the temperature variations of the QBO and EAO measurable SC modulations at polar latitudes near the tropopause, and such signatures have been reported in the literature. Contrary to conventional interpretation, however, we suggest that the effects are generated at least in part by the meridional circulation, and planetary waves presumably, which redistribute the energy from the equatorial region where wave forcing is very efficient and thereby amplifies the SC influence.
The relative motion problem is analyzed, as a linearized case, and as a numerically determined solution to provide a time history of the geometries representing the motion state. The displacement history and the hodographs for families of solutions are provided, analytically and graphically, to serve as an aid to understanding this problem area. Linearized solutions to relative motion problems of orbiting particles are presented for the impulsive and fixed thrust cases. Second order solutions are described to enhance the accuracy of prediction. A method was developed to obtain accurate, numerical solutions to the intercept and rendezvous problem; and, special situations are examined. A particular problem related to relative motions, where the motion traces develop a cusp, is examined in detail. This phenomenon is found to be dependent on a particular relationship between orbital eccentricity and the inclination between orbital planes. These conditions are determined, and, example situations are presented and discussed.
An analysis of linear wave modes associated with supersonic jets confined inside an acoustically lined rectangular duct is presented. Mathematical formulations are given for the vortex-sheet model and continuous mean flow model of the jet flow profiles. Detailed dispersion relations of these waves in a two-dimensional confined jet as well as an unconfined free jet are computed. Effects of the confining duct and the liners on the jet instability and acoustic waves are studied numerically. It is found that the effect of the liners is to attenuate waves that have supersonic phase velocities relative to the ambient flow. Numerical results also show that the growth rates of the instability waves could be reduced significantly by the use of liners. In addition, it is found that the upstream propagating neutral waves of an unconfined jet could become attenuated when the jet is confined.
A class of contoured wall fuel injectors was designed to enable shock-enhancement of hypervelocity mixing for supersonic combustion ramjet applications. Previous studies of these geometries left unresolved questions concerning the relative importance of various axial vorticity sources in mixing the injectant with the freestream. The present study is a numerical simulation of two generic fuel injectors which is aimed at elucidating the relative roles of axial vorticity sources including: baroclinic torque through shock-impingement, cross-stream shear, turning of boundary layer vorticity, shock curvature, and diffusive flux. Both the magnitude of the circulation, and the location of vorticity with respect to the mixing interface were considered. Baroclinic torque and cross-stream shear were found to be most important in convectively mixing the injectant with the freestream, with the former providing for deposition of vorticity directly on the fuel/air interface.
CGL relations including the effect of finite ion Larmor radius are used to consider a class of parametric instabilities of finite-amplitude, circularly polarized Alfven waves in high-beta plasmas. The disperison relation governing the instabilities is a sixth-order polynomial which is solved numerically. There are two types of instabilities: a modulational instability at k is less than k(0) and a relatively weak and narrow bandwidth instability at k is less than approximately k(0), where k and k(0) are the wavenumbers of the unstable density fluctuation and the 'pump' wave, respectively. It is shown that these instabilities can occur for left-handed pump waves and that the modulational instability is unstable over a very broad band in k with a maximum growth rate at finite k is not equal to 0.
The design, analysis and predicted performance of the Brilliant Eyes Ten-Kelvin Sorption Cryocooler Experiment (BETSCE) is described from a thermal perspective. BETSCE is a shuttle side-wall mounted cryogenic technology demonstration experiment planned for launch in November 1994. BETSCE uses a significant amount of power (about 500 W peak) and the resultant heat must be rejected passively with radiators, as BETSCE has no access to the active cooling capability of the shuttle. It was a major challenge to design and configure the individual hardware assemblies, with their relatively large radiators, to enable them to reject their heat while satisfying numerous severe shuttle-imposed constraints. This paper is a useful case study of a small shuttle payload that needs to reject relatively high heat loads passively in a highly constrained thermal environment. The design approach described is consistent with today's era of 'faster, better, cheaper' small-scale space missions.
The question of entropy stability for discrete approximations to hyperbolic systems of conservation laws is studied. The amount of numerical viscosity present in such schemes is quantified and related to their entropy stability by means of comparison. To this end, two main ingredients are used: entropy variables and the construction of certain entropy conservative schemes in terms of piecewise-linear finite element approximations. It is then shown that conservative schemes are entropy stable, if and (for three-point schemes) only if, they contain more numerical viscosity than the abovementioned entropy conservation ones.
Numerical results obtained in a simplified version of the one on one aerial combat problem are presented. The primary aim of the data is to specify the roles of pursuer and evader as functions of the relative geometry and of the significant physical parameters of the problem. Numerical results are given in a case in which the slower aircraft is more maneuverable than the faster aircraft. A third order dynamic model of the relative motion is described, for which the state variables are relative range, bearing, and heading. The ranges at termination are arbitary in the present version of the problem, so the weapon systems of both aircraft can be visualized as forward firing high velocity weapons, which must be aimed at the tail pipe of the evader. It was found that, for the great majority of the ralative geometries, each aircraft can evade the weapon system of the other.
In this paper, we introduce a third-order nodal-gradient cell-centered finite volume scheme applicable to mixed-element grids based on a single numerical flux per face. Previously, we considered a split-quadrilateral-face approach, where we split each quadrilateral face of a non-tetrahedral cell into two triangles and then apply second- and third-order schemes designed for tetrahedral grids. This approach is relatively simple to implement but expensive because it requires two numerical flux evaluations per quadrilateral face. In this study, we consider another approach, where the surface flux quadrature is performed with a single numerical flux per face with flux corrections. It requires a 3×3 correction matrix to be stored at each face but is more economical than the split-quadrilateral-face approach. The two approaches are compared for three-dimensional flows on irregular hexahedral grids, and their relative merits are discussed.
It has been estimated that aviation accidents are typically preceded by numerous minor incidents arising from the same causal factors that ultimately produced the accident. Accident databases provide in-depth information on a relatively small number of occurrences, however incident databases have the potential to provide insights into the human factors of Remotely Piloted Aircraft System (RPAS) operations based on a larger volume of less-detailed reports. Currently, there is a lack of incident data dealing with the human factors of unmanned aircraft systems. An exploratory study is being conducted to examine the feasibility of collecting voluntary critical incident reports from RPAS pilots. Twenty-three experienced RPAS pilots volunteered to participate in focus groups in which they described critical incidents from their own experience. Participants were asked to recall (1) incidents that revealed a system flaw, or (2) highlighted a case where the human operator contributed to system resilience or mission success. Participants were asked to only report incidents that could be included in a public document. During each focus group session, a note taker produced a de-identified written record of the incident narratives. At the end of the session, participants reviewed each written incident report, and made edits and corrections as necessary. The incidents were later analyzed to identify contributing factors, with a focus on design issues that either hindered or assisted the pilot during the events. A total of 90 incidents were reported. This presentation focuses on incidents that involved the management of the command and control (C2) link. The identified issues include loss of link, interference from undesired transmissions, voice latency, accidental control transfer, and the use of the lost link timer, or lost link OK features.
The effects of various approximations on the results of three-dimensional computation of turbulent compressible convection are examined and related to the phenomenon of bubble breaking in stellar convection. The problem of viscosity and the numerical issue of efficiency versus accuracy in such computations are discussed, and a model is formulated and tested with regard to turbulent viscosity and accuracy. Results on the convergence to a statistically stationary state, the effects of time accuracy on the statistics, the effects of the adjustable constant parameter c(mu), and the effects of aspect ratio, horizontal grid size, and vertical grid distribution are discussed. A preliminary look at the three-dimensional flow structure is given.
As a result of advances related to the combination of computer system technology and numerical modeling, computational aerodynamics has emerged as an essential element in aerospace vehicle design methodology. NASA has, therefore, initiated the Numerical Aerodynamic Simulation (NAS) Program with the objective to provide a basis for further advances in the modeling of aerodynamic flowfields. The Program is concerned with the development of a leading-edge, large-scale computer facility. This facility is to be made available to Government agencies, industry, and universities as a necessary element in ensuring continuing leadership in computational aerodynamics and related disciplines. Attention is given to the requirements for computational aerodynamics, the principal specific goals of the NAS Program, the high-speed processor subsystem, the workstation subsystem, the support processing subsystem, the graphics subsystem, the mass storage subsystem, the long-haul communication subsystem, the high-speed data-network subsystem, and software.
The purpose of this study is to determine if phase synchronization can be used to reduce the net radiated sound power from two rotors. Phase synchronization implies that the rotors have the same rotational speed with a fixed relative azimuthal blade position, or phase. The concept is evaluated both experimentally and numerically. Measurements of source directivity and thrust are initially compared with predictions to confirm that the model accurately captures the relevant trends. The model is then used to explore the design space and identify relevant parameters. Both experimental and numerical results show that the radiated sound power at the blade passage frequency can be reduced by appropriately controlling the relative azimuthal phase of the rotors. Vehicle level predictions are also provided for a notional octocopter, comparing two different modes of operation. Predictions show that phase synchronization can be used to achieve a 4-5 dB reduction of the sound pressure level at the blade passage frequency nearly everywhere on the ground plane beneath the vehicle.