Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Spreading models”

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 217 records · Page 12

Ignition and Wind Effects on the Transition to Flame Spread in a Microgravity Environment

The fundamental processes involved in ignition and flame spread have been extensively studied, they are generally studied separately without combining ignition and flame spread through the transition process. Moreover, the majority of the flame spread studies assume steady state flame spread. To study the transient aspects of ignition and the transition to flame spread, a time-dependent numerical model has been developed in which a thin strip is ignited along its width with either a pilot wire, radiative source, or both. For this configuration, the model is two-dimensional. Ignition is initiated away from either end of the sample, and two different flame fronts spread in opposite directions; one flows with and the other against a prescribed external flow. Usually, the flame spread is initiated by ignition at one end of the sample with or against a slow external flow, yielding one flame front. The present configuration is more realistic than the one flame front case because there is interaction between the two flames during the ignition and transition stages.

McGrattan, Kevin B.↗

Stagnation Flow in Streamer Boundaries

Flow tubes adjacent to closed magnetic field lines on the boundaries of streamers can have spreading factors which change rapidly with height. Numerical models in this thin layer are subject to uncertainties. Here we use an analytic model of magnetically closed and adjacent open regions to compute the spreading factor close to the closed field lines and the flow on the open field lines. We axe able to recover the shape of the closed field region, or helmet, and show the evolution of the helmet under slowly increasing temperature. We also show why and when stagnation flow can occur in the vicinity of the cusp at the top of the helmet.

Suess, Steven T.↗

3D Material Response of the MSL Heatshield Using NuSil-Coated PICA

The Mars Science Laboratory (MSL) was protected during its atmospheric entry by an instrumented heatshield that used NASA's Phenolic Impregnated Carbon Ablator (PICA) material [1]. PICA is a lightweight carbon fiber/polymeric resin material that offers outstanding performance for protecting probes during planetary entry. Data from the Mars Entry Descent and Landing Instrument (MEDLI) suite on MSL offers unique in-flight validation data for models of material response and atmospheric entry. MEDLI recorded, among other things, time-resolved in-depth temperature data of PICA using thermocouple sensors assembled in the MEDLI Integrated Sensor Plugs (MISP) [2]. A space-grade silicone-based coating commercially known as NuSil CV-1144-0 [3] was applied to the entire MSL heatshield, including the MEDLI plugs, to mitigate the spread of dust from PICA. Modeling the thermal response of PICA-NuSil (PICA-N) system is still an open challenge. Ground testing of PICA-N models exhibited surface temperature jumps of the order of 150 K due to oxide scale formation and sub-sequent NuSil burn-off. It is therefore critical to include a validated model for the material response of the coating in engineering codes. A test campaign has been conducted at the NASA’s Langley HyMETS [4] facility to screen the response of PICA-N and gather detailed data on its behavior [5]. A first model of PICA-N thermal response has been developed using the Hy-METS experiments [6]. The objective of this work is to analyze the material response of the latest PICA-N model compared to the engineering model used to simulate the entry of MSL. The environment and material response around the MSL aeroshell during Mars atmospheric entry is simulated using a collection of tools. The Direct Simulation Monte Carlo SPARTA code [7] is used in the rarefied regime, the Data Parallel Line Relaxation (DPLR) code [8] is used in the continuum regime and radiative heating conditions are provided by the Nonequilibrium air radiation (NEQAIR) code [9] to estimate the environmental conditions. The thermal response inside the material is computed using the Porous material Analysis Toolbox based on Open-FOAM (PATO) [10,11,12]. Thermodynamic and chemistry properties are estimated using the Mutation++ library [13]. The approach implemented in PATO as a first cut PICA-N thermal response model is outlined in Figure 1. While the recession is less than the coating thickness, the Surface mass and energy balance Boundary Condition (SBC) uses the NuSil B’ tables. Once the recession removes the coating, the usual PICA B’ tables are used for the SBC. The B’ tables are computed using an equilibrium solver implemented in Mutation++, given the temperature, pressure, blowing rate, composition of the pyrolysis and environment gases, and the condensed species at the surface. Preliminary results of the 3D material response of the MSL heat-shield at the peak heating (80 sec after Entry Interface) are shown in Figure 2. Current NASA’s mission to Mars, Mars 2020, used the spare heatshield of MSL for thermal protection during entry, descent, and landing. In preparation for Mars 2020 post-flight analysis, the PATO high-fidelity material response capability was benchmarked against flight data from MEDLI. This effort represents an important milestone toward the development of validated predictive capabilities for designing thermal protection systems for planetary probes. This bench-marking is awaiting the final release of the MEDLI-2 data.

Aerospace↗

Quantifying the Effect of Mixing on the Mean Age of Air in CCMVal-2 and CCMI-1 Models

The stratospheric age of air (AoA) is a useful measure of the overall capabilities of a general circulation model (GCM) to simulate stratospheric transport. Previous studies have reported a large spread in the simulation of AoA by GCMs and coupled chemistry-climate models (CCMs). Compared to observational estimates, simulated AoA is mostly too low. Here we attempt to untangle the processes that lead to the AoA differences between the models and between models and observations. AoA is influenced by both mean transport by the residual circulation and two-way mixing; we quantify the effects of these processes using data from the CCM inter-comparison projects CCMVal-2 (Chemistry-Climate Model Validation Activity 2) and CCMI-1 (Chemistry-Climate Model Initiative, phase 1). Transport along the residual circulation is measured by the residual circulation transit time (RCTT). We interpret the difference between AoA and RCTT as additional aging by mixing. Aging by mixing thus includes mixing on both the resolved and subgrid scale. We find that the spread in AoA between the models is primarily caused by differences in the effects of mixing and only to some extent by differences in residual circulation strength. These effects are quantified by the mixing efficiency, a measure of the relative increase in AoA by mixing. The mixing efficiency varies strongly between the models from 0.24 to 1.02. We show that the mixing efficiency is not only controlled by horizontal mixing, but by vertical mixing and vertical diffusion as well. Possible causes for the differences in the models' mixing efficiencies are discussed. Differences in subgrid-scale mixing (including differences in advection schemes and model resolutions) likely contribute to the differences in mixing efficiency. However, differences in the relative contribution of resolved versus parameterized wave forcing do not appear to be related to differences in mixing efficiency or AoA.

CCMI↗

Transition Heat Transfer Modeling Based on the Characteristics of Turbulent Spots

While turbulence models are being developed which show promise for simulating the transition region on a turbine blade or vane, it is believed that the best approach with the greatest potential for practical use is the use of models which incorporate the physics of turbulent spots present in the transition region. This type of modeling results in the prediction of transition region intermittency which when incorporated in turbulence models give a good to excellent prediction of the transition region heat transfer. Some models are presented which show how turbulent spot characteristics and behavior can be employed to predict the effect of pressure gradient and Mach number on the transition region. The models predict the spot formation rate which is needed, in addition to the transition onset location, in the Narasimha concentrated breakdown intermittency equation. A simplified approach is taken for modeling turbulent spot growth and interaction in the transition region which utilizes the turbulent spot variables governing transition length and spot generation rate. The models are expressed in terms of spot spreading angle, dimensionless spot velocity, dimensionless spot area, disturbance frequency and Mach number. The models are used in conjunction with a computer code to predict the effects of pressure gradient and Mach number on the transition region and compared with VKI experimental turbine data.

Simon, Fred↗

Are Solar Activity Variations Amplified by the QBO: A Modeling Study

Solar cycle activity effects (SCAE) in the lower and middle atmosphere, reported in several studies, are difficult to explain on the basis of the small changes in solar radiation that accompany the 11-year cycle. It is therefore natural to speculate that dynamical processes may come into play to produce a leverage. Such a leverage may be provided by the Quasi-Biennial Oscillation (QBO) in the zonal circulation of the stratosphere, which has been linked to solar activity variations. Driven primarily by wave mean flow interaction, the QBO period and its amplitude are variable but are also strongly influenced by the seasonal cycle in the solar radiation. This influence extends to low altitudes and is referred to as 'downward control'. Small changes in the solar radiative forcing may produce small changes in the period and phase of the QBO, but these in turn may produce measurable differences in the wind field. Thus, the QBO may be an amplifier of solar activity variations and a natural conduit of these variations to lower altitudes. To test this hypothesis, we conducted experiments with a 2D version of our Numerical Spectral Model that incorporates Hines' Doppler Spread Parameterization for small-scale gravity waves (GW). Solar cycle radiance variations (SCRV) are accounted for by changing the radiative heating rate on a logarithmic scale from 0.1% at the surface to 1% at 50 km to 10% at 100 km. With and without SCRV, but with the same GW flux, we then conduct numerical experiments to evaluate the magnitude of the SCAE in the zonal circulation. The numerical results indicate that, under certain conditions, the SCAE is significant and can extend to lower altitudes where the SCRV is small. For a modeled QBO period of 30 months, we find that the seasonal cycle in the solar forcing acts as a strong pacemaker to lock up the phase and period of the QBO. The SCAE then shows up primarily as a distinct but relatively weak amplitude modulation. But with a different QBO period between 30 and 34 (or less than 30, presumably) months, the seasonal phase lock is weak. Solar flux variations in the seasonal cycle then cause variations in the QBO period and phase. These amplify the SCAE to produce relatively large variations in the wind field. The SCAE in this case extends to mid-latitudes.

Mengel, J. G.↗

Solar Cycle Variations and Equatorial Oscillations: Modeling Study

Solar cycle activity effects (SCAE) in the lower and middle atmosphere, reported in several studies, are difficult to explain on the basis of the small changes in solar radiation that accompany the 11-year cycle, It is therefore natural to speculate that dynamical processes may come into play to produce a leverage. Such a leverage may be provided by the Quasi-Biennial Oscillation (QBO) in the zonal circulation of the stratosphere, which has been linked to solar activity variations. Driven primarily by wave mean flow interaction, the QBO period and its amplitude are variable but are also strongly influenced by the seasonal cycle in the solar radiation. This influence extends to low altitudes referred to as "downward control". Relatively small changes in solar radiative forcing can produce small changes in the period and phase of the QBO, but this in turn can produce measurable differences in the wind field. Thus, the QBO may be an amplifier of solar activity variations and a natural conduit of these variations to lower altitudes. To test this hypothesis, we conducted experiments with a 2D (two-dimensional) version of our Numerical Spectral Model that incorporates Hines' Doppler Spread Parameterization for small-scale gravity waves (GW). Solar cycle radiance variations (SCRV) are accounted for by changing the radiative heating rate on a logarithmic scale from 0.1 % at the surface to 1 % at 50 km to 10% at 100 km. With and without SCRV, but with the same GW flux, we then conduct numerical experiments to evaluate the magnitude of the SCAE in the zonal circulation. The numerical results indicate that, under certain conditions, the SCAE is significant and can extend to lower altitudes where the SCRV is inconsequential. At 20-km the differences in the modeled wind velocities are as large as 5 m/s. For a modeled QBO period of 30 months, we find that the seasonal cycle in the solar forcing (through the Semi-annual Oscillation (SAO)) acts as a strong pacemaker to lockup the phase and period of the QBO. The SCAE then shows up primarily as a distinct but relatively weak amplitude modulation. But with the QBO period between 30 and 34 (or less than 30, presumably) months, the seasonal phase lock is weak. Solar flux radiance variations in the seasonal cycle then cause variations in the QBO period and phase that amplify the SCAE to produce relatively large variations in the wind field. These variations also extend to mid latitudes.

Mayr, H. G.↗

Modelling of equatorial phase and amplitude scintillations from OGO-6 and AE irregularity data

The technique of modelling equatorial phase and amplitude scintillations from in-situ measurements of irregularity amplitude, ambient ionization density and irregularity power spectrum is outlined. The in-situ measurements do not directly provide information on irregularity layer thickness for which radar backscatter observations are utilized. The equatorial scintillation models developed from OGO-6 and AE-C in-situ measurements during the two solstices are presented. Considerable differences in the longitude variation of scintillations are predicted during the two solstices. The model estimates are in good agreement with the available ground based phase and amplitude measurements. Problems associated with the use of bottomside spread-F data for transionospheric propagation modelling at VHF/UHF are also discussed.

Basu, S.↗

Diurnal Forcing of Planetary Atmospheres

Much progress has been made on calculations of the Martian seasonal water cycle using the Mars Climate Model developed for this purpose. Two papers, documenting the model and the water transport results obtained with it have been published in the Journal of Geophysical Research - Planets. An additional paper describing results related to the evolution of the seasonal water cycle as a result of orbital changes was published in Advances in Space Research. Since that time, further studies have concentrated on the consequences of the soil adsorption required to match the observed water cycle and its relation to the stability of ground ice and other potential water reservoirs. Earth-related studies have concentrated on incorporating an efficient and realistic microphysical model into the Ames Stratospheric General Circulation Model used to simulate the spread of the ML Pinatubo and other volcanic clouds in the stratosphere. In addition, visualizations of the simulations are being incorporated into a video describing the UARS mission. A paper describing the new stratospheric aerosol microphysics package (and its consequences for volcanic cloud evolution) will be submitted in the near future. The paper will discuss the relative importance of condensation and coagulation to early particle growth and the separation of the cloud by sedimentation of the larger particles. A more general paper which highlights the observation that particle number densities did not increase dramatically after the ML Pinatubo eruption is planned. Simulations of atmospheric transport will be extended to include studies of terrestrial tropospheric tracers using the Fifth-Generation Penn State/NCAR Mesoscale Model.

Houben, Howard C.↗

TPSAS-NF1676L-29139-DND

The Arctic has been warming at a rate outpacing globally-averaged warming by 2-3 times, a phenomenon known as "Arctic Amplification" (AA). AA is evident in observations of Arctic surface temperature over the last century as well as in model predictions. The accuracy of model predictions is crucial to determine our response to climate change, as Arctic changes have consequences for the atmospheric circulation, sea level rise, and the carbon cycle. Unfortunately, models show a wide range of possible futures and intermodel spread in projections of Arctic temperature is larger than for any other region: a 2°C global temperature increase- the limit set by the Paris Climate Agreement- would lead to an Arctic temperature change between 3-7°C according to Coupled Model Intercomparison 5 (CMIP5) models. This intermodel spread represents a key deficiency in the understanding of Arctic feedback processes and how changes in various physical parameters (e.g. sea ice cover, clouds, heat transport and surface fluxes) affect the temperature response. Using a process-oriented decomposition of the simulate Arctic surface energy budget, we attribute the largest contribution to the CMIP5 intermodel spead in projected Arctic warming to the Barents-Kara Sea region. Our results reveal that the increasing seasonal fluxing of ocean heat content is a significant part of the intermodel spread in AA and relates strongly to change in surface turbulent fluxes. Our results indicate that if we are to significantly reduce the intermodel spread in projected AA we must focus on understanding the behavior and physical processes that drive surface energy flows in the B-K sea region.

Patrick Taylor↗

Shear flow beneath oceanic plates - Local nonsimilarity boundary layers for olivine rheology

The principle of local similarity, which has been used to model the two-dimensional boundary layers in the oceanic upper mantle, permits calculation of the temperature, velocity, and stress fields with essentially analytic techniques. Finite difference numerical methods are hard pressed to resolve the detail required by the large variation of viscosity between the lithosphere and the asthenosphere. In this paper the local similarity approximation has been justified by quantitatively evaluating the effect of nonsimilarity due to viscous heating, nonlinear temperature- and pressure-dependent rheology, buoyancy, adiabatic cooling, etc. Nonsimilar effects produce only small modifications of the locally similar boundary layers; important geophysical observables such as surface heat flux and ocean floor topography are given to better than 10 percent by the locally similar solution. A posteriori evaluations of the terms neglected in the boundary layer simplification of the complete equations have been conducted on the locally similar temperature and velocity profiles close to the spreading ridge. The boundary layer models are valid to depths of 100 km at 3 m.y. and 10 km at 0.3 m.y.

Yuen, D. A.↗

Better Constraining Supercooled Clouds Could Reduce Projected Warming Spread

The increase of climate sensitivity to rising greenhouse gas concentrations in the coupled model intercomparison project phase 6 (CMIP6) Earth system models (ESMs) compared to CMIP5 ESMs is primarily attributed to a larger extratropical cloud response to climate change, referred to as cloud feedback. The ratio of supercooled liquid cloud water relative to all cloud water, termed liquid phase ratio (LPR), which has also notably increased in many recent ESMs, is thought to be a primary driver of the extratropical cloud feedback increase. Unlike the preponderance of previous studies that compare native model LPR directly with observations, here we evaluate LPR over three ESM generations against Cloud–Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) using an instrument simulator approach, which mimics instrument limitations and uses consistent cloud definitions and resolutions. We find that current coupled model intercomparison project (CMIP) ESMs collectively simulate greater LPR than previous CMIP generations and overestimate LPR compared to observations, contrary to previous findings, likely driven by past inconsistent comparisons of ESM outputs with CALIPSO observations. We further show that greater LPR in ESM present-day climate is unexpectedly correlated with a smaller extratropical cloud feedback, attributable to a decrease of cloud amount exceeding the increase of cloud optical depth. Finally, our results suggest that improving constraints on model LPR using our evaluation framework would likely reduce the spread in CMIP6 climate sensitivities, owing to its effect on extratropical cloud feedback from supercooled clouds.

greenhouse gas concentrations↗

Premature Deaths in Africa Due to Particulate Matter Under High and Low Warming Scenarios

Sustainable development and climate change mitigation can provide enormous public health benefits via improved air quality, especially in polluted areas. We use the latest state-of-the-art composition-climate model simulations to contrast human exposure to fine particulate matter in Africa under a “baseline” scenario with high material consumption, population growth, and warming to that projected under a sustainability scenario with lower consumption, population growth, and warming. Evaluating the mortality impacts of these exposures, we find that under the low warming scenario annual premature deaths due to PM2.5 are reduced by roughly 515,000 by 2050 relative to the high warming scenario (100,000, 175,000, 55,000, 140,000, and 45,000 in Northern, West, Central, East, and Southern Africa, respectively). This reduction rises to ∼800,000 by the 2090s, though by that time much of the difference is attributable to the projected differences in population. By contrast, during the first half of the century benefits are driven predominantly by emissions changes. Depending on the region, we find large intermodel spreads of ∼25%–50% in projected future exposures owing to different physics across the ensemble of 6 global models. The spread of projected deaths attributable to exposure to fine particulate matter, including uncertainty in the exposure-response function, are reduced in every region to ∼20%–35% by the non-linear exposure-response function. Differences between the scenarios have an even narrower spread of ∼5%–25% and are highly statistically significant in all regions for all models. These results provide valuable information for policy-makers to consider when working toward climate change mitigation and sustainable development goals.

PM2.5↗

Shock Formation of Slow Magnetosonic Waves in Coronal Plumes

We investigate the height of shock formation in coroner plumes for slow magnetosonic waves. The models take into account plume geometric spreading, heat conduction and radiative damping. The wave parameters as well as the spreading functions of the plumes and the base magnetic field strengths are given by empirical constraints mostly from Solar and Heliospheric Observatory/Ultraviolet Coronagraph Spectrometer (SOHO/UVCS). Our models show that shock formation occurs at low coronal heights, i.e., within 1.3 solar radius, depending on the model parameters. The shock formation is calculated using the well-established wave breaking condition given by the intersection of C+ characteristics in the space-time plane. Our models show that shock heating by slow magnetosonic waves is expected to be relevant at most heights in solar coronal plumes, although slow magnetosonic waves are most likely not a solely operating energy supply mechanism.

Cuntz, Manfred↗

Solid surface combustion experiment flame spread in a quiescent, microgravity environment implications of spread rate and flame structure

A unique environment in which flame spreading, a phenomenon of fundamental, scientific interest, has importance to fire safety is that of spacecraft in which the gravitational acceleration is low compared with that of the Earth, i.e., microgravity. Experiments aboard eight Space Shuttle missions between October 1990 and February 1995 were conducted using the Solid Surface Combustion Experiment (SSCE) payload apparatus in an effort to determine the mechanisms of gas-phase flame spread over solid fuel surfaces in the absence of any buoyancy induced or externally imposed oxidizer flow. The overall SSCE effort began in December of 1984. The SSCE apparatus consists of a sealed container, approximately 0.039 cu m, that is filled with a specified O2/N2 mixture at a prescribed pressure. Five of the experiments used a thin cellulosic fuel, ashless filter paper, 3 cm wide x 10 cm long, 0.00825 cm half-thickness, ignited in five different ambient conditions. Three of the experiments, the most recent, used thick polymethylmethacrylate (PMMA) samples 0.635 cm wide x 2 cm long, 0.32 cm half-thickness. Three experiments, STS 41, 40 and 43, were designed to evaluate the effect of ambient pressure on flame spread over the thin cellulosic fuel while flights STS 50 and 47 were at the same pressure as two of the earlier flights but at a lower oxygen concentration in order to evaluate the effect of ambient oxygen level on the flame spread process at microgravity. For the PMMA flights, two experiments, STS 54 and 63, were at the same pressure but different oxygen concentrations while STS 64 was at the same oxygen concentration as STS 63 but at a higher pressure. Two orthogonal views of the experiments were recorded on 16 mm cine-cameras operating at 24 frames/s. In addition to filmed images of the side view of the flames and surface view of the burning samples, solid- and gas-phase temperatures were recorded using thermocouples. The experiment is battery powered and follows an automated sequence upon activation by the Shuttle Crew. In this study we separate the SSCE data into two groups according to the fuel type: (1) thin cellulose; and (2) thick PMMA. The experimental spread rates are compared with prediction from a number of models in an effort to uncover the important physics that characterize microgravity flame spread. Both steady and unsteady solutions are employed to explore the flame evolution, especially for thick fuels. Finally, the flame structure in downward spread is compared with the microgravity flame structure and modeling results to delineate the difference between the two configurations and the influence of normal gravity.

Bundy, Matthew↗

Modeling and Prediction of Wildfire Hazard in Southern California, Integration of Models with Imaging Spectrometry

Large urban wildfires throughout southern California have caused billions of dollars of damage and significant loss of life over the last few decades. Rapid urban growth along the wildland interface, high fuel loads and a potential increase in the frequency of large fires due to climatic change suggest that the problem will worsen in the future. Improved fire spread prediction and reduced uncertainty in assessing fire hazard would be significant, both economically and socially. Current problems in the modeling of fire spread include the role of plant community differences, spatial heterogeneity in fuels and spatio-temporal changes in fuels. In this research, we evaluated the potential of Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) and Airborne Synthetic Aperture Radar (AIRSAR) data for providing improved maps of wildfire fuel properties. Analysis concentrated in two areas of Southern California, the Santa Monica Mountains and Santa Barbara Front Range. Wildfire fuel information can be divided into four basic categories: fuel type, fuel load (live green and woody biomass), fuel moisture and fuel condition (live vs senesced fuels). To map fuel type, AVIRIS data were used to map vegetation species using Multiple Endmember Spectral Mixture Analysis (MESMA) and Binary Decision Trees. Green live biomass and canopy moisture were mapped using AVIRIS through analysis of the 980 nm liquid water absorption feature and compared to alternate measures of moisture and field measurements. Woody biomass was mapped using L and P band cross polarimetric data acquired in 1998 and 1999. Fuel condition was mapped using spectral mixture analysis to map green vegetation (green leaves), nonphotosynthetic vegetation (NPV; stems, wood and litter), shade and soil. Summaries describing the potential of hyperspectral and SAR data for fuel mapping are provided by Roberts et al. and Dennison et al. To utilize remotely sensed data to assess fire hazard, fuel-type maps were translated into standard fuel models accessible to the FARSITE fire spread simulator. The FARSITE model and BEHAVE are considered industry standards for fire behavior analysis. Anderson level fuels map, generated using a binary decision tree classifier are available for multiple dates in the Santa Monica Mountains and at least one date for Santa Barbara. Fuel maps that will fill in the areas between Santa Barbara and the Santa Monica Mountains study sites are in progress, as part of a NASA Regional Earth Science Application Center, the Southern California Wildfire Hazard Center. Species-level maps, were supplied to fire managing agencies (Los Angeles County Fire, California Department of Forestry). Research results were published extensively in the refereed and non-refereed literature. Educational outreach included funding of several graduate students, undergraduate intern training and an article featured in the California Alliance for Minorities Program (CAMP) Quarterly Journal.

Roberts, Dar A.↗

Stagnation Flow in Streamer Boundaries

It is possible for flow tubes adjacent to closed magnetic field lines on the boundaries of streamers to have spreading factors that change rapidly according to height. Mathematical models in this thin layer are prone to uncertainties. In this paper we use an analytical model of magnetically closed and adjacent open regions to calculate the spreading factor close to the closed field lines and near the flow on the open field lines. We recovered the shape of the closed field region, or helmet, and show the evolution of the helmet under increasing temperature. We also show why and when stagnation flow can occur in the area of the cusp at the top of the helmet.

Suess, S. T.↗