Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrodynamics modeling”

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 127 records · Page 7

Passage of Bolides Through the Atmosphere

Different fragmentation models are applied to a number of events, including the entry of TC3 2008 asteroid in order to reproduce existing observational data. Keywords meteoroid entry fragmentation modeling 1 Introduction Fragmentation is a very important phenomenon which occurs during the meteoroid entry into the atmosphere and adds more drastic effects than mere deceleration and ablation. Modeling of bolide fragmentation (100 106 kg in mass) may be divided into several approaches. Detail fitting of observational data (deceleration and/or light curves) allows the determination of some meteoroid parameters (ablation and shape-density coefficients, fragmentation points, amount of mass loss) (Ceplecha et al. 1993; Ceplecha and ReVelle 2005). Observational data with high accuracy are needed for the gross-fragmentation model (Ceplecha et al. 1993), which is used for the analysis of European and Desert bolide networks data. Hydrodynamical models, which describe the entry of the meteoroid including evolution of its material, are applied mainly for large bodies (>106 kg) (Boslough et al. 1994; Svetsov et al. 1995; Shuvalov and Artemieva 2002, and others). Numerous papers were devoted to the application of standard equations for large meteoroid entry in the attempts to reproduce dynamics and/or radiation for different bolides and to predict meteorite falls. These modeling efforts are often supplemented by different fragmentation models (Baldwin and Sheaffer, 1971; Borovi.ka et al. 1998; Artemieva and Shuvalov, 2001; Bland and Artemieva, 2006, and others). The fragmentation may occur in different ways. For example, few large fragments are formed. These pieces initially interact through their shock waves and then continue their flight independently. The progressive fragmentation model suggests that meteoroids are disrupted into fragments, which continue their flight as independent bodies and may be disrupted further. Similar models were suggested in numerous papers, beginning with Levin (1956) and initial interaction of fragments started to be taken into account after the paper by Passey and Melosh (1980). The progressive fragmentation model with lateral spreading of formed fragments is widely used (Artemieva and Shuvalov, 1996; Nemtchinov and Popova, 1997; Borovi.ka et al. 1998; Bland and Artemieva, 2006).

Popova, O.↗

3D Modeling of Solar-Type Stars to Characterize Stellar Jitter

Detection of Earth-mass planets requires measurements of radial velocity with extreme precision. To capture the tiny disturbances caused by a planet's motion, it is necessary to understand and characterize the host star’s turbulent dynamics in order to apply proper filtering to the observational data. We take advantage of current computational and technological capabilities to develop 3D realistic models of the stellar subsurface convection and atmospheres and thereby estimate the photospheric jitter. We have identified an initial set of target stars, obtained initial conditions using the MESA code, and obtained initial 3D radiative models of the stellar surfaces and atmospheres with a spatial resolution of 50km. We present initial 3D radiative hydrodynamic model results of the planet-hosting star HD209458.

SMD↗

Scanning laser Doppler velocimeter system simulation for sensing aircraft wake vortices

A model that simulates the interaction of a laser-Doppler velocimeter with an aircraft wake flowfield is described. A hydrodynamic model is developed which represents the trailing vortex sheet and wind shear as discrete free vortices distributed over a two-dimensional grid. A sensor model is formulated for scanning both in range and in angle to produce a fan beam configuration without frequency translation. Output of this model is a frequency spectrum vs both range and angle. Once the spectrum is evaluated, simulations of the data analysis procedure are carried out. Patterns of various features of the signature are presented in range-elevation angle plots. The problem of locating the vortex centers is discussed as a pattern recognition problem and as a point target problem.

Thomson, J. A. L.↗

Cometary activity, discrete outgassing areas, and dust-jet formation

Conceptual models for various types of features observed in cometary comae (jets, spirals, halos, fans, etc.), their computer simulation, and the hydrodynamic models for jet formation are critically reviewed, and evidence for anisotropic, strongly collimated flows of ejecta emanating from discrete active regions (vents) on the rotating cometary nuclei is presented. Techniques employed to generate synthetic comet images that simulate the features observed are described, and their relevance to the primary objects of coma-morphology studies is discussed. Modeling of temporal variations in the water emission from discrete active regions suggests that production curves asymmetric with respect to perihelion should be commonplace. Critical comparisons with the activity profiles of Enke's comet and with light curves of disappearing comets and comets that undergo outbursts are presented. Recent developments in the understanding of the processes that cause the nongravitational perturbations of cometary motions are reviewed, and the observed discontinuities are identified with the birth of new sources and/or deactivation of old vents.

Sekanina, Z.↗

Nanoflare Heating During the Evolution of an AR

Nanoflares are thought to be prime candidates to heat the solar active regions. It is challenging to detect individual nanoflares with our present instrumentations. Determining the frequency and magnitude of the nanoflares are crucial to understand their contribution in coronal heating. In order to comprehend the role of nanoflares in coronal heating, researchers often combine the observed plasma emission with the coronal model. Here, by assuming nanoflare heating scenarios, we will study the evolution of an AR using the field-aligned hydrodynamic model. The frequency of heating will subsequently be determined by comparing the stimulated emission to the EUV and X-ray observation. We will continue to investigate how the heating frequency varies with the evolution of the AR. Further we will discuss the effect of observational parameters (e.g., exposure time, instrument energy range etc.) on the estimation of the nanoflare properties.

nano↗

Comparison of hydrodynamic and semi-kinetic treatments for plasma flow along closed field lines

Hydrodynamic and semi-kinetic treatments of plasma flow along closed geomagnetic field lines are compared. The hydrodynamic treatment is based on a simplified 16-moment set of transport equations as the equations for the heat flows are not solved; the heat flows are treated heuristically. The semi-kinetic treatment is based on a particle code. The comparison deals with the distributions of the plasma density, flow velocity, and parallel and perpendicular temperatures as obtained from the two treatments during the various stages of the flow. In the kinetic treatment, the appropriate boundary condition is the prescription of the velocity distribution functions for the particles entering the flux tubes at the ionospheric boundaries; those particles leaving the system are determined by the processes occurring in the flux tube. The prescribed distributions are half-Maxwellian with temperature T(sub 0) and density n(sub 0). In the hydrodynamic model, the prescribed boundary conditions are on density (n(sub 0)), flow velocity (V(sub 0)) and temperature (T(sub 0). It was found that results from the hydrodynamic treatment critically depend on V(sub 0); for early stages of the flow this treatment yields results in good agreement with those from the kinetic treatment, when V(sub 0) = square root of (kT(sub 0)/2 (pi)m), which is the average velocity of particles moving in a given direction for a Maxwellian distribution. During this early stage, the flows developing form the conjugate ionospheres show some distinct transitions. For the first hour or so, the flows are highly supersonic and penetrate deep into the opposite hemispheres, and both hydrodynamics and kinetic treatments yield almost similar features. It is found that during this period heatflow effects are negligibly small. When a flow penetrates deep into the opposite hemisphere, the kinetic treatment predicts reflection and setting up of counterstreaming. In contrast, the hydrodynamic treatment yields a shock in the flow. The reasons for this difference in the two treatments is discussed, showing that in view of the relatively warm ions, the coupling of ion beams and the consequent shock formation in the offequatorial region are not likely due to the enhancements in the beam temperatures. The counterstreaming in the kinetic treatment and the shock in the hydrodynamic treatment first advance upward to the equator and then downward to the ionospheric boundary from where the flow originated. The transit time for this advancement is found to be about 1 hour for the respective models. After 2 hours or so, both models predict that the flows from the ionospheric boundaries are generally subsonic with respect to the local ion-sound speed. At late stages of the flow, when a substantial fraction of ions entering the flux tube begin to return back in the kinetic treatment, the hydrodynamic treatment with the boundary condition V(sub 0) = square root of (kT(sub 0)/2(pi)m) yields an over-refilling, and the choice of V(sub 0) becomes uncertain.

Singh, Nagendra↗

Ionosphere and atmosphere of the moon in the geomagnetic tail

The paper presents calculations of the densities and energies of the various constituents of the lunar ionosphere during the time that the moon is in the geomagnetic tail; the surface concentrations of neon and argon are calculated from a theoretical model to be 3,900 and 1,700, respectively. It is found that a hydrostatic model of the ionospheric plasma is inadequate because the gravitational potential energy of the plasma is considerably smaller than its thermal energy. A hydrodynamic model, comparable to that used to describe the solar wind, is developed to obtain plasma densities and flow velocities as functions of altitude. The electromagnetic properties of the quiescent ionosphere are then investigated, and it is concluded that plasma effects on lunar induction can be neglected for quiescent conditions in the geomagnetic tail lobes.

Daily, W. D.↗

A Global Ocean Tide Model From TOPEX/POSEIDON Altimetry: GOT99.2

Goddard Ocean Tide model GOT99.2 is a new solution for the amplitudes and phases of the global oceanic tides, based on over six years of sea-surface height measurements by the TOPEX/POSEIDON satellite altimeter. Comparison with deep-ocean tide-gauge measurements show that this new tidal solution is an improvement over previous global models, with accuracies for the main semidiurnal lunar constituent M2 now below 1.5 cm (deep water only). The new solution benefits from use of prior hydrodynamic models, several in shallow and inland seas as well as the global finite-element model FES94.1. This report describes some of the data processing details involved in handling the altimetry, and it provides a comprehensive set of global cotidal charts of the resulting solutions. Various derived tidal charts are also provided, including tidal loading deformation charts, tidal gravimetric charts, and tidal current velocity (or transport) charts. Finally, low-degree spherical harmonic coefficients are computed by numerical quadrature and are tabulated for the major short-period tides; these are useful for a variety of geodetic and geophysical purposes, especially in combination with similar estimates from satellite laser ranging.

Ray, Richard D.↗

Chemical-dynamical models of the Venus mesosphere based upon diurnal microwave CO variations

An attempt is made to explain the fact that the large variation in Venus CO abundance with planetary phase using a combination of photochemical and kinematical models. Present knowledge about the Venus mesosphere is first summarized, emphasizing the more completely measured lower and upper atmosphere. The predictions of Dickinson and Ridley's (1977) hydrodynamic modeling are compared to the results of microwave observations of CO. It is noted that the nightside CO bulge above 90-95 km altitude indicated by microwave measurements is a primary characteristic of these models, whereas the opposite phase behavior of CO between 80 and 90 km is not predicted by them. The result of diurnal photochemical models for the Venus mesosphere are presented and an attempt is made to reproduce the phase behavior of CO between 80 and 90 km in the Venus mesosphere. Possible diurnal variations due to chemistry and vertical eddy diffusion are considered.

Clancy, R. T.↗

Absolute wind velocities in the lower thermosphere of Venus using infrared heterodyne spectroscopy

NASA's IR Telescope Facility and the McMath Solar Telescope have yielded absolute wind velocities in the Venus thermosphere for December 1985 to March 1987 with sufficient spatial resolution for circulation model discrimination. A qualitative analysis of beam-integrated winds indicates subsolar-to-antisolar circulation in the lower thermosphere; horizontal wind velocity was derived from a two-parameter model wind field of subsolar-antisolar and zonal components. A unique model fit common to all observing periods possessed 120 m/sec subsolar-antisolar and 25 m/sec zonal retrograde components, consistent with the Bougher et al. (1986, 1988) hydrodynamical models for 110 km.

Goldstein, Jeffrey J.↗

O VI Emission from the Interstellar Medium

This project's primary goal was to examine the Local Bubble, a large hot bubble surrounding the solar neighborhood. In order to do this, we observed the 1032 and 1038 A resonance line emission from O VI in the bubble and used the,results to comment on models for the Local Bubble and its embedded clouds. In order to maximize the signal to noise of our spectrum, we combined the awarded guest investigator observation with unpublished FUSE In Orbit Checkout observations. The resulting spectrum was sufficiently good as to enable us to place tight 2 sigma upper limits on the intensities of the 1032 and the 1038 A resonance lines. We also measured or placed upper limits on the other cosmic lines in the bandpass, including C III and C II. These are the first known ultraviolet emission line measurements and/or upper limits for the gas in the Local Bubble (as opposed to gas anywhere along long lines of sight). With the O VI upper limits, we were able to quantitatively evaluate competing theories for the origins of the Local Bubble. The upper limits are well below those expected in the Breitschwerdt model (which proposes that during its its early development, the Local Bubble rapidly expanded beyond its nascent cloud and, as a result, is now vastly underionized). The upper limits on the O VI resonance line doublet intensity and the measurement of the C III intensity, garnered from this project, combined with measurements of the O VI column density, garnered from another project, are so far below the predictions, that they make a good case for eliminating the Breitschwerdt model from the field of possibilities. Thus, instead of being vastly underionized, the Local Bubble is near ionizational equilibrium. In addition, the upper limits challenge the other well-known model for the Local Bubble. In that model, the Local Bubble was blown by a series of supernova explosions and winds and contains a myriad of evaporating clouds. The intensity of the O VI resonance line doublet predicted from models exceeds the observationally determined 2 sigma upper limits. The discrepancy prompts the exciting suggestion that there may be some astrophysics yet to be considered in standard hydrodynamic models. These measurements performed a second useful function. By subtract- ing them from the O VI intensity observed on long sight-lines through the Galactic halo and Local Bubble, we were able to estimate the intensity originating the Galactic halo.

Shelton, Robin L.↗

An Ellipsoidal Particle-Finite Element Method for Hypervelocity Impact Simulation

A number of coupled particle-element and hybrid particle-element methods have been developed for the simulation of hypervelocity impact problems, to avoid certain disadvantages associated with the use of pure continuum based or pure particle based methods. To date these methods have employed spherical particles. In recent work a hybrid formulation has been extended to the ellipsoidal particle case. A model formulation approach based on Lagrange's equations, with particles entropies serving as generalized coordinates, avoids the angular momentum conservation problems which have been reported with ellipsoidal smooth particle hydrodynamics models.

Shivarama, Ravishankar↗

Pressure-confined Lyman-alpha clouds

Results are presented of numerical models of pressure-confined spherical gas clouds which produce absorption resembling the low to intermediate atomic column density lines found in high-redshift QSO spectra. One-dimensional hydrodynamical models including electron conduction are described, and the rate equations are solved to find ionization and excitation states. Results are presented for both static and adiabatically expanding confining media covering a range of initial pressures. It is found that Ly-alpha lines are very similar over a wide range of conditions and that the most promising diagnostic of pressure is to compare the column density in H I to that in He I and He II. No single-pressure model can explain the wide range of observed H I column densities.

Baron, E.↗

Physics of heavy ions (1989-1990)

The results from studies on polar wind ion heating due to kinetic ion beam instabilities and the effects of such ion heating on the outflow of O(+) in the polar wind are presented and discussed. First, the linear instabilities associated with an O(+) and H(+) polar wind plasma in the presence of O(+) and H(+) beams for a range of O(+)/H(+) beam densities, T(sub e)/T(sub i), and ion beam speeds were examined. Then, nonlinear heating of the polar wind ions was studied, using numerical simulations. The O(+) and H(+) polar wind ions were modeled by isotropic Maxwellian distributions, and the electrons, O(+) beams, and H(+) beams were modeled by drifting Maxwellian distributions. The effects of the kinetic ion heating on the outflow of the polar wind ions were examined from the ionosphere, using a time-dependent hydrodynamic model. A numerical code to solve the O(+) and H(+) continuity and momentum equations in a flux tube from ionospheric to magnetospheric altitudes were developed. The effects of ion heating were included by allowing for the altitudinal variation of the ion temperatures in the momentum equation. The ion temperature profiles were specified based on the ion heating characteristics found from previous kinetic simulations. It was assumed that heating occurred above 1500 km and increased to a saturated value of temperature that was obtained directly from the kinetic simulation study. The characteristics of the dynamical polar wind without ion heating were studied, and a flux tube on closed field lines that suddenly became open at t = 0 was simulated. Then, the effects of ion heating were included. To gain some physical insight, two limiting cases were considered: preferential H(+) heating and preferential O(+) heating. How O(+) heating can lead to enhanced polar wind O(+) fluxes in the polar magnetosphere is shown.

Ashour-Abdalla, M.↗

Shock waves in a beta Doradus model

The accumulation of observational evidence and interpretation supports the conclusion that beta Doradus has running shock waves in its atmosphere. It is the work by Hutchinson with its predictions for the phases of the shock waves that prompted this author to construct a numerical hydrodynamical model for the atmosphere of beta Doradus. A short description of the procedure of the calculation follows. The results for the beta Doradus model contain a comparison with the observations. It is of particular interest to verify that the atmosphere contains multiple shock waves as suggested by Hutchinson for the ultraviolet observations.

Hill, S. J.↗

A High-Resolution X-Ray and Optical Study of SN1006: Asymmetric Expansion and Small-Scale Structure in a Type Ia Supernova Remnant

We introduce a deep (670 ks) X-ray survey of the entire SN 1006 remnant from the Chandra X-Ray Observatory, together with a deep Hα image of SN 1006 from the 4 m Blanco telescope at CTIO. Comparison with Chandra images from 2003 gives the first measurement of the X-ray proper motions around the entire periphery, carried out over a 9 yr baseline. We find that the expansion velocity varies significantly with azimuth. The highest velocity of approx.7400 km/s (almost 2.5 times that in the northwest (NW)) is found along the southeast (SE) periphery, where both the kinematics and the spectra indicate that most of the X-ray emission stems from ejecta that have been decelerated little, if at all. Asymmetries in the distribution of ejecta are seen on a variety of spatial scales. Si-rich ejecta are especially prominent in the SE quadrant, while O and Mg are more uniformly distributed, indicating large-scale asymmetries arising from the explosion itself. Neon emission is strongest in a sharp filament just behind the primary shock along the NWrim, where the pre-shock density is highest. Here the Ne is likely interstellar, while Ne within the shell may include a contribution from ejecta. Within the interior of the projected shell we find a few isolated "bullets" of what appear to be supernova ejecta that are immediately preceded by bowshocks seen in Hα, features that we interpret as ejecta knots that have reached relatively dense regions of the surrounding interstellar medium, but that appear in the interior in projection. Recent three-dimensional hydrodynamic models for Type Ia supernovae display small-scale features that strongly resemble the ones seen in X-rays in SN 1006; an origin in the explosion itself or from subsequent hydrodynamic instabilities both remain viable options. We have expanded the search for precursor X-ray emission ahead of a synchrotron-dominated shock front, as expected from diffusive shock acceleration theory, to numerous regions along both the northeast and southwest rims of the shell. Our data require that a precursor be thinner than about 3, and fainter than about 5% of the post-shock peak. These limits suggest that the magnetic field is amplified by a factor of seven or more in a narrow precursor region, promoting diffusive particle acceleration.

galaxies: active↗

Modeling of dense reactive granular flows

Because of the lack of fundamental knowledge of both the fluid mechanical and thermal behavior of multiphase flow systems, the design of these processes has been necessarily based on empirical correlations and experiments performed in laboratory or pilot scale units. Detailed models and simulations of these systems would allow for the optimization of the processes and their scaling, without the need for expensive and time consuming testing. In this paper, we present the proposed hydrodynamic model in detail.

granular↗

Constraining Nanoflare Frequency Through Observations and Models

Nanoflares are considered prime candidates for heating the non-flaring solar corona. However, direct detection of individual nanoflares with present-generation instruments remains challenging. Understanding the frequency and magnitude of nanoflares is crucial for comprehending their role in coronal heating. In this study, we employ field-aligned hydrodynamic model to simulate the emission from a non-flaring X-ray bright point (XBP) observed by the Marshall Grazing Incidence X-ray Spectrometer (MaGIXS). The length and magnetic field strength of the coronal loops associated with this XBP are derived from magnetic field extrapolation of observed photospheric magnetograms by SDO/HMI. Each loop is assumed to be heated by random nanoflares, whose frequency and magnitude are determined by the loop length and their field strength. The simulation results are then compared and matched against the observation to determine the nanoflare heating frequency. We will discuss the findings of this study. Furthermore, we will discuss the evolution of nanoflare heating frequency during the formation and decay of an active region (AR) on the solar disk.

coronalheating↗