Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase field model”

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 109 records · Page 6

Ring Current Development During Storm Main Phase

The development of the ring current ions in the inner magnetosphere during the main phase of a magnetic storm is studied. The temporal and spatial evolution of the ion phase space densities in a dipole field are calculated using a three dimensional ring current model, considering charge exchange and Coulomb losses along drift paths. The simulation starts with a quiet time distribution. The model is tested by comparing calculated ion fluxes with Active Magnetospheric Particle Tracer Explorers/CCE measurement during the storm main phase on May 2, 1986. Most of the calculated omnidirectional fluxes are in good agreement with the data except on the dayside inner edge (L less than 2.5) of the ring current, where the ion fluxes are underestimated. The model also reproduces the measured pitch angle distributions of ions with energies below 10 keV. At higher energy, an additional diffusion in pitch angle is necessary in order to fit the data. The role of the induced electric field on the ring current dynamics is also examined by simulating a series of substorm activities represented by stretching and collapsing the magnetic field lines. In response to the impulsively changing fields, the calculated ion energy content fluctuates about a mean value that grows steadily with the enhanced quiescent field.

Fok, Mei-Ching↗

Ozone Sensitivity to Varying Greenhouse Gases and Ozone-Depleting Substances in CCMI-1 Simulations

Ozone fields simulated for the first phase of the Chemistry-Climate Model Initiative (CCMI-1) will be used as forcing data in the 6th Coupled Model Intercomparison Project. Here we assess, using reference and sensitivity simulations produced for CCMI-1, the suitability of CCMI-1 model results for this process, investigating the degree of consistency amongst models regarding their responses to variations in individual forcings. We consider the influences of methane, nitrous oxide, a combination of chlorinated or brominated ozone-depleting substances, and a combination of carbon dioxide and other greenhouse gases. We find varying degrees of consistency in the models' responses in ozone to these individual forcings, including some considerable disagreement. In particular, the response of total-column ozone to these forcings is less consistent across the multi-model ensemble than profile comparisons. We analyse how stratospheric age of air, a commonly used diagnostic of stratospheric transport, responds to the forcings. For this diagnostic we find some salient differences in model behaviour, which may explain some of the findings for ozone. The findings imply that the ozone fields derived from CCMI-1 are subject to considerable uncertainties regarding the impacts of these anthropogenic forcings. We offer some thoughts on how to best approach the problem of generating a consensus ozone database from a multi-model ensemble such as CCMI-1.

CCMI↗

Magnetic Field Configuration of Active Region NOAA 6555 at the Time of a Long Duration Flare on 23 March 1991: An Exception to Standard Flare Reconnection Model

The high-resolution H(sub alpha) images observed during the decay phase of a long duration flare on 23 March 1991 are used to study the three-dimensional magnetic field configuration of the active region NOAA 6555. Whereas, all the large flares in NOAA 6555 occurred at the location of high magnetic shear and flux emergence, this long duration flare was observed in the region of low magnetic shear at the photosphere. The H(sub alpha) loop activity started soon after the maximum phase of the flare. There were few long loop at the initial phase of the activity. Some of these were sheared in the chromosphere at an angle of about 45 deg with the east-west axis. Gradually, increasing number of shorter loops, oriented along the east-west axis, started appearing. The chromospheric Dopplergrams show blue-shifts at the end points of the loops. By using different magnetic field models, we have extrapolated the photospheric magnetograms to the chromospheric heights. The magnetic field lines computed by using the potential field model correspond to most of the observed H(sub alpha) loops. The height of the H(sub alpha) loops were derived by comparing them with the computed field lines. From the temporal evolution of the H(sub alpha) loop activity, we derive the negative rate of appearance of H(sub alpha) features as a function of height. It is found that the field lines oriented along one of the neutral lines was sheared and low lying. The higher field lines were mostly potential. The paper also outlines a possible scenario for describing the post-flare stage of the observed long duration flare.

Choudhary, Debi Prasad↗

The 'convection surge' mechanism of ion acceleration during substorms

To explain the presence of the so-called 'bounce-phase-bunched' ion distributions (eV range to about 20 keV) observed in the earth's geosynchronous magnetosphere in association with substorm expansion phases, Quinn and Southwood proposed the 'convection surge' mechanism of ion energization. This mechanism is associated with a sudden earthward displacement of curved field lines resulting from the short-lived application of an intense, east-west electric field. A quantitative computer model of the mechanism has been constructed here. It is found that the mechanism easily generates the bounce-phase-bunched ion distributions in question. Additionally, the mechanism can generate dramatically field-aligned distributions. It is proposed that a key general characteristic of geosynchronous ion pitch angle distributions, field-aligned low energies, and field-perpendicular high energies is a signature of this energization mechanism. The presence of field-aligned ions cannot be presumed to be a signature of recent ionospheric extraction.

Mauk, B. H.↗

Consistent large-eddy simulation of a temporal mixing layer laden with evaporating drops : Part 1. Direct numerical simulation formulation and a priori analysis

This paper is organized as follows. In § 2, we present the DNS formulation. This formulation consists of a set of coupled conservation equations for the gas phase and the liquid phase (i.e. the drops). The coupling describing the interaction between the two phases is provided in the gas-phase conservation equations by source terms. In § 3, we analyse the DNS results; in particular, we evaluate the adopted form of the heat and mass fluxes. The gas-phase LES equations are developed in § 4 and assumptions that may simplify these equations are examined using the DNS database. Thereafter, the contribution of the source terms is examined by considering the budget of the LES equations. The final form of the LES equations are then presented, in which only the SGS fluxes and the FSTs need be modelled. A parallel perspective on the importance of the source terms is presented in § 5, wherein the irreversible entropy production (i.e. the dissipation) expression for two-phase flow with phase change is derived and evaluated for transitional DNS and the filtered flow fields. Models for the FSTs appearing in the LES equations are investigated in § 6, while several constantcoefficient SGS models are proposed and calibrated in § 7. The present study aims to create simple models so that the modelling approach can be validated by means of an SGS models are proposed and calibrated in § 7. The present study aims to create simple models so that the modelling approach can be validated by means of an a posteriori study.

"Okong'o, Nora A.", 'Bellan, Josette↗

High-Fidelity Micromechanics Model Enhanced for Multiphase Particulate Materials

This 3-year effort involves the development of a comprehensive micromechanics model and a related computer code, capable of accurately estimating both the average response and the local stress and strain fields in the individual phases, assuming both elastic and inelastic behavior. During the first year (fiscal year 2001) of the investigation, a version of the model called the High-Fidelity Generalized Method of Cells (HFGMC) was successfully completed for the thermo-inelastic response of continuously reinforced multiphased materials with arbitrary periodic microstructures (refs. 1 and 2). The model s excellent predictive capability for both the macroscopic response and the microlevel stress and strain fields was demonstrated through comparison with exact analytical and finite element solutions. This year, HFGMC was further extended in two technologically significant ways. The first enhancement entailed the incorporation of fiber/matrix debonding capability into the two-dimensional version of HFGMC for modeling the response of unidirectionally reinforced composites such as titanium matrix composites, which exhibit poor fiber/matrix bond. Comparison with experimental data validated the model s predictive capability. The second enhancement entailed further generalization of HFGMC to three dimensions to enable modeling the response of particulate-reinforced (discontinuous) composites in the elastic material behavior domain. Next year, the three-dimensional version will be generalized to encompass inelastic effects due to plasticity, viscoplasticity, and damage, as well as coupled electromagnetothermomechanical (including piezoelectric) effects.

Pindera, Marek-Jerzy↗

Progress in Turbulence Detection via GNSS Occultation Data

The increased availability of radio occultation (RO) data offers the ability to detect and study turbulence in the Earth's atmosphere. An analysis of how RO data can be used to determine the strength and location of turbulent regions is presented. This includes the derivation of a model for the power spectrum of the log-amplitude and phase fluctuations of the permittivity (or index of refraction) field. The bulk of the paper is then concerned with the estimation of the model parameters. Parameter estimators are introduced and some of their statistical properties are studied. These estimators are then applied to simulated log-amplitude RO signals. This includes the analysis of global statistics derived from a large number of realizations, as well as case studies that illustrate various specific aspects of the problem. Improvements to the basic estimation methods are discussed, and their beneficial properties are illustrated. The estimation techniques are then applied to real occultation data. Only two cases are presented, but they illustrate some of the salient features inherent in real data.

Cornman, L. B.↗

The propagation of solar cosmic-ray bursts.

Description of a model showing analytically the three phases of anisotropy which occur during solar cosmic-ray events in the 7.5 to 21 MeV kinetic-energy interval and reported by McCracken et al. (1971): (1) a highly anisotropic, near field-aligned, initial phase, (2) a convective phase, and (3) a late-time phase in which the anisotropy is approximately perpendicular to the mean interplanetary magnetic field. The model is based on the cosmic-ray particles being convectively transported out from the sun, undergoing anisotropic diffusion along the interplanetary magnetic-field lines, and losing energy by adiabatic deceleration or by collision processes. The event is seen simply as a pulse moving outward from the sun after a cosmic-ray burst with a negative density-gradient in front of it and a positive gradient behind.

Ng, C. K.↗

Beaming of Uranian kilometric radiation and implications for the planetary magnetic field

Radio waves from Uranus at 481.2 kHz detected by Voyager 2 show a strong beaming effect. When the data are analyzed according to the magnetic latitude of the observer in an offset tilted dipole field model, a sharply defined emission beam extending from 25 to 55 deg latitude is evident consistent with a source over the south magnetic pole. The lack of prominent emission from the north magnetic pole can be attributed to the large spatial offset of the magnetic dipole resulting in much lower source frequencies. For proper alignment of the beam, an advance of the dipole axis longitude by 23 deg from the standard model is required. This phase shift is interpreted as the result of higher order magnetic moments producing a skew of the field lines at low altitudes in the polar source region. This conclusion further substantiates models of the interior of Uranus which predict a complex intrinsic field with significant nondipolar contributions to the surface magnetic field.

Barbosa, D. D.↗

Production of flickering aurora and field-aligned electron flux by electromagnetic ion cyclotron waves

Recent observations have suggested that flickering aurora is produced by a modulation of the field-aligned component of the electron flux within an auroral arc. It is proposed that a portion of the field-aligned electrons are of ionospheric origin and that these electrons are accelerated and their flux modulated by electromagnetic ion cyclotron waves that occur below the main acceleration region on auroral arc field lines. A model of the electromagnetic ion cyclotron wave shows that the parallel phase velocity of the wave increase as the wave propagates toward the ionosphere. A test particle calculation shows that ionospheric electrons trapped or reflected by the wave are accelerated to energies of several keV and that their flux is modulated at the wave frequency. The relative amplitudes of the model wave electric fields are consistent with the observations of small-scale low-frequency ionospheric and magnetospheric electric fields near auroral arcs of approximately 10 mV/m and 100 mV/m, respectively. The large-amplitude ion cyclotron waves also produce a ponderomotive force and a self-consistent ambipolar electric field. Energy considerations show that the downward energy flux in the electromagnetic ion cyclotron wave can be several percent of the total downward auroral electron energy flux.

Temerin, M.↗

Magnitude of B(sub Z) in the neutral sheet of the magnetotail

Statistical estimates of the average value of B(sub Z) in the magnetotail neutral sheet between 10 and 22 R(sub E) in X, and -10 to 10 R(sub E) in Y (GSM coordinates) are given for different phases of geomagnetic activity. With few exceptions the average value of B(sub Z) is found to be between 5 and 8 nT, but generally near 7 nT, irrespective of the type or phase of activity. This is considerably higher than is predicted in current magnetic field models of the magnetotail. The magnetic field within 2 R(sub E) of the neutral sheet is examined to determine whether B(sub Z) remains approximately constant. We do not find this to be true, B(sub Z) decreasing significantly over this distance from Delta-Z = 0. The structure of the cross-tail current and variations in current density with substorm phase appear to be more complicated than is commonly accepted.

Huang, C. Y.↗

Dynamical capture of field stars by interstellar clouds

Evidence suggesting the existence of old stars in young galactic clusters, has motivated an investigation of a model in which such stars are captured from the field during the gas cloud collapse phase of star formation. A hyperbolic field star cannot be gravitationally captured by a static cloud since the star's total energy is conserved. However, if the cloud's potential is an explicit function of time, the energy of a penetrating star can be decreased. As a first step in modeling the dynamical capture of field stars, capture by freely falling, homogeneous spherical cloud is considered. This simplified model is solved, and a prediction is obtained for number of captured field stars in terms of the cloud parameters and dispersion velocities of the field stars.

Whitman, P. G.↗

Forecasting Lightning Threat using Cloud-resolving Model Simulations

As numerical forecasts capable of resolving individual convective clouds become more common, it is of interest to see if quantitative forecasts of lightning flash rate density are possible, based on fields computed by the numerical model. Previous observational research has shown robust relationships between observed lightning flash rates and inferred updraft and large precipitation ice fields in the mixed phase regions of storms, and that these relationships might allow simulated fields to serve as proxies for lightning flash rate density. It is shown in this paper that two simple proxy fields do indeed provide reasonable and cost-effective bases for creating time-evolving maps of predicted lightning flash rate density, judging from a series of diverse simulation case study events in North Alabama for which Lightning Mapping Array data provide ground truth. One method is based on the product of upward velocity and the mixing ratio of precipitating ice hydrometeors, modeled as graupel only, in the mixed phase region of storms at the -15\dgc\ level, while the second method is based on the vertically integrated amounts of ice hydrometeors in each model grid column. Each method can be calibrated by comparing domainwide statistics of the peak values of simulated flash rate proxy fields against domainwide peak total lightning flash rate density data from observations. Tests show that the first method is able to capture much of the temporal variability of the lightning threat, while the second method does a better job of depicting the areal coverage of the threat. A blended solution is designed to retain most of the temporal sensitivity of the first method, while adding the improved spatial coverage of the second. Weather Research and Forecast Model simulations of selected North Alabama cases show that this model can distinguish the general character and intensity of most convective events, and that the proposed methods show promise as a means of generating quantitatively realistic fields of lightning threat. However, because models tend to have more difficulty in correctly predicting the instantaneous placement of storms, forecasts of the detailed location of the lightning threat based on single simulations can be in error. Although these model shortcomings presently limit the precision of lightning threat forecasts from individual runs of current generation models, the techniques proposed herein should continue to be applicable as newer and more accurate physically-based model versions, physical parameterizations, initialization techniques and ensembles of cloud-allowing forecasts become available.

McCaul, E. W., Jr.↗

Modelling of SAR polarisation phase difference from trees

The data for polarization phase difference Delta Phi between the HH- and VV-polarized backscattered waves from tree-covered fields were obtained with an airborne synthetic aperture radar at 1.225 GHz. The mean values over tree-covered fields were derived from the images of the phase difference and were examined as a function of incident beam angle from 15 to 55 deg. A theoretical model for simulating these data, based on the electromagnetic wave scatterings from the tree trunk and its branches, both of which are assumed as very long dielectric cylinders was developed. The radius and direction of a tree branch are taken as random variables and are chosen by a Monte Carlo method to encounter the incident waves in producing the scattering events. The Monte Carlo simulated results are in good agreement with the observations within experimental uncertainty.

Mo, Tsan↗

Gamma-Ray Light Curves from Pulsar Magnetospheres with Finite Conductivity

The Fermi Large Area Telescope has provided an unprecedented database for pulsar emission studies that includes gamma-ray light curves for over 100 pulsars. Modeling these light curves can reveal and constrain the geometry of the particle accelerator, as well as the pulsar magnetic field structure. We have constructed 3D magnetosphere models with finite conductivity, that bridge the extreme vacuum and force-free solutions used in previous light curves modeling. We are investigating the shapes of pulsar gamma-ray light curves using these dissipative solutions with two different approaches: (l) assuming geometric emission patterns of the slot gap and outer gap, and (2) using the parallel electric field provided by the resistive models to compute the trajectories and . emission of the radiating particles. The light curves using geometric emission patterns show a systematic increase in gamma-ray peak phase with increasing conductivity, introducing a new diagnostic of these solutions. The light curves using the model electric fields are very sensitive to the conductivity but do not resemble the observed Fermi light curves, suggesting that some screening of the parallel electric field, by pair cascades not included in the models, is necessary

Harding, A. K.↗

Reply to Comments by Tsurutani et al. on "Modeling Extreme 'Carrington-Type' Space Weather Events Using Three-Dimensional Global MHD Simulations"

In this response, we address the three main comments by Tsurutani et al. (2018, http://doi.org/10.1002/2017JA024779) namely, unusually high plasma density, interplanetary magnetic field intensity, and fast storm recovery phase. The authors agree that there is room to improve the modeling by taking into account these comments and other aspects that were not fully explored during our initial work. We are already in the process of undertaking a more comprehensive modeling project.

Ngwira, Chigomezyo M.↗

Heliomagnetic latitude dependence of the heliospheric magnetic field

Previous studies have revealed systematic variations of the interplanetary magnetic field with heliographic latitude. Luhmann et al. (1987) modeled Pioneer Venus (PVO) and ISEE-3 observations by assuming an asymmetric dependence on heliolatitude with stronger fields in the northern hemisphere. In a subsequent study, using data from ISEE-3/ICE and IMP-8, Burton et al. (1990) found evidence for a similar asymmetry. However, neither model has been completely successful. The model derived from PVO/ICE observations agrees quite well near solar maximum but shows significant discrepancies during the descending phase of the solar cycle. The model derived from the ICE/IMP-8 comparison suffers from significant phase delays between the difference in field magnitude at the two spacecraft and their latitude difference. In an attempt to account for these phase shifts, the IMP-8 and ICE data have been reexamined in heliomagnetic coordinates which are defined by the orientation of the solar magnetic dipole. The latitude and longitude of the dipole inferred from the data have then been compared with those implicit in source surface calculations. The IMP/ICE correlations have been extended into the recent solar maximum and descending phase. Comparisons have also been carried out between IMP-8 and Ulysses as it traveled to -30 deg south heliographic latitude.

Burton, M. E.↗

Hot gas ingestion test results of a two-poster vectored thrust concept with flow visualization in the NASA Lewis 9- x 15-foot low speed wind tunnel

A 9.2 percent scale Short Takeoff and Vertical Landing (STOVL) hot gas ingestion model was designed and built by McDonnell Douglas Corporation (MCAIR) and tested in the Lewis Research Center 9 x 15 foot Low Speed Wind Tunnel (LSWT). Hot gas ingestion, the entrainment of heated engine exhaust into the inlet flow field, is a key development issure for advanced short takeoff and vertical landing aircraft. Flow visualization from the Phase 1 test program, which evaluated the hot ingestion phenomena and control techniques, is covered. The Phase 2 test program evaluated the hot gas ingestion phenomena at higher temperatures and used a laser sheet to investigate the flow field. Hot gas ingestion levels were measured for the several forward nozzle splay configurations and with flow control/life improvement devices (LIDs) which reduced the hot gas ingestion. The model support system had four degrees of freedom - pitch, roll, yaw, and vertical height variation. The model support system also provided heated high-pressure air for nozzle flow and a suction system exhaust for inlet flow. The test was conducted at full scale nozzle pressure ratios and inlet Mach numbers. Test and data analysis results from Phase 2 and flow visualization from both Phase 1 and 2 are documented. A description of the model and facility modifications is also provided. Headwind velocity was varied from 10 to 23 kn. Results are presented over a range of nozzle pressure ratios at a 10 kn headwind velocity. The Phase 2 program was conducted at exhaust nozzle temperatures up to 1460 R and utilized a sheet laser system for flow visualization of the model flow field in and out of ground effects. The results reported are for nozzle exhaust temperatures up to 1160 R. These results will contain the compressor face pressure and temperature distortions, the total pressure recovery, the inlet temperature rise, and the environmental effects of the hot gas. The environmental effects include the ground plane contours, the model airframe heating, and the location of the ground flow separation.

Johns, Albert L.↗