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

Three-dimensional solidification and melting using magnetic field control

A new two-fluid mathematical model for fully three dimensional steady solidification under the influence of an arbitrary acceleration vector and with or without an arbitrary externally applied steady magnetic field have been formulated and integrated numerically. The model includes Joule heating and allows for separate temperature dependent physical properties within the melt and the solid. Latent heat of phase change during melting/solidification was incorporated using an enthalpy method. Mushy region was automatically captured by varying viscosity orders of magnitude between liquidus and solidus temperature. Computational results were obtained for silicon melt solidification in a parallelepiped container cooled from above and from a side. The results confirm that the magnetic field has a profound influence on the solidifying melt flow field thus changing convective heat transfer through the boundaries and the amount and shape of the solid accrued. This suggests that development of a quick-response algorithm for active control of three dimensional solidification is feasible since it would require low strength magnetic fields.

Dulikravich, George S.↗

Interferometric Tomographic Measurement of an Instataneous Flow Field Under Adverse Environments

Measurement of an instantaneous flow field by interferometric tomography, that is, reconstruction of a three-dimensional refractive-index field from multi-directional projection data, has been conducted. In order to simulate the expected experimental arrangement at a wind tunnel, reconstructions are made from a restricted view angle less than 40 degrees and incomplete projections. In addition, appreciable ambient air and experimental setup disturbances are present. A new phase-stepping technique, based on a generalized phase-stepping approach of a four-bucket model, is applied for expeditious and accurate phase information extraction from projection interferograms under the harsh environments. Phase errors caused by the various disturbances, which can include ambient refractive-index change, optical component disturbance, hologram repositioning error, etc., are partially compensated with a linear corrective model. A new computational tomographic technique based on a series expansion approach was also utilized to efficiently deal with arbitrary boundary shapes and the continuous flow fields in reconstruction. The results of the preliminary investigation are encouraging; however, the technique needs to be further developed in the future through refinement of the approaches reported here and through hybridization with previously developed techniques. Keywords: interferometry, tomography, phase stepping

Yu, En-Xi↗

Global Simulation of Proton Precipitation Due to Field Line Curvature During Substorms

The low latitude boundary of the proton aurora (known as the Isotropy Boundary or IB) marks an important boundary between empty and full downgoing loss cones. There is significant evidence that the IB maps to a region in the magnetosphere where the ion gyroradius becomes comparable to the local field line curvature. However, the location of the IB in the magnetosphere remains in question. In this paper, we show simulated proton precipitation derived from the Field Line Curvature (FLC) model of proton scattering and a global magnetohydrodynamic simulation during two substorms. The simulated proton precipitation drifts equatorward during the growth phase, intensifies at onset and reproduces the azimuthal splitting published in previous studies. In the simulation, the pre-onset IB maps to 7-8 RE for the substorms presented and the azimuthal splitting is caused by the development of the substorm current wedge. The simulation also demonstrates that the central plasma sheet temperature can significantly influence when and where the azimuthal splitting takes place.

Field Line Curvature↗

A generic multi-flex-body dynamics, controls simulation tool for space station

An order (n) multiflex body Space Station simulation tool is introduced. The flex multibody modeling is generic enough to model all phases of Space Station from build up through to Assembly Complete configuration and beyond. Multibody subsystems such as the Mobile Servicing System (MSS) undergoing a prescribed translation and rotation are also allowed. The software includes aerodynamic, gravity gradient, and magnetic field models. User defined controllers can be discrete or continuous. Extensive preprocessing of 'body by body' NASTRAN flex data is built in. A significant aspect, too, is the integrated controls design capability which includes model reduction and analytic linearization.

London, Ken W.↗

Lunar gravitational field estimation and the effects of mismodeling upon lunar satellite orbit prediction

Lunar spherical harmonic gravity coefficients are estimated from simulated observations of a near-circular low altitude polar orbiter disturbed by lunar mascons. Lunar gravity sensing missions using earth-based nearside observations with and without satellite-based far-side observations are simulated and least squares maximum likelihood estimates are developed for spherical harmonic expansion fit models. Simulations and parameter estimations are performed by a modified version of the Smithsonian Astrophysical Observatory's Planetary Ephemeris Program. Two different lunar spacecraft mission phases are simulated to evaluate the estimated fit models. Results for predicting state covariances one orbit ahead are presented along with the state errors resulting from the mismodeled gravity field. The position errors from planning a lunar landing maneuver with a mismodeled gravity field are also presented. These simulations clearly demonstrate the need to include observations of satellite motion over the far side in estimating the lunar gravity field. The simulations also illustrate that the eighth degree and order expansions used in the simulated fits were unable to adequately model lunar mascons.

Davis, John H.↗

A proposed production model of rapid subauroral ion drifts and their relationship to substorm evolution

The temporal relationship between subauroral ion drifts (SAIDs) and the phases of an auroral substorm is examined on the basis of multisatellite data. The time of expansive phase onset is identified and the time at which recovery begins is estimated. SAIDs are found to typically occur well after substorm onset (more than 30 min), during the substorm recovery phase. Substantial westward ion drifts and field-aligned currents are observed well equatorward of the auroral oval during the expansion phase of a substorm, but the drifts lack the narrow spike signature associated with SAIDs. A phenomenological model of SAID production that qualitatively agrees with the observed ionospheric signatures and substorm temporal relationship is proposed.

Anderson, P. C.↗

Use of coherence and phase data between two receivers in evaluation of noise environments

For certain types of noise control problems, where transducers cannot be mounted on suspected sources, valuable information can often be obtained by comparing the coherence and phase data measured between two closely spaced microphones with analytical models deduced from the physics of the problem. However, the application of such analysis techniques must be pursued with care, particularly when the measurements are made in a reverberant area. A simple illustration is presented where the acoustic field in the test section of a wind tunnel is evaluated by modelling the field as a combination of diffuse noise due to the boundary layer turbulence in the test section and propagating noise generated by the tunnel fan and possible flow disturbances outside the test section. The coherence and phase between two closely spaced microphones in the tunnel test section are predicted for various ratios of diffuse to propagating noise contributions and compared to actual measurements under several different tunnel operating conditions.

Piersol, A. G.↗

A finite element model of conduction, convection, and phase change near a solid/melt interface

Detailed understanding of heat transfer and fluid flow is required for many aerospace thermal systems. These systems often include phase change and operate over a range of accelerations or effective gravitational fields. An approach to analyzing such systems is presented which requires the simultaneous solution of the conservation laws of energy, momentum, and mass, as well as an equation of state. The variable property form of the governing equations are developed in two-dimensional Cartesian coordinates for a Newtonian fluid. A numerical procedure for solving the governing equations is presented and implemented in a computer program. The Galerkin form of the finite element method is used to solve the spatial variation of the field variables, along with the implicit Crank-Nicolson time marching algorithm. Quadratic Langrangian elements are used for the internal energy and the two components of velocity. Linear Lagrangian elements are used for the pressure. The location of the solid/liquid interface as well as the temperatures are determined form the calculated internal energy and pressure. This approach is quite general in that it can describe heat transfer without phase change, phase change with a sharp interface, and phase change without an interface. Analytical results from this model are compared to those of other researchers studying transient conduction, convection, and phase change and are found to be in good agreement. The numerical procedure presented requires significant computer resources, but this is not unusual when compared to similar studies by other researchers. Several methods are suggested to reduce the computational times.

Viterna, Larry A.↗

Extreme energetic particle decreases near geostationary orbit - A manifestation of current diversion within the inner plasma sheet

A qualitative model of magnetic field reconfiguration as might result from neutral line formation in the central plasma sheet late in a substorm growth phase is considered. It is suggested that magnetic reconnection probably begins before the substorm expansion phase and that cross-tail current is enhanced across the plasma sheet both earthward and tailward of a limited region near the neutral line. Such an enhanced cross-tail current earthward of the original X line region may contribute to thinning the plasma sheet substantially, and this would in turn affect the drift currents in that location, thus enhancing the current even closer toward the earth. In this way a redistribution and progressive diversion of normal cross-tail current throughout much of the inner portion of the plasma sheet could occur. The resulting intensified current, localized at the inner edge of the plasma sheet, would lead to a very thin plasma confinement region. This would explain the very taillike field and extreme particle dropouts often seen late in substorm growth phases.

Baker, D. N.↗

Three-dimensional aspects of radiative transfer in remote sensing of precipitation: Application to the 1986 COHMEX storm

Several multifrequency techniques for passive microwave estimation of precipitation based on the absorption and scattering properties of hydrometers have been proposed in the literature. In the present study, plane-parallel limitations are overcome by using a model based on the discrete-ordinates method to solve the radiative transfer equation in three-dimensional rectangular domains. This effectively accounts for the complexity and variety of radiation problems encountered in the atmosphere. This investigation presents result for plane-parallel and three-dimensional radiative transfer for a precipitating system, discusses differences between these results, and suggests possible explanations for these differences. Microphysical properties were obtained from the Colorado State University Regional Atmospehric Modeling System and represent a hailstorm observed during the 1986 Cooperative Huntsville Meteorological Experiment. These properties are used as input to a three-dimensional radiative transfer model in order to simulate satellite observation of the storm. The model output consists of upwelling brightness temperatures at several of the frequencies on the Special Sensor Microwave/Imager. The radiative transfer model accounts for scattering and emission of atmospheric gases and hydrometers in liquid and ice phases. Brightness temperatures obtained from the three-dimensional model of this investigation indicate that horizontal inhomogeneities give rise to brightness temperature fields that can be quite different from fields obtained using plane-parallel radiative transfer theory. These differences are examined for various resolutions of the satellite sensor field of view. In adddition, the issue of boundary conditions for three-dimensional atmospheric radiative transfer is addressed.

Haferman, J. L.↗

Advanced Propulsion Physics Lab: Eagleworks Investigations

Eagleworks Laboratory is an advanced propulsions physics laboratory with two primary investigations currently underway. The first is a Quantum Vacuum Plasma Thruster (QVPT or Q-thrusters), an advanced electric propulsion technology in the development and demonstration phase. The second investigation is in Warp Field Interferometry (WFI). This is an investigation of Dr. Harold "Sonny" White's theoretical physics models for warp field equations using optical experiments in the Electro Optical laboratory (EOL) at Johnson Space Center. These investigations are pursuing technology necessary to enable human exploration of the solar system and beyond.

Scogin, Tyler↗

Delay and Disruption Tolerant Networking MACHETE Model

To verify satisfaction of communication requirements imposed by unique missions, as early as 2000, the Communications Networking Group at the Jet Propulsion Laboratory (JPL) saw the need for an environment to support interplanetary communication protocol design, validation, and characterization. JPL's Multi-mission Advanced Communications Hybrid Environment for Test and Evaluation (MACHETE), described in Simulator of Space Communication Networks (NPO-41373) NASA Tech Briefs, Vol. 29, No. 8 (August 2005), p. 44, combines various commercial, non-commercial, and in-house custom tools for simulation and performance analysis of space networks. The MACHETE environment supports orbital analysis, link budget analysis, communications network simulations, and hardware-in-the-loop testing. As NASA is expanding its Space Communications and Navigation (SCaN) capabilities to support planned and future missions, building infrastructure to maintain services and developing enabling technologies, an important and broader role is seen for MACHETE in design-phase evaluation of future SCaN architectures. To support evaluation of the developing Delay Tolerant Networking (DTN) field and its applicability for space networks, JPL developed MACHETE models for DTN Bundle Protocol (BP) and Licklider/Long-haul Transmission Protocol (LTP). DTN is an Internet Research Task Force (IRTF) architecture providing communication in and/or through highly stressed networking environments such as space exploration and battlefield networks. Stressed networking environments include those with intermittent (predictable and unknown) connectivity, large and/or variable delays, and high bit error rates. To provide its services over existing domain specific protocols, the DTN protocols reside at the application layer of the TCP/IP stack, forming a store-and-forward overlay network. The key capabilities of the Bundle Protocol include custody-based reliability, the ability to cope with intermittent connectivity, the ability to take advantage of scheduled and opportunistic connectivity, and late binding of names to addresses.

Segui, John S.↗

Modeling the Superstorm in November 2003

The superstorm on 20.21 November 2003 was the largest geomagnetic storm in solar cycle 23 as measured by Dst, which attained a minimum value of .422 nT. We have simulated this storm to understand how particles originating from the solar wind and ionosphere get access to the magnetosphere and how the subsequent transport and energization processes contribute to the buildup of the ring current. The global electromagnetic configuration and the solar wind H+ distribution are specified by the Lyon-Fedder-Mobarry (LFM) magnetohydrodynamics model. The outflow of H+ and O+ ions from the ionosphere are also considered. Their trajectories in the magnetosphere are followed by a test-particle code. The particle distributions at the inner plasma sheet established by the LFM model and test-particle calculations are then used as boundary conditions for a ring current model. Our simulations reproduce the rapid decrease of Dst during the storm main phase and the fast initial phase of recovery. Shielding in the inner magnetosphere is established at early main phase. This shielding field lasts several hours and then breaks down at late main phase. At the peak of the storm, strong penetration of ions earthward to L shell of 1.5 is revealed in the simulation. It is surprising that O+ is significant but not the dominant species in the ring current in our calculation for this major storm. It is very likely that substorm effects are not well represented in the models and O+ energization is underestimated. Ring current simulation with O+ energy density at the boundary set comparable to Geotail observations produces excellent agreement with the observed symH. As expected in superstorms, ring current O+ is the dominant species over H+ during the main to mid-recovery phase of the storm.

Fok, Mei-Ching↗

Application of Powder Diffraction Methods to the Analysis of Short- and Long-Range Atomic Order in Nanocrystalline Diamond and SiC: The Concept of the Apparent Lattice Parameter (alp)

Two methods of the analysis of powder diffraction patterns of diamond and SiC nanocrystals are presented: (a) examination of changes of the lattice parameters with diffraction vector Q ('apparent lattice parameter', alp) which refers to Bragg scattering, and (b), examination of changes of inter-atomic distances based on the analysis of the atomic Pair Distribution Function, PDF. Application of these methods was studied based on the theoretical diffraction patterns computed for models of nanocrystals having (i) a perfect crystal lattice, and (ii), a core-shell structure, i.e. constituting a two-phase system. The models are defined by the lattice parameter of the grain core, thickness of the surface shell, and the magnitude and distribution of the strain field in the shell. X-ray and neutron experimental diffraction data of nanocrystalline SiC and diamond powders of the grain diameter from 4 nm up to micrometers were used. The effects of the internal pressure and strain at the grain surface on the structure are discussed based on the experimentally determined dependence of the alp values on the Q-vector, and changes of the interatomic distances with the grain size determined experimentally by the atomic Pair Distribution Function (PDF) analysis. The experimental results lend a strong support to the concept of a two-phase, core and the surface shell structure of nanocrystalline diamond and SiC.

Palosz, B.↗

The solar wind cycle, the sunspot cycle, and the corona.

Recent theories of the solar cycle and of coronal heating strongly suggest that solar cycle variations of different quantities (i.e., sunspots, coronal green line, etc.) ought not to be expected to be in phase with one another. In agreement with this notion it is noted that the shape of the corona typical of a 'maximum' eclipse occurs 1.5 yr before sunspot maximum, compared with 2 yr as might be expected from Leighton's 'standard' model. Further, it is argued that the phase of the solar wind cycle can be determined from geomagnetic observations. Using this phase, a solar cycle variation of 100 km/sec in the solar wind velocity and 1 gamma in the magnetic field intensity becomes apparent. In general, the solar wind cycle lags behind the coronal-eclipse-form cycle by 3 yr, compared with the 2 yr that might be expected from model calculations.

Hirshberg, J.↗

Three Dimensional Atmospheric Radiative Transfer-Applications and Methods Comparison

We review applications of 3D radiative transfer in the atmosphere, emphasizing the wide spectrum of scales important to remote sensing and modeling of cloud fields, and the characteristic scales introduced into observed radiances and fluxes by the distribution of photon pathlengths at conservative and absorbing wavelengths. We define the "plane-parallel bias", which is a measure of the importance of 3D cloud structure in large-scale models, and the "independent pixel errors" that quantify the significance of 3D effects in remote sensing, and emphasize their relative magnitude and scale dependence. A variety of approaches in current use in 3D radiative transfer, and issues of speed, accuracy, and flexibility are summarized. We also describe a recently initiated "International Intercomparison of 3-Dimensional Radiation Codes", or I3RC. I3RC is a 3-phase effort that has as its goals to: (1) understand the errors and limits of 3D methods; (2) provide "baseline" cases for future 3D code development; (3) promote sharing of 3D tools; (4) derive guidelines for 3D tool selection; and (5) improve atmospheric science education in 3D radiative transfer. Selected results from Phases 1 and 2 of I3RC are discussed. These are taken from five cloud fields: a 1D field of bar clouds, a 2D radar-derived field, a 3D Landsat-derived field, a stratiform cloud from the model of C. Moeng, and a convective cloud from the model of B. Stevens. Computations have been carried out for three monochromatic wavelengths (one conservative, one absorptive, and one thermal) and two solar zenith angles (0, 60 degrees).

Cahalan, Robert F.↗

Precipitation Processes Derived from TRMM Satellite Data, Cloud Resolving Model and Field Campaigns

Rainfall is a key link in the hydrologic cycle and is a primary heat source for the atmosphere. The vertical distribution of latent-heat release, which is accompanied by rainfall, modulates the large-scale circulations of the tropics and in turn can impact midlatitude weather. This latent heat release is a consequence of phase changes between vapor, liquid. and solid water. Present large-scale weather and climate models can simulate cloud latent heat release only crudely thus reducing their confidence in predictions on both global and regional scales. In this paper, NASA Tropical Rainfall Measuring (TRMM) precipitation radar (PR) derived rainfall information and the Goddard Convective and Stratiform Heating (CSH) algorithm used to estimate the four-dimensional structure of global monthly latent heating and rainfall profiles over the global tropics from December 1997 to October 2000. Rainfall latent heating and radar reflectively structure between ENSO (1997-1998 winter) and non-ENSO (1998-1999 winter) periods are examined and compared. The seasonal variation of heating over various geographic locations (i.e. Indian ocean vs west Pacific; Africa vs S. America) are also analyzed. In addition, the relationship between rainfall latent heating maximum heating level), radar reflectively and SST are examined.

Tao, W.-K.↗