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.
Eight different experiments are carried out, ranging from 8 to 20 days, with a primitive equation model consisting of five atmospheric layers and one oceanic layer. The purpose is to investigate the relative importance of radiative heating and deep cumulus condensation, orography, and initial conditions for the development of the mean monsoon circulation from June to July in the 0-180 deg E, 25 deg S - 55 deg N region. Two alternative initial states are used, one based on the observed monthly mean June pressure distribution, the other on the mean June zonal average pressure distribution. Whereas large-scale condensation and actual land and sea distributions are included in every experiment, deep cumulus cone condensation and radiative heating are always taken together. It is found that the means sea level pressure distribution and low-level flow pattern produced by these simulations are in the main determined by the diabatic heating distribution and are influenced somewhat by orography; they are, however, almost independent of the initial state. The low pressure systems are found to develop faster when diurnal variation of solar radiation is allowed for; this is thought to be due mainly to the more vigorous cumulus activity during the day.
In order to determine the degree to which a coupled dynamic thermodynamic model can reproduce the seasonal cycle of the Antarctic Sea ice in the Weddell Sea, a series of two-year simulations of the Weddell Sea ice pack were carried out. The simulations employed Hibler's (1979) model and used a one-day time step on an 18 x 15 grid with a resolution of 222 km. Daily atmospheric data from 1979 were used to drive the simulations. The simulations yielded a seasonal cycle of ice with maximum and minimum extents close to that observed. Except for portions of the western Weddell, the advance of the ice is found to be primarily thermodynamic in nature, while the rapid decay depends critically on the presence of both leads and lateral ice advection. In early summer the ice motion causes a residual tongue of ice to extend eastward from the Antarctic Peninsula in agreement with observations. Mean ice thicknesses are consistent with observations, as are the mean drift rates of about 5 km/day.
A mathematical model of glaze and rime ice accretion on a two-dimensional airfoil is presented. The model employs standard methods for calculating the flow field and cloud droplet trajectories. This data is used as input to a thermodynamic analysis of the ice accretion process which includes liquid runback. The structure of the model allows for dynamic updating of the droplet collection efficiency and heat and mass transfer processes as the accreted ice shape changes. This results in the improved shape prediction over previously described methods. Simulations made with the model are compared to experimental results obtained for a NACA 0012 airfoil in the NASA/Lewis Icing Research Tunnel. The agreement with experiment is found to be generally satisfactory for some cases, but a need for improvements is indicated by others. Suggestions for improvements are made.
In the present investigation, the 'thin layer' Navier-Stokes equations are used to compute the flow about a delta wing and a strake-delta wing configuration. Both configurations possess blunt noses and rounded leading edges. The computational grid about these configurations is generated using a newly developed three-dimensional grid-generation code which solves a set of Poisson equations with spherical coordinate variables using an alternating direction implicit (ADI) scheme. Computational results are obtained for an isolated wing with 60 deg sweep and a strake-wing configuration with 80-60 deg sweep. Computations for angles of attack in the range from 6 to 30 deg are presented. Attention is given to the effect of the angle of attack, the Mach number, and the Reynolds number. Theoretical results are compared with experimental data.
A one-dimensional, nonperiodic hybrid code in which ion dynamics are treated exactly, while those of electrons are omitted by neglecting electron inertia and pressure, is employed in the simulation of quasi-parallel bow shock structures. It is found that for an Alfven Mach number smaller than 3, the magnetic field profile of the shock is laminar or quasi-laminar, and that the upstream waves are right-hand polarized whistlers. Downstream waves are absent or low in amplitudes. For Alfven Mach numbers greater than 3, the magnetic field profile of the shock is turbulent, and the upstream waves are again right-hand polarized whistlers. The transition from laminar-subcritical to turbulent-supercritical shock structures is shown to be due to the firehose instability which occurs at Alfven Mach numbers greater than about 3.
There exists observational and modeling evidence to the effect that sea surface temperature (SST) anomalies in the equatorial Pacific have a significant effect upon wintertime climate in the Pacific/North American sector of the hemisphere. In connection with the emergence of a more definitive observational basis for designing general circulation model (GCM) experiments to investigate the atmospheric response to equatorial SST anomalies, it was recommended to conduct a series of GCM experiments. The present study involves a series of extended (1-2 month) integrations.Attention is given to a review of previous GCM experiments, the plan for the present study, the general circulation model and its simulation of the winter and summer circulation described by Shukla et al. (1981), tropical precipitation, upper level flow, low-level circulation and temperature, and the tropical tropospheric temperature.
Preliminary results of four runs with a three-dimensional model of the effects of vertical wind shear on cloud top height/temperature structure and the internal properties of isolate midlatitude thunderstorms are reported. The model is being developed as an aid to analyses of GEO remote sensing satellite data. The grid is a 27 x 27 x 20 mesh with 2 km horizontal resolution and 0.9 vertical resolution. The total grid is 54 km on a side and 18 km deep. A second-order Crowley scheme for advecting momentum is extended with a third-order correction for spatial truncation error, and the earth-relative horizontal surface wind components are decreased to 50 percent of their values at 0.45 km. A temperature increase with height is included, together with an initial impulse consisting of a nonrotating cylindrical weak buoyant updraft 10 km in radius. The results of the runs are discussed in terms of the time variation of the vertical velocity extrema, the effects of strong and weak shear on a storm, the cloud top height, the Lagrangian dynamics of a thermal couplet, and data from a real storm.
The vortex-filament method was applied to the simulation of vortex breakdown. The principal vortex region was represented by multiple filaments, and an axial velocity component was induced by a spiral winding of the filaments. An accuracy check was performed for a cylindrical swirling flow with simple analytical expressions for the axial and theta velocities. The result suggests that the flow field can be simulated to any accuracy by increasing the number of filaments. An axisymmetric-type vortex breakdown was simulated, with experimental data serving as upstream conditions. The calculated axial- and theta-velocity contours show the breakdown of the vortex, including a rapid change in the vortex core, followed axially by a recovery zone and then a second breakdown. When three dimensional initial data are used the second breakdown appears to be of the spiral type in correspondence with experimental observations. The present method can easily be used to simulate other types of vortex breakdown or other vortex flows with axial velocity.
The incompressible, Lagrangian, triangular grid code, SPLISH, was converted for the study of flows in and around fuel droplets. This involved developing, testing and incorporating algorithms for surface tension and viscosity. The major features of the Lagrangian method and the algorithms are described. Benchmarks of the algorithms are given. Several calculations are presented for kerosene droplets in air. Finally, extensions which make the code compressible and three dimensional are discussed.
A computer program was developed to predict the trajectory, ground deposition, and drift of liquid sprays injected into the wake of an agricultural aircraft in ground effect. The program uses a horseshoe vortex wake model and includes the effects of liquid droplet evaporation, crosswind, the propeller slipstream, ground effect, and tunnel walls on small scale models. This user's guide includes several case examples demonstrating user options. A complete listing of the FORTRAN program is provided.
The formulation of the problem of finite elastic-viscoplastic deformation of thin beams is described. The finite element analysis of a representative problem is presented and various solution methods are surveyed.
Computer simulations have become an increasingly popular, important and insightful tool for studying space plasmas. This review describes MHD and particle simulations, both of which treat the plasma and the electromagnetic field in which it moves in a self consistent fashion but on drastically different spatial and temporal scales. The complementary roles of simulation, observations and theory are stressed. Several examples of simulations being carried out in the area of magnetospheric plasma physics are described to illustrate the power, potential and limitations of the approach.
Inlet flowfields for airbreathing missiles are calculated by a two-dimensional computational method. A supersonic freestream is assumed to allow the forebody calculation to be uncoupled from the inlet calculation. The inlet calculation employs an implicit, time-marching, finite difference procedure to solve the Euler equations formulated in body-fitted coordinates. The method can be used for a flowfield with both subsonic and supersonic regions and is found to converge rapidly for supercritical inlet operation. For subcritical inlet operation, however, convergence to steady state is slow.
Two-dimensional Euler and Navier-Stokes solutions of the flow through three inlet/diffuser configurations with terminal shock systems are reported. Calculations without bleed indicate that the terminal shock location is very sensitive to the outflow back pressure. For cases where there are no available experimental results, it becomes difficult to estimate the back pressure that will result in a terminal shock. Estimates based on quasi-one-dimensional analysis are not found adequate for complex two-dimensional flows. It is found that since the flow downstream of the terminal shock is subsonic, and what happens at the outflow boundary affects the flow inside the inlet, enough of the subsonic diffuser must be modeled to accurately predict the terminal shock region. The diffuser portion should be fairly long with the outflow boundary occurring in a region of more or less uniform flow to be able to prescribe a uniform back pressure. The third configuration studied was investigated with and without incorporating bleed in the code. It is found that the use of bleed stabilizes the shock location and allows solutions which without bleed result in unstarting of the inlet. Comparisons are made with available experimental data.
A 'thin-layer' Navier-Stokes code capable of predicting steady state viscous flows is applied to complicated three-dimensional flow fields. The code is written in a generalized coordinate system, and a recently developed grid generation procedure is used for the flow-field discretization. Application is made to the vortical flow over a delta wing at high angle of attack, and the computed results are compared with experimental results. The results indicate that the present method can capture the physical phenomenon well. It is observed that a leading-edge separation vortex is formed over the wing, as is a secondary separation vortex near the leading edge. The flow field behind the trailing edge is also well described. Application is also made to the transonic flow over the Shuttle configuration. The result appears to be reasonable even though no experimental data are available for comparison. These results indicate that the present approach is capable of computing complicated three-dimensional flow fields.
Eleven sensitivity experiments that were made with general circulation models to see how land-surface boundary conditions can influence the rainfall, temperature, and motion fields of the atmosphere are discussed. In one group of experiments, different soil moistures or albedos are prescribed as time-invariant boundary conditions. In a second group, different soil moistures or different albedos are initially prescribed, and the soil moisture (but not the albedo) is allowed to change with time according to the governing equations for soil moisture. In a third group, the results of constant versus time-dependent soil moistures are compared. Previously announced in STAR as N83-27536
Previously cited in issue 15, p. 2346, Accession no. A82-31959