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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

A fast operator perturbation method for the solution of the special relativistic equation of radiative transfer in spherical symmetry

A fast method for the solution of the radiative transfer equation in rapidly moving spherical media, based on an approximate Lambda-operator iteration, is described. The method uses the short characteristic method and a tridiagonal approximate Lambda-operator to achieve fast convergence. The convergence properties and the CPU time requirements of the method are discussed for the test problem of a two-level atom with background continuum absorption and Thomson scattering. Details of the actual implementation for fast vector and parallel computers are given. The method is accurate and fast enough to be incorporated in radiation-hydrodynamic calculations.

Hauschildt, P. H.↗
NASA NTRS2018

Modeling Fiber Kinking at the Microscale and Mesoscale

A computational micromechanics (CMM) model is employed to interrogate the assumptions of a recently developed mesoscale continuum damage mechanics (CDM) model for fiber kinking. The CMM model considers an individually discretized three dimensional fiber and surrounding matrix accounting for nonlinearity in the fiber, matrix plasticity, fiber/matrix interface debonding, and geometric nonlinearity. Key parameters of the CMM model were measured through experiments. In particular, a novel experimental technique to characterize the in situ longitudinal compressive strength of carbon fibers through indentation of micropillars is presented. The CDM model is formulated on the basis of Budiansky's fiber kinking theory (FKT) with a constitutive deformation-decomposition approach to alleviate mesh size sensitivity. In contrast to conventional mesoscale CDM models that prescribe a constitutive response directly, the response of the proposed model is an outcome of material nonlinearity and large rotations of the fiber direction following FKT. Comparison of the predictions from the CMM and CDM models shows remarkable correlation in strength, post-peak residual stress, and fiber rotation, with less than 10% difference between the two models in most cases. Additional comparisons are made with several fiber kinking models proposed in the literature to highlight the efficacy of the two models. Finally, the CMM model is exercised in parametric studies to explore opportunities to improve the longitudinal compression strength of a ply through the use of nonconventional microstructures.

Herraez, Miguel↗
NASA NTRS2007

Convergent Close-Coupling Approach to Electron-Atom Collisions

It was with great pleasure and honour to accept the invitation to make a presentation at the symposium celebrating the life-long work of Aaron Temkin and Richard Drachman. The work of Aaron Temkin was particularly influential on our own during the development of the CCC method for electron-atom collisions. There are a number of key problems that need to be dealt with when developing a general computational approach to such collisions. Traditionally, the electron energy range was subdivided into the low, intermediate, and high energies. At the low energies only a finite number of channels are open and variational or close-coupling techniques could be used to obtain accurate results. At high energies an infinite number of discrete channels and the target continuum are open, but perturbative techniques are able to yield accurate results. However, at the intermediate energies perturbative techniques fail and computational approaches need to be found for treating the infinite number of open channels. In addition, there are also problems associated with the identical nature of electrons and the difficulty of implementing the boundary conditions for ionization processes. The beauty of the Temkin-Poet model of electron-hydrogen scattering is that it simplifies the full computational problem by neglecting any non-zero orbital angular momenta in the partial-wave expansion, without loosing the complexity associated with the above-mentioned problems. The unique nature of the problem allowed for accurate solution leading to benchmark results which could then be used to test the much more general approaches to electron-atom collision problems. The immense value of the Temkin-Poet model is readily summarised by the fact that the initial papers of Temkin and Poet have been collectively cited around 250 times to date and are still being cited in present times. Many of the citations came from our own work during the course of the development of the CCC method, which we now describe.

Bray, Igor↗
NASA NTRS2012

Computational Implementation of a Thermodynamically Based Work Potential Model For Progressive Microdamage and Transverse Cracking in Fiber-Reinforced Laminates

A continuum-level, dual internal state variable, thermodynamically based, work potential model, Schapery Theory, is used capture the effects of two matrix damage mechanisms in a fiber-reinforced laminated composite: microdamage and transverse cracking. Matrix microdamage accrues primarily in the form of shear microcracks between the fibers of the composite. Whereas, larger transverse matrix cracks typically span the thickness of a lamina and run parallel to the fibers. Schapery Theory uses the energy potential required to advance structural changes, associated with the damage mechanisms, to govern damage growth through a set of internal state variables. These state variables are used to quantify the stiffness degradation resulting from damage growth. The transverse and shear stiffness of the lamina are related to the internal state variables through a set of measurable damage functions. Additionally, the damage variables for a given strain state can be calculated from a set of evolution equations. These evolution equations and damage functions are implemented into the finite element method and used to govern the constitutive response of the material points in the model. Additionally, an axial failure criterion is included in the model. The response of a center-notched, buffer strip-stiffened panel subjected to uniaxial tension is investigated and results are compared to experiment.

Pineda, Evan J.↗
NASA NTRS1978

Computation of three-dimensional nozzle-exhaust flow fields with the GIM code

A methodology is introduced for constructing numerical analogs of the partial differential equations of continuum mechanics. A general formulation is provided which permits classical finite element and many of the finite difference methods to be derived directly. The approach, termed the General Interpolants Method (GIM), can combined the best features of finite element and finite difference methods. A quasi-variational procedure is used to formulate the element equations, to introduce boundary conditions into the method and to provide a natural assembly sequence. A derivation is given in terms of general interpolation functions from this procedure. Example computations for transonic and supersonic flows in two and three dimensions are given to illustrate the utility of GIM. A three-dimensional nozzle-exhaust flow field is solved including interaction with the freestream and a coupled treatment of the shear layer. Potential applications of the GIM code to a variety of computational fluid dynamics problems is then discussed in terms of existing capability or by extension of the methodology.

Spradley, L. W.↗
NASA NTRS1982

Quiet sun observations of the Al I autoionization lines lambda 1932 and lambda 1936

Quet sun observations obtained during a rocket flight of the Al I autoionization lines 1932 A and 1936 A at solar pointings ranging from mu = 0.73 out to the visible limb. Absolute intensities are estimated to be accurate to about + or - 20%. These lines weaken progressively with decreasing mu but never go into emission before finally disappearing with the continuum just beyond the visible solar limb. The observations are then compared with LTE line profiles computed through the quiet sun atmosphere of Vernazza et al. (1976). Several areas of disagreement between the synthetic and observed profiles are discussed.

Cook, J. W.↗
NASA NTRS1993

Effects of inclination angle on the spectra of X-ray binaries

Iron emission and absorption features from nine X-ray binaries are compared with a model that predicts characteristics of iron K-shell features corresponding to different viewing angles for the accretion disk. In this model, emission and absorption features arise from an accretion disk corona whose structure is computed under the assumptions of ionization, thermal, and hydrostatic balance of gas illuminated by the central continuum source. Hence the equivalent widths of the emission and absorption features are predicted without the introduction of any free parameters to control the physical conditions in the disk corona. Our data/model comparison suggests an explanation for the equivalent widths of the iron K-alpha emission and the depth of the absorption edge in terms of (1) the angle of inclination between the observer's line of sight and the plane of the accretion disk and (2) the shape of the incident source spectrum.

Vrtilek, S. D.↗
NASA NTRS2003

Observational Investigation of Solar Interior and Atmosphere

The Imaging Vector Magnetograph (IVM) has been modified to make it easier to observe at more than one spectral line. The cell holding the blocking filter has been replaced by a four-position filter wheel, so that changing to a different line is a matter of a few minutes rather than the several hours it used to take to disassemble the cell and install a new filter. Three new filters have been obtained, for Na 1589.6 nm, Fe 1630.25 nm, and H 1656.3 nm. The new filters have better bandpass profiles than the ones they replaced: somewhat wider, with flatter tops and steeper wings. This results in a reduction of parasitic light coming from adjacent Fabry-Perot orders, from seven percent to about two percent, and flattens the apparent continuum. The Mees CCD Imaging Spectrograph (MCCD) was upgraded under this grant, with a new control computer and data system. The camera was replaced with a faster, larger-format frame-transfer camera. Final integration of the upgrades is not yet complete, but tests indicate that the system cadence will be improved by a factor of five to ten, while increasing the spatial coverage by a factor of two (depending on observation options). Synoptic observations with the IVM and MCCD continue to be conducted daily, to the extent permitted by the fact that we have a single observer responsible for the observations. The older Haleakala Stokes Polarimeter is also used to make a daily vector magnetogram, normally of the region selected by the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) duty scientists. This instrument, however, is showing its age to the extent that its maintenance is becoming something of a challenge. We also run a white light full-disk imager and a video H alpha prominence camera, continuously during times of observations. Of particular interest, we obtained rapid-cadence observations of the 2003 July 15 white light flare with both the IVM and MCCD. The vector magnetograms show no obvious difference between the preflare and postflare configurations. The photospheric magnetogram sequence and H alpha spectra we obtained can be combined with spacecraft observations to construct a rather detailed picture of the events of this flare. Several papers are in preparation.

Kuhn, Jeffrey R.↗
NASA NTRS2018

Numerically Simulating an Expanding Continuum Jet into a Surrounding Non-Continuum Region

A numerical technique that simulates the interaction between a continuum jet and a free-molecular, rarefied or transitional flow field is presented and applied. The technique is based on selecting a boundary between the expanding jet plume and interacting flow field defined by the Bird breakdown parameter, P. By properly choosing the boundary, the jet plume computed by computational fluid dynamics (CFD) is uncoupled from the interaction region that is simulated by direct simulation Monte Carlo (DSMC). The mechanics of uncoupling the continuum and non-continuum regions is discussed first. Also, a discussion of expanding continuum plume breakdown and the formulation of the Bird breakdown parameter is presented. To show the usefulness, examples of applying the CFD-DSMC technique are given. Many of the examples were taken from projects worked by the staff of the Aerothermodynamics Branch at the NASA Langley Research Center. These examples show a variety of applications of the technique and allow those examples not previously documented to be formally presented.

Jet interaction↗
NASA NTRS1981

Theory, computation, and application of exponential splines

A generalization of the semiclassical cubic spline known in the literature as the exponential spline is discussed. In actuality, the exponential spline represents a continuum of interpolants ranging from the cubic spline to the linear spline. A particular member of this family is uniquely specified by the choice of certain tension parameters. The theoretical underpinnings of the exponential spline are outlined. This development roughly parallels the existing theory for cubic splines. The primary extension lies in the ability of the exponential spline to preserve convexity and monotonicity present in the data. Next, the numerical computation of the exponential spline is discussed. A variety of numerical devices are employed to produce a stable and robust algorithm. An algorithm for the selection of tension parameters that will produce a shape preserving approximant is developed. A sequence of selected curve-fitting examples are presented which clearly demonstrate the advantages of exponential splines over cubic splines.

Mccartin, B. J.↗
NASA NTRS1990

A computational procedure for multibody systems including flexible beam dynamics

A computational procedure suitable for the solution of equations of motions for flexible multibody systems has been developed. A fully nonlinear continuum approach capable of accounting for both finite rotations and large deformations has been used to model a flexible beam component. The beam kinematics are referred directly to an inertial reference frame such that the degrees of freedom embody both the rigid and flexible deformation motions. As such, the beam inertia expression is identical to that of rigid body dynamics. The nonlinear coupling between gross body motion and elastic deformation is contained in the internal force expression. Numerical solution procedures for the integration of spatial kinematic systems can be directily applied to the generalized coordinates of both the rigid and flexible components. An accurate computation of the internal force term which is invariant to rigid motions is incorporated into the general solution procedure.

Downer, J. D.↗
NASA NTRS1978

Two-dimensional computer code for design assessment of SRM aft closure environment

A general interpolants method for constructing numerical analogs of the partial differential equations of continuum mechanics and which combines the best features of both finite element and finite difference methods was used in a two dimensional axisymmetric analysis of the space shuttle solid motor aft closure region flowfield. The numerical technique used is discussed as well as the results of the analysis. A steady state solution for the submerged nozzle nose region indicates the development of an area gas flow recirculation near the propellant burning surface boundary and in the region of the underside of the nozzle nose. This recirculation region remains spatially fixed for successive iterations (greater than 800) of the solution.

Freeman, J. A.↗
NASA NTRS1972

Total radiative intensity calculations for 100% H2 and 87% H2-13% He.

Isothermal radiative intensity calculations for 100% H2 and 87% H2-13% He are presented for temperatures of 10,000-25,000 K, density ratios of .0001 to .1, and path lengths of 1.0-30.0 cm. The actual spectral details of the absorption coefficient were computed for 16,000 points from 240 to 30,000 A by summing the various line and continuum radiative processes at each point. This method should result in a very accurate calculation of radiative emission, including an accurate accounting for reabsorption due to overlapping lines.

Stickford, G. H., Jr.↗
NASA NTRS1992

Solution of the Burnett equations for hypersonic flows near the continuum limit

The INCA code, a three-dimensional Navier-Stokes code for analysis of hypersonic flowfields, was modified to analyze the lower reaches of the continuum transition regime, where the Navier-Stokes equations become inaccurate and Monte Carlo methods become too computationally expensive. The two-dimensional Burnett equations and the three-dimensional rotational energy transport equation were added to the code and one- and two-dimensional calculations were performed. For the structure of normal shock waves, the Burnett equations give consistently better results than Navier-Stokes equations and compare reasonably well with Monte Carlo methods. For two-dimensional flow of Nitrogen past a circular cylinder the Burnett equations predict the total drag reasonably well. Care must be taken, however, not to exceed the range of validity of the Burnett equations.

Imlay, Scott T.↗
NASA NTRS1997

Soot Formation in Freely-Propagating Laminar Premixed Flames

Soot formation within hydrocarbon-fueled flames is an important unresolved problem of combustion science. Thus, the present study is considering soot formation in freely-propagating laminar premixed flames, exploiting the microgravity environment to simplify measurements at the high-pressure conditions of interest for many practical applications. The findings of the investigation are relevant to reducing emissions of soot and continuum radiation from combustion processes, to improving terrestrial and spacecraft fire safety, and to developing methods of computational combustion, among others. Laminar premixed flames are attractive for studying soot formation because they are simple one-dimensional flows that are computationally tractable for detailed numerical simulations. Nevertheless, studying soot-containing burner-stabilized laminar premixed flames is problematical: spatial resolution and residence times are limited at the pressures of interest for practical applications, flame structure is sensitive to minor burner construction details so that experimental reproducibility is not very good, consistent burner behavior over the lengthy test programs needed to measure soot formation properties is hard to achieve, and burners have poor durability. Fortunately, many of these problems are mitigated for soot-containing, freely-propagating laminar premixed flames. The present investigation seeks to extend work in this laboratory for various soot processes in flames by observing soot formation in freely-propagating laminar premixed flames. Measurements are being made at both Normal Gravity (NG) and MicroGravity (MG), using a short-drop free-fall facility to provide MG conditions.

Lin, K.-C.↗
NASA NTRS1998

Computational Study of Surface Tension and Wall Adhesion Effects on an Oil Film Flow Underneath an Air Boundary Layer

The fringe-imaging skin friction (FISF) technique, which was originally developed by D. J. Monson and G. G. Mateer at Ames Research Center and recently extended to 3-D flows, is the most accurate skin friction measurement technique currently available. The principle of this technique is that the skin friction at a point on an aerodynamic surface can be determined by measuring the time-rate-of-change of the thickness of an oil drop placed on the surface under the influence of the external air boundary layer. Lubrication theory is used to relate the oil-patch thickness variation to shear stress. The uncertainty of FISF measurements is estimated to be as low as 4 percent, yet little is known about the effects of surface tension and wall adhesion forces on the measured results. A modified version of the free-surface Navier-Stokes solver RIPPLE, developed at Los Alamos National Laboratories, was used to compute the time development of an oil drop on a surface under a simulated air boundary layer. RIPPLE uses the volume of fluid method to track the surface and the continuum surface force approach to model surface tension and wall adhesion effects. The development of an oil drop, over a time period of approximately 4 seconds, was studied. Under the influence of shear imposed by an air boundary layer, the computed profile of the drop rapidly changes from its initial circular-arc shape to a wedge-like shape. Comparison of the time-varying oil-thickness distributions computed using RIPPLE and also computed using a greatly simplified numerical model of an oil drop equation which does not include surface tension and wall adhesion effects) was used to evaluate the effects of surface tension on FISF measurement results. The effects of surface tension were found to be small but not necessarily negligible in some cases.

Celic, Alan↗