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 307 records · Page 17

Specialists Meeting on Rotorcraft Dynamics, Moffett Field, Calif., February 13-15, 1974, Proceedings

Analysis of specific problems in rotorcraft dynamics. Topics include hingeless rotor theory, dynamic stall modelling, periodic systems identification, analysis of complex systems with phasing matrices, flapping stability, flap-lag dynamics at high advance ratios, finite element analysis and fuselage free-vibration characteristics, coupled rotor/frame vibration methods, gust response characteristics with unsteady stall effects, antiresonance theory, cyclic feathering motions and dynamic loads, control load envelope shaping, rotor aeroelasticity, use of Floquet theory, theory of proprotors and tilt-rotors, two-bladed teetering rotors, stability of air and ground resonance, vertical-plane pendulum absorbers, multicyclic jet-flap control, engine/frame interface dynamics, and others. The minutes of the question and answer periods following the presentations are presented in the supplement. Individual items are announced in this issue.

Source record↗

Preliminary analysis of University of North Dakota aircraft data from the FIRE Cirrus IFO-2

The stated goals of the First ISCCP (International Satellite Cloud Climatology Project) Regional Experiment (FIRE) are 'to promote the development of improved cloud and radiation parameterization for use in climate models, and to provide for assessment and improvement of ISCCP projects'. FIRE Phase 2 has focused on the formation, maintenance and dissipation of cirrus and marine stratocumulus cloud systems. These objectives have been approached through a combination of modeling, extended-time observations and intensive field observation (IFO) periods. The work under this grant was associated with the FIRE Cirrus IFO 2. This field measurement program was conducted to obtain observations of cirrus cloud systems on a range of scales from the synoptic to the microscale, utilizing simultaneous measurements from a variety of ground-based, satellite and airborne platforms. By combining these remote and in situ measurements a more complete picture of cirrus systems can be obtained. The role of the University of North Dakota in Phase 2 was three-fold: to collect in situ microphysical data during the Cirrus IFO 2; to process and archive these data; and to collaborate in analyses of IFO data. This report will summarize the activities and findings of the work performed under this grant; detailed description of the data sets available and of the analyses are contained in the Semi-annual Status Reports submitted to NASA.

Poellot, Michael R.↗

Solar Cycle Variations and Equatorial Oscillations: Modeling Study

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

Mayr, H. G.↗

Flow interaction experiment. Volume 1: Aerothermal modeling, phase 2

An experimental and computational study is reported for the flow of a turbulent jet discharging into a rectangular enclosure. The experimental configurations consisting of primary jets only, annular jets only, and a combination of annular and primary jets are investigated to provide a better understanding of the flow field in an annular combustor. A laser Doppler velocimeter is used to measure mean velocity and Reynolds stress components. Major features of the flow field include recirculation, primary and annular jet interaction, and high turbulence. A significant result from this study is the effect the primary jets have on the flow field. The primary jets are seen to create statistically larger recirculation zones and higher turbulence levels. In addition, a technique called marker nephelometry is used to provide mean concentration values in the model combustor. Computations are performed using three levels of turbulence closures, namely k-epsilon model, algebraic second moment (ASM), and differential second moment (DSM) closure. Two different numerical schemes are applied. One is the lower-order power-law differencing scheme (PLDS) and the other is the higher-order flux-spline differencing scheme (FSDS). A comparison is made of the performance of these schemes. The numerical results are compared with experimental data. For the cases considered in this study, the FSDS is more accurate than the PLDS. For a prescribed accuracy, the flux-spline scheme requires a far fewer number of grid points. Thus, it has the potential for providing a numerical error-free solution, especially for three-dimensional flows, without requiring an excessively fine grid. Although qualitatively good comparison with data was obtained, the deficiencies regarding the modeled dissipation rate (epsilon) equation, pressure-strain correlation model, and the inlet epsilon profile and other critical closure issues need to be resolved before one can achieve the degree of accuracy required to analytically design combustion systems.

Nikjooy, M.↗

Flow interaction experiment. Volume 2: Aerothermal modeling, phase 2

An experimental and computational study is reported for the flow of a turbulent jet discharging into a rectangular enclosure. The experimental configurations consisting of primary jets only, annular jets only, and a combination of annular and primary jets are investigated to provide a better understanding of the flow field in an annular combustor. A laser Doppler velocimeter is used to measure mean velocity and Reynolds stress components. Major features of the flow field include recirculation, primary and annular jet interaction, and high turbulence. A significant result from this study is the effect the primary jets have on the flow field. The primary jets are seen to create statistically larger recirculation zones and higher turbulence levels. In addition, a technique called marker nephelometry is used to provide mean concentration values in the model combustor. Computations are performed using three levels of turbulence closures, namely k-epsilon model, algebraic second moment (ASM), and differential second moment (DSM) closure. Two different numerical schemes are applied. One is the lower-order power-law differencing scheme (PLDS) and the other is the higher-order flux-spline differencing scheme (FSDS). A comparison is made of the performance of these schemes. The numerical results are compared with experimental data. For the cases considered in this study, the FSDS is more accurate than the PLDS. For a prescribed accuracy, the flux-spline scheme requires a far fewer number of grid points. Thus, it has the potential for providing a numerical error-free solution, especially for three-dimensional flows, without requiring an excessively fine grid. Although qualitatively good comparison with data was obtained, the deficiencies regarding the modeled dissipation rate (epsilon) equation, pressure-strain correlation model, and the inlet epsilon profile and other critical closure issues need to be resolved before one can achieve the degree of accuracy required to analytically design combustion systems.

Nikjooy, M.↗

An analysis of the annual cycle in upper stratospheric ozone

The mid-latitude upper stratospheric ozone profiles obtained by the solar backscatter ultraviolet instrument on the Nimbus 7 satellite show a clear annual cycle both in the absolute ozone amounts between 0.98 and 15.6 mbar and in the magnitude of disturbances that reveal themselves as longitudinal structure. At the lowest pressure analyzed a winter maximum in ozone exists, but as one progresses downward in altitude a shift in the temporal phase of the annual cycle occurs in the vicinity of 3 to 4 mbar. Comparison of the observed behavior with the predictions of a one-dimensional photochemical model shows a systematic tendency for calculated ozone amount to be 20-27 percent below the data for pressures less than 7.8 mbar. The chemical model successfully predicts the change in phase of the annual cycle, although at a pressure greater than observed. Diagnosis of model results shows the observed shift to be closely coupled to the magnitude of the ozone column density near 3-4 mbar. The wavelength-dependent attenuation of the solar radiation field by ozone alters the relative magnitude of the molecular oxygen and ozone dissociation rates, leading to a change in the temporal phase of the annual cycle.

Frederick, J. E.↗

Longwave radiation effects of the Harmattan haze

Infrared in situ radiance observations at 8.0 to 14 microns and 9.5 to 11.5 microns of the West African Harmattan haze during the 1974 Global Atmospheric Research Project Atlantic Tropical Experiment field phase made possible the determination of some of the radiative properties of this tropospheric phenomenon. This in turn permitted development of a simple calculation model for radiative transfer through the haze. Radiometric observations of the dust haze, reaching from 600 m to 6.25 km, were analyzed for haze IR transmission. A transfer model incorporating these transmission properties gave an average calculated IR cooling rate of 0.09 C/h for the entire haze layer as compared to a haze-free cloudless-troposphere cooling rate of 0.06 C/h for the same levels. The haze volume-absorption coefficient was approximately 0.042/km for layers of all depths. This uniformity of the haze transmission was further evident in the direct correlation of its transmission and optical depth.

Kuhn, P. M.↗

Effect of particle presence on the incompressible inviscid flow through a two dimensional compressor cascade

A computational model for calculating the flow field, in the presence of solid particles, through a two dimensional compressor cascade is presented. The method treats the particle phase in the Lagrangian system and the fluid phase in the Eulerian system. The equations of momentum of the fluid phase are modified to account for the momentum exchange between the two phases. The resulting modified momentum and continuity equations are reduced to the conventional stream function-vorticity formulation. The analysis yields the change in the blade surface pressure distribution, the total pressure and velocity in the flow field. Based on this, it is possible to predict the decrease in compressor performance under the presence of particles.

Balan, C.↗

Radial diffusion and losses of energetic protons in the 5 to 12 R(S) region of Saturn's magnetosphere

Phase space density profiles for nearly equatorial mirroring energetic protons, derived from Pioneer 11 measurements during the 1979 Saturn encounter, were analyzed using solutions of the time-averaged radial diffusion equation in a dipolar planetary magnetic field. Then, several loss models were considered, ranging from a model that takes into account satellite absorption only to the one that accounts for satellite plus ring-E absorption plus added distributed losses, and the corresponding form of the time-averaged radial diffusion coefficient was determined for each case by a minimum-variance fit to the data-derived phase space densiy profiles. The model comparisons clearly indicated the need for distributed losses, in addition to satellite absorption, in order to explain the smooth decline of the experimental energetic proton phase space density profiles in the 5-12-Saturn radius region of Saturn's magnetosphere.

Hood, L. L.↗

First evaluation of the main parameters in the dynamics of the Small Expendable Deployer System (SEDS) for a tethered satellite

The dynamics of the motion of the Small Expendable Deployer System (SEDS) is studied by using a simplified model in which no external forces (except the gravity gradient field) are applied on the tethered body and the tether is assumed massless. The dynamics of SEDS operation is modeled as a sequence of two phases: the deployment phase and the swing phase. For the first one the velocity dependent forces are found to force the tether forward from the local vertical. When the deployment ends, Coriolis effects vanish and the swing phase begins, which is characterized by a wide free libration. The time duration as well as velocity, acceleration and tension of the tethered body are estimated for both deployment and swing phases.

Deluca, Luigi↗

Summary of completed project

The primary goal of the work in this proposal has been the determination of the requirement for acceleration of particles, especially electrons, during solar flares. This is one of the important and yet unsolved problems in solar flares and other astronomical objects. This is done by detailed theoretical investigation of interaction of accelerated particles with ambient plasma including interaction with particles, waves, and large scale magnetic fields. Comparison with observed impulsive phase radiation such as X-rays, gamma-rays, microwaves, and type III radiation is used for constraining the models and for determination of the characteristics of the accelerated particles. Steady state and time-dependent numerical codes have been developed based on the Fokker-Planck method for study of the particle acceleration and transport during a flare. Using these, the correlation between characteristics of radiation at various wavelengths have been investigated and some constraints have been set on the model parameters. During the last three year period of this grant, dissertation research of three graduate students have been partially supported by this grant. Each of these theses has dealt with various observed properties of flare emission with particular emphasis on the interrelationship between radiation in at least two energy bands: X-rays and gamma-rays, X-rays and microwaves, and X-rays and type III radiation. This has allowed considerable tightening of constraints on models, and inclusion of wave particle interaction has allowed direct calculations of acceleration processes.

Petrosian, Vahe↗

Upward Flame Spread Over Thin Solids in Partial Gravity

The effects of partial-gravity, reduced pressure, and sample width on upward flame spread over a thin cellulose fuel were studied experimentally and the results were compared to a numerical flame spread simulation. Fuel samples 1-cm, 2-cm, and 4-cm wide were burned in air at reduced pressures of 0.2 to 0.4 atmospheres in simulated gravity environments of 0.1-G, 0.16-G (Lunar), and 0.38-G (Martian) onboard the NASA KC-135 aircraft and in normal-gravity tests. Observed steady flame propagation speeds and pyrolysis lengths were approximately proportional to the gravity level. Flames spread more quickly and were longer with the wider samples and the variations with gravity and pressure increased with sample width. A numerical simulation of upward flame spread was developed including three-dimensional Navier-Stokes equations, one-step Arrhenius kinetics for the gas phase flame and for the solid surface decomposition, and a fuel-surface radiative loss. The model provides detailed structure of flame temperatures, the flow field interactions with the flame, and the solid fuel mass disappearance. The simulation agrees with experimental flame spread rates and their dependence on gravity level but predicts a wider flammable region than found by experiment. Some unique three-dimensional flame features are demonstrated in the model results.

Feier, I. I.↗

Highly Variable Cycle Nozzle Concept: Validation of Flow and Noise Predictions

Results from experimental and numerical studies of highly Variable Cycle (HVC) exhaust model were presented. The model was designed and fabricated under a Supersonics NRA awarded to Rolls-Royce. The model had a lobed mixer for the core stream nozzle, and elliptic fan stream nozzle, and an ejector. Experiments included far-field acoustic array, phased array, and Particle Image Velocimetry (PIV) measurements. Numerical studies included flow simulations using the WIND-US code and far-field acoustic solutions using an acoustic analogy developed by Goldstein (2003) and Leib and Goldstein (2011). Far-field acoustic measurements showed increased noise levels over the round baseline nozzle when using non-static forward flight conditions. Phased array measurements showed noise sources near the ejector doors when tones were produced for small ejector door positions. Ejector door separation identified in the experiments was reproduced in the numerical flow simulations. Acoustic solutions were unable to match levels measured in the peak jet noise direction indicating additional development work is needed to predict noise from highly three-dimensional flows.

Halbig, Michael C.↗

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.

Gilson, M. L.↗

Data Assimilation and Uncertainties in Early Solar Cycle Predictions

Stochastic nature of solar activity variations together with our limited knowledge of the dynamo mechanism and subsurface dynamics causes uncertainty in predictions of the solar cycle. For improving the physics-­‐based predictions we can take advantage of the mathematical data assimilation approach that allows us to take into account both, observational errors and model uncertainties, and provide estimates of the next solar cycle along with prediction uncertainties. In this study we use the Parker's migratory dynamo model together with the equation of magnetic helicity balance, which reproduces main properties of the sunspot cycles and allow us to minimize discrepancies between the observed global activity variations and the model solution. The test simulation runs show that a reliable prediction can be obtained for two phases of preceding solar cycle: 1) if the polar field reversals shortly after the solar maxima (strong toroidal field and weak poloidal field), and 2) during the solar minima (strongest poloidal and weak toroidal fields). The early estimate of Cycle 25 obtained by this method shows that this cycle will start in 2019 – 2020, reach the maximum in 2023 – 2024, and that the mean sunspot number at the maximum will be about 90 (for the v2.0 sunspot number series).

Data Assimilation↗

Analytical and experimental studies of an optimum multisegment phased liner noise suppression concept

Results are presented from detailed analytical studies made to define methods for obtaining improved multisegment lining performance by taking advantage of relative placement of each lining segment. Properly phased liner segments reflect and spatially redistribute the incident acoustic energy and thus provide additional attenuation. A mathematical model was developed for rectangular ducts with uniform mean flow. Segmented acoustic fields were represented by duct eigenfunction expansions, and mode-matching was used to ensure continuity of the total field. Parametric studies were performed to identify attenuation mechanisms and define preliminary liner configurations. An optimization procedure was used to determine optimum liner impedance values for a given total lining length, Mach number, and incident modal distribution. Optimal segmented liners are presented and it is shown that, provided the sound source is well-defined and flow environment is known, conventional infinite duct optimum attenuation rates can be improved. To confirm these results, an experimental program was conducted in a laboratory test facility. The measured data are presented in the form of analytical-experimental correlations. Excellent agreement between theory and experiment verifies and substantiates the analytical prediction techniques. The results indicate that phased liners may be of immediate benefit in the development of improved aircraft exhaust duct noise suppressors.

Sawdy, D. T.↗

Quantitative modeling of the 'convection surge' mechanism of ion acceleration

To explore the properties of the 'convection surge' mechanism postulated by Quinn and Southwood (1982), a quantitative computer model of the mechanism is constructed. The model incorporates the dipolarization of the field line shapes as well as the results of recent transient electric field measurements by the AST-6 satellite. The computer model numerically integrates the equations of motion of ensembles of ions populating a single flux tube that goes through the convection surge or 'dipolarization' transformation. The numerical results confirm that the mechanism easily generates the bounce-phase-bunched ion distributions in question, and show that the convection surge can generate dramatically field-aligned ion distributions. It is hypothesized that the convection surge mechanism is fundamentally associated with the processes responsible for transporting tail ion populations to the geosynchronous regions of the earth's magnetosphere.

Mauk, B. H.↗

A Three-Dimensional MHD Simulation of the Solar Wind for a Tilted-Dipole Magnetic Field on the Sun

Using a three-dimensional MHD model, we simulate the global steady-state structure of the solar corona and solar wind for a dipole magnetic field on the Sun inclined by 30 degrees to the solar rotation axis. This represents the solar conditions typical for a declining phase of solar cycle. The computations can extend from the coronal base out to 100-AU and at large heliospheric distances includes the effects of interstellar neutral hydrogen and their interaction with solar wind protons. The simulations can model the formation of corotating interaction regions and the heliospheric current sheet. The simulations are also capable of describing very strong rarefaction regions that include embedded sub-Alfvenic regions that form on the trailing edge of a fast flows.

Goldstein, Melvyn L.↗