SEARCH · Search NASA

Results for “atmospheric drag modeling”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 325 records · Page 18

All sources NASA NTRS DOE OSTI

NASA NTRS2006

Structural Analysis and Testing of the Inflatable Re-entry Vehicle Experiment (IRVE)

The Inflatable Re-entry Vehicle Experiment (IRVE) is a 3.0 meter, 60 degree half-angle sphere cone, inflatable aeroshell experiment designed to demonstrate various aspects of inflatable technology during Earth re-entry. IRVE will be launched on a Terrier-Improved Orion sounding rocket from NASA s Wallops Flight Facility in the fall of 2006 to an altitude of approximately 164 kilometers and re-enter the Earth s atmosphere. The experiment will demonstrate exo-atmospheric inflation, inflatable structure leak performance throughout the flight regime, structural integrity under aerodynamic pressure and associated deceleration loads, thermal protection system performance, and aerodynamic stability. Structural integrity and dynamic response of the inflatable will be monitored with photogrammetric measurements of the leeward side of the aeroshell during flight. Aerodynamic stability and drag performance will be verified with on-board inertial measurements and radar tracking from multiple ground radar stations. In addition to demonstrating inflatable technology, IRVE will help validate structural, aerothermal, and trajectory modeling and analysis techniques for the inflatable aeroshell system. This paper discusses the structural analysis and testing of the IRVE inflatable structure. Equations are presented for calculating fabric loads in sphere cone aeroshells, and finite element results are presented which validate the equations. Fabric material properties and testing are discussed along with aeroshell fabrication techniques. Stiffness and dynamics tests conducted on a small-scale development unit and a full-scale prototype unit are presented along with correlated finite element models to predict the in-flight fundamental mod

Lindell, Michael C.↗

NASA NTRS1995

Coupled Ablation, Heat Conduction, Pyrolysis, Shape Change and Spallation of the Galileo Probe

The Galileo probe enters the atmosphere of Jupiter in December 1995. This paper presents numerical methodology and detailed results of our final pre-impact calculations for the heat shield response. The calculations are performed using a highly modified version of a viscous shock layer code with massive radiation coupled with a surface thermochemical ablation and spallation model and with the transient in-depth thermal response of the charring and ablating heat shield. The flowfield is quasi-steady along the trajectory, but the heat shield thermal response is dynamic. Each surface node of the VSL grid is coupled with a one-dimensional thermal response calculation. The thermal solver includes heat conduction, pyrolysis, and grid movement owing to surface recession. Initial conditions for the heat shield temperature and density were obtained from the high altitude rarefied-flow calculations of Haas and Milos. Galileo probe surface temperature, shape, mass flux, and element flux are all determined as functions of time along the trajectory with spallation varied parametrically. The calculations also estimate the in-depth density and temperature profiles for the heat shield. All this information is required to determine the time-dependent vehicle mass and drag coefficient which are necessary inputs for the atmospheric reconstruction experiment on board the probe.

Milos, Frank S.↗

DOE OSTI2025

Improving the QBO Forcing by Resolved Waves With Vertical Grid Refinement in E3SMv2

Abstract The quasi‐biennial oscillation (QBO) is the dominate mode of variability in the tropical stratosphere and plays an important role in stratospheric dynamics and chemistry. The QBO is notably deficient in many climate models, including the Energy Exascale Earth System Model (E3SM) developed by the US Department of Energy. In this work, we refine the lower stratospheric vertical grid spacing from roughly 1 km to 500 m to facilitate more realistic equatorial wave activity in the lower stratosphere in E3SM version 2. The refinement results in a simulated QBO with a reasonable amplitude and easterly‐westerly transition in both directions, but still has a longer period than observed, slower easterly downward propagation speed, and shallower vertical depth. Similar refinement in the multi‐scale modeling framework configuration of E3SM yields similar improvements. By analyzing the forcing contributions from different wave types, we find that most of the QBO forcing still comes from parameterized gravity wave drag from convection. The improved QBO forcing contributions from resolved waves, especially equatorial Kelvin waves and resolved small scale waves, can be attributed to the grid refinement.

E3SM↗

NASA NTRS2011

Stratospheric Collection of Dust from Comet 73P/Schwassmann-Wachmann 3

Interplanetary dust particles (IDPs) collected in the stratosphere are unique materials that are compositionally distinct from meteorites. Astronomical observations and dynamical models indicate that both asteroids and short-period comets are significant sources of IDPs. IDPs having fragile, porous structures, unequilibrated, anhydrous mineralogy, and high atmospheric entry velocities are thought to derive from comets, whereas asteroidal IDPs are identified by their compact structure, hydrated mineralogy and low atmospheric entry velocities. Uncertainty remains in the classification of asteroidal and cometary IDPs owing to our limited sampling of comets and the asteroid belt and the complex dynamical histories of most IDPs in space. Most IDPs spend thousands of years in space prior to being accreted by the Earth. During this time, dust particles undergo orbital evolution, including gradual reduction in their perihelion and eccentricity as a result of Poynting-Robertson drag. Planetary encounters may also significantly change their orbital parameters. Consequently, it is generally not possible to identify the specific parent body of a given IDP. However, it has been proposed that it is possible to identify dust from comets that have formed Earth-crossing dust trails. In this case, the dust particles have been in space for such a short period of time (a few decades or less) that their orbits have not significantly changed. Furthermore, these fresh IDPs could be identified in the laboratory from their short space-exposure histories (low solar noble gas abundance and lack of solar flare tracks). NASA flew several dedicated IDP collection missions attempting to collect dust from comet 26P/Grigg-Skjellerup, the best candidate identified. Remarkably, many particles from those collectors exhibit unusual properties, including low abundances of solar noble gases and high abundances of presolar grains. These observations are consistent with the dust particles originating from comet Grigg-Skjellerup (hereafter G-S). This study considers the prospects for collection of dust from comet 73P/Schwassmann-Wachmann 3 (hereafter SW3). SW3 is a small (2 km diameter) Jupiter family comet whose perihelion is close to and just inside the Earth's orbit. The orbit of SW3 is suitable for producing a low-velocity Earth-crossing dust stream and is the likely parent of the Tau Herculid meteor stream. This study complements a previously published model of the SW3 meteor stream that predicted a very low level of activity for grains 100 micron -- 100 mm in size.

Messenger, Scott R.↗

NASA NTRS2000

Feeling Gravity's Pull: Gravity Modeling. The Gravity Field of Mars

Most people take the constant presence of gravitys pull for granted. However, the Earth's gravitational strength actually varies from location to location. This variation occurs because mass, which influences an object's gravitational pull, is not evenly distributed within the planet. Changes in topography, such as glacial movement, an earthquake, or a rise in the ocean level, can subtly affect the gravity field. An accurate measurement of the Earth's gravity field helps us understand the distribution of mass beneath the surface. This insight can assist us in locating petroleum, mineral deposits, ground water, and other valuable substances. Gravity mapping can also help notice or verify changes in sea surface height and other ocean characteristics. Such changes may indicate climate change from polar ice melting and other phenomena. In addition, gravity mapping can indicate how land moves under the surface after earthquakes and other plate tectonic processes. Finally, changes in the Earth's gravity field might indicate a shift in water distribution that could affect agriculture, water supplies for population centers, and long-term weather prediction. Scientists can map out the Earth's gravity field by watching satellite orbits. When a satellite shifts in vertical position, it might be passing over an area where gravity changes in strength. Gravity is only one factor that may shape a satellite's orbital path. To derive a gravity measurement from satellite movement, scientists must remove other factors that might affect a satellite's position: 1. Drag from atmospheric friction. 2. Pressure from solar radiation as it heads toward Earth and. as it is reflected off the surface of the Earth 3. Gravitational pull from the Sun, the Moon, and other planets in the Solar System. 4. The effect of tides. 5. Relativistic effects. Scientists must also correct for the satellite tracking process. For example, the tracking signal must be corrected for refraction through the atmosphere of the Earth. Supercomputers can calculate the effect of gravity for specific locations in space following a mathematical process known as spherical harmonics, which quantifies the gravity field of a planetary body. The process is based on Laplace's fundamental differential equation of gravity. The accuracy of a spherical harmonic solution is rated by its degree and order. Minute variations in gravity are measured against the geoid, a surface of constant gravity acceleration at mean sea level. The geoid reference gravity model strength includes the central body gravitational attraction (9.8 m/sq s) and a geopotential variation in latitude partially caused by the rotation of the Earth. The rotational effect modifies the shape of the geoid to be more like an ellipsoid, rather than a perfect, circle. Variations of gravity strength from the ellipsoidal reference model are measured in units called milli-Galileos (mGals). One mGal equals 10(exp -5) m/sq s. Research projects have also measured the gravity fields of other planetary bodies, as noted in the user profile that follows. From this information, we may make inferences about our own planet's internal structure and evolution. Moreover, mapping the gravity fields of other planets can help scientists plot the most fuel-efficient course for spacecraft expeditions to those planets.

Lemoine, Frank↗

DOE OSTI2024

X-ray imaging and electron temperature evolution in laser-driven magnetic reconnection experiments at the national ignition facility

We present results from x-ray imaging of high-aspect-ratio magnetic reconnection experiments driven at the National Ignition Facility. Two parallel, self-magnetized, elongated laser-driven plumes are produced by tiling 40 laser beams. A magnetic reconnection layer is formed by the collision of the plumes. A gated x-ray framing pinhole camera with micro-channel plate detector produces multiple images through various filters of the formation and evolution of both the plumes and current sheet. As the diagnostic integrates plasma self-emission along the line of sight, two-dimensional electron temperature maps ⟨Te⟩Y are constructed by taking the ratio of intensity of these images obtained with different filters. The plumes have a characteristic temperature ⟨Te⟩Y=240 ± 20 eV at 2 ns after the initial laser irradiation and exhibit a slow cooling up to 4 ns. The reconnection layer forms at 3 ns with a temperature ⟨Te⟩Y=280 ± 50 eV as the result of the collision of the plumes. The error bars of the plumes and current sheet temperatures separate at 4 ns, showing the heating of the current sheet from colder inflows. Using a semi-analytical model, we survey various heating mechanisms in the current sheet. We find that reconnection energy conversion would dominate at low density (ne≲7×1018 cm−3) and electron-ion collisional drag at high-density (≳1019 cm−3).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1972

Convection electric fields and polar thermospheric winds.

Use of the qualitative ideas of convection electric fields over the earth's polar regions to demonstrate the importance of ion drag in establishing a thermospheric wind system. Recent measurements indicate that uniform electric fields of 10 to 40 mV/m are a regular feature of the polar-cap ionosphere. Calculations of the neutral thermospheric wind, using these measured fields in a simple ionospheric model, have been made. The time scale for motion of the neutral gas ranges from less than 1 hour at F-region heights to about 2 hours in the dynamo region of the ionosphere. It has been found that the viscosity of the atmosphere is important in determining the winds in the dynamo region. Results are given that show ion-temperature enhancements of hundreds of degrees that are due to ion-neutral frictional effects. In addition, the total deposition rate of convection energy in the polar thermosphere is shown to be of the same order of magnitude as that due to absorption of solar EUV radiation. The implications of these results for the dynamics and energetics of the thermosphere are discussed.

Fedder, J. A.↗

NASA NTRS1999

Climatological Structures of the GRIPS Models: Mean States and Forcing

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.

Pawson, Steven↗

NASA NTRS1977

The low-latitude circulation of Mars

A steady-state two-layer atmospheric model is used in an attempt to simulate the circulation of the equatorial regions of Mars, with consideration of both thermal and mechanical orographic effects. A steady-state solution is sought with the aim of producing mean seasonal circulations upon which tides can be superimposed, and the relative importance of orographic forcing at low latitudes is evaluated. It is assumed that the atmospheres over a topographic feature and over flat terrain possess the same radiative and convective properties, simple dissipation laws are introduced which represent a surface drag effect, and six first-order difference equations are obtained and then reduced to two linear complex second-order coupled difference equations. Analysis of the response of the equinoctial basic wind field to a Gaussian orographic feature situated at the equator and zero deg longitude reveals a number of similarities and differences between the Martian and terrestrial responses. Most of the similarities are explained in terms of the dynamics of a rotating atmosphere, while the major differences are at least partially explained by the different optical properties, lower boundaries, and smaller mass of the Martian atmosphere.

Webster, P. J.↗

NASA NTRS1990

Clast sizes in terrestrial and Martian ignimbrite lag deposits

A theoretical model for predicting the sizes of clasts which may be expected in proximal ignimbrite deposits on the earth and Mars is developed. The model can provide an explanation of the eruption mechanism. It takes into account the inside and outside conditions of a vent from which a large-scale pyroclastic flow takes place. The drag forces exerted by the dusty gas on relatively large clasts are also taken into account. The sizes of the largest clasts that contribute to a proximal lag breccia deposit can thus be deduced. The analysis has shown that the pressure at the surface vent level in an ignimbrite-forming eruption is independent of the planetary atmospheric pressure; the gravity and vent size then become the major factors controlling the size of the largest clast which can be erupted.

Wilson, Lionel↗

NASA NTRS2014

Elastic Model Transitions: A Hybrid Approach Utilizing Quadratic Inequality Constrained Least Squares (LSQI) and Direct Shape Mapping (DSM)

A method for transitioning linear time invariant (LTI) models in time varying simulation is proposed that utilizes a hybrid approach for determining physical displacements by augmenting the original quadratically constrained least squares (LSQI) algorithm with Direct Shape Mapping (DSM) and modifying the energy constraints. The approach presented is applicable to simulation of the elastic behavior of launch vehicles and other structures that utilize discrete LTI finite element model (FEM) derived mode sets (eigenvalues and eigenvectors) that are propagated throughout time. The time invariant nature of the elastic data presents a problem of how to properly transition elastic states from the prior to the new model while preserving motion across the transition and ensuring there is no truncation or excitation of the system. A previous approach utilizes a LSQI algorithm with an energy constraint to effect smooth transitions between eigenvector sets with no requirement that the models be of similar dimension or have any correlation. This approach assumes energy is conserved across the transition, which results in significant non-physical transients due to changing quasi-steady state energy between mode sets, a phenomenon seen when utilizing a truncated mode set. The computational burden of simulating a full mode set is significant so a subset of modes is often selected to reduce run time. As a result of this truncation, energy between mode sets may not be constant and solutions across transitions could produce non-physical transients. In an effort to abate these transients an improved methodology was developed based on the aforementioned approach, but this new approach can handle significant changes in energy across mode set transitions. It is proposed that physical velocities due to elastic behavior be solved for using the LSQI algorithm, but solve for displacements using a two-step process that independently addresses the quasi-steady-state and non-steady-state contributions to the elastic displacement. For structures subject to large external forces, such as thrust or atmospheric drag, it is imperative to capture these forces when solving for elastic displacement. To simplify the mathematical formulation, assumptions are made regarding mass matrix normalization, constant external forcing, and constant viscous damping. These simplifications allow for direct solutions to the quasi-steady-state displacements through a process titled Direct Shape Mapping. DSM solves for the displacements using the eigenvalues of the elastic modes and the external forcing and returns a set of elastic displacements dictated by the eigenvectors of the post-transition mode set. For the non-steady-state contributions to displacement we formulate a LSQI problem that is constrained by energy of the non-steady state terms. The contributions from the quasi-steady-state and non-steady state solutions are then combined to obtain the physical displacements associated with the new set of eigenvectors. Results for the LSQI-DSM approach show significant reduction/complete removal of transients across mode set transitions while maintaining elastic motion from the prior state. For time propagation applications employing discrete elastic models that need to be transitioned in time and where running with full a full mode set is not feasible, the method developed offers a practical solution to simulating vehicle elasticity.

Hannan, Mike R.↗

NASA NTRS2012

Simulation Tools Model Icing for Aircraft Design

Here s a simple science experiment to try: Place an unopened bottle of distilled water in your freezer. After 2-3 hours, if the water is pure enough, you will notice that it has not frozen. Carefully pour the water into a bowl with a piece of ice in it. When it strikes the ice, the water will instantly freeze. One of the most basic and commonly known scientific facts is that water freezes at around 32 F. But this is not always the case. Water lacking any impurities for ice crystals to form around can be supercooled to even lower temperatures without freezing. High in the atmosphere, water droplets can achieve this delicate, supercooled state. When a plane flies through clouds containing these droplets, the water can strike the airframe and, like the supercooled water hitting the ice in the experiment above, freeze instantly. The ice buildup alters the aerodynamics of the plane - reducing lift and increasing drag - affecting its performance and presenting a safety issue if the plane can no longer fly effectively. In certain circumstances, ice can form inside aircraft engines, another potential hazard. NASA has long studied ways of detecting and countering atmospheric icing conditions as part of the Agency s efforts to enhance aviation safety. To do this, the Icing Branch at Glenn Research Center utilizes a number of world-class tools, including the Center s Icing Research Tunnel and the NASA 607 icing research aircraft, a "flying laboratory" for studying icing conditions. The branch has also developed a suite of software programs to help aircraft and icing protection system designers understand the behavior of ice accumulation on various surfaces and in various conditions. One of these innovations is the LEWICE ice accretion simulation software. Initially developed in the 1980s (when Glenn was known as Lewis Research Center), LEWICE has become one of the most widely used tools in icing research and aircraft design and certification. LEWICE has been transformed over the years from strictly a research tool to one used routinely by industry and other government agencies. Glenn contractor William Wright has been the architect of this development, supported by a team of researchers investigating icing physics, creating validation data, and ensuring development according to standard software engineering practices. The program provides a virtual simulation environment for determining where water droplets strike an airfoil in flight, what kind of ice would result, and what shape that ice would take. Users can enter geometries for specific, two-dimensional cross sections of an airfoil or other airframe surface and then apply a range of inputs - different droplet sizes, temperatures, airspeeds, and more - to model how ice would build up on the surface in various conditions. The program s versatility, ease of use, and speed - LEWICE can run through complex icing simulations in only a few minutes - have contributed to it becoming a popular resource in the aviation industry.

Source record↗

NASA NTRS1982

Comparison of evolving photovoltaic and nuclear power systems for earth orbital applications

Photovoltaic and fission reactor orbital power systems are compared in terms of the end-to-end system power-to-mass ratios. Three PV systems are examined, i.e., a solid substrate with a cell array and a NiCd battery, a modified SEP array and an NiH2 battery, and a 62-micron Si cell array and a fuel cell. All arrays were modeled to be 13.5% efficient and to produce 25 kW dc. The SP-100 reactor consists of the heat source, radiation shield, heat pipes to transfer thermal energy from the reactor to thermoelectric elements, and a waste heat radiator. Consideration is given to system applications in orbits ranging from LEO to GEO, and to mission durations of 1, 5, and 10 yr. PV systems are concluded to be flight-proven, useful out of radiation belts, and best for low to moderate power levels. Limitations exist for operations where atmospheric drag may become a factor and due to the size of a large PV power supply. Space nuclear reactors will continue under development and uses at high power levels and in low altitude orbits are foreseen.

Rockey, D. E.↗

NASA NTRSDate not supplied

Small Satellite-sized Hypersonic Inflatable Aerodynamic Decelerators for Interplanetary Science Missions

To make the most of ridesharing opportunities, small satellite (SmallSat) mission designers endeavor to pack as much payload into a SmallSat-class form factor as possible. The mass and volume constraints of this smaller vehicle class present a challenge for interplanetary mission sets that require a means of achieving orbit insertion at their destination of interest. For a fully propulsive orbit insertion design, this may translate to the propellant mass being a significant fraction of the overall vehicle mass and prolonged insertion time. Aerocapture is a single quick maneuver that can significantly reduce the required propellant mass for orbit insertion. Because aerocapture uses a planet’s atmosphere to achieve the necessary change in velocity, a protective aeroshell is needed. The constraints imposed on secondary payloads render traditional rigid aeroshells mass and space prohibitive for the SmallSat class of vehicles; thus, warranting consideration of deployable designs that can be stowed compactly until needed for atmospheric entry. The Hypersonic Inflatable Aerodynamic Decelerator (HIAD) is a deployable aeroshell that leverages inflatable toroids to achieve the large drag area needed for aerodynamic deceleration. While the technology is currently being analyzed for Mars human-scale missions, it has the potential applicability for interplanetary SmallSat-scale missions as well. This paper highlights a study conducted during an internship at NASA Langley Research Center to investigate the feasibility of using a scaled-down HIAD design in SmallSat aerocapture missions. Several scaling methodologies are investigated including use of parametric models and direct computer-aided design (CAD) model scaling. Candidate HIAD configurations that conform to secondary payload adapter requirements are identified. The Program to Optimize Simulated Trajectories II (POST2) is utilized to conduct orbit insertion performance and trajectory sensitivity studies using the candidate configurations at Earth, Venus, and Mars. The results of the study indicate that multiple SmallSat-sized HIAD designs, targeting a range of SmallSat payload classes, are feasible for planetary aerocapture missions to Mars and Venus as well as Earth-based aerocapture missions.

Shelly C. Mann↗

DOE OSTI2023

Terrain-Influenced Winds and Fire-Fire Interactions in Wildland Fire Simulations [Dissertation]

Ensemble-based approaches to prescribed fire planning cannot be supported by computational fluid dynamics based models like FIRETEC and the Wildland-Urban Interface Fire Dynamics Simulator (WFDS) because they are too computationally expensive and cannot leverage large eddy simulation approaches like CAWFE and WRF-SFIRE because they have too coarse of resolution. QUIC-Fire was developed to fill this gap but it cannot currently address complex terrain, that is typical for instance in the Western United States. This dissertation describes a variety of improvements made to QUIC-Fire and its various incorporated algorithms in an effort to make it a viable tool in simulating wildland and prescribed fires on terrain. The modifications made to QUIC-Fire are described in three chapters. The first chapter describes the extension of the diagnostic wind model QUIC-URB, the wind engine of QUIC-Fire, to a terrain-following coordinate system. The terraininfluenced winds it generates are analyzed and compared. In particular, this chapter presents the mathematical derivation of the wind solver leading to a linear system of equations that are solved through the successive over-relaxation method. The model is validated against a standard test used in previous works (the Askervein Hill) and against a new dataset from measurements in the Socorro Mountains, New Mexico. The terrain-following implementation captures the correct phenomenology for the isolated Askervein Hill, with a wind speed up at the top of the hill. The model agrees well with measurements on the upwind side of the peak, but overestimates speed-up on the downwind side of the hill. This is due to the inability of the model to generate flow separation and wake-eddy dynamics. On a common laptop, the divergence-free wind field is obtained in 6 s, making the solver appealing for coupled fire-atmosphere simulations. The Socorro Mountain is highly complex, with many cliff faces, peaks, and valleys. Although the model captures the magnitude and direction of inlet and outlet areas of the domain, it performs rather poorly in the valley region and in the regions near the steep cliffs. Hence, the model shows good agreement with data in areas of open sloped terrain but lacks in areas where flow separation and thermally driven effects may be present (neither effect is addressed in this work). In the second chapter the implementation of the terrain-following version of QUIC-URB into QUIC-Fire, and the necessary changes needed to include terrain are described. No changes to the underlying fire spread algorithm are made other than what is required to correctly account for the inclusion of terrain. Previously published FIRETEC results that use five different topographies that share the same centerline profile are compared to simulation results from the modified QUIC-Fire that use the same topographies and fuels. QUIC-Fire results show overall similar behaviors in terms of how the topographies affect fire shapes and trends in spread rates. Due to the terrain-following version of QUIC-URB being unable to generate flow separations at the crest of hills, fire spread rates in these regions across all topographies are over-predicted when compared to FIRETEC. Lateral fire growth shows similar trends with FIRETEC between topographies but does not capture the increase in spread due to a diagonal interface between grassland and forested fuel region of the domain. These results suggest that there are three algorithms within QUIC-Fire that could use improvement: how flame tilt angle is accounted for, the incorporation of non-local drag effects, and the inclusion of the wake-eddy parameterizations that are used in QUIC-URB. Lastly, the third chapter describes a modification to the initial guess used for the calculation of the QUIC-URB mass-conserved wind solution during fire simulations. The modification is aimed at improving fire-fire interactions in QUIC-Fire simulations. The modification consists of using the solution from the previous timestep as the starting point for the calculation of the solution for the next timestep. Fire-fire interactions is greatly improved by the change but a new source of error is introduced. Due to how plumes are modelled in QUIC-Fire the new solution contains errors where gaps in the plume structure are present. However, these errors are mostly limited to the upper atmosphere, where they do not affect fire behavior at the surface, and their magnitude isn’t significant enough to discount the amount of new fire phenomenology now captured in QUIC-Fire with the change.

58 GEOSCIENCES↗

NASA NTRS1999

Use of On-Line Tracers as a Diagnostic Tool in General Circulation Model Development: Transport Between the Troposphere and Stratosphere - 2

A key component of climate/chemistry modeling is how to handle the influx into (and egress from) the troposphere. This is especially important when considering tropospheric ozone, and its precursors (e.g., NO(x) from aircraft). A study has been conducted with various GISS models to determine the minimum requirements necessary for producing realistic troposphere-stratosphere exchange. Four on-line tracers are employed: CFC-11 and SF6 for mixing from the troposphere into the stratosphere, Rn222 for vertical mixing within the troposphere, and 14C for mixing from the stratosphere into the troposphere. Four standard models are tested, with varying vertical resolution, gravity wave drag and location of the model top, and additional subsidiary models are employed to examine specific features. The results show that proper vertical transport between the troposphere and stratosphere in the GISS models requires lifting the top of the model considerably out of the stratosphere, and including gravity wave drag in the lower stratosphere. Increased vertical resolution without these aspects does not improve troposphere-stratosphere exchange. The transport appears to be driven largely by the residual circulation within the stratosphere; associated E-P flux convergences require both realistic upward propagating energy from the troposphere, and realistic pass-through possibilities. A 23 layer version with a top at the mesopause and incorporating gravity wave drag appears to have reasonable stratospheric-tropospheric exchange, in terms of both the resulting tracer distributions and atmospheric mass fluxes.

Rind, David H.↗

NASA NTRS2002

Implementation of the NCAR Community Land Model (CLM) in the NASA/NCAR finite-volume Global Climate Model (fvGCM)

In this study, the NCAR CLM version 2.0 land-surface model was integrated into the NASA/NCAR fvGCM. The CLM was developed collaboratively by an open interagency/university group of scientists and based on well-proven physical parameterizations and numerical schemes that combine the best features of BATS, NCAR-LSM, and IAP94. The CLM design is a one-dimensional point model with 1 vegetation layer, along with sub-grid scale tiles. The features of the CLM include 10-uneven soil layers with water, ice, and temperature states in each soil layer, and five snow layers, with water flow, refreezing, compaction, and aging allowed. In addition, the CLM utilizes two-stream canopy radiative transfer, the Bonan lake model and topographic enhanced streamflow based on TOPMODEL. The DAO fvGCM uses a genuinely conservative Flux-Form Semi-Lagrangian transport algorithm along with terrain- following Lagrangian control-volume vertical coordinates. The physical parameterizations are based on the NCAR Community Atmosphere Model (CAM-2). For our purposes, the fvGCM was run at 2 deg x 2.5 deg horizontal resolution with 55 vertical levels. The 10-year climate from the fvGCM with CLM2 was intercompared with the climate from fvGCM with LSM, ECMWF and NCEP. We concluded that the incorporation of CLM2 did not significantly impact the fvGCM climate from that of LSM. The most striking difference was the warm bias in the CLM2 surface skin temperature over desert regions. We determined that the warm bias can be partially attributed to the value of the drag coefficient for the soil under the canopy, which was too small resulting in a decoupling between the ground surface and the canopy. We also discovered that the canopy interception was high compared to observations in the Amazon region. A number of experiments were then performed focused on implementing model improvements. In order to correct the warm bias, the drag coefficient for the soil under the canopy was considered a function of LAI (Leaf Area Index). Analysis of the results revealed that there was a substantial impact, and the warm and dry bias in the CLM2 was significantly reduced. For the interception scheme, the canopy throughfall was increased to allow for more infiltration of precipitation into the soil, resulting in increased low-level moisture and a decrease in the interception loss ratio (canopy evaporation to precipitation).

Radakovich, Jon D.↗

NASA NTRS2016

Tropical Waves and the Quasi-Biennial Oscillation in a 7-km Global Climate Simulation

This study investigates tropical waves and their role in driving a quasi-biennial oscillation (QBO)-like signal in stratospheric winds in a global 7-km-horizontal-resolution atmospheric general circulation model. The Nature Run (NR) is a 2-year global mesoscale simulation of the Goddard Earth Observing System Model, version 5 (GEOS-5). In the tropics, there is evidence that the NR supports a broad range of convectively generated waves. The NR precipitation spectrum resembles the observed spectrum in many aspects, including the preference for westward-propagating waves. However, even with very high horizontal resolution and a healthy population of resolved waves, the zonal force provided by the resolved waves is still too low in the QBO region and parameterized gravity wave drag is the main driver of the NR QBO-like oscillation (NRQBO). The authors suggest that causes include coarse vertical resolution and excessive dissipation. Nevertheless, the very-high-resolution NR provides an opportunity to analyze the resolved wave forcing of the NR-QBO. In agreement with previous studies, large-scale Kelvin and small-scale waves contribute to the NRQBO driving in eastward shear zones and small-scale waves dominate the NR-QBO driving in westward shear zones. Waves with zonal wavelength,1000 km account for up to half of the small-scale (,3300 km) resolved wave forcing in eastward shear zones and up to 70% of the small-scale resolved wave forcing in westward shear zones of the NR-QBO.

NR-QBO↗