Search NASA⌕ Search

SEARCH · Search NASA

Results for “volume integral equation”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

Optimal Observations for Variational Data Assimilation

An important aspect of Ocean state estimation is the design of an observing system that allows the efficient study of climate aspects in the ocean. A solution of the design problem is presented here in terms of optimal observations that emerge as nondimensionalized singular vectors of the modified data resolution matrix. The actual computation is feasible only for scalar quantities in the limit of large observational errors. In the framework of a lo resolution North Atlantic primitive equation model it is demonstrated that such optimal observations when applied to determining the strength of the volume and heat transport across the Greenland-Scotland ridge, perform significantly better than traditional section data. On seasonal to inter-annual time-scales optimal observations are located primarily along the continental shelf and information about heat-transport, wind stress and stratification is being communicated via boundary waves and advective processes. On time-scales of about a month, sea surface height observations appear to be more efficient in reconstructing the cross-ridge heat transport than hydrographic observations. Optimal observations also provide a tool for understanding how the ocean state is effected by anomalies of integral quantities such as meridional heat transport.

Koehl, Armin↗
NASA NTRS1987

Theory versus experiment for the rotordynamic coefficients of labyrinth gas seals. I - A two control volume model

The basic equations are derived for a two-control-volume model for compressible flow in a labyrinth seal. The recirculation velocity in the cavity is incorporated into the model for the first time. The flow is assumed to be completely turbulent and isoenergetic. The wall friction factors are determined using the Blasius formula. Jet flow theory is used for the calculation of the recirculation velocity in the cavity. Linearized zeroth- and first-order perturbation equations are developed for small motion about a centered position by an expansion in the eccentricity ratio. The zeroth-order pressure distribution is found by satisfying the leakage equation while the circumferential velocity distribution is determined by satisfying the momentum equations. The first-order equations are solved by a separation of variable solution. Integration of the resultant pressure distribution along and around the seal defines the reaction force developed by the seal and the corresponding dynamic coefficients.

Scharrer, J. K.↗
NASA NTRS2016

Modeling and Simulation of Radiative Compressible Flows in Aerodynamic Heating Arc-Jet Facility

Numerical simulations of an arc heated flow inside NASA's 20 [MW] Aerodynamics heating facility (AHF) are performed in order to investigate the three-dimensional swirling flow and the current distribution inside the wind tunnel. The plasma is considered in Local Thermodynamics Equilibrium(LTE) and is composed of Air-Argon gas mixture. The governing equations are the Navier-Stokes equations that include source terms corresponding to Joule heating and radiative cooling. The former is obtained by solving an electric potential equation, while the latter is calculated using an innovative massively parallel ray-tracing algorithm. The fully coupled system is closed by the thermodynamics relations and transport properties which are obtained from Chapman-Enskog method. A novel strategy was developed in order to enable the flow solver and the radiation calculation to be preformed independently and simultaneously using a different number of processors. Drastic reduction in the computational cost was achieved using this strategy. Details on the numerical methods used for space discretization, time integration and ray-tracing algorithm will be presented. The effect of the radiative cooling on the dynamics of the flow will be investigated. The complete set of equations were implemented within the COOLFluiD Framework. Fig. 1 shows the geometry of the Anode and part of the constrictor of the Aerodynamics heating facility (AHF). Fig. 2 shows the velocity field distribution along (x-y) plane and the streamline in (z-y) plane.

Pasma flows↗
NASA NTRS1995

Flow simulations about steady-complex and unsteady moving configurations using structured-overlapped and unstructured grids

The limiting factor in simulating flows past realistic configurations of interest has been the discretization of the physical domain on which the governing equations of fluid flow may be solved. In an attempt to circumvent this problem, many Computational Fluid Dynamic (CFD) methodologies that are based on different grid generation and domain decomposition techniques have been developed. However, due to the costs involved and expertise required, very few comparative studies between these methods have been performed. In the present work, the two CFD methodologies which show the most promise for treating complex three-dimensional configurations as well as unsteady moving boundary problems are evaluated. These are namely the structured-overlapped and the unstructured grid schemes. Both methods use a cell centered, finite volume, upwind approach. The structured-overlapped algorithm uses an approximately factored, alternating direction implicit scheme to perform the time integration, whereas, the unstructured algorithm uses an explicit Runge-Kutta method. To examine the accuracy, efficiency, and limitations of each scheme, they are applied to the same steady complex multicomponent configurations and unsteady moving boundary problems. The steady complex cases consist of computing the subsonic flow about a two-dimensional high-lift multielement airfoil and the transonic flow about a three-dimensional wing/pylon/finned store assembly. The unsteady moving boundary problems are a forced pitching oscillation of an airfoil in a transonic freestream and a two-dimensional, subsonic airfoil/store separation sequence. Accuracy was accessed through the comparison of computed and experimentally measured pressure coefficient data on several of the wing/pylon/finned store assembly's components and at numerous angles-of-attack for the pitching airfoil. From this study, it was found that both the structured-overlapped and the unstructured grid schemes yielded flow solutions of comparable accuracy for these simulations. This study also indicated that, overall, the structured-overlapped scheme was slightly more CPU efficient than the unstructured approach.

Newman, James C., III↗
NASA NTRS1998

Nuclear modules for space electric propulsion

Analysis of interplanetary cargo and piloted missions requires calculations of the performances and masses of subsystems to be integrated in a final design. In a preliminary and scoping stage the designer needs to evaluate options iteratively by using fast computer simulations. The Oak Ridge National Laboratory (ORNL) has been involved in the development of models and calculational procedures for the analysis (neutronic and thermal hydraulic) of power sources for nuclear electric propulsion. The nuclear modules will be integrated into the whole simulation of the nuclear electric propulsion system. The vehicles use either a Brayton direct-conversion cycle, using the heated helium from a NERVA-type reactor, or a potassium Rankine cycle, with the working fluid heated on the secondary side of a heat exchanger and lithium on the primary side coming from a fast reactor. Given a set of input conditions, the codes calculate composition. dimensions, volumes, and masses of the core, reflector, control system, pressure vessel, neutron and gamma shields, as well as the thermal hydraulic conditions of the coolant, clad and fuel. Input conditions are power, core life, pressure and temperature of the coolant at the inlet of the core, either the temperature of the coolant at the outlet of the core or the coolant mass flow and the fluences and integrated doses at the cargo area. Using state-of-the-art neutron cross sections and transport codes, a database was created for the neutronic performance of both reactor designs. The free parameters of the models are the moderator/fuel mass ratio for the NERVA reactor and the enrichment and the pitch of the lattice for the fast reactor. Reactivity and energy balance equations are simultaneously solved to find the reactor design. Thermalhydraulic conditions are calculated by solving the one-dimensional versions of the equations of conservation of mass, energy, and momentum with compressible flow.

Difilippo, F. C.↗
NASA NTRS2023

An Evaluation of A Physics-Based Firn Model and A Semi-Empirical Firn Model Across the Greenland Ice Sheet (1980–2020)

The Greenland Ice Sheet's (GrIS) firn layer buffers the ice sheet's contribution to sea level rise by storing meltwater in its pore space. However, available pore space and meltwater retention capability is lost due to ablation of the firn layer and refreezing of meltwater as near-surface ice slabs in the firn. Understanding how firn properties respond to climate is important for constraining the GrIS's future contribution to sea level rise in a warming climate. Observations of firn density provide detailed information about firn properties, but they are spatially and temporally limited. Here we use two firn models, the physics-based SNOWPACK model and the Community Firn Model configured with a semi-empirical densification equation (CFM-GSFC), to quantify firn properties across the GrIS from 1980 through 2020. We use an identical forcing (Modern-Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2) atmospheric reanalysis) for SNOWPACK and the CFM-GSFC in order to isolate firn model differences. To evaluate the models, we compare simulated firn properties, including firn air content (FAC), to measurements from the Surface Mass Balance and Snow on Sea Ice Working Group (SUMup) dataset of snow and firn density. Both models perform well (mean absolute percentage errors of 14 % in SNOWPACK and 16 % in the CFM-GSFC), though their performance is hindered by the spatial resolution of the atmospheric forcing. In the ice-sheet-wide simulations, the 1980–1995 average spatially integrated FAC (i.e., air volume in the firn) for the upper 100 m is 34 645 km 3 from SNOWPACK and 28 581 km 3 from the CFM-GSFC. The discrepancy in the magnitude of the modeled FAC stems from differences in densification with depth and variations in the sensitivity of the models to atmospheric forcing. In more recent years (2005–2020), both models simulate substantial depletion of pore space. During this period, the spatially integrated FAC across the entire GrIS decreases by 3.2 % (−66.6 km 3 yr −1 ) in SNOWPACK and 1.5 % (−17.4 km 3 yr −1 ) in the CFM-GSFC. These differing magnitudes demonstrate how model differences propagate throughout the FAC record. Over the full modeled record (1980–2020), SNOWPACK simulates a loss of pore space equivalent to 3 mm of sea level rise buffering, while the CFM-GSFC simulates a loss of 1 mm. The greatest depletion in FAC is along the margins and especially along the western margin where observations and models show the formation of near-surface, low-permeability ice slabs that may inhibit meltwater storage.

Megan Thompson-Munson↗
NASA NTRS2018

Thin-Film Solar Arrays for Small Spacecraft

The use of thin-film based solar arrays for space applications has long been recognized as an advantageous power generation option. Thinner materials yield a mass savings, equating to lighter launch loads, larger arrays, and/or more payload allocation. Further, their mechanical flexibility lends itself well to stowage and deployment schemes, allowing for a vast improvement to both specific power (W/kg - watts per kilogram) as well as stowed power density (W/m3 - watts per cubic meter). A key application of thin-film space solar arrays is in the small satellite community; where spacecraft are largely power-starved. Their need for higher power generation coupled with the extreme mass and volume restrictions of the small spacecraft bus is driving the requirement for advanced solar arrays. And, when coupled with the relatively short operational requirements (e.g. less than 2 years) of small spacecraft, yields the perfect application for these thin-film arrays. In this presentation, NASA’s recent work on developing the Lightweight Integrated Solar Array and anTenna (LISA-T - also known as Lightweight Integrated Solar Array and Transceiver)), a thin-film solar array for small spacecraft, will be discussed. The array configuration will be shown alongside bench top deployment videos, a summary of environmental testing to date, and forward plans to a flight test. Advanced, robotic additive manufacturing of the thin-film web will also be discussed.

Carr, John A.↗
NASA NTRS1991

The MHOST finite element program: 3-D inelastic analysis methods for hot section components. Volume 1: Theoretical manual

Formulations and algorithms implemented in the MHOST finite element program are discussed. The code uses a novel concept of the mixed iterative solution technique for the efficient 3-D computations of turbine engine hot section components. The general framework of variational formulation and solution algorithms are discussed which were derived from the mixed three field Hu-Washizu principle. This formulation enables the use of nodal interpolation for coordinates, displacements, strains, and stresses. Algorithmic description of the mixed iterative method includes variations for the quasi static, transient dynamic and buckling analyses. The global-local analysis procedure referred to as the subelement refinement is developed in the framework of the mixed iterative solution, of which the detail is presented. The numerically integrated isoparametric elements implemented in the framework is discussed. Methods to filter certain parts of strain and project the element discontinuous quantities to the nodes are developed for a family of linear elements. Integration algorithms are described for linear and nonlinear equations included in MHOST program.

Nakazawa, Shohei↗
NASA NTRS1992

Catalytic methods using molecular oxygen for treatment of PMMS and ECLSS waste streams, volume 2

Catalytic oxidation has proven to be an effective addition to the baseline sorption, ion exchange water reclamation technology which will be used on Space Station Freedom (SSF). Low molecular weight, polar organics such as alcohols, aldehydes, ketones, amides, and thiocarbamides which are poorly removed by the baseline multifiltration (MF) technology can be oxidized to carbon dioxide at low temperature (121 C). The catalytic oxidation process by itself can reduce the Total Organic Carbon (TOC) to below 500 ppb for solutions designed to model these waste waters. Individual challenges by selected contaminants have shown only moderate selectivity towards particular organic species. The combined technology is applicable to the more complex waste water generated in the Process Materials Management System (PMMS) and Environmental Control and Life Support System (ECLSS) aboard SSF. During the phase 3 Core Module Integrated Facility (CMIF) water recovery tests at NASA MSFC, real hygiene waste water and humidity condensate were processed to meet potable specifications by the combined technology. A kinetic study of catalytic oxidation demonstrates that the Langmuir-Hinshelwood rate equation for heterogeneous catalysts accurately represent the kinetic behavior. From this relationship, activation energy and rate constants for acetone were determined.

Akse, James R.↗
NASA NTRS2007

Using MODIS and GLAS Data to Develop Timber Volume Estimates in Central Siberia

The boreal forest is the Earth's largest terrestrial biome, covering some 12 million km2 and accounting for about one third of this planet's total forest area. Mapping of boreal forest's type, structure parameters and biomass are critical for understanding the boreal forest's significance in the carbon cycle, its response to and impact on global climate change. Ground based forest inventories, have much uncertainty in the inventory data, particularly in remote areas of Siberia where sampling is sparse and/or lacking. In addition, many of the forest inventories that do exist for Siberia are now a decade or more old. Thus, available forest inventories fail to capture the current conditions. Changes in forest structure in a particular forest-type and region can change significantly due to changing environment conditions, and natural and anthropogenic disturbance. Remote sensing methods can potentially overcome these problems. Multispectral sensors can be used to provide vegetation cover maps that show a timely and accurate geographic distribution of vegetation types rather than decade old ground based maps. Lidar sensors can be used to directly obtain measurements that can be used to derive critical forest structure information (e.g., height, density, and volume). These in turn can used to estimate biomass components using allometric equations without having to use out dated forest inventory. Finally, remote sensing data is ideally suited to provide a sampling basis for a rigorous statistical estimate of the variance and error bound on forest structure measures. In this study, new remote sensing methods were applied to develop estimates timber volume using NASA's MODerate resolution Imaging Spectroradiometer (MODIS) and unique waveform data of the geoscience laser altimeter system (GLAS) for a 10 deg x 10 deg area in central Siberia. Using MODIS and GLAS data, maps were produced for cover type and timber volume for 2003, and a realistic variance (error bound) for timber volume was calculated for the study area. In this 'study we used only GLAS footprints that had a slope value of less than 10 deg. This was done to avoid large errors due to the effect of slope on the GLAS models. The method requires the integration of new remote sensing methods with available ground studies of forest timber volume conducted in Russian forests. The results were compared to traditional ground forest inventory methods reported in the literature and to ground truth collected in the study area.

Ranson, K. Jon↗
NASA NTRS2020

Apparent dust size discrepancy in aerosol reanalysis in north African dust after long-range transport

North African dust reaches the southeastern United States every summer. Size-resolved dust mass measurements taken in Miami, Florida, indicate that more than one-half of the surface dust mass concentrations reside in particles with geometric diameters less than 2.1 µm, while vertical profiles of micropulse lidar depolarization ratios show dust reaching above 4 km during pronounced events. These observations are compared to the representation of dust in the Modern-Era Retrospective analysis for Research and Applications, version 2 (MERRA-2) aerosol reanalysis and closely related Goddard Earth Observing System model version 5 (GEOS-5) Forward Processing (FP) aerosol product, both of which assimilate satellite-derived aerosol optical depths using a similar protocol and inputs. These capture the day-to-day variability in aerosol optical depth well, in a comparison to an independent sun-photometer-derived aerosol optical depth dataset. Most of the modeled dust mass resides in diameters between 2 and 6 µm, in contrast to the measurements. Model-specified mass extinction efficiencies equate light extinction with approximately 3 times as much aerosol mass, in this size range, compared to the measured dust sizes. GEOS-5 FP surface-layer sea salt mass concentrations greatly exceed observed values, despite realistic winds and relative humidities. In combination, these observations help explain why, despite realistic total aerosol optical depths, (1) free-tropospheric model volume extinction coefficients are lower than those retrieved from the micro-pulse lidar, suggesting too-low model dust loadings in the free troposphere, and (2) model dust mass concentrations near the surface can be higher than those measured. The modeled vertical distribution of dust, when captured, is reasonable. Large, aspherical particles exceeding the modeled dust sizes are also occasionally present, but dust particles with diameters exceeding 10 µm contribute little to the measured total dust mass concentrations after such long-range transport. Remaining uncertainties warrant a further integrated assessment to confirm this study's interpretations.

MEERA-2 aerosols reanalysis↗
NASA NTRS2002

Digital Distortion Caused by Traveling- Wave-Tube Amplifiers Simulated

Future NASA missions demand increased data rates in satellite communications for near real-time transmission of large volumes of remote data. Increased data rates necessitate higher order digital modulation schemes and larger system bandwidth, which place stricter requirements on the allowable distortion caused by the high-power amplifier, or the traveling-wave-tube amplifier (TWTA). In particular, intersymbol interference caused by the TWTA becomes a major consideration for accurate data detection at the receiver. Experimentally investigating the effects of the physical TWTA on intersymbol interference would be prohibitively expensive, as it would require manufacturing numerous amplifiers in addition to acquiring the required digital hardware. Thus, an accurate computational model is essential to predict the effects of the TWTA on system-level performance when a communication system is being designed with adequate digital integrity for high data rates. A fully three-dimensional, time-dependent, TWT interaction model has been developed using the electromagnetic particle-in-cell code MAFIA (Solution of Maxwell's equations by the Finite-Integration-Algorithm). It comprehensively takes into account the effects of frequency-dependent AM (amplitude modulation)/AM and AM/PM (phase modulation) conversion, gain and phase ripple due to reflections, drive-induced oscillations, harmonic generation, intermodulation products, and backward waves. This physics-based TWT model can be used to give a direct description of the effects of the nonlinear TWT on the operational signal as a function of the physical device. Users can define arbitrary excitation functions so that higher order modulated digital signals can be used as input and that computations can directly correlate intersymbol interference with TWT parameters. Standard practice involves using communication-system-level software packages, such as SPW, to predict if adequate signal detection will be achieved. These models use a nonlinear, black-box model to represent the TWTA. The models vary in complexity, but most make several assumptions regarding the operation of the high-power amplifier. When the MAFIA TWT interaction model was used, these assumptions were found to be in significant error. In addition, digital signal performance, including intersymbol interference, was compared using direct data input into the MAFIA model and using the system-level analysis tool SPW for several higher order modulation schemes. Results show significant differences in predicted degradation between SPW and MAFIA simulations, demonstrating the significance of the TWTA approximations made in the SPW model on digital signal performance. For example, a comparison of the SPW and MAFIA output constellation diagrams for a 16-ary quadrature amplitude modulation (16-QAM) signal (data shown only for second and fourth quadrants) is shown. The upper-bound degradation was calculated from the corresponding eye diagrams. In comparison to SPW simulations, the MAFIA data resulted in a 3.6-dB larger degradation.

Kory, Carol L.↗
NASA NTRS1993

Fuel injector: Air swirl characterization aerothermal modeling, phase 2, volume 1

A well integrated experimental/analytical investigation was conducted to provide benchmark quality relevant to a prefilming type airblast fuel nozzle and its interaction with the combustor dome air swirler. The experimental investigation included a systematic study of both single-phase flows that involved single and twin co-axial jets with and without swirl. A two-component Phase Doppler Particle Analyzer (PDPA) was used to document the interaction of single and co-axial air jets with glass beads that simulate nonevaporating spray and simultaneously avoid the complexities associated with fuel atomization processes and attendant issues about the specification of relevant boundary conditions. The interaction of jets with methanol spray produced by practical airblast nozzle was also documented in the spatial domain of practical interest. Model assessment activities included the use of three turbulence models (k-epsilon, algebraic second moment (ASM), and differential second moment (DSM)) for the carrier phase, deterministic or stochastic Lagrangian treatment of the dispersed phase, and advanced numerical schemes. Although qualitatively good comparison with data was obtained for most of the cases investigated, the model deficiencies in regard to modeled dissipation rate transport equation, single length scale, pressure-strain correlation, and other critical closure issues need to be resolved before one can achieve the degree of accuracy required to analytically design combustion systems.

Nikjooy, M.↗
NASA NTRS1993

Fuel Injector: Air swirl characterization aerothermal modeling, phase 2, volume 2

A well integrated experimental/analytical investigation was conducted to provide benchmark quality data relevant to prefilming type airblast fuel nozzle and its interaction with combustor dome air swirler. The experimental investigation included a systematic study of both single-phase flows that involved single and twin co-axial jets with and without swirl. A two-component Phase Doppler Particle Analyzer (PDPA) equipment was used to document the interaction of single and co-axial air jets with glass beads that simulate nonevaporating spray and simultaneously avoid the complexities associated with fuel atomization processes and attendant issues about the specification of relevant boundary conditions. The interaction of jets with methanol spray produced by practical airblast nozzle was also documented in the spatial domain of practical interest. Model assessment activities included the use of three turbulence models (k-epsilon, algebraic second moment (ASM) and differential second moment (DSM)) for the carrier phase, deterministic or stochastic Lagrangian treatment of the dispersed phase, and advanced numerical schemes. Although qualitatively good comparison with data was obtained for most of the cases investigated, the model deficiencies in regard to modeled dissipation rate transport equation, single length scale, pressure-strain correlation, and other critical closure issues need to be resolved before one can achieve the degree of accuracy required to analytically design combustion systems.

Nikjooy, M.↗
NASA NTRS2009

Method for Measuring the Volume-Scattering Function of Water

The volume scattering function (VSF) of seawater affects visibility, remote sensing properties, in-water light propagation, lidar performance, and the like. Currently, it s possible to measure only small forward angles of VSF, or to use cumbersome, large, and non-autonomous systems. This innovation is a method of measuring the full range of VSF using a portable instrument. A single rapid-sensing photosensor is used to scan a green laser beam, which delivers the desired measurement. By using a single sensor, inter-calibration is avoided. A compact design is achieved by using drift-free detector electronics, fiber optics, and a new type of photomultiplier. This provides a high angular resolution of 1 or better, as well as the ability to focus in on a VSF region of particular interest. Currently, the total scattering of light is measured as a difference from the other two parts of the light budget equation. This innovation will allow the direct calculation of the total scattering of light by taking an integral of the VSF over all angles. This directly provides one of the three components of the light budget equation, allowing greater versatility in its calculation.

Agrawal, Yogesh C.↗
NASA NTRSDate not supplied

Explicit Discontinuous Galerkin Methods for Conservation Laws

The two explicit DG methods in this study are based on a ‘predictor-corrector’ formulation, the first introduced by Lörcher, Gassner, and Munz (2007, 2008) called space–time expansion discontinuous Galerkin or STE-DG scheme, and the second, introduced independently by the author (Huynh 2006, 2013) called the upwind moment scheme. The predictor step of the two methods is essentially identical using a Cauchy-Kovalevsky (CK) procedure, which involves no interaction of the data among neighboring cells. The corrector step also shares the same space-time integration formulation and is where interaction of the data among neighboring cells takes place; the difference, however, is in how the resulting space-time volume integral is estimated. As a consequence of the different estimates, for the case of advection in one spatial dimension (1D), the moment scheme has a CFL (Courant-Friedrichs-Lewy) condition of 1 for all p and is accurate to order 2p+1, i.e., it possesses the super accuracy property, whereas the STE-DG method has a more restrictive CFL condition and is accurate to the expected order of p+1. For 1D advection, compared with the CFL conditions of 1/(2p+1) of standard RK-DG (Runge-Kutta) scheme where space and time discretization are of the same order, the moment scheme allows a significantly larger time step size. It also turns out that the scheme yields a result identical to Van Leer’s scheme III (1977), which amounts to shifting the data a distance of advection corresponding to the time step and projecting the result onto the space of polynomial solutions. Contrary to Van Leer’s approach, however, the space-time ‘predictor-corrector’ formulation facilitates extensions to the case of systems of equations. Concerning 2D extensions, in the case of advection, when the flow is along the diagonal direction, the CFL conditions for the moment schemes become restrictive as will be shown by Fourier (Von Neumann) stability and accuracy analyses. Since the moment scheme employs the right Radau points as collocation points in time, the method is closely related to the implicit Radau IIA scheme, which is stable for any time step size. The role of Radau IIA in relieving stability restriction for these explicit DG schemes remains to be explored

Discontinuous Galerkin↗
NASA NTRS2000

Model Study of the Sources of Upper Tropospheric Ozone Over The Atlantic, with Emphasis on Interhemispheric Asymmetries

Measurements of ozone and aerosols from ships transecting the Atlantic Ocean reveal considerable interhemispheric asymmetries in the vertical structure of the ozone distribution. There is much more upper tropospheric ozone to the South of the Equator than in the northern hemisphere. These profiles suggest that the air in the southern hemispheric tropics and subtropics has a source in the stratosphere, while that in the northern hemisphere originates in the boundary layer. A model experiment designed to investigate this problem will be discussed. They are based on a new general circulation model, based on a Finite-Volume dynamical core and the physical parametrizations from the Community Climate Model at the National Center for Atmospheric Research (NCAR). The model simulates a climatological state in reasonable agreement with observations, comparable to other contemporary climate models, but has a much more accurate transport capability than most other present-day models. The experiment involves using a suite of tracers, designed to highlight air from different source regions. Some of these (water vapor and ozone) represent real atmospheric quantities, while others are markers which indicate whether the air originates in the stratosphere, the boundary layer, or has traversed regions of convection (for example). The experiments were run for periods of up to one month, initialized from snapshots from a long climate integration of the model. A sample of results will be presented, chosen to illustrate the different sources of air in the upper troposphere and the physical mechanisms by which these air parcels arrived there.

Pawson, Steven↗
NASA NTRS1998

Effect of Apex Flap Deflection on Vertical Tail Buffeting

A computational study of the effect of vortex breakdown location on vertical tail buffeting is conducted. The position of the breakdown is modified by employing an apex flap deflected by an experimentally determined optimal angle. The delayed breakdown flow and buffeting response is then compared to the nominal undeflected case. This multidisciplinary problem is solved sequentially for the fluid flow, the elastic tail deformations and the grid displacements. The fluid flow is simulated by time accurately solving the unsteady, compressible, Reynolds-averaged Navier-Stokes equations using an implicit, upwind, flux-difference splitting finite volume scheme. The elastic vibrations of the tails are modeled by uncoupled bending and torsion beam equations. These equations are solved accurately in time using the Galerkin method and a five-stage Runge-Kutta-Verner scheme. The grid for the fluid dynamics calculations is continuously deformed using interpolation functions to disperse the displacements smoothly throughout the computational domain. An angle-of-attack of 35 deg.is chosen such that the wing primary-vortex cores experience vortex breakdown and the resulting turbulent wake flow impinges on tile vertical tails. The dimensions and material properties of the vertical tails are chosen such that the deflections are large enough to insure interaction with the flow, and the natural frequencies are high enough to facilitate a practical computational solution. Results are presented for a baseline uncontrolled buffeting case and a delayed breakdown case in which the apex flap has been deflected 15 deg. The flap was found to be very effective in delaying the breakdown, increasing the location from 50%c to 94%c, which resulted in a 6% increase in lift coefficient and pitching moment. However, the integrated buffet loads and tip responses were roughly equivalent for the two cases.

Massey, Steven J.↗