Search NASA⌕ Search

SEARCH · Search NASA

Results for “phase field simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 271 records · Page 15

Nonlinear evolution of radiation-driven thermally unstable fluids

The nonlinear evolution of a radiation-driven thermally unstable planar fluid is simulated numerically using a semiimplicit finite-difference algorithm. When the equilibrium state of the fluid is perturbed by random initial excitation of the velocity field, dense, cool, two-dimensional structures are found to form in a rarer, warmer surrounding medium. The nonlinear phase of evolution is characterized by the turbulent contraction of the condensed region, accompanied by a significant increase in the amount of energy radiated. It is found that, if the random velocity perturbation has a sufficiently large amplitude, the fluid will not form condensed structures. Finally, the relationship of these results to observations of the solar chromosphere, transition region, and corona is discussed.

Dahlburg, R. B.↗

Simulation of ion-cyclotron-like modes in a magnetoplasma with transverse inhomogeneous electric field

The new kinetic ion-cyclotronlike (IC-like) instability driven by a localized electric field perpendicular to the ambient magnetic field is investigated using particle simulation. The existence of the new mode is demonstrated and the mode is distinguished from the well-known Kelvin-Helmholtz branch and the kinetic modification of this branch. The issue of the initial loading of particles in phase space in the simulation is clarified. The simulation results show that the growing ion waves are associated with the small vortices in the linear state, which evolve into a nonlinear stage dominated by large vortices with lower frequencies.

Nishikawa, K.-I.↗

Direct numerical simulations of stably-stratified sheared turbulence: Implications for oceanic mixing

Direct numerical simulations of the time evolution of homogeneous stably stratified turbulent shear flows have been performed for several Richardson numbers Ri and Reynolds numbers R(sub lambda) in earlier works. The results show excellent agreement with length scale models developed from laboratory experiments to characterize oceanic turbulence. When the Richardson number Ri is less than the stationary value Ri(sub s), the turbulence intensity grows at all scales, and the growth rate appears to be a function of Ri. The size of the vertical density inversions also increases. On the other hand, when Ri is greater than or equal to Ri(sub s) the largest turbulent eddies become vertically constrained by buoyancy when the Ellison (turbulence) scale L(sub E) and the Ozmidov (buoyancy) scale L(sub O) are equal. At this point, the mixing efficiency is maximal and corresponds to a flux Richardson number R(sub f) = 0.20. The vertical mass flux becomes counter-gradient when epsilon = 19(nu)N(exp 2) and vertical density overturns are suppressed in less than half a Brunt-Vaisala period. The results of the simulations were also recast in terms of the Hydrodynamic Phase Diagram introduced in fossil turbulence models. The so-called point of fossilization occurs when epsilon = 4DCN(exp 2); Gibson proposed 13DCN(exp 2). This value is in agreement with indirect laboratory observations and field observations. Finally, the validity of the steady-state models to estimate vertical eddy diffusivities in the oceanic thermocline is discussed.

Itsweire, E. C.↗

Gyrophase-bunched Electrons: Cluster Observations and Simulations

The interaction of the solar wind with the Earth's bow shock is responsible for a number of phenomena both upstream and downstream. One of the least understood and studied of these is the generation of gyrophase-bunched electrons. We will describe initial work searching for and analyzing such events in (Plasma Electron And Current Experiment (PEACE) data from CLUSTER. Plasma simulations will complement the data analysis. Recent analyses suggest that gyrophase-bunched electrons are present rather frequently in the upstream region, due either to energization in the reflection off the shock front or, more likely, to phase trapping in locally produced whistler wave fields.

Gurgiolo, C.↗

Effect of Aft Rotor on the Inter-Rotor Flow of an Open Rotor Propulsion System

The effects of the aft rotor on the inter-rotor flow field of an open rotor propulsion rig were examined. A Particle Image Velocimetry (PIV) dataset that was acquired phase locked to the front rotor position has been phase averaged based on the relative phase angle between the forward and aft rotors. The aft rotor phase was determined by feature tracking in raw PIV images through an image processing algorithm. The effects of the aft rotor potential field on the inter-rotor flow were analyzed and shown to be in reasonably good agreement with Computational Fluid Dynamics (CFD) simulations. The aft rotor position was shown to have a significant upstream effect, with implications for front rotor interaction noise. It was found that the aft rotor had no substantial effect on the position of the forward rotor tip vortex but did have a small effect on the circulation strength of the vortex when the rotors were highly loaded.

Counter Rotating↗

Radar Differential Phase Signatures of Ice Orientation for the Prediction of Lightning Initiation and Cessation

The majority of lightning-related casualties typically occur during thunderstorm initiation (e.g., first flash) or dissipation (e.g., last flash). The physics of electrification and lightning production during thunderstorm initiation is fairly well understood. As such, the literature includes a number of studies presenting various radar techniques (using reflectivity and, if available, other dual-polarimetric parameters) for the anticipation of initial electrification and first lightning flash. These radar techniques have shown considerable skill at forecasting first flash. On the other hand, electrical processes and lightning production during thunderstorm dissipation are not nearly as well understood and few, if any, successful techniques have been developed to anticipate the last flash and subsequent cessation of lightning. One promising approach involves the use of dual-polarimetric radar variables to infer the presence of oriented ice crystals in lightning producing storms. In the absence of strong vertical electric fields, ice crystals fall with their largest (semi-major) axis in the horizontal associated with gravitational and aerodynamic forces. In thunderstorms, strong vertical electric fields (100-200 kV m(sup -1)) have been shown to orient small (less than 2 mm) ice crystals such that their semi-major axis is vertical (or nearly vertical). After a lightning flash, the electric field is typically relaxed and prior radar research suggests that ice crystals rapidly resume their preferred horizontal orientation. In active thunderstorms, the vertical electric field quickly recovers and the ice crystals repeat this cycle of orientation for each nearby flash. This change in ice crystal orientation from primarily horizontal to vertical during the development of strong vertical electric fields prior to a lightning flash forms the physical basis for anticipating lightning initiation and, potentially, cessation. Research has shown that radar reflectivity (Z) and other co-polar back-scattering radar measurements like differential reflectivity (Z(sub dr)) typically measured by operational dual-polarimetric radars are not sensitive to these changes in ice crystal orientation. However, prior research has demonstrated that oriented ice crystals cause significant propagation effects that can be routinely measured by most dual-polarimetric radars from X-band (3 cm) to S-band (10 cm) wavelengths using the differential propagation phase shift (often just called differential phase, phi(sub dp)) or its range derivative, the specific differential phase (K(sub dp)). Advantages of the differential phase include independence from absolute or relative power calibration, attenuation, differential attenuation and relative insensitivity to ground clutter and partial beam occultation effects (as long as the signal remains above noise). In research mode, these sorts of techniques have been used to anticipate initial cloud electrification, lightning initiation, and cessation. In this study, we develop a simplified model of ice crystal size, shape, orientation, dielectric, and associated radar scattering and propagation effects in order to simulate various idealized scenarios of ice crystals responding to a hypothetical electric field and their dual-polarimetric radar signatures leading up to lightning initiation and particularly cessation. The sensitivity of the K(sub dp) ice orientation signature to various ice properties and radar wavelength will be explored. Since K(sub dp) is proportional to frequency in the Rayleigh- Gans scattering regime, the ice orientation signatures should be more obvious at higher (lower) frequencies (wavelengths). As a result, simulations at radar wavelengths from 10 cm down to 1 cm (Ka-band) will be conducted. Resonance effects will be considered using the T-matrix method. Since most K(sub dp) Vbased observations have been shown at S-band, we will present ice orientation signatures from C-band (UAH/NASA ARMOR) and X-bd (UAH MAX) dual-polarimetric radars located in Northern Alabama. Issues related to optimal radar scanning for the detection of oriented ice will be discussed. Preliminary suggestions on how these differential phase signatures of oriented ice could contribute to lightning initiation and cessation algorithms will be presented.

Carey, L.D.↗

Short hyperdynamic profiles influence primate temperature regulation

Primates have been shown to be sensitive to hyperdynamic fields. That is, when exposed to + 2Gz, body temperature falls. The purpose of this study was to examine the relative sensitivity of these animals to short centrifugation profiles which mimic the gravitational envelope seen on the Space Shuttle during launch (8 minutes, 2.9 Gz max) and re-entry (19 min, 1.7 Gz max). Four loosely restrained squirrel monkeys, isolated from additional external stimuli, were exposed to these profiles. During launch simulation, the temperatures never fell markedly below control levels. However, subsequent to return to 1G, the recovery phase showed decreases in body temperature in all four animals averaging 0.4 C over the next 10 to 15 minutes. The two animals exposed to the reentry profile showed decreases in body temperature within five minutes of the onset of centrifugation. Maximum fall in body temperature was reached by the end of the centrifugation phase and averaged 0.7 C. Thus, the temperature regulation system of this primate is sensitive to short hyperdynamic field exposures.

Fuller, C. A.↗

Consistent Evaluation of ACOS-GOSAT, BESD-SCIAMACHY, CarbonTracker, and MACC Through Comparisons to TCCON

Consistent validation of satellite CO2 estimates is a prerequisite for using multiple satellite CO2 measurements for joint flux inversion, and for establishing an accurate long-term atmospheric CO2 data record. Harmonizing satellite CO2 measurements is particularly important since the differences in instruments, observing geometries, sampling strategies, etc. imbue different measurement characteristics in the various satellite CO2 data products. We focus on validating model and satellite observation attributes that impact flux estimates and CO2 assimilation, including accurate error estimates, correlated and random errors, overall biases, biases by season and latitude, the impact of coincidence criteria, validation of seasonal cycle phase and amplitude, yearly growth, and daily variability. We evaluate dry-air mole fraction (X(sub CO2)) for Greenhouse gases Observing SATellite (GOSAT) (Atmospheric CO2 Observations from Space, ACOS b3.5) and SCanning Imaging Absorption spectroMeter for Atmospheric CHartographY (SCIAMACHY) (Bremen Optimal Estimation DOAS, BESD v2.00.08) as well as the CarbonTracker (CT2013b) simulated CO2 mole fraction fields and the Monitoring Atmospheric Composition and Climate (MACC) CO2 inversion system (v13.1) and compare these to Total Carbon Column Observing Network (TCCON) observations (GGG2012/2014). We find standard deviations of 0.9, 0.9, 1.7, and 2.1 parts per million vs. TCCON for CT2013b, MACC, GOSAT, and SCIAMACHY, respectively, with the single observation errors 1.9 and 0.9 times the predicted errors for GOSAT and SCIAMACHY, respectively. We quantify how satellite error drops with data averaging by interpreting according to (error(sup 2) equals a(sup 2) plus b(sup 2) divided by n (with n being the number of observations averaged, a the systematic (correlated) errors, and b the random (uncorrelated) errors). a and b are estimated by satellites, coincidence criteria, and hemisphere. Biases at individual stations have year-to-year variability of 0.3 parts per million, with biases larger than the TCCON predicted bias uncertainty of 0.4 parts per million at many stations. We find that GOSAT and CT2013b under-predict the seasonal cycle amplitude in the Northern Hemisphere (NH) between 46 and 53 degrees North latitude, MACC over-predicts between 26 and 37 degrees North latitude, and CT2013b under-predicts the seasonal cycle amplitude in the Southern Hemisphere (SH). The seasonal cycle phase indicates whether a data set or model lags another data set in time. We find that the GOSAT measurements improve the seasonal cycle phase substantially over the prior while SCIAMACHY measurements improve the phase significantly for just two of seven sites. The models reproduce the measured seasonal cycle phase well except for at Lauder_125HR (CT2013b) and Darwin (MACC). We compare the variability within 1 day between TCCON and models in June-July-August; there is correlation between 0.2 and 0.8 in the NH, with models showing 10-50 percent the variability of TCCON at different stations and CT2013b showing more variability than MACC. This paper highlights findings that provide inputs to estimate flux errors in model assimilations, and places where models and satellites need further investigation, e.g., the SH for models and 45-67 degrees North latitude for GOSAT and CT2013b.

ACOS-GOSAT↗

Sources, Sinks, and Transatlantic Transport of North African Dust Aerosol: A Multimodel Analysis and Comparison With Remote Sensing Data

This study evaluates model-simulated dust aerosols over North Africa and the North Atlantic from five global models that participated in the Aerosol Comparison between Observations and Models phase II model experiments. The model results are compared with satellite aerosol optical depth (AOD) data from Moderate Resolution Imaging Spectroradiometer (MODIS), Multiangle Imaging Spectroradiometer (MISR), and Sea-viewing Wide Field-of-view Sensor, dust optical depth (DOD) derived from MODIS and MISR, AOD and coarse-mode AOD (as a proxy of DOD) from ground-based Aerosol Robotic Network Sun photometer measurements, and dust vertical distributions/centroid height from Cloud Aerosol Lidar with Orthogonal Polarization and Atmospheric Infrared Sounder satellite AOD retrievals. We examine the following quantities of AOD and DOD: (1) the magnitudes over land and over ocean in our study domain, (2) the longitudinal gradient from the dust source region over North Africa to the western North Atlantic, (3) seasonal variations at different locations, and (4) the dust vertical profile shape and the AOD centroid height (altitude above or below which half of the AOD is located). The different satellite data show consistent features in most of these aspects; however, the models display large diversity in all of them, with significant differences among the models and between models and observations. By examining dust emission, removal, and mass extinction efficiency in the five models, we also find remarkable differences among the models that all contribute to the discrepancies of model-simulated dust amount and distribution. This study highlights the challenges in simulating the dust physical and optical processes, even in the best known dust environment, and stresses the need for observable quantities to constrain the model processes.

Africa↗

One-dimensional numerical study of charged particle trajectories in turbulent electrostatic wave fields

The paper describes a numerical simulation of electron trajectories in weak random electric fields under conditions that are approximately true for Langmuir waves whose wavelength is much longer than the Debye length. Two types of trajectory calculations were made: (1) the initial particle velocity was made equal to the mean phase velocity of the waves, or (2) it was equal to 0.7419 times the mean velocity of the waves, so that the initial velocity differed substantially from all phase velocities of the wave spectrum. When the autocorrelation time is much greater than the trapping time, the particle motion can change virtually instantaneously from one of three states - high-velocity, low-velocity, or trapped state - to another. The probability of instantaneous transition from a high- or low-velocity state becomes small when the difference between the particle velocity and the mean phase velocity of the waves becomes high in comparison to the trapping velocity. Diffusive motion becomes negligible under these conditions also.

Graham, K. N.↗

Ultrasonic techniques for repair of aircraft structures with bonded composite patches

This is a paper on a research and development project to demonstrate a novel ultrasonic process for the field application of boron/epoxy (B/Ep) patches for repair of aircraft structures. The first phase of the project was on process optimization and testing to develop the most practical ultrasonic processing techniques. Accelerated testing and aging behavior of precured B/Ep patches, which were ultrasonically bonded to simulated B-52 wing panel assemblies, were performed by conducting flight-by-flight spectrum loading fatigue tests. The spectrum represented 2340 missions/flights or 30 years of service. The effects of steady-state applied temperature and prior exposure of the B/Ep composite patches were evaluated. Representative experimental results of this phase of the project are presented.

Smith, S. H.↗

Comparison of Measurements and FluorMOD Simulations for Solar Induced Chlorophyll Fluorescence and Reflectance of a Corn Crop under Nitrogen Treatments [SIF and Reflectance for Corn]

The FLuorescence Explorer (FLEX) satellite concept is one of six semifinalist mission proposals selected in 2006 for pre-Phase studies by the European Space Agency (ESA). The FLEX concept proposes to measure passive solar induced chlorophyll fluorescence (SIF) of terrestrial ecosystems. A new spectral vegetation Fluorescence Model (FluorMOD) was developed to include the effects of steady state SIF on canopy reflectance. We used our laboratory and field measurements previously acquired from foliage and canopies of corn (Zea mays L.) under controlled nitrogen (N) fertilization to parameterize and evaluate FluorMOD. Our data included biophysical properties, fluorescence (F) and reflectance spectra for leaves; reflectance spectra of canopies and soil; solar irradiance; plot-level leaf area index; and canopy SIF emissions determined using the Fraunhofer Line Depth principal for the atmospheric telluric oxygen absorption features at 688 nm (O2-beta) and 760 nm (O2-alpha). FluorMOD simulations implemented in the default "look-up-table" mode did not reproduce the observed magnitudes of leaf F, canopy SIF, or canopy reflectance. However, simulations for all of these parameters agreed with observations when the default FluorMOD information was replaced with measurements, although N treatment responses were underestimated. Recommendations were provided to enhance FluorMOD's potential utility in support of SIF field experiments and studies of agriculture and ecosystems.

Middleton, Elizabeth M.↗

Simulating "Mars on Earth"

By now, everyone who's heard of the Haughton-Mars Project knows that we travel to Devon Island to learn how people will live and work on Mars. But how do we learn about Mars operations from what happens in the Arctic? We must document our experience--traverses, life in the hab, instrument deployment, communications, and so on. Then we must analyze and formally model what happens. In short, while most scientists are studying the crater, other scientists must be studying the expedition itself. That's what I have done in the past four field seasons. I study field science, both as it naturally occurs at Haughton (unconstrained by a "Mars Sam") and as a constrained experiment using the Flashline Mars Arctic Research Station. During the second week of July 2001, I lived and worked in the hab as part of the Phase 2 crew of six. Besides participating in all activities, I took many photographs and time lapse video. The result of my work will be a computer simulation of how we lived and worked in the hab. It won't be a model of particular people or even my own phase per se, but a pastiche that demonstrates (a proof of concept) that we have appropriate tools for simulating the layout of the hab and daily routines followed by the group and individual scientists. Activities-how people spend their time-are the focus of my observations for building such a simulation model.

Clancey, William J.↗

Double-diffusive convection with sidewalls

Stommel et al. (1956) have first described an instability, known as thermosolutal convection, thermohaline convection, or double-diffusive convection. This instability may occur in the case of a fluid in a gravitational field with two diffusing components present. The present study is concerned with the effect of sidewalls on flow in the fingering regime in the absence of applied horizontal gradients. The work was motivated by numerical results obtained on the basis of a simulation of thermosolutal convection occurring during the unidirectional solidification of a binary alloy. In this case, the unperturbed solute field in the liquid ahead of the solidifying planar interface has an exponential vertical profile because of the rejection or preferential incorporation of solute by the solid phase.

Mcfadden, G. B.↗

Solar Wind Ablation of Terrestrial Planet Atmospheres

Internal plasma sources usually arise in planetary magnetospheres as a product of stellar ablation processes. With the ignition of a new star and the onset of its ultraviolet and stellar wind emissions, much of the volatiles in the stellar system undergo a phase transition from gas to plasma. Condensation and accretion into a disk is replaced by radiation and stellar wind ablation of volatile materials from the system- Planets or smaller bodies that harbor intrinsic magnetic fields develop an apparent shield against direct stellar wind impact, but UV radiation still ionizes their gas phases, and the resulting internal plasmas serve to conduct currents to and from the central body along reconnected magnetic field linkages. Photoionization and thermalization of electrons warms the ionospheric topside, enhancing Jeans' escape of super-thermal particles, with ambipolar diffusion and acceleration. Moreover, observations and simulations of auroral processes at Earth indicate that solar wind energy dissipation is concentrated by the geomagnetic field by a factor of 10-100, enhancing heavy species plasma and gas escape from gravity, and providing more current carrying capacity. Thus internal plasmas enable coupling with the plasma, neutral gas and by extension, the entire body. The stellar wind is locally loaded and slowed to develop the required power. The internal source plasma is accelerated and heated, inflating the magnetosphere as it seeks escape, and is ultimately blown away in the stellar wind. Bodies with little sensible atmosphere may still produce an exosphere of sputtered matter when exposed to direct solar wind impact. Bodies with a magnetosphere and internal sources of plasma interact more strongly with the stellar wind owing to the magnetic linkage between the two created by reconnection.

Moore, Thomas Earle↗

Space vacuum investigations

Inaccessible-Non Repairable- Costly: The meaning of these words is absolute in the event of a complete failure of an orbiting spacecraft. Adequate design, or rather, minimizing chance of failure, challenges all scientists and engineers. Where experience has not been developed, the risk involved in a new application must be weighed against the consequence of a design failure. In space work, the spectrum of consequence is extremely wide. It ranges from loss of scientific data on one hand, to loss of life on the other.. There is little question that knowledge of the space environment and its effects must be determined to produce reliable spacecraft. Tests to prove the design must be performed with equipment capable of simulating the effects of the environment in which the design is to operate. This does not necessarily mean that the natural environment must be duplicated. The degree of environmental simulation that should be achieved depends primarily upon the significance of the physical effects it causes. The purpose of space vacuum simulation in an environmental test program should be accurately conceived and the testing phased to obtain maximum information from the limited samples available for test. This discussion is confined to some of the factors involved in simulating the effects of space vacuum environment. This survey is by no means complete, but will touch on some of the highlights in this field: 1. A description of the environment 2. The level of simulation necessary to evaluate its effects 3. Test equipment capability 4. Experience in test 5. Experience in orbit 6. Observations & future emphasis

VACUUM EFFECT↗

Modeling of convection phenomena in Bridgman-Stockbarger crystal growth

Thermal convection phenomena in a vertically oriented Bridgman-Stockbarger apparatus were modeled by computer simulations for different gravity conditions, ranging from earth conditions to extremely low gravity, approximate space conditions. The modeling results were obtained by the application of a state-of-the art, transient, multi-dimensional, completely densimetrically coupled, discrete-element computational model which was specifically developed for the simulation of flow, temperature, and species concentration conditions in two-phase (solid-liquid) systems. The computational model was applied to the simulation of the flow and the thermal conditions associated with the convection phenomena in a modified Germanium-Silicon charge enclosed in a stationary fused-silica ampoule. The results clearly indicated that the gravitational field strength influences the characteristics of the coherent vortical flow patterns, interface shape and position, maximum melt velocity, and interfacial normal temperature gradient.

Carlson, F. M.↗

Energy buildup in coronal magnetic flux tubes

A time-dependent two-dimensional MHD simulation is used to study the response of the magnetic field in coronal loops to photospheric motion. From an initially uniform field, circular sections of the ends of the loop are slowly rotated to represent the photospheric motion. The evolution of the field and flow is characterized by three phases: (1) a phase of negligible kinetic energy where the current and field are predominantly parallel; (2) a phase where the field twist increases, the axial field at and near the axis increases, and the axial field decreases in two cylindrical regions away from the axis; and (3) a phase in which a significant portion of the field makes several rotations at large radii, with a corresponding reducton in the axial field to a few percent of the initial value.

Steinolfson, R. S.↗