Search NASA⌕ Search

SEARCH · Search NASA

Results for “time-scale”

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 253 records · Page 14

Thermal-Acoustic Analysis of a Metallic Integrated Thermal Protection System Structure

A study is undertaken to investigate the response of a representative integrated thermal protection system structure under combined thermal, aerodynamic pressure, and acoustic loadings. A two-step procedure is offered and consists of a heat transfer analysis followed by a nonlinear dynamic analysis under a combined loading environment. Both analyses are carried out in physical degrees-of-freedom using implicit and explicit solution techniques available in the Abaqus commercial finite-element code. The initial study is conducted on a reduced-size structure to keep the computational effort contained while validating the procedure and exploring the effects of individual loadings. An analysis of a full size integrated thermal protection system structure, which is of ultimate interest, is subsequently presented. The procedure is demonstrated to be a viable approach for analysis of spacecraft and hypersonic vehicle structures under a typical mission cycle with combined loadings characterized by largely different time-scales.

Behnke, Marlana N.↗

Quantifying VOC-Reaction Tracers, Ozone Production, and Continuing Aerosol Production Rates in Urban and Far-Downwind Atmospheres

We have found a surprisingly informative decomposition of the complex question of smoggy ozone production (basically, [HO2] in a more locally determined field of [NO]) in the process of linked investigations of modestly smoggy Eastern North America (by NASA aircraft, July 2004) and rather polluted Flushing, NYC (Queens College, July, 2001). In both rural and very polluted situations, we find that a simple contour graph parameterization of the local principal ozone production rate can be estimated using only the variables [NO] and j(sub rads) [HCHO]: Po(O3) = c (j(sub rads) [HCHO])(sup a) [HCHO](sup b). Here j(sub rads) is the photolysis of HCHO to radicals, presumably capturing many harder-UV photolytic processes and the principle ozone production is that due to HO2; mechanisms suggest that ozone production due to RO2 is closely correlated, often suggesting a limited range of different proportionality factors. The method immediately suggests a local interpretation for concepts of VOC limitation and NOx limitation. We believe that the product j(sub rads) [HCHO] guages the oxidation rate of observed VOC mixtures in a way that also provides [HO2] useful for the principle ozone production rate k [HO2] [NO], and indeed, all ozone chemical production. The success of the method suggests that dominant urban primary-HCHO sources may transition to secondary plume-HCHO sources in a convenient way. Are there other, simple, near-terminal oxidized VOC's which help guage ozone production and aerosol particle formation? Regarding particles, we report on, to the extent NASA Research resources allow, on appealing relationships between far-downwind (Atlantic PBL) HCHO and very fine aerosol (including sulfate. Since j(sub rads) [HCHO] provides a time-scale, we may understand distant-plume particle production in a more quantitative manner. Additionally we report on a statistical search in the nearer field for relationships between glyoxals (important near-terminal aromatic and isoprene reaction products) and aerosol production, looking for VOC's that might be most implicated. All three variables j(sub rads), [HCHO], and [NO] are relatively easily measured in widespread air pollution monitoring networks, and all are deducible form space-borne observations, though estimation of [NO] from [NO2] (the species observable from space) may require care. We report also on airborne and surface observations of HCHO, suggesting that concentrated (urban) and more diffuse (forest) sources may be distinguishable from space. The use of the 3.58 micron microwindow for HCHO remote sensing should allow much sharper resolution of HCHO than the UV. UV sensing requires large and expensive instruments, but even these seem justified since formaldehyde is so informative.

Chatfield, Robert↗

Fast X-Ray Timing: A Window into the Strong-Field Regime

The dynamical time-scales in the vicinity of neutron star surfaces and black hole horizons are in the millisecond range. Over the past decade, timing signatures on such scales, either quasi-periodic oscillations (QPOs) and/or coherent pulsations, have been discovered and studied from both neutron stars and black holes with NASA's Rossi X-ray Timing Explorer, Although theoretical interpretations are still hotly debated, these timing properties almost certainly reflect the dynamics of matter in regions dominated by relativistic gravity. I will survey our current understanding of these timing properties, with a focus on how they might he used as probes of fundamental physics.

Strohmayer, Tod↗

A Singular Perturbation Approach for Time-Domain Assessment of Phase Margin

This paper considers the problem of time-domain assessment of the Phase Margin (PM) of a Single Input Single Output (SISO) Linear Time-Invariant (LTI) system using a singular perturbation approach, where a SISO LTI fast loop system, whose phase lag increases monotonically with frequency, is introduced into the loop as a singular perturbation with a singular perturbation (time-scale separation) parameter Epsilon. First, a bijective relationship between the Singular Perturbation Margin (SPM) max and the PM of the nominal (slow) system is established with an approximation error on the order of Epsilon(exp 2). In proving this result, relationships between the singular perturbation parameter Epsilon, PM of the perturbed system, PM and SPM of the nominal system, and the (monotonically increasing) phase of the fast system are also revealed. These results make it possible to assess the PM of the nominal system in the time-domain for SISO LTI systems using the SPM with a standardized testing system called "PM-gauge," as demonstrated by examples. PM is a widely used stability margin for LTI control system design and certification. Unfortunately, it is not applicable to Linear Time-Varying (LTV) and Nonlinear Time-Varying (NLTV) systems. The approach developed here can be used to establish a theoretical as well as practical metric of stability margin for LTV and NLTV systems using a standardized SPM that is backward compatible with PM.

Zhu, J. Jim↗

Coherent Eigenmodes in Homogeneous MHD Turbulence

The statistical mechanics of Fourier models of ideal, homogeneous, incompressible magnetohydrodynamic (MHD) turbulence is discussed, along with their relevance for dissipative magnetofluids. Although statistical theory predicts that Fourier coefficients of fluid velocity and magnetic field are zero-mean random variables, numerical simulations clearly show that certain coefficients have a non-zero mean value that can be very large compared to the associated standard deviation, i.e., we have coherent structure. We use eigenanalysis of the modal covariance matrices in the probability density function to explain this phenomena in terms of `broken ergodicity', which is defined to occur when dynamical behavior does not match ensemble predictions on very long time-scales. We provide examples from 2-D and 3-D magnetohydrodynamic simulations of homogeneous turbulence, and show new results from long-time simulations of MHD turbulence with and without a mean magnetic field

Shebalin, John V.↗

Broken Ergodicity in Ideal, Homogeneous, Incompressible Turbulence

We discuss the statistical mechanics of numerical models of ideal homogeneous, incompressible turbulence and their relevance for dissipative fluids and magnetofluids. These numerical models are based on Fourier series and the relevant statistical theory predicts that Fourier coefficients of fluid velocity and magnetic fields (if present) are zero-mean random variables. However, numerical simulations clearly show that certain coefficients have a non-zero mean value that can be very large compared to the associated standard deviation. We explain this phenomena in terms of broken ergodicity', which is defined to occur when dynamical behavior does not match ensemble predictions on very long time-scales. We review the theoretical basis of broken ergodicity, apply it to 2-D and 3-D fluid and magnetohydrodynamic simulations of homogeneous turbulence, and show new results from simulations using GPU (graphical processing unit) computers.

Morin, Lee↗

Spatially and Temporally-Resolved Multi-Parameter Interferometric Rayleigh Scattering

A novel approach to simultaneously measure the translational temperature, bulk velocity, and density in gases by collecting, referencing, and analyzing nanosecond time-scale Rayleigh scattered light from molecules is described. A narrow-band pulsed laser source is used to probe two largely separated measurement locations, one of which is used for reference. The elastically scattered photons containing information from both measurement locations are collected at the same time and analyzed spectrally using a planar Fabry - Perot interferometer. A practical means of referencing the measurement of velocity using the laser frequency, and the density and temperature using the information from the reference measurement location maintained at constant properties is described. To demonstrate the technique single-shot spectra of elastic scattered light are obtained in a near zero velocity H2-air Hencken burner flame and simultaneously in an N2-filled gas cell. A simplified Gaussian distribution model to the scattered light spectra is used to obtain the flame properties. Corrections to this model are applied at lower gas temperatures when the simplified Gaussian approximation is no longer suitable. The near-zero measured velocity as a function of the measured flame temperature, and a comparison of the measured flame density and temperature with the perfect gas law are presented.

Bivolaru, Daniel↗

Analysis of Ice Plains of Filchner/Ronne Ice Shelf Using ICESat Data

We use repeat-track laser altimeter data from the Ice, Cloud, and land Elevation Satellite (ICESat) to map the grounding zone of Filchner/Ronne Ice Shelf (FRIS), Antarctica. Repeated passes of ICESat reveal ice flexure in the grounding zone occurs as the ice shelf responds to ocean height changes due primarily to tides. In the course of our mapping, we have confirmed or identified three major "ice plains", regions of low surface slope near the GZ where the ice is close to hydrostatic equilibrium: one on Institute Ice Stream, another to its east, and another west of Foundation Ice Stream. The vertical information from repeated ICESat tracks enables us to study the topography and flexure characteristics across these three ice plains, and we use this to develop a classification scheme for ice plains based on their surface topography and their state of flotation. We show that one of these ice plains indicates changes in lateral extent on short time-scales, depending on the state of the ocean tide. Understanding the location and nature of ice plains is important for ice sheet modeling, since they add uncertainty to the absolute boundary between floating and grounded ice.

Brunt, Kelly M.↗

Adaptive Numerical Algorithms in Space Weather Modeling

Space weather describes the various processes in the Sun-Earth system that present danger to human health and technology. The goal of space weather forecasting is to provide an opportunity to mitigate these negative effects. Physics-based space weather modeling is characterized by disparate temporal and spatial scales as well as by different physics in different domains. A multi-physics system can be modeled by a software framework comprising of several components. Each component corresponds to a physics domain, and each component is represented by one or more numerical models. The publicly available Space Weather Modeling Framework (SWMF) can execute and couple together several components distributed over a parallel machine in a flexible and efficient manner. The framework also allows resolving disparate spatial and temporal scales with independent spatial and temporal discretizations in the various models. Several of the computationally most expensive domains of the framework are modeled by the Block-Adaptive Tree Solar wind Roe Upwind Scheme (BATS-R-US) code that can solve various forms of the magnetohydrodynamics (MHD) equations, including Hall, semi-relativistic, multi-species and multi-fluid MHD, anisotropic pressure, radiative transport and heat conduction. Modeling disparate scales within BATS-R-US is achieved by a block-adaptive mesh both in Cartesian and generalized coordinates. Most recently we have created a new core for BATS-R-US: the Block-Adaptive Tree Library (BATL) that provides a general toolkit for creating, load balancing and message passing in a 1, 2 or 3 dimensional block-adaptive grid. We describe the algorithms of BATL and demonstrate its efficiency and scaling properties for various problems. BATS-R-US uses several time-integration schemes to address multiple time-scales: explicit time stepping with fixed or local time steps, partially steady-state evolution, point-implicit, semi-implicit, explicit/implicit, and fully implicit numerical schemes. Depending on the application, we find that different time stepping methods are optimal. Several of the time integration schemes exploit the block-based granularity of the grid structure. The framework and the adaptive algorithms enable physics based space weather modeling and even forecasting.

Toth, Gabor↗

Modeling the Frozen-In Anticyclone in the 2005 Arctic Summer Stratosphere

Immediately following the breakup of the 2005 Arctic spring stratospheric vortex, a tropical air mass, characterized by low potential vorticity (PV) and high nitrous oxide (N2O), was advected poleward and became trapped in the easterly summer polar vortex. This feature, known as a "Frozen-In Anticyclone (FrIAC)", was observed in Earth Observing System (EOS) Aura Microwave Limb Sounder (MLS) data to span the potential temperature range from approximately 580 to 1100 K (approximately 25 to 40 km altitude) and to persist from late March to late August 2005. This study compares MLS N2O observations with simulations from the Global Modeling Initiative (GMI) chemistry and transport model, the GEOS-5/MERRA Replay model, and the VanLeer Icosahedral Triangular Advection isentropic transport model to elucidate the processes involved in the lifecycle of the FrIAC which is here divided into three distinct phases. During the "spin-up phase" (March to early April), strong poleward flow resulted in a tight isolated anticyclonic vortex at approximately 70-90 deg N, marked with elevated N2O. GMI, Replay, and VITA all reliably simulted the spin-up of the FrIAC, although the GMI and Replay peak N2O values were too low. The FrIAC became trapped in the developing summer easterly flow and circulated around the polar region during the "anticyclonic phase" (early April to the end of May). During this phase, the FrIAC crossed directly over the pole between the 7th and 14th of April. The VITA and Replay simulations transported the N2O anomaly intact during this crossing, in agreement with MLS, but unrealistic dispersion of the anomaly occurred in the GMI simulation due to excessive numerical mixing of the polar cap. The vortex associated with the FrIAC was apparently resistant to the weak vertical hear during the anticyclonic phase, and it thereby protected the embedded N20 anomaly from stretching. The vortex decayed in late May due to diabatic processes, leaving the N2O anomaly exposed to horizontal and vertical wind shears during the "shearing phase" (June to August). The observed lifetime of the FrIAC during this phase is consistent with time-scales calculated from the ambient horizontal and vertical wind shear. Replay maintained the horizontal structure of the N2O anomaly similar to NILS well into August. The VITA simulation also captured the horizontal structure of the FrIAC during this phase, but VITA eventually developed fine-scale N2O structure not observed in MLS data.

Allen, D. R.↗

Constraints on CME Evolution from in situ Observations of Ionic Charge States

We present a novel procedure for deriving the physical properties of Coronal Mass Ejections (CMES) in the corona. Our methodology uses in-situ measurements of ionic charge states of C, O, Si and Fe in the heliosphere and interprets them in the context of a model for the early evolution of ICME plasma, between 2 - 5 R-solar. We find that the data can be fit only by an evolution that consists of an initial heating of the plasma, followed by an expansion that ultimately results in cooling. The heating profile is consistent with a compression of coronal plasma due to flare reconnect ion jets and an expansion cooling due to the ejection, as expected from the standard CME/flare model. The observed frozen-in ionic charge states reflect this time-history and, therefore, provide important constraints for the heating and expansion time-scales, as well as the maximum temperature the CME plasma is heated to during its eruption. Furthermore, our analysis places severe limits on the possible density of CME plasma in the corona. We discuss the implications of our results for CME models and for future analysis of ICME plasma composition.

Gruesbeck, Jacob R.↗

The Sun, Its Extended Corona, the Interplanetary Space, the Earth's Magnetosphere, Ionosphere, Middle and Low Atmosphere, are All Parts of a Complex System - the Heliosphere

Various manifestations of solar activity cause disturbances known as space weather effects in the interplanetary space, near-Earth environment, and all the Earth's "spheres. Longterm variations in the frequency, intensity and relative importance of the manifestations of solar activity are due to the slow changes in the output of the solar dynamo, and they define space climate. Space climate governs long-term variations in geomagnetic activity and is the primary natural driver of terrestrial climate. To understand how the variable solar activity affects the Earth's environment, geomagnetic activity and climate on both short and long time scales, we need to understand the origins of solar activity itself and its different manifestations, as well as the sequence of coupling processes linking various parts of the system. This session provides a forum to discuss the chain of processes and relations from the Sun to the Earth's surface: the origin and long-term and short-term evolution of solar activity, initiation and temporal variations in solar flares, CMEs, coronal holes, the solar wind and its interaction with the terrestrial magnetosphere, the ionosphere and its connection to the neutral dominated regions below and the plasma dominated regions above, the stratosphere, its variations due to the changing solar activity and its interactions with the underlying troposphere, and the mechanisms of solar influences on the lower atmosphere on different time-scales. Particularly welcome are papers highlighting the coupling processes between the different domains in this complex system.

Gopalswamy, Natchimuthuk↗

Searching for Soft Relativistic Jets in Core-Collapse Supernovae with the IceCube Optical Follow-up Program

Context. Transient neutrino sources such as Gamma-Ray Bursts (GRBs) and Supernovae (SNe) are hypothesized to emit bursts of high-energy neutrinos on a time-scale of < or approx.100 s. While GRB neutrinos would be produced in high relativistic jets, core-collapse SNe might host soft-relativistic jets, which become stalled in the outer layers of the progenitor star leading to an efficient production of high-energy neutrinos. Aims. To increase the sensitivity to these neutrinos and identify their sources, a low-threshold optical follow-up program for neutrino multiplets detected with the IceCube observatory has been implemented. Methods. If a neutrino multiplet, i.e. two or more neutrinos from the same direction within 100 s, is found by IceCube a trigger is sent to the Robotic Optical Transient Search Experiment, ROTSE. The 4 ROTSE telescopes immediately start an observation program of the corresponding region of the sky in order to detect an optical counterpart to the neutrino events. Results. No statistically significant excess in the rate of neutrino multiplets has been observed and furthermore no coincidence with an optical counterpart was found. Conclusions. The search allows, for the first time, to set stringent limits on current models predicting a high-energy neutrino flux from soft relativistic hadronic jets in core-collapse SNe. We conclude that a sub-population of SNe with typical Lorentz boost factor and jet energy of 10 and 3 x 10(exp 51) erg, respectively, does not exceed 4:2% at 90% confidence.

Abbasi, R.↗

New Insights on the Accretion Disk-Winds Connection in Radio-Loud AGNs from Suzaku

From the spectral analysis of long Suzaku observations of five radio-loud AGNs we have been able to discover the presence of ultra-fast outflows with velocities ,,approx.0.1 c in three of them, namely 3C III, 3C 120 and 3C 390.3. They are consistent with being accretion disk winds/outflows. We also performed a follow-up on 3C III to monitor its outflow on approx.7 days time-scales and detected an anti-correlated variability of a possible relativistic emission line with respect to blue-shifted Fe K features, following a flux increase. This provides the first direct evidence for an accretion disc-wind connection in an AGN. The mass outflow rate of these outflows can be comparable to the accretion rate and their mechanical power can correspond to a significant fraction of the bolometric luminosity and is comparable to their typical jet power. Therefore, they can possibly play a significant role in the expected feedback from AGNs and can give us further clues on the relation between the accretion disk and the formation of winds/jets.

Tombesi, F.↗

Variability of Jupiter's Five-Micron Hot Spot Inventory

Global upheavals on Jupiter involve changes in the albedo of entire axisymmetric regions, lasting several years, with the last two occurring in 1989 and 2006. Against this backdrop of planetary-scale changes, discrete features such as the Great Red Spot (GRS), and other vortices exhibit changes on shorter spatial- and time-scales. We track the variability of the discrete equatorial 5-micron hot spots, semi-evenly spaced in longitude and confined to a narrow latitude band centered at 6.5degN (southern edge of the North Equatorial Belt, NEB), abundant in Voyager images. Tantalizingly similar patterns were observed in the visible (bright plumes and blue-gray regions), where reflectivity in the red is anti-correlated with 5-microns thermal radiance. Ortiz et al. (1998, GRL, 103) characterized the latitude and drift rates of the hot spots, including the descent of the Galileo probe at the southern edge of a 5-micron hot spot, as the superposition of equatorial Rossby waves, with phase speeds between 99 - 103m/s, relative to System III. We note that the high 5-micron radiances correlate well but not perfectly with high 8.57-micron radiances. Because the latter are modulated primarily by changes in the upper ammonia (NH3) ice cloud opacity, this correlation implies that changes in the ammonia ice cloud field may be responsible for the variability seen in the 5-~m maps. During the NEB fade (2011 - early 2012), however, these otherwise ubiquitous features were absent, an atmospheric state not seen in decades. The ongoing NEB revival indicates nascent 5-~m hot spots as early as April 2012, with corresponding visible dark spots. Their continuing growth through July 2012 indicates the possit.le re-establishment of Rossby waves. The South Equatorial Belt (SEB) and NEB revivals began similarly with an instability that developed into a major outbreak, and many similarities in the observed propagation of clear regions.

Yanamandra-Fisher, Padma A.↗

Comparative Analysis of Upper Ocean Heat Content Variability from Ensemble Operational Ocean Analyses

Upper ocean heat content (HC) is one of the key indicators of climate variability on many time-scales extending from seasonal to interannual to long-term climate trends. For example, HC in the tropical Pacific provides information on thermocline anomalies that is critical for the longlead forecast skill of ENSO. Since HC variability is also associated with SST variability, a better understanding and monitoring of HC variability can help us understand and forecast SST variability associated with ENSO and other modes such as Indian Ocean Dipole (IOD), Pacific Decadal Oscillation (PDO), Tropical Atlantic Variability (TAV) and Atlantic Multidecadal Oscillation (AMO). An accurate ocean initialization of HC anomalies in coupled climate models could also contribute to skill in decadal climate prediction. Errors, and/or uncertainties, in the estimation of HC variability can be affected by many factors including uncertainties in surface forcings, ocean model biases, and deficiencies in data assimilation schemes. Changes in observing systems can also leave an imprint on the estimated variability. The availability of multiple operational ocean analyses (ORA) that are routinely produced by operational and research centers around the world provides an opportunity to assess uncertainties in HC analyses, to help identify gaps in observing systems as they impact the quality of ORAs and therefore climate model forecasts. A comparison of ORAs also gives an opportunity to identify deficiencies in data assimilation schemes, and can be used as a basis for development of real-time multi-model ensemble HC monitoring products. The OceanObs09 Conference called for an intercomparison of ORAs and use of ORAs for global ocean monitoring. As a follow up, we intercompared HC variations from ten ORAs -- two objective analyses based on in-situ data only and eight model analyses based on ocean data assimilation systems. The mean, annual cycle, interannual variability and longterm trend of HC have been analyzed

Xue, Yan↗

X-Ray Evidence for the Accretion Disc-Outflow Connection in 3C 111

We present the spectral analysis of three Suzaku X-ray Imaging Spectrometer observations of 3C III requested to monitor the predicted variability of its ultrafast outflow on approximately 7 d time-scales. We detect an ionized iron emission line in the first observation and a blueshifted absorption line in the second, when the flux is approximately 30 per cent higher. The location of the material is constrained at less than 0.006 pc from the variability. Detailed modelling supports an identification with ionized reflection off the accretion disc at approximately 20-100rg from the black hole and a highly ionized and massive ultrafast outflow with velocity approximately 0.1c, respectively. The outflow is most probably accelerated by radiation pressure, but additional magnetic thrust cannot be excluded. The measured high outflow rate and mechanical energy support the claims that disc outflows may have a significant feedback role. This work provides the first direct evidence for an accretion disc-outflow connection in a radio-loud active galactic nucleus, possibly linked also to the jet activity.

Tombesi, Frank↗

Current Flow and Pair Creation at Low Altitude in Rotation-Powered Pulsars' Force-Free Magnetospheres: Space Charge Limited Flow

We report the results of an investigation of particle acceleration and electron-positron plasma generation at low altitude in the polar magnetic flux tubes of rotation-powered pulsars, when the stellar surface is free to emit whatever charges and currents are demanded by the force-free magnetosphere. We apply a new 1D hybrid plasma simulation code to the dynamical problem, using Particle-in-Cell methods for the dynamics of the charged particles, including a determination of the collective electrostatic fluctuations in the plasma, combined with a Monte Carlo treatment of the high-energy gamma-rays that mediate the formation of the electron-positron pairs.We assume the electric current flowing through the pair creation zone is fixed by the much higher inductance magnetosphere, and adopt the results of force-free magnetosphere models to provide the currents which must be carried by the accelerator. The models are spatially one dimensional, and designed to explore the physics, although of practical relevance to young, high-voltage pulsars. We observe novel behaviour (a) When the current density j is less than the Goldreich-Julian value (0 < j/j(sub GJ) < 1), space charge limited acceleration of the current carrying beam is mild, with the full Goldreich-Julian charge density comprising the charge densities of the beam and a cloud of electrically trapped particles with the same sign of charge as the beam. The voltage drops are of the order of mc(sup 2)/e, and pair creation is absent. (b) When the current density exceeds the Goldreich-Julian value (j/j(sub GJ) > 1), the system develops high voltage drops (TV or greater), causing emission of curvature gamma-rays and intense bursts of pair creation. The bursts exhibit limit cycle behaviour, with characteristic time-scales somewhat longer than the relativistic fly-by time over distances comparable to the polar cap diameter (microseconds). (c) In return current regions, where j/j(sub GJ) < 0, the system develops similar bursts of pair creation. These discharges are similar to those encountered in previous calculations by Timokhin of pair creation when the surface has a high work function and cannot freely emit charge. In cases (b) and (c), the intermittently generated pairs allow the system to simultaneously carry the magnetospherically prescribed currents and adjust the charge density and average electric field to force-free conditions. We also elucidate the conditions for pair creating beam flow to be steady (stationary with small fluctuations in the rotating frame), finding that such steady flows can occupy only a small fraction of the current density parameter space exhibited by the force-free magnetospheric model. The generic polar flow dynamics and pair creation are strongly time dependent. The model has an essential difference from almost all previous quantitative studies, in that we sought the accelerating voltage (with pair creation, when the voltage drops are sufficiently large; without, when they are small) as a function of the applied current.

Timokhin, A. N.↗