Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Differentiable physics”

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 577 records · Page 32

Geophysical data and the interior of the moon

Geophysical data and physical properties of the lunar interior are considered, giving attention to density, gravity field, viscosity, the strength of lunar materials, electrical conductivity, and magnetic properties. Seismic data and structure of the lunar interior are discussed, taking into account moonquakes and lunar tectonism, the velocity structure, and the compositional implications of the velocity structure. Questions regarding the thermal state and the evolution of the moon are also explored. The data and the models presented characterize the moon as a differentiated body which evolved relatively early in its history.

Toksoz, M. N.↗

Chemoviscosity modeling for thermosetting resins - I

A new analytical model for chemoviscosity variation during cure of thermosetting resins was developed. This model is derived by modifying the widely used WLF (Williams-Landel-Ferry) Theory in polymer rheology. Major assumptions involved are that the rate of reaction is diffusion controlled and is linearly inversely proportional to the viscosity of the medium over the entire cure cycle. The resultant first order nonlinear differential equation is solved numerically, and the model predictions compare favorably with experimental data of EPON 828/Agent U obtained on a Rheometrics System 4 Rheometer. The model describes chemoviscosity up to a range of six orders of magnitude under isothermal curing conditions. The extremely non-linear chemoviscosity profile for a dynamic heating cure cycle is predicted as well. The model is also shown to predict changes of glass transition temperature for the thermosetting resin during cure. The physical significance of this prediction is unclear at the present time, however, and further research is required. From the chemoviscosity simulation point of view, the technique of establishing an analytical model as described here is easily applied to any thermosetting resin. The model thus obtained is used in real-time process controls for fabricating composite materials.

Hou, T. H.↗

Research highlights: June 1990 - May 1991

Linear instability calculations at MSFC have suggested that the Geophysical Fluid Flow Cell (GFFC) should exhibit classic baroclinic instability at accessible parameter settings. Interest was in the mechanisms of transition to temporal chaos and the evolution of spatio-temporal chaos. In order to understand more about such transitions, high resolution numerical experiments for the physically simplest model of two layer baroclinic instability were conducted. This model has the advantage that the numerical code is exponentially convergent and can be efficiently run for very long times, enabling the study of chaotic attractors without the often devastating effects of low-order trunction found in many previous studies. Numerical algorithms for implementing an empirical orthogonal function (EOF) analysis of the high resolution numerical results were completed. Under conditions of rapid rotation and relatively low differential heating, convection in a spherical shell takes place as columnar banana cells wrapped around the annular gap, but with axes oriented along the axis of rotation; these were clearly evident in the GFFC experiments. The results of recent numerical simulations of columnar convection and future research plans are presented.

Source record↗

A Probabilistic Framework for the Validation and Certification of Computer Simulations

The paper presents a methodology for quantifying, propagating, and managing the uncertainty in the data required to initialize computer simulations of complex phenomena. The purpose of the methodology is to permit the quantitative assessment of a certification level to be associated with the predictions from the simulations, as well as the design of a data acquisition strategy to achieve a target level of certification. The value of a methodology that can address the above issues is obvious, specially in light of the trend in the availability of computational resources, as well as the trend in sensor technology. These two trends make it possible to probe physical phenomena both with physical sensors, as well as with complex models, at previously inconceivable levels. With these new abilities arises the need to develop the knowledge to integrate the information from sensors and computer simulations. This is achieved in the present work by tracing both activities back to a level of abstraction that highlights their commonalities, thus allowing them to be manipulated in a mathematically consistent fashion. In particular, the mathematical theory underlying computer simulations has long been associated with partial differential equations and functional analysis concepts such as Hilbert spares and orthogonal projections. By relying on a probabilistic framework for the modeling of data, a Hilbert space framework emerges that permits the modeling of coefficients in the governing equations as random variables, or equivalently, as elements in a Hilbert space. This permits the development of an approximation theory for probabilistic problems that parallels that of deterministic approximation theory. According to this formalism, the solution of the problem is identified by its projection on a basis in the Hilbert space of random variables, as opposed to more traditional techniques where the solution is approximated by its first or second-order statistics. The present representation, in addition to capturing significantly more information than the traditional approach, facilitates the linkage between different interacting stochastic systems as is typically observed in real-life situations.

Ghanem, Roger↗

The Oscillations and Stability of Differentially Rotating Neutron Stars

Studies of the oscillations and stability of neutron stars are motivated by the fact that vibrating neutron stars are a promising source of gravitational waves. One important factor is the influence of differential rotation, which is likely to arise in a neutron star at times, such as the immediate aftermath of the supernova, when we expect strong vibrations. I will discuss two phenomena unique to differentially rotating systems: dynamical shear instabilities, and the existence of a co-rotation band (a frequency band in which mode pattern speed matches the local angular velocity). Using a simple model, we have found dynamical shear instabilities that arise where modes cross into the co-rotation band, if the degree of differential rotation exceeds a certain threshold. We are currently investigating whether this mechanism operates in more realistic stellar models, and whether it is responsible for the dynamical instabilities occurring at low ratios of kinetic to potential energy that have been observed by several authors. I will present the latest results of these studies. Another topic of investigation is the nature of oscillations within the co-rotation band. The band gives rise to a continuous spectrum whose collective physical perturbation exhibits complicated temporal behaviour. I will also discuss the existence of modes within the continuous spectrum that appear physically indistinguishable from the discrete modes outside the band, despite the apparently singular nature of their eigenfunctions.

Watts, A.↗

Intelligent Engine Systems Work Element 1.3: Sub System Health Management

The objectives of this program were to develop health monitoring systems and physics-based fault detection models for engine sub-systems including the start, lubrication, and fuel. These models will ultimately be used to provide more effective sub-system fault identification and isolation to reduce engine maintenance costs and engine down-time. Additionally, the bearing sub-system health is addressed in this program through identification of sensing requirements, a review of available technologies and a demonstration of a demonstration of a conceptual monitoring system for a differential roller bearing. This report is divided into four sections; one for each of the subtasks. The start system subtask is documented in section 2.0, the oil system is covered in section 3.0, bearing in section 4.0, and the fuel system is presented in section 5.0.

Ashby, Malcolm↗

Process modifications for improved carbon fiber composites: Alleviation of the electrical hazards problem

Attempts to alleviate carbon-fiber-composite electrical hazards during airplane crash fires through fiber gasification are described. Thermogravimetric and differential scanning calorimetric experiments found several catalysts that caused fibers to combust when composites were exposed to test fires. Composites were tested in the 'Burn-Bang' apparatus and in high voltage electrical detection grid apparatus. In a standard three minute burn test modified composites released no fibers, while state-of-the-art composites released several hundred fiber fragments. Expected service life with and without catalytic modification was studied and electron microscopy and X-ray microanalysis furnished physical appearance and chemical composition data. An acrylic acid polymer fiber coating was developed that wet the carbon fiber surface uniformly with the catalyst, providing a marked contrast with the uneven coats obtained by solution-dipping.

Ramohalli, K.↗

Grid and design variables sensitivity analyses for NACA four-digit wing-sections

Two distinct parameterization procedures are developed for investigating the grid sensitivity with respect to design parameters of a wing-section example. The first procedure is based on traditional (physical) relations defining NACA four-digit wing-sections. The second is advocating a novel (geometrical) parameterization using spline functions such as NURBS (Non-Uniform Rational B-Splines) for defining the wing-section geometry. An interactive algebraic grid generation technique, known as Hermite Cubic Interpolation, is employed to generate C-type grids around wing-sections. The grid sensitivity of the domain with respect to design and grid parameters has been obtained by direct differentiation of the grid equations. A hybrid approach is proposed for more geometrically complex configurations. A comparison of the sensitivity coefficients with those obtained using a finite-difference approach has been made to verify the feasibility of the approach. The aerodynamic sensitivity coefficients are obtained using the compressible two-dimensional thin-layer Navier-Stokes equations.

Sadrehaghighi, Ideen↗

A mathematical model for simulating noise suppression of lined ejectors

A mathematical model containing the essential features embodied in the noise suppression of lined ejectors is presented. Although some simplification of the physics is necessary to render the model mathematically tractable, the current model is the most versatile and technologically advanced at the current time. A system of linearized equations and the boundary conditions governing the sound field are derived starting from the equations of fluid dynamics. A nonreflecting boundary condition is developed. In view of the complex nature of the equations, a parametric study requires the use of numerical techniques and modern computers. A finite element algorithm that solves the differential equations coupled with the boundary condition is then introduced. The numerical method results in a matrix equation with several hundred thousand degrees of freedom that is solved efficiently on a supercomputer. The model is validated by comparing results either with exact solutions or with approximate solutions from other works. In each case, excellent correlations are obtained. The usefulness of the model as an optimization tool and the importance of variable impedance liners as a mechanism for achieving broadband suppression within a lined ejector are demonstrated.

Watson, Willie R.↗

A physical and chemical model of early lunar history

A 'cool' thermal model of the moon's early history is discussed in terms of lunar petrology. Heat from the totally molten outer half of the moon's volume was, according to the model, lost to space and to the lunar interior, so that, barring additions of heat from external sources, all petrogenesis operating exclusively on material of the initially totally molten zone must have occured in an environment of decreasing temperatures. Mare basalts would result from hybridization by migration, mixing, and reequilibration of a variety of intercumulus liquids. Evidence is considered for the layered structure and a significant structural boundary that should result from differentiation of the approximately 350-km-thick initially totally molten zone. Magnetization of lunar rocks is considered.

Hubbard, N. J.↗

Computer simulation of the mathematical modeling involved in constitutive equation development: Via symbolic computations

Development of new material models for describing the high temperature constitutive behavior of real materials represents an important area of research in engineering disciplines. Derivation of mathematical expressions (constitutive equations) which describe this high temperature material behavior can be quite time consuming, involved and error prone; thus intelligent application of symbolic systems to facilitate this tedious process can be of significant benefit. A computerized procedure (SDICE) capable of efficiently deriving potential based constitutive models, in analytical form is presented. This package, running under MACSYMA, has the following features: partial differentiation, tensor computations, automatic grouping and labeling of common factors, expression substitution and simplification, back substitution of invariant and tensorial relations and a relational data base. Also limited aspects of invariant theory were incorporated into SDICE due to the utilization of potentials as a starting point and the desire for these potentials to be frame invariant (objective). Finally not only calculation of flow and/or evolutionary laws were accomplished but also the determination of history independent nonphysical coefficients in terms of physically measurable parameters, e.g., Young's modulus, was achieved. The uniqueness of SDICE resides in its ability to manipulate expressions in a general yet predefined order and simplify expressions so as to limit expression growth. Results are displayed when applicable utilizing index notation.

Arnold, S. M.↗

The 1992/93 eclipse of 31 Cygni

Extensive new ultraviolet and optical spectra of an atmospheric eclipse define the physical properties throughout the wind and chromosphere of 31 Cyg. These data require mass loss of approximately 3 x 10(exp -8) solar masses/yr in a wind that may merge smoothly into the chromosphere. Considerations of how energy is injected into the wind, however, suggest that the chromosphere and wind are separate structures. Most, if not all, of the velocity structure in metallic lines, which we have heretofore simulated with Doppler widths in the range 15-25 km/s, results from differential expansion of the atmosphere. Electron densities in the inner R(sub star) of the chromosphere are in the range 1.5 x 10(exp 9) - 2 x 10(exp 8)/cu cm, which implies clumping of the gas. The ionization of oxygen and nitrogen is consistent with clumping by factors of 3-30 in the outer chromosphere, roughly the amount required to give enough gas pressure to support the chromosphere. Chromospheric gas in 31 Cyg becomes hotter with increasing height, thus with decreasing optical depth, in a way that seems similar for all the zeta Aur binaries. Excitation temperature for Fe II in 1992/93 rises from about 5000 K at the deepest points sampled to about 12,500 K high in the wind. Strengths of violet Balmer lines give an excitation temperature for hydrogen of 6200-6500 K above a radial mass column density of approximately 0.01 g/sq cm. This amount of excitation implies that Lyman-alpha is thermalized beneath about 4 x 10(exp -3) g/sq cm. The outer atmosphere was symmetrical to within a factor of 2 in 1992/93, although it was clearly variable at this level, and it had similar mass column densities as in 1982. One manifestation of the variability was a flow toward the B star at phases 0.013-0.022 spanning velocities 45-100 km/s. Several lines of evidence point to a complicated and variable ionization in the wind: At large distance from the K star, measured mass column densities are less by up to a factor of 3 than required by a smoothly flowing wind. Also, much of the gas beyond r = 350 solar radii = 1.75 R(sub star) has a negative radial velocity. Radial velocities of the shell lines imply the outer atmosphere rotates in the direction of orbital motion, possibly through the deflection of wind flow lines in this direction.

Eaton, Joel A.↗

An analysis of turbulent diffusion flame in axisymmetric jet

The kinetic theory of turbulent flow was employed to study the mixing limited combustion of hydrogen in axisymmetric jets. The integro-differential equations in two spatial and three velocity coordinates describing the combustion were reduced to a set of hyperbolic partial differential equations in the two spatial coordinates by a binodal approximation. The MacCormick's finite difference method was then employed for solution. The flame length was longer than that predicted by the flame-sheet analysis, and was found to be in general agreement with a recent experimental result. Increase of the turbulence energy and scale resulted in an enhancement of the combustion rate and, hence, in a shorter flame length. Details of the numerical method as well as of the physical findings are discussed.

Chung, P. M.↗

Evolution of inclined galactic gas disks. I - A cloud-fluid approach

In this paper, numerical models are developed to follow the time-dependent behavior of inclined gas disks in nonspherical galaxies. Differential equations are developed which describe the time rate of change of the orbits of mass elements composing a thin, smoothly warping and twisting gas disk. These equations include multiple moments due to the nonsphericity of the galactic mass distribution, Coriolis effects arising from determination of orbits in a reference frame which tumbles with the surface figure of the galaxy, and viscous forces arising from cloud-cloud collisions. A cloud-fluid approach is used in determining the viscous forces. Physical arguments are used to determine the maximum permissible local value of the coefficient of kinematic viscosity. For small initial inclinations, analytic expressions are obtained for the inclination as a function of time and for the settling time, and the features of the solution are discussed.

Steiman-Cameron, Thomas Y.↗

An analysis of laminar free-convection flow and heat transfer about a flat plate paralled to the direction of the generating body force

The free-convection flow and heat transfer (generated by a body force) about a flat plate parallel to the direction of the body force are formally analyzed and the type of flow is found to be dependent on the Grashof number alone. For large Grashof numbers (which are of interest in aeronautics), the flow is of the boundary-layer type and the problem is reduced in a formal manner, which is analogous to Prandtl's forced-flow boundary-layer theory, to the simultaneous solution of two ordinary differential equations subject to the proper boundary conditions. Velocity and temperature distributions for Prandtl numbers of 0.01, 0.72, 0.733, 1, 1, 10, 100, and 1000 are computed, and it is shown that velocities and Nusselt numbers of the order of magnitude of those encountered in forced-convection flows may be obtained in free-convection flows. The theoretical and experimental velocity and temperature distributions are in good agreement. A flow and a heat-transfer parameter, from which the important physical quantities such as shear stress and heat-transfer rate can be computed, are derived as functions of Prandtl number alone.

Ostrach, Simon↗

Spherical means of solutions of partial differential equations in a conical region

The spherical means of the solutions of a linear partial differential equation Lu = f in a conical region are studied. The conical region is bounded by a surface generated by curvilinear xi lines and by two truncating xi surfaces. The spherical mean is the average of u over a constant xi surface. Conditions on the linear differential operator, L, and on the orthogonal coordinates xi, eta, and zeta are established so that the problem for the determination of the spherical mean of the solution subjected to the appropriate boundary and initial conditions can be reduced to a problem with only one space variable. Conditions are then established so that the spherical mean of the solution in one conical region will be proportional to that of a known solution in another conical region. Applications to various problems of mathematical physics and their physical interpretations are presented.

Ting, L.↗

Finite Element Modeling of the NASA Langley Aluminum Testbed Cylinder

The NASA Langley Aluminum Testbed Cylinder (ATC) was designed to serve as a universal structure for evaluating structural acoustic codes, modeling techniques and optimization methods used in the prediction of aircraft interior noise. Finite element models were developed for the components of the ATC based on the geometric, structural and material properties of the physical test structure. Numerically predicted modal frequencies for the longitudinal stringer, ring frame and dome component models, and six assembled ATC configurations were compared with experimental modal survey data. The finite element models were updated and refined, using physical parameters, to increase correlation with the measured modal data. Excellent agreement, within an average 1.5% to 2.9%, was obtained between the predicted and measured modal frequencies of the stringer, frame and dome components. The predictions for the modal frequencies of the assembled component Configurations I through V were within an average 2.9% and 9.1%. Finite element modal analyses were performed for comparison with 3 psi and 6 psi internal pressurization conditions in Configuration VI. The modal frequencies were predicted by applying differential stiffness to the elements with pressure loading and creating reduced matrices for beam elements with offsets inside external superelements. The average disagreement between the measured and predicted differences for the 0 psi and 6 psi internal pressure conditions was less than 0.5%. Comparably good agreement was obtained for the differences between the 0 psi and 3 psi measured and predicted internal pressure conditions.

Grosveld, Ferdinand W.↗

Ionosphere-exosphere coupling through charge exchange and momentum transfer in hydrogen-proton collisions

The implications of a traditional assumption of exospheric physics, that collisions of hydrogen atoms and protons preferentially result in charge exchange with negligible momentum transfer are examined. Initially adopted as a necessary convenience to accommodate limited computer resources in exosphere model calculations, this approximation results in a direct transformation of the proton velocity distribution into a hot component of neutral hydrogen. With expanding computational facilities, the need for the approximation has passed. As the first step toward its replacement with a realistic, quantum mechanical model of the H - H(+) collision process, differential and cumulative cross sections were calculated for quantum elastic scattering of indistinguishable nuclei for a fine grid of encounter energies and scattering angles. These data are used to study the nature of ionosphere-exosphere coupling through H - H(+) collisions, and to demonstrate that the distribution of velocities of scattered H produced in the traditional exospheric charge exchange approximation, as well as that arising from an alternative, fluid dynamic approach, leads to unacceptable abundances of coronal atoms in long-term, highly elliptic trajectories.

Hodges, R. R., Jr.↗