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

Statistical discrete particle simulation studies

A particle simulation code was developed that is suited to Multiple Instruction Multiple Data multiprocessor computers. The resulting code is designed to be portable between a large class of multiprocessor computer architectures, including hypercubes (e.g., Intel iPSC/860), shared memory machines (e.g., Cray, SGI, DASH), and simple uniprocessors (e.g., SUN, VAX, IBM PC). This code provides a powerful general framework that can be adapted to specific applications, such as thermomechanical modeling, general 3-D geometry support, or visualization support via a codeveloped CPlot data reduction and visualization system. Chemistry models were implemented and validated with comparison to continuum solutions of thermochemically relaxing gas mixtures and with experimental results from the high speed flow about a circular cylinder. Furthermore, the exchange of energy between translational and internal modes was modeled. New models were developed, combining greater efficiency than earlier phenomenological models. They offer greater physical detail by addressing vibrational energy as a discretely distributed quantity.

Mcdonald, Jeffrey D.↗
NASA NTRS2012

Plume-Free Stream Interaction Heating Effects During Orion Crew Module Reentry

During reentry of the Orion Crew Module (CM), vehicle attitude control will be performed by firing reaction control system (RCS) thrusters. Simulation of RCS plumes and their interaction with the oncoming flow has been difficult for the analysis community due to the large scarf angles of the RCS thrusters and the unsteady nature of the Orion capsule backshell environments. The model for the aerothermal database has thus relied on wind tunnel test data to capture the heating effects of thruster plume interactions with the freestream. These data are only valid for the continuum flow regime of the reentry trajectory. A Direct Simulation Monte Carlo (DSMC) analysis was performed to study the vehicle heating effects that result from the RCS thruster plume interaction with the oncoming freestream flow at high altitudes during Orion CM reentry. The study was performed with the DSMC Analysis Code (DAC). The inflow boundary conditions for the jets were obtained from Data Parallel Line Relaxation (DPLR) computational fluid dynamics (CFD) solutions. Simulations were performed for the roll, yaw, pitch-up and pitch-down jets at altitudes of 105 km, 125 km and 160 km as well as vacuum conditions. For comparison purposes (see Figure 1), the freestream conditions were based on previous DAC simulations performed without active RCS to populate the aerodynamic database for the Orion CM. Other inputs to the analysis included a constant Orbital reentry velocity of 7.5 km/s and angle of attack of 160 degrees. The results of the study showed that the interaction effects decrease quickly with increasing altitude. Also, jets with highly scarfed nozzles cause more severe heating compared to the nozzles with lower scarf angles. The difficulty of performing these simulations was based on the maximum number density and the ratio of number densities between the freestream and the plume for each simulation. The lowest altitude solutions required a substantial amount of computational resources (up to 1800 processors) to simulate approximately 2 billion molecules for the refined (adapted) solutions.

Marichalar, J.↗
NASA NTRS2004

A Hierarchical Approach to Fracture Mechanics

Recent research conducted under NASA LaRC's Creativity and Innovation Program has led to the development of an initial approach for a hierarchical fracture mechanics. This methodology unites failure mechanisms occurring at different length scales and provides a framework for a physics-based theory of fracture. At the nanoscale, parametric molecular dynamic simulations are used to compute the energy associated with atomic level failure mechanisms. This information is used in a mesoscale percolation model of defect coalescence to obtain statistics of fracture paths and energies through Monte Carlo simulations. The mathematical structure of predicted crack paths is described using concepts of fractal geometry. The non-integer fractal dimension relates geometric and energy measures between meso- and macroscales. For illustration, a fractal-based continuum strain energy release rate is derived for inter- and transgranular fracture in polycrystalline metals.

Saether, Erik↗
NASA NTRS2022

Recommended Direct Simulation Monte Carlo Collision Model Parameters for Reacting Methane Flows

While hydrocarbon combustion and gas-surface interactions have been simulated for a variety of applications in continuum solvers, such processes have been rarely investigated in Direct Simulation Monte Carlo (DSMC). Thus Variable Hard/Soft Sphere (VHS/VSS)collision parameters are not available for most molecules of interest in the decomposition of methane, a significant pyrolysis gas product and potential atmospheric species for Titan re-entries. As these properties are essential for accurate modeling of gas transport and thermochemical effects, a method must be devised to compute them for use in DSMC. In this study, the methodology outlined in Stephani et al is utilized to compute the VHS/VSS parameters by fitting the Ω(1,1)and Ω(2,2) collision integrals [1]. These collision integrals are typically computed from intermolecular potentials and are used in CFD simulations to calculate transport properties such as viscosity, thermal conductivity, and diffusion. This approach ensures consistency between micro-scale collisions and macro-scale transport properties. To produce the necessary collision integrals needed for the fitting routine, we utilize the method described by Laricchiuta et al, where collision integrals are computed based on the constituent particles’ polarizability and effective electron number [2].Preliminary computations of pure-species Ω(1,1)and Ω(2,2)collision integrals produced from the coupled approach of Stephani et al and Laricchiuta et al for methane show good agreement with those computed via available ab-initio Potential Energy Surface Computations, as seen in Figure 1. Above 300K, no more than 12% deviation is observed between the two results over the temperature range of the fitted PES data [3].Current results of this study include full collision-specificVSS transport properties for a complex-species mixture involving major byproducts of pure methane decomposition. Targeted final results will include a comprehensive list of VSS parameters relevant toDSMC simulations of Titan atmospheric entry, based on the reduced kinetic model of Savajano et al [4].

M R Gosma↗
NASA NTRS1993

Hypersonic blunt body wake computations using DSMC and Navier-Stokes solvers

Numerical results obtained with direct simulation Monte Carlo (DSMC) and Navier-Stokes methods are presented for Mach 20 nitrogen flow about a 70-deg blunted cone. The flow conditions simulated are those that can be obtained in existing low-density hypersonic wind tunnels. Three sets of flow conditions are considered with freestream Knudsen numbers ranging from 0.03 to 0.001. The focus is on the wake structure: how does the wake structure change as a function of rarefaction, what are the afterbody levels of heating, and to what limits are continuum models realistic as rarefaction in the wake is progressively increased. Calculations are made with and without an afterbody sting. Results for the afterbody sting are emphasized in anticipation of an experimental study for the current flow conditions and model configuration. The Navier-Stokes calculations were made with and without slip boundary conditions. Comparisons of the results obtained with the two simulation methodologies are made for both flowfield structure and surface quantities.

Moss, James N.↗
NASA NTRS1973

Adhesive-bonded scarf and stepped-lap joints

Continuum mechanics solutions are derived for the static load-carrying capacity of scarf and stepped-lap adhesive-bonded joints. The analyses account for adhesive plasticity and adherend stiffness imbalance and thermal mismatch. The scarf joint solutions include a simple algebraic formula which serves as a close lower bound, within a small fraction of a per cent of the true answer for most practical geometries and materials. Digital computer programs were developed and, for the stepped-lap joints, the critical adherend and adhesive stresses are computed for each step. The scarf joint solutions exhibit grossly different behavior from that for double-lap joints for long overlaps inasmuch as that the potential bond shear strength continues to increase with indefinitely long overlaps on the scarf joints. The stepped-lap joint solutions exhibit some characteristics of both the scarf and double-lap joints. The stepped-lap computer program handles arbitrary (different) step lengths and thickness and the solutions obtained have clarified potentially weak design details and the remedies. The program has been used effectively to optimize the joint proportions.

Hart-Smith, L. J.↗
NASA NTRS1997

Predictive Modeling in Plasma Reactor and Process Design

Research continues toward the improvement and increased understanding of high-density plasma tools. Such reactor systems are lauded for their independent control of ion flux and energy enabling high etch rates with low ion damage and for their improved ion velocity anisotropy resulting from thin collisionless sheaths and low neutral pressures. Still, with the transition to 300 mm processing, achieving etch uniformity and high etch rates concurrently may be a formidable task for such large diameter wafers for which computational modeling can play an important role in successful reactor and process design. The inductively coupled plasma (ICP) reactor is the focus of the present investigation. The present work attempts to understand the fundamental physical phenomena of such systems through computational modeling. Simulations will be presented using both computational fluid dynamics (CFD) techniques and the direct simulation Monte Carlo (DSMC) method for argon and chlorine discharges. ICP reactors generally operate at pressures on the order of 1 to 10 mTorr. At such low pressures, rarefaction can be significant to the degree that the constitutive relations used in typical CFD techniques become invalid and a particle simulation must be employed. This work will assess the extent to which CFD can be applied and evaluate the degree to which accuracy is lost in prediction of the phenomenon of interest; i.e., etch rate. If the CFD approach is found reasonably accurate and bench-marked with DSMC and experimental results, it has the potential to serve as a design tool due to the rapid time relative to DSMC. The continuum CFD simulation solves the governing equations for plasma flow using a finite difference technique with an implicit Gauss-Seidel Line Relaxation method for time marching toward a converged solution. The equation set consists of mass conservation for each species, separate energy equations for the electrons and heavy species, and momentum equations for the gas. The sheath is modeled by imposing the Bohm velocity to the ions near the walls. The DSMC method simulates each constituent of the gas as a separate species which would be analogous in CFD to employing separate species mass, momentum, and energy equations. All particles including electrons are moved and allowed to collide with one another with the stipulation that the electrons remain tied to the ions consistent with the concept of ambipolar diffusion. The velocities of the electrons are allowed to be modified during collisions and are not confined to a Maxwellian distribution. These benefits come at a price in terms of computational time and memory. The DSMC and CFD are made as consistent as possible by using similar chemistry and power deposition models. Although the comparison of CFD and DSMC is interesting, the main goal of this work is the increased understanding of high-density plasma flowfields that can then direct improvements in both techniques. This work is unique in the level of the physical models employed in both the DSMC and CFD for high-density plasma reactor applications. For example, the electrons are simulated in the present DSMC work which has not been done before for low temperature plasma processing problems. In the CFD approach, for the first time, the charged particle transport (discharge physics) has been self-consistently coupled to the gas flow and heat transfer.

Hash, D. B.↗
NASA NTRS1973

Mathematical modeling of spinning elastic bodies for modal analysis.

The problem of modal analysis of an elastic appendage on a rotating base is examined to establish the relative advantages of various mathematical models of elastic structures and to extract general inferences concerning the magnitude and character of the influence of spin on the natural frequencies and mode shapes of rotating structures. In realization of the first objective, it is concluded that except for a small class of very special cases the elastic continuum model is devoid of useful results, while for constant nominal spin rate the distributed-mass finite-element model is quite generally tractable, since in the latter case the governing equations are always linear, constant-coefficient, ordinary differential equations. Although with both of these alternatives the details of the formulation generally obscure the essence of the problem and permit very little engineering insight to be gained without extensive computation, this difficulty is not encountered when dealing with simple concentrated mass models.

Likins, P. W.↗
NASA NTRS1990

One-dimensional hybrid continuum/particle simulation approach for rarefied hypersonic flows

A hybrid numerical technique, based on local application of continuum and particle methods, is applied to the one-dimensional normal shock structure in a monatomic gas. The methods are coupled iteratively via a flux boundary condition in a region where both are applicable. The validity of each method through the shock structure is quantified. The physical and numerical aspects of coupling are discussed and a variety of coupling methods are assessed. Results are presented to show the improvement in computational efficiency available via the hybrid approach.

Wadsworth, D. C.↗
NASA NTRS1994

X-ray refelection from photoionized media in active galactic nuclei

We calculate the spectrum of X-ray radiation and reprocessed by a partly ionized optically thick medium in an active galactic nucleus. We self-consistently calculate the ionization balance and thermal balance in the medium along with the distribution of X-ray intensity with optical depth. In addition to absorption or scattering of the incident X-rays, we also compute the spectrum of X-rays emitted by the material, including lines, edges, and bremsstrahlung. The albedo of the medium depends primarily on the X-ray ionization parameter (ratio of incident flux to gas density, zeta(sub Chi), and secondarily on the UV flux generated by dissipation inside the disk; we locate the critical range of zeta(sub Chi) over which the albedo increases from small to nearly unity. While the continuum reflection is very weak below 10 keV when zeta(sub Chi) is small, significnat fluxes are emitted in atomic lines and edges in this energy range. In the limit of large zeta(sub Chi), the albedo below 10 keV increases, but reflection in this band is never gray: some photoelectric absorption remains up to rather large values of zeta(sub Chi), while at still higher values, inverse Compton scattering amplifies the soft X-ray flux. These features are sufficiently sharp that current and near-future X-ray experiments should permit diagnostic measures of zeta(sub Chi).

Zycki, Piotr T.↗
NASA NTRS2001

Higher-Order Theory for Functionally Graded Materials

Functionally graded materials (FGM's) are a new generation of engineered materials wherein the microstructural details are spatially varied through nonuniform distribution of the reinforcement phase(s). Engineers accomplish this by using reinforcements with different properties, sizes, and shapes, as well as by interchanging the roles of the reinforcement and matrix phases in a continuous manner (ref. 1). The result is a microstructure that produces continuously or discretely changing thermal and mechanical properties at the macroscopic or continuum scale. This new concept of engineering the material's microstructure marks the beginning of a revolution both in the materials science and mechanics of materials areas since it allows one, for the first time, to fully integrate the material and structural considerations into the final design of structural components. Functionally graded materials are ideal candidates for applications involving severe thermal gradients, ranging from thermal structures in advanced aircraft and aerospace engines to computer circuit boards. Owing to the many variables that control the design of functionally graded microstructures, full exploitation of the FGM's potential requires the development of appropriate modeling strategies for their response to combined thermomechanical loads. Previously, most computational strategies for the response of FGM's did not explicitly couple the material's heterogeneous microstructure with the structural global analysis. Rather, local effective or macroscopic properties at a given point within the FGM were first obtained through homogenization based on a chosen micromechanics scheme and then subsequently used in a global thermomechanical analysis.

Aboudi, J.↗
NASA NTRS1987

Anticipation of the landing shock phenomenon in flight simulation

An aircraft landing may be described as a controlled crash because a runway surface is intercepted. In a simulation model the transition from aerodynamic flight to weight on wheels involves a single computational cycle during which stiff differential equations are activated; with a significant probability these initial conditions are unrealistic. This occurs because of the finite cycle time, during which large restorative forces will accompany unrealistic initial oleo compressions. This problem was recognized a few years ago at Ames Research Center during simulation studies of a supersonic transport. The mathematical model of this vehicle severely taxed computational resources, and required a large cycle time. The ground strike problem was solved by a described technique called anticipation equations. This extensively used technique has not been previously reported. The technique of anticipating a significant event is a useful tool in the general field of discrete flight simulation. For the differential equations representing a landing gear model stiffness, rate of interception and cycle time may combine to produce an unrealistic simulation of the continuum.

Mcfarland, Richard E.↗
NASA NTRS1972

X-ray study of solar plage regions and a small flare.

The O VII and Ne IX resonance line and 5 keV continuum emission of a small solar flare are interpreted on the basis of a model in which the emission measure of the region is an exponentially decreasing function of temperature. Good agreement with all three sets of data is obtained - a result not possible on the basis of an isothermal approximation. The change in the emission of the flare over a two minute interval appears to be consistent with a simple cooling of the region. Computation of the contribution function versus temperature for each emission indicates that each predominantly comes from a physically separate volume within the feature. The emission measure distribution for the flare is distinctly flatter than that of five plage regions which tend to have similar slopes and which are in turn flatter than the general corona.

Acton, L. W.↗
NASA NTRS1996

Spray Combustion Modeling with VOF and Finite-Rate Chemistry

A spray atomization and combustion model is developed based on the volume-of-fluid (VOF) transport equation with finite-rate chemistry model. The gas-liquid interface mass, momentum and energy conservation laws are modeled by continuum surface force mechanisms. A new solution method is developed such that the present VOF model can be applied for all-speed range flows. The objectives of the present study are: (1) to develop and verify the fractional volume-of-fluid (VOF) cell partitioning approach into a predictor-corrector algorithm to deal with multiphase (gas-liquid) free surface flow problems; (2) to implement the developed unified algorithm in a general purpose computational fluid dynamics (CFD) code, Finite Difference Navier-Stokes (FDNS), with droplet dynamics and finite-rate chemistry models; and (3) to demonstrate the effectiveness of the present approach by simulating benchmark problems of jet breakup/spray atomization and combustion. Modeling multiphase fluid flows poses a significant challenge because a required boundary must be applied to a transient, irregular surface that is discontinuous, and the flow regimes considered can range from incompressible to highspeed compressible flows. The flow-process modeling is further complicated by surface tension, interfacial heat and mass transfer, spray formation and turbulence, and their interactions. The major contribution of the present method is to combine the novel feature of the Volume of Fluid (VOF) method and the Eulerian/Lagrangian method into a unified algorithm for efficient noniterative, time-accurate calculations of multiphase free surface flows valid at all speeds. The proposed method reformulated the VOF equation to strongly couple two distinct phases (liquid and gas), and tracks droplets on a Lagrangian frame when spray model is required, using a unified predictor-corrector technique to account for the non-linear linkages through the convective contributions of VOF. The discontinuities within the sharp interface will be modeled as a volume force to avoid stiffness. Formations of droplets, tracking of droplet dynamics and modeling of the droplet breakup/evaporation, are handled through the same unified predictor-corrector procedure. Thus the new algorithm is non-iterative and is flexible for general geometries with arbitrarily complex topology in free surfaces. The FDNS finite-difference Navier-Stokes code is employed as the baseline of the current development. Benchmark test cases of shear coaxial LOX/H2 liquid jet with atomization/combustion and impinging jet test cases are investigated in the present work. Preliminary data comparisons show good qualitative agreement between data and the present analysis. It is indicative from these results that the present method has great potential to become a general engineering design analysis and diagnostics tool for problems involving spray combustion.

Chen, Yen-Sen↗
NASA NTRS2004

Anomolous Fatigue Crack Growth Phenomena in High-Strength Steel

The growth of a fatigue crack through a material is the result of a complex interaction between the applied loading, component geometry, three-dimensional constraint, load history, environment, material microstructure and several other factors. Previous studies have developed experimental and computational methods to relate the fatigue crack growth rate to many of the above conditions, with the intent of discovering some fundamental material response, i.e. crack growth rate as a function of something. Currently, the technical community uses the stress intensity factor solution as a simplistic means to relate fatigue crack growth rate to loading, geometry and all other variables. The stress intensity factor solution is a very simple linear-elastic representation of the continuum mechanics portion of crack growth. In this paper, the authors present fatigue crack growth rate data for two different high strength steel alloys generated using standard methods. The steels exhibit behaviour that appears unexplainable, compared to an aluminium alloy presented as a baseline for comparison, using the stress intensity factor solution.

Forth, Scott C.↗
NASA NTRS1992

A continuum method for modeling surface tension

In the novel method presented for modeling the effects of surface tension on fluid motion, the interfaces between fluids with different, color-represented properties are finite-thickness transition regions across which the color varies continuously. A force density proportional to the surface curvature of constant color is defined at each point in the transition region; this force-density is normalized in such a way that the conventional description of surface tension on an interface is recovered when the ratio of local transition-reion thickness to local curvature radius approaches zero. The properties of the method are illustrated by computational results for 2D flows.

Brackbill, J. U.↗
NASA NTRS1996

Atomic and Molecular Data for Interstellar Studies: A Status Report

Most interstellar species have a large fraction of their electronic transitions at far ultraviolet wavelengths. Observations at these wavelengths reveal spectra rich in absorption lines seen against the continuum of a background source, such as a hot star in our Galaxy, a supernova in a nearby galaxy, or even a bright nucleus in an active galaxy. Most of the observations continue to be made with space-borne instruments, but recent work includes measurements of extragalactic material at large redshifts obtained at high resolution with large ground-based telescopes (e.g., the Keck Telescope). The combination of precise experimental oscillator strengths, large-scale computations, and astronomical spectra with high signal-to-noise ratios are yielding a set of self-consistent-values that span a range in strength in excess of 100 for more and more species. The large range is important for studies involving the different environments probed by the various background sources. This review highlights recent work on the atomic species. Si II, S I, and Fe II, and on the molecules, CO and C2.

Federman, Steven R.↗
NASA NTRS1990

Galerkin approximations for dissipative magnetohydrodynamics

A Galerkin approximation scheme is proposed for voltage-driven, dissipative magnetohydrodynamics. The trial functions are exact eigenfunctions of the linearized continuum equations and represent helical deformations of the axisymmetric, zero-flow, driven steady state. The lowest nontrivial truncation is explored: one axisymmetric trial function and one helical trial function each for the magnetic and velocity fields. The system resembles the Lorenz approximation to Benard convection, but in the region of believed applicability, its dynamical behavior is rather different, including relaxation to a helically deformed state similar to those that have emerged in the much higher resolution computations of Dahlburg et al.

Chen, Hudong↗