Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuum 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 1,171 records · Page 65

Physical Models for Accreting Pulsars at High Luminosity

A new window for better understanding the accretion onto strongly magnetized neutron stars in X-ray binaries is opening. In these systems the accreted material follows the magnetic field lines as it approaches the neutron star, forming accretion columns above the magnetic poles. The plasma falls toward the neutron star surface at near-relativistic speeds, losing energy by emitting X-rays. The X-ray spectral continua are commonly described using phenomenological models, i.e., power laws with different types of curved cut-offs at higher energies. Here we consider high luminosity pulsars. In these systems the mass transfer rate is high enough that the accreting plasma is thought to be decelerated in a radiation-dominated radiative shock in the accretion columns. While the theory of the emission from such shocks had already been developed by 2007, a model for direct comparison with X-ray continuum spectra in xspec or isis has only recently become available. Here we analyze the broadband X-ray spectra of the accreting pulsars Centaurus X-3 and 4U1626-67 obtained withNuSTAR. We present results from traditional empirical modeling as well as successfully apply the radiation-dominated radiative shock model. We also fit the energy-dependent pulse profiles of 4U 1626-67 using a new relativistic light bending model.

Pottschmidt, Katja↗

Detection of Photospheric Features in the Near-Infrared Spectrum of a Class 0 Protostar

We present a near-infrared K-band R ≃ 1500 Keck spectrum of S68N, a Class 0 protostar in the Serpens molecular cloud. The spectrum shows a very red continuum, CO absorption bands, weak or non-existent atomic metal absorptions, and H2 emission lines. The near-IR H2 emission is consistent with excitation in shocks or by X-rays but not by UV radiation. We model the absorption component as a stellar photosphere plus circumstellar continuum emission with wavelength-dependent extinction. A Markov Chain Monte Carlo analysis shows that the most likely model parameters are consistent with a low-temperature, low-gravity photosphere with signifi cant extinction and no more than modest continuum veiling. Its T(eff) ≃ 3260 K effective temperature is similar to that of older, more evolved pre-main-sequence stars, but its surface gravity log g ≃ 2:4 cm/sq.s is approximately 1 dex lower. This implies that the radius of this protostar is a factor of ~ 3 larger than that of 106 yr old T Tauri stars. Its low veiling is consistent with a circumstellar disk having intrinsic near-IR emission that is less than or equal to that of more evolved Class I protostars. Along with the high extinction, this suggests that most of the circumstellar material is in a cold envelope, as expected for a Class 0 protostar. This is the first known detection and analysis of a Class 0 protostar absorption spectrum.

Thomas P Greene↗

Beam-generated upper hybrid noise in Jupiter's outer magnetosphere

A model for generation of upper hybrid waves in Jupiter's outer magnetosphere is presented. Energetic electrons accelerated on high-latitude auroral field lines at 1 R(J) altitudes are assumed to stream outwards to the distant magnetosphere. The large decrease in the local magnetic field strength results in a highly collimated, field-aligned electron beam due to conservation of the particle's magnetic moment. As the beam travels outward it passes through the plasma sheet boundary layer region located on the edges of the plasma sheet. Because of the large field-aligned anisotropy of the beam, the distribution is unstable to the excitation of the upper hybrid waves which amplify and then undergo mode conversion at the f(p) layer to generate the Jovian continuum radiation. A novel feature of the model is the beam-anisotropic heat flux instability which drives UH waves unstable without the need for any positive slope on the electron distribution function.

Barbosa, D. D.↗

Molecular Hydrogen in the Quiescent Disk of SW UMA

The data have been acquired and re-reduced by Alex Lobel. We do not detect the molecular hydrogen transitions expected from models of a molecular disk illuminated by Ly alpha photons. We are currently working on an interpretation of the O VI emission lines to determine whether they arise in the disk or the boundary layer. The upper limit on continuum flux will provide an important constraint on models of white dwarf heating in cataclysmic variables.

Raymond, John C.↗

Peridynamics and surrogate modeling of pressure-driven well stimulation

In this work we use the peridynamics theory of solid mechanics to simulate fracture in an annular rock domain subject to an in-situ stress and create surrogate models that predict the area of the resulting cracks. Peridynamics is a non-local formulation of continuum mechanics that naturally accommodates material discontinuities. Furthermore, unlike other fracture modeling techniques there is no need to provide information about the crack path. We utilize the peridynamics code Peridigm and take a two-stage approach to fracture modeling. First an implicit solve is performed to compute the in-situ stress state. We then execute an explicit solve where a pressure loading designed to emulate fluid-driven hydraulic fracture is applied at the borehole and transmitted to the pre-stressed rock. We present results from polynomial and single and multi-level Gaussian process surrogate models constructed from a sampling study of the peridynamics model. The surrogates predict crack area given a measure of the in-situ stress anisotropy and rise time and amplitude of the pressure loading. These surrogates take a minuscule fraction of peridynamics model's running time to evaluate and are a step towards enabling advanced optimization and uncertainty quantification workflows that require many model evaluations.

42 ENGINEERING↗

Gravitational Acceleration Effects on Macrosegregation: Experiment and Computational Modeling

Experiments were performed under terrestrial gravity (1g) and during parabolic flights (10-2 g) to study the solidification and macrosegregation patterns of Al-Cu alloys. Alloys having 2% and 5% Cu were solidified against a chill at two different cooling rates. Microscopic and Electron Microprobe characterization was used to produce microstructural and macrosegregation maps. In all cases positive segregation occurred next to the chill because shrinkage flow, as expected. This positive segregation was higher in the low-g samples, apparently because of the higher heat transfer coefficient. A 2-D computational model was used to explain the experimental results. The continuum formulation was employed to describe the macroscopic transports of mass, energy, and momentum, associated with the solidification phenomena, for a two-phase system. The model considers that liquid flow is driven by thermal and solutal buoyancy, and by solidification shrinkage. The solidification event was divided into two stages. In the first one, the liquid containing freely moving equiaxed grains was described through the relative viscosity concept. In the second stage, when a fixed dendritic network was formed after dendritic coherency, the mushy zone was treated as a porous medium. The macrosegregation maps and the cooling curves obtained during experiments were used for validation of the solidification and segregation model. The model can explain the solidification and macrosegregation patterns and the differences between low- and high-gravity results.

Leon-Torres, J.↗

The X-ray spectrum of SN 1006

The supernova (SN) of A.D. 1006 represents the brightest in recorded history. The X-ray observations of its remnant appear to be full of puzzles. The present paper has the objective to reevaluate published X-ray data on SN 1006, using detailed nonequilibrium spectral computations and a self-consistent hydrodynamic SNR model. In a summary of the main conclusions, it is pointed out that an observed power-law continuum can be understood as mainly free-free emission from a reverse shock into SN ejecta with an initially uniform density profile. The spectral index is correctly predicted by such a model. Most of the continuum is produced by fully ionized gas in the dense outer layers of ejecta. The X-ray spectrum is not consistent with a solar composition or a uniform composition but is consistent with a layered composition of ejecta containing carbon on the outside, and processed heavy elements on the inside. There is no difficulty in including a large mass of iron.

Hamilton, A. J. S.↗

Spatially resolved methane band photometry of Jupiter. III - Cloud vertical structures for several axisymmetric bands and the Great Red Spot

The paper presents cloud structure models for Jupiter's Great Red Spot, Equatorial and North Tropical Zones, North and South Temperate Zones, and North and South Polar Regions. The models are based on images of Jupiter in three methane bands and nearby continuum; radiative transfer calculations include multiple scattering and absorption from three aerosol layers. The model results include the transition in the upper-cloud altitude to 3 km lower altitude from the tropical zones to temperate zones and polar regions, a N/S asymmetry in cloud thickness in the tropical and temperature zones, and the presence of aerosols up to about 0.3 bar in the Great Red Spot and Equatorial Zone. It is concluded that polarization data are sensitive to aerosols in and above the upper cloud layer but insensitive to deeper cloud structure.

West, R. A.↗

Hybrid MD-Nernst Planck Model of Alpha-hemolysin Conductance Properties

Motivated by experiments in which an applied electric field translocates polynucleotides through an alpha-hemolysin protein channel causing ionic current transient blockade, a hybrid simulation model is proposed to predict the conductance properties of the open channel. Time scales corresponding to ion permeation processes are reached using the Poisson-Nemst-Planck (PNP) electro-diffusion model in which both solvent and local ion concentrations are represented as a continuum. The diffusion coefficients of the ions (K(+) and Cl(-)) input in the PNP model are, however, calculated from all-atom molecular dynamics (MD). In the MD simulations, a reduced representation of the channel is used. The channel is solvated in a 1 M KCI solution, and an external electric field is applied. The pore specific diffusion coefficients for both ionic species are reduced 5-7 times in comparison to bulk values. Significant statistical variations (17-45%) of the pore-ions diffusivities are observed. Within the statistics, the ionic diffusivities remain invariable for a range of external applied voltages between 30 and 240mV. In the 2D-PNP calculations, the pore stem is approximated by a smooth cylinder of radius approx. 9A with two constriction blocks where the radius is reduced to approx. 6A. The electrostatic potential includes the contribution from the atomistic charges. The MD-PNP model shows that the atomic charges are responsible for the rectifying behaviour and for the slight anion selectivity of the a-hemolysin pore. Independent of the hierarchy between the anion and cation diffusivities, the anionic contribution to the total ionic current will dominate. The predictions of the MD-PNP model are in good agreement with experimental data and give confidence in the present approach of bridging time scales by combining a microscopic and macroscopic model.

Cozmuta, Ioana↗

Mesoscale Computational Investigation of Shocked Heterogeneous Materials with Application to Large Impact Craters

The propagation of shock waves through target materials is strongly influenced by the presence of small-scale structure, fractures, physical and chemical heterogeneities. Pre-existing fractures often create craters that appear square in outline (e.g. Meteor Crater). Reverberations behind the shock from the presence of physical heterogeneity have been proposed as a mechanism for transient weakening of target materials. Pre-existing fractures can also affect melt generation. In this study, we are attempting to bridge the gap in numerical modeling between the micro-scale and the continuum, the so-called meso-scale. To accomplish this, we are developing a methodology to be used in the shock physics hydrocode (CTH) using Monte-Carlo-type methods to investigate the shock properties of heterogeneous materials. By comparing the results of numerical experiments at the micro-scale with experimental results and by using statistical techniques to evaluate the performance of simple constitutive models, we hope to embed the effect of physical heterogeneity into the field variables (pressure, stress, density, velocity) allowing us to directly imprint the effects of micro-scale heterogeneity at the continuum level without incurring high computational cost.

Crawford, D. A.↗

Effect of tropospheric aerosols upon atmospheric infrared cooling rates

An investigation has been made of the impact of wind-blown dust particles upon local climate of arid regions. The case of Northwest India is specifically considered, where a dense layer of dust persists for several months during the summer. In order to examine the effect of this dust layer on the infrared radiative flux and cooling rates, a method is presented for calculating the IR flux within a dusty atmosphere which allows the use of gaseous band models and is applicable in the limit of small single scattering albedo and pronounced forward scattering. The participating components of the atmosphere are assumed to be water vapor and spherical quartz particles only. The atmospheric window is partially filled by including the water vapor continuum bands for which empirically obtained transmission functions have been used. It is shown that radically different conclusions may be drawn on dust effects if the continuum absorption is not considered. The radiative transfer model, when applied to a dusty atmosphere, indicates that there is a moderate enhancement in the atmospheric greenhouse and a 10% increase in the mean IR radiative cooling rate, relative to the dust free case, within the lower troposphere. These results have been compared with previous work by other authors in the context of the possibility of dust layers inhibiting local precipitation.

Harshvardhan, MR.↗

HPTAM, a two-dimensional Heat Pipe Transient Analysis Model, including the startup from a frozen state

A two-dimensional Heat Pipe Transient Analysis Model, 'HPTAM,' was developed to simulate the transient operation of fully-thawed heat pipes and the startup of heat pipes from a frozen state. The model incorporates: (a) sublimation and resolidification of working fluid; (b) melting and freezing of the working fluid in the porous wick; (c) evaporation of thawed working fluid and condensation as a thin liquid film on a frozen substrate; (d) free-molecule, transition, and continuum vapor flow regimes, using the Dusty Gas Model; (e) liquid flow and heat transfer in the porous wick; and (f) thermal and hydrodynamic couplings of phases at their respective interfaces. HPTAM predicts the radius of curvature of the liquid meniscus at the liquid-vapor interface and the radial location of the working fluid level (liquid or solid) in the wick. It also includes the transverse momentum jump condition (capillary relationship of Pascal) at the liquid-vapor interface and geometrically relates the radius of curvature of the liquid meniscus to the volume fraction of vapor in the wick. The present model predicts the capillary limit and partial liquid recess (dryout) in the evaporator wick, and incorporates a liquid pooling submodel, which simulates accumulation of the excess liquid in the vapor core at the condenser end.

Tournier, Jean-Michel↗

Comparisons of Air Radiation Model with Shock Tube Measurements

This paper presents an assessment of the predictive capability of shock layer radiation model appropriate for NASA s Orion Crew Exploration Vehicle lunar return entry. A detailed set of spectrally resolved radiation intensity comparisons are made with recently conducted tests in the Electric Arc Shock Tube (EAST) facility at NASA Ames Research Center. The spectral range spanned from vacuum ultraviolet wavelength of 115 nm to infrared wavelength of 1400 nm. The analysis is done for 9.5-10.5 km/s shock passing through room temperature synthetic air at 0.2, 0.3 and 0.7 Torr. The comparisons between model and measurements show discrepancies in the level of background continuum radiation and intensities of atomic lines. Impurities in the EAST facility in the form of carbon bearing species are also modeled to estimate the level of contaminants and their impact on the comparisons. The discrepancies, although large is some cases, exhibit order and consistency. A set of tests and analyses improvements are proposed as forward work plan in order to confirm or reject various proposed reasons for the observed discrepancies.

Bose, Deepak↗

A radiative transfer model for surface radiation budget studies

A high resolution spectral radiative transfer model of the troposphere for computing downwelling radiance and flux density at the surface in the 5-200 micron region has been developed. The model is fast and accurate and takes into consideration all major and minor constituents active in this region, namely, water vapor, carbon dioxide, ozone, nitrous oxide, and methane. It uses the quasi-random band model of absorption for computing atmospheric transmittances due to various bands. Continuum absorption by water vapor in the 8-14 micron region has also been taken into account. The model has been used to examine the sensitivity of downwelling flux to variations of surface water vapor pressure (or surface relative humidity), water vapor scale-height and concentrations of other gaseous constituents of the atmosphere. In addition, sensitivity to changes in a number of cloud parameters, namely, fractional cloud cover, cloud height, and high-cloud emissivity, has also been examined. Results obtained with the present model compare well with those from the empirical and semi-empirical correlations. The present work demonstrates, for the first time, the feasibility of using a high-resolution narrow-band model for computing this flux for a meteorological application.

Gupta, S. K.↗

SUSIM/UARS observations of the 120 to 300 nm flux variations during the maximum of the solar cycle - Inferences for the 11-year cycle

Since October 1991, the Solar Ultraviolet Spectral Irradiance Monitor on board the Upper Atmosphere Research Satellite has been measuring the solar spectral irradiance from 120 to 400 nm with low (5 nm) and moderate (1.1 nm) resolution on a daily basis. By scaling the modulation measured over four solar rotations to proxy indicators of solar activity, we estimate the amplitude of the solar ultraviolet irradiance variation during the 11-year activity cycle, extending the wavelength coverage of existing empirical variability models to 120 nm and distinguishing the variability of the emission lines from that of the underlying solar continuum. These estimates are compared with results from another empirical variability model and from direct measurements.

Lean, Judith↗

Cometary compact H II regions are stellar-wind bow shocks

Comet-shaped H II regions, like G34.3 + 0.2, are easily explained as bow shocks created by wind-blowing massive stars moving supersonically through molecular clouds. The required velocities of the stars through dense clumps are less than about 10 km/s, comparable to the velocity dispersion of stars in OB associations. An analytic model of bow shocks matches the gross characteristics seen in the radio continuum and the velocity structure inferred from hydrogen recombination and molecular line observations. The champagne flow model cannot account for these structures. VLBI observations of masers associated with the shells of cometary compact H II regions should reveal tailward proper motions predominantly parallel to the shell, rather than perpendicular. It is predicted that over a decade baseline, high signal-to-noise VLA observations of this class of objects will show headward pattern motion in the direction of the symmetry axis, but not expansion. Finally, shock-generated and coronal infrared lines are also predicted.

Van Buren, Dave↗

Modeling of Wall-Bounded Complex Flows and Free Shear Flows

Various wall-bounded flows with complex geometries and free shear flows have been studied with a newly developed realizable Reynolds stress algebraic equation model. The model development is based on the invariant theory in continuum mechanics. This theory enables us to formulate a general constitutive relation for the Reynolds stresses. Pope was the first to introduce this kind of constitutive relation to turbulence modeling. In our study, realizability is imposed on the truncated constitutive relation to determine the coefficients so that, unlike the standard k-E eddy viscosity model, the present model will not produce negative normal stresses in any situations of rapid distortion. The calculations based on the present model have shown an encouraging success in modeling complex turbulent flows.

Shih, Tsan-Hsing↗

Extended quantitative rescattering model for simulating high-order harmonic streaking spectra by synchronization of an intense IR laser and a time-delayed attosecond XUV pulse

In this work, we theoretically investigate the modulated high-harmonic generation (HHG) driven by an intense few-cycle infrared (IR) laser field and a weak extreme-ultraviolet (XUV) pulse at a delayed time. We establish an extended quantitative rescattering (EQRS) model to simulate the HHG streaking spectra, with the ideas of correcting the IR ionization and the transition from the ground to continuum states in the strong-field approximation. The EQRS model has an accuracy comparable to that from “exactly” solving the time-dependent Schrödinger equation (TDSE). We reveal that the fringes in the streaking spectra are caused by the interference between the attosecond XUV pulse and harmonics resulting from different recombination pathways under the intense IR laser. We then demonstrate that the XUV pulse can be accurately retrieved by treating the single-atom TDSE results or macroscopic propagation results as the “input” data. This work provides with a tool for efficiently simulating and thoroughly analyzing the XUV-assisted HHG, which could also enhance its capability for tracing the electron dynamics involved in the strong-field phenomena.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗