Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computing continuum”

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 361 records · Page 20

Experimental investigation of a wall stablized arc with comparisons to theoretical predictions

A wall-stabilized, constricted-arc heater was investigated by extensive measurements of pressure, voltage, and wall heat flux at different axial locations along the constrictor. Radial temperature profiles were measured at one axial location by the radiance of a single line method and also by the continuum radiance at a wavelength of 4955 A. The arc heater with an inside diameter of 2.37 cm and a constrictor length of 42 cm was operated with nitrogen flow rate from about 2 to 10 g/sec, and are currents from 400 to 1200 amperes. Based on a literature search, two different digital computer programs were selected which were used to predict the pertinent parameters by numerical solutions to the governing differential equations. The program selections were based on completeness and availability. In general the predictions were adequate when laminar flow was assumed, and the greatest difficulty was encountered in the predictions of wall heat flux. The results of the measurements are shown in graphical form, along with the theoretical predictions from both computer programs.

Wells, W. L.↗

Boundary formulations for sensitivity analysis without matrix derivatives

A new hybrid approach to continuum structural shape sensitivity analysis employing boundary element analysis (BEA) is presented. The approach uses iterative reanalysis to obviate the need to factor perturbed matrices in the determination of surface displacement and traction sensitivities via a univariate perturbation/finite difference (UPFD) step. The UPFD approach makes it possible to immediately reuse existing subroutines for computation of BEA matrix coefficients in the design sensitivity analysis process. The reanalysis technique computes economical response of univariately perturbed models without factoring perturbed matrices. The approach provides substantial computational economy without the burden of a large-scale reprogramming effort.

Kane, J. H.↗

The Micromechanics of the Moving Contact Line

The proposed research is divided into three components concerned with molecular structure, molecular orientation, and continuum averages of discrete systems. In the experimental program, we propose exploring how changes in interfacial molecular structure generate contact line motion. Rather than rely on the electrostatic and electrokinetic fields arising from the molecules themselves, we augment their interactions by an imposed field at the solid/liquid interface. By controling the field, we can manipulate the molecular structure at the solid/liquid interface. In response to controlled changes in molecular structure, we observe the resultant contact line motion. In the analytical portion of the proposed research we seek to formulate a system of equations governing fluid motion which accounts for the orientation of fluid molecules. In preliminary work, we have focused on describing how molecular orientation affects the forces generated at the moving contact line. Ideally, as assumed above, the discrete behavior of molecules can be averaged into a continuum theory. In the numerical portion of the proposed research, we inquire whether the contact line region is, in fact, large enough to possess a well-defined average. Additionally, we ask what types of behavior distinguish discrete systems from continuum systems. Might the smallness of the contact line region, in itself, lead to behavior different from that in the bulk? Taken together, our proposed research seeks to identify and accurately account for some of the molecular dynamics of the moving contact line, and attempts to formulate a description from which one can compute the forces at the moving contact line.

Han, Minsub↗

Structure and Soot Formation Properties of Laminar Flames

Soot formation within hydrocarbon-fueled flames is an important unresolved problem of combustion science for several reasons: soot emissions are responsible for more deaths than any other combustion-generated pollutant, thermal loads due to continuum radiation from soot limit the durability of combustors, thermal radiation from soot is mainly responsible for the growth and spread of unwanted fires, carbon monoxide emissions associated with soot emissions are responsible for most fire deaths, and limited understanding of soot processes in flames is a major impediment to the development of computational combustion. Motivated by these observations, soot processes within laminar premixed and nonpremixed (diffusion) flames are being studied during this investigation. The study is limited to laminar flames due to their experimental and computational tractability, noting the relevance of these results to practical flames through laminar flamelet concepts. Nonbuoyant flames are emphasized because buoyancy affects soot processes in laminar diffusion flames whereas effects of buoyancy are small for most practical flames. This study involves both ground- and space-based experiments, however, the following discussion will be limited to ground-based experiments because no space-based experiments were carried out during the report period. The objective of this work was to complete measurements in both premixed and nonpremixed flames in order to gain a better understanding of the structure of the soot-containing region and processes of soot nucleation and surface growth in these environments, with the latter information to be used to develop reliable ways of predicting soot properties in practical flames. The present discussion is brief, more details about the portions of the investigation considered here can be found in refs. 8-13.

El-Leathy, A. M.↗

Mechanics of composite materials: Recent advances; Proceedings of the Symposium, Virginia Polytechnic Institute and State University, Blacksburg, VA, August 16-19, 1982

The present conference on the mechanics of composites discusses microstructure's influence on particulate and short fiber composites' thermoelastic and transport properties, the elastoplastic deformation of composites, constitutive equations for viscoplastic composites, the plasticity and fatigue of metal matrix composites, laminate damping mechanisms, the micromechanical modeling of Kevlar/epoxy composites' time-dependent failure, the variational characterization of waves in composites, and computational methods for eigenvalue problems in composite design. Also discussed are the elastic response of laminates, elastic coupling nonlinear effects in unsymmetrical laminates, elasticity solutions for laminate problems having stress singularities, the mechanics of bimodular composite structures, the optimization of laminated plates and shells, NDE for laminates, the role of matrix cracking in the continuum constitutive behavior of a damaged composite ply, and the energy release rates of various microcracks in short fiber composites.

Hashin, Z.↗

Large Eddy Simulation (LES) of Particle-Laden Temporal Mixing Layers

High-fidelity models of plume-regolith interaction are difficult to develop because of the widely disparate flow conditions that exist in this process. The gas in the core of a rocket plume can often be modeled as a time-dependent, high-temperature, turbulent, reacting continuum flow. However, due to the vacuum conditions on the lunar surface, the mean molecular path in the outer parts of the plume is too long for the continuum assumption to remain valid. Molecular methods are better suited to model this region of the flow. Finally, granular and multiphase flow models must be employed to describe the dust and debris that are displaced from the surface, as well as how a crater is formed in the regolith. At present, standard commercial CFD (computational fluid dynamics) software is not capable of coupling each of these flow regimes to provide an accurate representation of this flow process, necessitating the development of custom software. This software solves the fluid-flow-governing equations in an Eulerian framework, coupled with the particle transport equations that are solved in a Lagrangian framework. It uses a fourth-order explicit Runge-Kutta scheme for temporal integration, an eighth-order central finite differencing scheme for spatial discretization. The non-linear terms in the governing equations are recast in cubic skew symmetric form to reduce aliasing error. The second derivative viscous terms are computed using eighth-order narrow stencils that provide better diffusion for the highest resolved wave numbers. A fourth-order Lagrange interpolation procedure is used to obtain gas-phase variable values at the particle locations.

Bellan, Josette↗

Effects of plasma mass flow on Alfven wave phase mixing in coronal loops

Ground-based and space observations suggest that flow phenomena are intrinsic to the solar corona loops. The present analysis addresses flow effects on Alfven wave behavior in coronal loops. The effect of flow on the shear Alfven wave continuum and the associated damping process resulting from phase-mixing is investigated. Flow Doppler-shifts the Alfven wave continuum, thus displacing the range of frequencies for optimum heating. The present study explores Alfven resonance effects on surface waves propagating along coronal loops that are under the influence of flow. Coupling to the Alfven resonance leads to damping of the surface waves. The damping rate in the presence of sub-Alfvenic flows is computed. Flow increases the damping rate of certain modes, thus enhancing the heating effectiveness of those modes. Therefore, these calculations suggest that flow is potentially an important ingredient for an accurate evaluation of heating the solar corona by dissipation of Alfven waves via phase-mixing.

Peredo, M.↗

Molecular Hydrogen in the Quiescent Disk of SW UMa

The FUSE observation has been reduced and a paper has been submitted to ApJ. The analysis has been slow because of the very noisy quality of the data, but we have derived line profile information for O VI and limits to the continuum brightness which place an interesting limit on the white dwarf temperature. The primary results are that a narrow O VI emission component seems to arise from the accretion flow onto the white dwarf itself, in agreement with cooling flow models for the X-ray spectra of low accretion rate dwarf novae. The broad component of the O VI lines is weaker than the observed C IV emission, suggesting that the UV line emission from the disk comes from photoionized plasma. A secondary result is that there is no H-2 fluorescent emission. The upper limits indicate that if molecular gas is present in the disk, it is shielded from Ly alpha photons by a layer of atomic hydrogen on the disk surface. We also derive an upper limit to the continuum level is below that observed by IUE. The limits are compatible with the lower end of the WD temperature range derived from IUE measurements, and they appear to agree with unpublished analysis of HST spectra. The grant has provided partial support for a data aide (Matt Povich) and a postdoc (Alex Lobel). It purchased a computer for M. Menou.

Raymond, John C.↗

A Light-Weight Inflatable Hypersonic Drag Device for Planetary Entry

The author has analyzed the use of a light-weight inflatable hypersonic drag device, called a ballute, (balloon + parachute) for flight in planetary atmospheres, for entry, aerocapture, and aerobraking. Studies to date include missions to Mars, Venus, Earth, Saturn, Titan, Neptune and Pluto. Data on a Pluto lander and a Mars orbiter will be presented to illustrate the concept. The main advantage of using a ballute is that aero deceleration and heating in atmospheric entry occurs at much smaller atmospheric density with a ballute than without it. For example, if a ballute has a diameter 10 times as large as the spacecraft, for unchanged total mass, entry speed and entry angle,the atmospheric density at peak convective heating is reduced by a factor of 100, reducing the peak heating by a factor of 10 for the spacecraft, and a factor of about 30 for the ballute. Consequently the entry payload (lander, orbiter, etc) is subject to much less heating, requires a much reduced thermal protection system (possibly only an MLI blanket), and the spacecraft design is therefore relatively unchanged from its vacuum counterpart. The heat flux on the ballute is small enough to be radiated at temperatures below 800 K or so. Also, the heating may be reduced further because the ballute enters at a more shallow angle, even allowing for the increased delivery angle error. Added advantages are a smaller mass ratio of entry system to total entry mass, and freedom from the low-density and transonic instability problems that conventional rigid entry bodies suffer, since the vehicle attitude is determined by the ballute, usually released at continuum conditions (hypersonic for an orbiter, and subsonic for a lander). Also, for a lander the range from entry to touchdown is less, offering a smaller footprint. The ballute derives an entry corridor for aerocapture by entering on a path that would lead to landing, and releasing the ballute adaptively, responding to measured deceleration, at a speed computed to achieve the desired orbiter exit conditions. For a lander an accurate landing point could be achieved by providing the lander with a small gliding capacity, using the large potential energy available from being subsonic at high altitude. Alternatively the ballute can be retained to act as a parachute or soft-landing device, or to float the payload as a buoyant aerobot. As expected, the ballute has smaller size for relatively small entry speeds, such as for Mars, or for the extensive atmosphere of a low-gravity planet such as Pluto. The author will discuss presently available ballute materials and a development program of aerodynamic tests and materials that would be required for ballutes to achieve their full potential.

McRonald, Angus D.↗

A Light-Weight Inflatable Hypersonic Drag Device for Planetary Entry

The author has analyzed the use of a light-weight inflatable hypersonic drag device, called a ballute, for flight in planetary atmospheres, for entry, aerocapture, and aerobraking. Studies to date include Mars, Venus, Earth, Saturn, Titan, Neptune and Pluto, and data on a Pluto lander and a Mars orbiter will be presented to illustrate the concept. The main advantage of using a ballute is that aero, deceleration and heating in atmospheric entry occurs at much smaller atmospheric density with a ballute than without it. For example, if a ballute has a diameter 10 times as large as the spacecraft, for unchanged total mass, entry speed and entry angle,the atmospheric density at peak convective heating is reduced by a factor of 100, reducing the heating by a factor of 10 for the spacecraft and a factor of 30 for the ballute. Consequently the entry payload (lander, orbiter, etc) is subject to much less heating, requires a much reduced thermal. protection system (possibly only an MLI blanket), and the spacecraft design is therefore relatively unchanged from its vacuum counterpart. The heat flux on the ballute is small enough to be radiated at temperatures below 800 K or so. Also, the heating may be reduced further because the ballute enters at a more shallow angle, even allowing for the increased delivery angle error. Added advantages are less mass ratio of entry system to total entry mass, and freedom from the low-density and transonic instability problems that conventional rigid entry bodies suffer, since the vehicle attitude is determined by the ballute, usually released at continuum conditions (hypersonic for an orbiter, and subsonic for a lander). Also, for a lander the range from entry to touchdown is less, offering a smaller footprint. The ballute derives an entry corridor for aerocapture by entering on a path that would lead to landing, and releasing the ballute adaptively, responding to measured deceleration, at a speed computed to achieve the desired orbiter exit conditions. For a lander an accurate landing point could be achieved by providing the lander with a small gliding capacity, using the large potential energy available from being subsonic at high altitude. Alternatively the ballute can be retained to act as a parachute or soft-landing device, or to float the payload as a buoyant aerobot. As expected, the ballute has smaller size for relatively small entry speeds, such as for Mars and Titan, or for the extensive atmosphere of a low-gravity planet such as Pluto. Details of a ballute to place a small Mars orbiter and a small Pluto lander will be given to illustrate the concept. The author will discuss presently available ballute materials and a development program of aerodynamic tests and materials that would be required for ballutes to achieve their full potential.

McRonald, Angus D.↗

Optical Measurement System for Strain Field Ahead of a Crack Tip for Lattice Structures

The aim of the ARAMADAS project is to automate the construction of cuboctahedral lattice structures. Lattice materials are appealing for aerospace applications due to their strength and stiffness at ultra-light densities. However, in order for any material to be realistically considered for such environments, it must also be damage tolerant. The ability of a material to absorb damage is characterized by its fracture toughness, which remains poorly characterized for lattice materials. Consequently, the objective of this research is to develop an optical measurement system to experimentally validate the strain field ahead of a crack tip in architecture lattice materials. Although the ability to predict the strain field ahead of a crack tip has been investigated for continuum materials, such behaviour of three-dimensional architectures is under-investigated. As such, we will use a custom optical measurement system to track deformation of the voxels in a side-cracked plate fracture specimen. The system shall use 3D pose estimation, stereo imaging, and possibly color tracking, in combination with optical flow algorithms, to compute information regarding the three-dimensional movement movement of the lattice nodes during mechanical testing. Bench top experiments will validate the optical measurement system and characterize precision. Additionally, the effect of multiple cameras on precision, as well as system scalability will be investigated. Final results will compare measured lattice deformation with finite element predictions.

Lennig, Miriam↗

The optical properties of platinum and gold in the vacuum ultraviolet

The optical constants of platinum and gold thin films have been determined in the spectral region of 40 to 200 nm by reflection measurements. The highly polarized continuum of synchrotron radiation emitted by the 240-MeV electron storage ring at the Physical Sciences Laboratory of the University of Wisconsin was used as a light source for the spectrum below 120 nm, while a windowless discharge lamp coupled to a normal incidence monochromator provided a source for the longer wavelengths. Optical constants were determined by a computer program based on iterative solutions to the Fresnel equations for reflection as a function of the angle of incidence.

Linton, R. C.↗

High spectral resolution ground-based observations of Venus in the 450- to 1250 kayser region

Ground-based observations of Venus were made with a 5-cm drive Michelson interferometer during December 1970 and December 1973. The thermal emission spectrum of the central portion of the apparent disk was recorded from 450-1250 kaysers with an apodized spectral resolution of 0.25 kaysers. All statistically significant sharp line-absorption features in the spectrum have been identified with gaseous CO2. Comparison between the observed spectrum and a synthetic spectrum computed from a model atmosphere, assuming gaseous CO2 and a sulfuric acid haze as opacity sources, indicates good agreement. A broad diffuse absorption feature associated with the sulfuric acid haze is evident in the 870- to 930 kayser region. With the exception of the rotational lines of the 927 kayser CO2 band, the above feature appears as a continuum down to 0.25 kayser resolution. In the 750- to 1250 kayser range, the spectrum exhibits moderate thermal contrast with maximum brightness temperatures of 234-238 K occurring near 825 kaysers. These temperatures are in general agreement with previous measurements.

Kunde, V. G.↗

Assessment of three numerical methods for the computation of a low-density plume flow

Results from three numerical methods including one based on the Navier-Stokes equations, one based on kinetic theory using the DSMC method, and one based on the Boltzmann equation with a Krook-type collision term are compared to each other and to experimental data for a model problem of heated nitrogen flow in a conical nozzle expanding into a vacuum. The problem simulates flow in a resistojet, a low-thrust, electrothermal rocket. The continuum method is applied to both the internal flow and near-field plume. The DSMC and Boltzmann methods are applied primarily to the plume. Experimental measurements of Pitot pressure and flow angle, taken with an apparatus that duplicates the model nozzle flow, are used in the comparisons.

Penko, Paul F.↗

Alaska Crustal Deformation: Finite Element Modeling Constrained by Geologic and Very Long Baseline Interferometry Data

We compute crustal motions in Alaska by calculating the finite element solution for an elastic spherical shell problem. The method we use allows the finite element mesh to include faults and very long baseline interferometry (VLBI) baseline rates of change. Boundary conditions include Pacific-North American (PA-NA) plate motions. The solution is constrained by the oblique orientation of the Fairweather-Queen Charlotte strike-slip faults relative to the PA-NA relative motion direction and the oblique orientation from normal convergence of the eastern Aleutian trench fault systems, as well as strike-shp motion along the Denali and Totschunda fault systems. We explore the effects that a range of fault slip constraints and weighting of VLBI rates of change has on the solution. This allows us to test the motion on faults, such as the Denali fault, where there are conflicting reports on its present-day slip rate. We find a pattern of displacements which produce fault motions generally consistent with geologic observations. The motion of the continuum has the general pattern of radial movement of crust to the NE away from the Fairweather-Queen Charlotte fault systems in SE Alaska and Canada. This pattern of crustal motion is absorbed across the Mackenzie Mountains in NW Canada, with strike-slip motion constrained along the Denali and Tintina fault systems. In south central Alaska and the Alaska forearc oblique convergence at the eastern Aleutian trench and the strike-shp motion of the Denali fault system produce a counterclockwise pattern of motion which is partially absorbed along the Contact and related fault systems in southern Alaska and is partially extruded into the Bering Sea and into the forearc parallel the Aleutian trench from the Alaska Peninsula westward. Rates of motion and fault slip are small in western and northern Alaska, but the motions we compute are consistent with the senses of strike-slip motion inferred geologically along the Kaltag, Kobuk Trench, and Thompson Creek faults and with the normal faulting observed in NW Alaska near Nome. The nonrigid behavior of our finite element solution produces patterns of motion that would not have been expected from rigid block models: strike-slip faults can exist in a continuum that has motion mostly perpendicular to their strikes, and faults can exhibit along-strike differences in magnitudes and directions.

Lundgren, Paul↗

The Nimbus 4 infrared spectroscopy experiment. II - Comparison of observed and theoretical radiances from 425-1450 wavelengths/cm

The Nimbus 4 infrared interferometer spectrometer (IRIS) measured the thermal emission of the earth's atmosphere and surface from 400-1600 wavelengths/cm with an apodized spectral resolution of 2.8 wavelengths/cm. A comparison of theoretical radiances, computed from in situ measurements and using a direct integration slant path atmospheric transmittance model, with the observed IRIS radiances has been made to verify the radiometric and spectral performance of the instrument and to assess the accuracy of the atmospheric transmittances. The radiance comparison has indicated a relatively constant difference of less than 5% in the water vapor continuum in the 425 to 550 wavelengths/cm and 750 to 1200 wavelengths/cm atmospheric 'window' regions, whereas in the 667 wavelengths/cm CO2 band the difference was 5-10%.

Kunde, V. G.↗

Acceleration efficiency in line-driven flows

We reexamine the physics of flow driving by line scattering of a continuum radiation source to determine the degree to which such line scattering can heat as well as accelerate the flow. Within the framework of the Sobolev theory for line transfer, we argue that the finite thermal width of the line scattering profile can lead to a significant 'Doppler heating' via photon frequency redistribution within a Sobolev resonance layer. Quantitative computation of this heating shows, however, that it is largely canceled by a corresponding cooling by the diffuse radiation. The resulting reduction in net Doppler heating or cooling means that the overall effect is only of limited importance in the energy balance of line-driven stellar winds. Through simple scaling relations, we compare the effect to other competing heating or cooling terms, including the ion-drag frictional heating recently discussed by Springmann and Pauldrach. We also provide a physical explanation of the unexpected cooling effect, and comment that its near cancellation of the anticipated heating provides another example of the tendency for ideal Sobolev theory to apply to a higher order than expected.

Gayley, Kenneth G.↗

Catastrophic Disruption of Comet ISON

We report submillimeter 450 and 850 microns dust continuum observations for comet C/2012 S1 (ISON) obtained at heliocentric distances 0.31-0.08 au prior to perihelion on 2013 November 28 (rh=0.0125 au). These observations reveal a rapidly varying dust environment in which the dust emission was initially point-like. As ISON approached perihelion, the continuum emission became an elongated dust column spread out over as much as 60″ (greater than 10(exp 5) km in the anti-solar direction. Deconvolution of the November 28.04 850 microns image reveals numerous distinct clumps consistent with the catastrophic disruption of comet ISON, producing approximately 5.2×10(exp 10) kg of submillimeter-sized dust. Orbital computations suggest that the SCUBA-2 emission peak coincides with the cometʼs residual nucleus.

comets: individual↗