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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

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.↗
NASA NTRS1996

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↗
NASA NTRS2001

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.↗
NASA NTRS1983

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.↗
NASA NTRS2012

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↗
DOE OSTI2023

Predictive scale-bridging simulations through active learning

Throughout computational science, there is a growing need to utilize the continual improvements in raw computational horsepower to achieve greater physical fidelity through scale-bridging over brute-force increases in the number of mesh elements. For instance, quantitative predictions of transport in nanoporous media, critical to hydrocarbon extraction from tight shale formations, are impossible without accounting for molecular-level interactions. Similarly, inertial confinement fusion simulations rely on numerical diffusion to simulate molecular effects such as non-local transport and mixing without truly accounting for molecular interactions. With these two disparate applications in mind, we develop a novel capability which uses an active learning approach to optimize the use of local fine-scale simulations for informing coarse-scale hydrodynamics. Our approach addresses three challenges: forecasting continuum coarse-scale trajectory to speculatively execute new fine-scale molecular dynamics calculations, dynamically updating coarse-scale from fine-scale calculations, and quantifying uncertainty in neural network models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

From classical thermodynamics to phase-field method

Phase-field method is a density-based computational method at the mesoscale for modeling and predicting the temporal microstructure and property evolution during materials processes. The focus of this article is on connecting the most common phase-field equations to the very basic first and second laws of classical thermodynamics through rudimentary irreversible thermodynamics. It briefly discusses the relations of the continuum phase-field equations to their counter parts at the microscopic and atomic levels. It attempts to clarify the contributions of long-range elastic, electrostatic, and magnetic interactions to domain structure evolution during structural, ferroelectric, and ferromagnetic phase transformations by separating order parameter changes due to the presence of quasi-static fields and those arising from phase transformations. A few examples are presented to demonstrate the possibility of employing the phase-field method to provide guidance to designing materials for optimum properties or discovering novel mesoscale phenomena or new materials functionalities. Here, the article ends with a brief perspective on a number of potential future directions on the development and applications of phase-field method beyond its traditional applications to structural alloys.

36 MATERIALS SCIENCE↗
NASA NTRS1990

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.↗
DOE OSTI2022

Uncertainty quantification in breakup reactions

Breakup reactions are one of the favored probes to study loosely bound nuclei, particularly those in the limit of stability forming a halo. In order to interpret such breakup experiments, the continuum discretized coupled channel method is typically used. In this study, the first Bayesian analysis of a breakup reaction model is performed. Herein we use a combination of statistical methods together with a three-body reaction model (the continuum discretized coupled channel method) to quantify the uncertainties on the breakup observables due to the parameters in the effective potential describing the loosely bound projectile of interest. The combination of tools we develop opens the path for a Bayesian analysis of not only breakup processes, but also a wide array of complex processes that require computationally intensive reaction models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2004

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.↗
NASA NTRS1995

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.↗
NASA NTRS2000

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.↗
NASA NTRS2018

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↗
DOE OSTI2024

Assessing Critical Conditions for Scour Near Obstructions using Bed Shear, Particle Onset of Motion Balances, and CFD-DEM Modeling of Granular Beds

Computational Fluid Dynamics combined with a Discrete Element Method is one of the computational methods that can be used to model multiphase flows. In this method various phases, gas and liquid or solid, are present in the same computational domain. The local averaged Navier–Stokes equations determine the flow of the continuous phase fluid and are solved using the traditional CFD finite volume approach. DEM is based on a Lagrangian formulation, which solves the equations of motion, expressed in ordinary differential equations, for representative particles as they move in space and time. The interactions between the continuous fluid phase and discrete solid phase are modeled with the use of Newton’s laws of motion via drag force. The particles interact with each other and with the boundaries of the fluid continuum, and the resulting contact forces are included in the equations of motion. The properties of solid particles and boundaries are treated as elastic bodies, with specified density, elastic modulus, and Poisson’s ratio. Particle shapes may vary from single spherical particles to more complex-shaped composite particles. The particles may be introduced into the domain by random or structured injection at a point, surface, or volume, depending on the application. More details on the formulation can be found in the Simcenter STAR-CCM+ User’s Manual and OpenFOAM website.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1972

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.↗
NASA NTRS1977

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.↗
NASA NTRS1993

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.↗
NASA NTRS1995

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↗