SEARCH · Engineering Papers
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.
The equilibrium strategy for N-person differential games can be found by studying a min-max problem subject to differential systems constraints. The differential constraints are penalized and finite elements are used to compute numerical solutions. Convergence proof and error estimates are given. Numerical results are also included and compared with those obtained by the dual method.
Asymptotic flows inside curved ducts of rectangular as well as polar cross sections have been analyzed using the Navier-Stokes equations in terms of the axial velocity and vorticity and the cross-flow stream function. Numerical solutions of the three second-order coupled elliptic partial differential equations governing this flow have been obtained using the multi-grid-SI scheme, to determine fine-grid results for configurations with strong secondary flow.
The fully differential cross section for photon-electron pair production is integrated numerically over phase space. Results are obtained for the astrophysically interesting case in which the interaction between an ultrarelativistic electron and a soft photon results in electron-positron pair production. The positron spectrum is a function of the energies of both the photon and the electron, as well as the angle of interaction. It is found that the energy at which the positron distribution peaks is inversely proportional to the photon energy and independent of the electron energy. The positron spectrum is integrated once more over initial electron energies for a power-law energy distribution of primary electrons. The same procedure is repeated for the recoil particle; it is shown that the peak of the recoil energy distribution depends linearly on the energy of the primary electron. Finally, semianalytical expressions are obtained for the energy losses of the primary electrons.
Conjugate residual methods for the solution of spectral equations are described. An inexact finite-difference operator is introduced as a preconditioner in the iterative procedures. Application of these techniques is limited to problems for which the symmetric part of the coefficient matrix is positive definite. Although the spectral equation is a very ill-conditioned and full matrix problem, the computational effort of the present iterative methods for solving such a system is comparable to that for the sparse matrix equations obtained from the application of either finite-difference or finite-element methods to the same problems. Numerical experiments are shown for a self-adjoint elliptic partial differential equation with Dirichlet boundary conditions, and comparison with other solution procedures for spectral equations is presented.
To predict small-particle diffusional mass transfer (deposition), including particle thermophoresis, transpiration cooling, and variable properties, the coupled ordinary differential equations governing self-similar laminar boundary layers are solved numerically. Under typical combustion turbine conditions, although diffusional deposition rates can be dramatically reduced by transpiration cooling (e.g., by some 5-decades for mainstream submicron particles corresponding to a Schmidt number of about 100 and a wall transpiration-cooled to Tw/Te = 0.8), actual deposition rate reductions will be smaller than previously expected (by about 1 decade for particles with Sc of about 100), owing to thermophoretic particle drift caused by the colder wall. Such microdroplets, small enough to behave like heavy molecules in combustion systems, are often important because they can cause adherence of the much larger ash particles which inertially impact on the same surface.
The issue of adopting the velocity components as dependent velocity variables for the Navier-Stokes flow computations is investigated. The viewpoint advocated is that a numerical algorithm should preferably honor both the physical conservation law in differential form and the geometric conservation law in discrete form. With the use of Cartesian velocity vector, the momentum equations in curvilinear coordinates can retain the full conservation-law form and satisfy the physical conservation laws. With the curvilinear velocity components, source terms appear in differential equations and hence the full conservation law form can not be retained. In discrete expressions, algorithms based on the Cartesian components can satisfy the geometric conservation-law form for convection terms but not for viscous terms; those based on the curvilinear components, on the other hand, cannot satisfy the geometric conservation-law form for either convection or viscous terms. Several flow solutions for domain with 90 and 360 degree turnings are presented to illustrate the issues of using the Cartesian velocity components and the staggered grid arrangement.
This paper discusses 'exotic' rotationally symmetric containers that admit an entire continuum of distinct equilibrium capillary free surfaces. The paper extends earlier work to a larger class of parameters and clarifies and simplifies the governing differential equations, while expressing them in a parametric form appropriate for numerical integration. A unified presentation suitable for both zero and nonzero gravity is given. Solutions for the container shapes are depicted graphically along with members of the free-surface continuum, and comments are given concerning possible physical experiments.
An algorithm for the computer simulation of large space structures under active control is considered. Linked lists are used in a matrix data structure to implement the trapezoidal rule on system differential equations. The use of the trapezoidal rule ensures that the numerical stability is equivalent to the system stability, which is essential for this type of simulation. The sparsity of the system matrices is exploited by the linked lists, and the algorithm efficiently steps through the lists in an orderly fashion. Results of simulations on a NASA large space structure experiment are reported.
An algorithm for the computer simulation of large space structures under active control is considered. Linked lists are used in a matrix data structure to implement the trapezoidal rule on the system differential equations. The use of the trapezoidal rule ensures that the numerical stability is equivalent to the system stability, which is essential for this type of simulation. The sparsity of the system matrices is exploited by the linked lists, and the algorithm efficiently steps through the lists in an orderly fashion. Results of simulations on a NASA large space structure experiment are reported.
Important characteristics of the aeroacoustic wave propagation are mostly encoded in their dispersion relations. Hence, a computational aeroacoustic (CAA) algorithm, which reasonably preserves these relations, was investigated. It was derived using an optimization procedure to ensure, that the numerical derivatives preserved the wave number and angular frequency of the differential terms in the linearized, 2-D Euler equations. Then, simulations were performed to validate the scheme and a compatible set of discretized boundary conditions. The computational results were found to agree favorably with the exact solutions. The boundary conditions were transparent to the outgoing waves, except when the disturbance source was close to a boundary. The time-domain data generated by such CAA solutions were often intractable until their spectra was analyzed. Therefore, the relative merits of three different methods were included in the study. For simple, periodic waves, the periodogram method produced better estimates of the steep-sloped spectra than the Blackman-Tukey method. Also, for this problem, the Hanning window was more effective when used with the weighted-overlapped-segment-averaging and Blackman-Tukey methods gave better results than the periodogram method. Finally, it was demonstrated that the representation of time domain-data was significantly dependent on the particular spectral analysis method employed.
The objectives are: Conduct a parametric theoretical and numerical investigation of vibro-convective buoyancy-driven flow in differentially heated cylindrical containers. Investigate buoyant vibro-convective transport regimes in Bridgman-type systems with a focus on the use of vibration to suppress, or control, convection in order to achieve transport control during crystal growth. Assess the feasibility of vibro-convective control as a means of offsetting "g-jitter" effects under microgravity conditions, Exchange information with the experimental group at the General Physics Institute (GPI) of the Russian Academy of Science who are undertaking a complementary experimental program.
One-dimensional water-hammer modeling involves the solution of two coupled non-linear hyperbolic partial differential equations (PDEs). These equations result from applying the principles of conservation of mass and momentum to flow through a pipe, and usually the assumption that the speed at which pressure waves propagate through the pipe is constant. In order to solve these equations for the interested quantities (i.e. pressures and flow rates), they must first be converted to a system of ordinary differential equations (ODEs) by either approximating the spatial derivative terms with numerical techniques or using the Method of Characteristics (MOC). The MOC approach is ideal in that no numerical approximation errors are introduced in converting the original system of PDEs into an equivalent system of ODEs. Unfortunately this resulting system of ODEs is bound by a time step constraint so that when integrating the equations the solution can only be obtained at fixed time intervals. If the fluid system to be modeled also contains dynamic components (i.e. components that are best modeled by a system of ODEs), it may be necessary to take extremely small time steps during certain points of the model simulation in order to achieve stability and/or accuracy in the solution. Coupled together, the fixed time step constraint invoked by the MOC, and the occasional need for extremely small time steps in order to obtain stability and/or accuracy, can greatly increase simulation run times. As one solution to this problem, a method for combining variable step integration (VSI) algorithms with the MOC was developed for modeling water-hammer in systems with highly dynamic components. A case study is presented in which reverse flow through a dual-flapper check valve introduces a water-hammer event. The predicted pressure responses upstream of the check-valve are compared with test data.
Simulations of airspace operational concepts can play a significant role in determining future paradigms that would allow for a safe increase in airspace density. In particular, airspace simulations which are capable of handling large numbers of aircraft act as an enabling capability for the testing of proposed airspace operational concepts. Simulations allowing for gradient based optimization methods are particularly attractive, since they would potentially allow for an efficient and empirical means to derive best operational practices. These could also allow for vehicle multidisciplinary design and optimization studies to include air traffic management considerations as a discipline. But any large scale simulation of airspace operations must include some methodology for addressing airspace separation requirements, which in the most direct sense would be tracked in a manner that computationally grows as a quadratic function of the number of simulated aircraft. Efficient indirect methods have been developed in certain contexts to address this limitation. However, any means of addressing separation requirements in a gradient based optimization context should be implemented by functions which provide analytic derivative information to maximize numerical precision and computational efficiency. In this paper, a fast and differentiable separation metric is described in application to gradient based optimization of airspace operations. Rather than computing the separation distance between every pair of aircraft in a simulation, this method effectively reduces the problem to a smaller relevant set using a geometric decomposition. This method guarantees that the smallest distance at all points in simulated time is determined exactly. When used in an optimization constraint context, this guarantees that a minimum separation is maintained between all pairs of aircraft. The presented metric has logarithmic computational growth with respect to the number of simulated aircraft, and is shown to perform well in a series of notional 2D airspace optimization problems when used to enforce specified airborne separation constraints. Results show that this is notably faster than a direct pairwise distance computing metric for optimizations involving both small and large numbers of aircraft, yet enforce separation requirements to the same tolerance. It is shown that this favorable scalability is an enabling capability for more sophisticated air traffic management conceptual studies.
Geodetic and geophysical investigations of the Galilean moon Callisto address fundamental questions regarding the formation and evolution of the Jovian system. Callisto's evolution and internal structure appear to significantly differ from the other Jovian satellites. Similarly-sized Ganymede is a highly evolved ice-rock moon with a differentiated interior, intrinsic magnetic field, and abundant surface evidence of internal activity. In contrast, Callisto's surface is ancient, and Galileo spacecraft data suggest its interior is only incompletely differentiated, despite the presumed presence of a sub-surface ocean. These properties make Callisto uniquely able to constrain the timing and nature of the Jovian system formation. The Magnetics, Altimetry, Gravity, and Imaging of Callisto (MAGIC) mission concept is conceived to fully characterize the properties of this enigmatic moon from its deep interior to the icy shell. Three main instruments are included as a scientific payload. Highly accurate measurements of Callisto's topography, magnetic field, and morphology are obtained by the onboard laser altimeter, magnetometer, and camera, respectively. The telecommunication system supports an additional gravity and radio science investigation. Long- and short-wavelength gravity anomalies afford powerful constraints on internal differentiation and the properties of the hydrosphere (water and ice). Comprehensive numerical simulations and covariance analyses of MAGIC mission scenarios presented in this paper show that the gravitational degree-2 normalized coefficients and the pole obliquity enable the determination of the moment of inertia with an accuracy better than 0.015%. The combination of gravity and altimetry measurements acquired by MAGIC are essential to the characterization of Callisto's interior if – as is likely – the degree-2 gravity includes non-hydrostatic terms. MAGIC's radio science data yield the estimation of Callisto's gravity field with spatial resolutions of <100 km. The combination of gravitational and deformation tides that are retrieved by the radio science and altimetry investigations, respectively, leads to the recovery of the rigid ice shell thickness to within ∼3 km. Together these datasets would resolve ambiguities inherent in Galileo flyby data, revealing Callisto's interior structure as well as the existence and properties of its postulated internal ocean.
During the past six months we have conducted significant research in several domains in order to clarify and understanding the aging and chemical failure mechanism of thermoplastics (PVDF or Tefzel) for pipes. We organized numerous analytical studies with methods including Fourier Transform Infrared Spectroscopy, Dynamic Mechanical Analysis, Differential Scanning Calorimetry, and Stress Relaxation experiments. In addition we have reanalyzed previous thermogravimetric data concerning the rate of deplasticization of Coflon pipe. We investigated a number of aged samples of both Tefzel and Coflon that were forwarded from MERL. We conducted stress relaxation experiments of Coflon pipe at several temperatures and determined an activation energy. We also examined the dynamic mechanical response PVDF during deplasticization and during methanol plasticization. We performed numerous DSC analyses to research the changing crystalline morphology. We have noted significant changes in crystallinity upon aging for both PVDF and Tefzel. Little variation in elemental composition was noted for many of the aged Coflon and Tefzel samples tested.
Lebesgue integration of derivatives of strongly-oscillatory functions is a recurring challenge in computational science and engineering. Integration by parts is an effective remedy for huge computational costs associated with Monte Carlo integration schemes. In case of Lebesgue integrals over a smooth manifold, however, integration by parts gives rise to a derivative of the density implied by charts describing the domain manifold. This paper focuses on the computation of that derivative, which we call the density gradient function, on general smooth manifolds. We analytically derive formulas for the density gradient and present examples of manifolds determined by popular differential equation-driven systems. Furthermore, we highlight the significance of the density gradient by demonstrating a numerical example of Monte Carlo integration involving oscillatory integrands.
The derivatives of the spectra are commonly used for quantification in Auger Electron Spectroscopy (AES) spectra, while the derivative of the KLL C Auger line has proven to be valuable in obtaining a measure of the relative proportions of sp 2 ‐ and sp 3 ‐hybridization using the D‐parameter in both AES and X‐ray Photoelectron Spectroscopy (XPS). Differentiation of X‐ray Photoelectron Spectroscopy (XPS) and Auger Electron Spectroscopy (AES) spectra by numerical means is presented and illustrated for polymeric, such as PEEK and Nylon, as well as for graphitic materials including highly ordered pyrolytic graphite and graphene oxide. The most commonly available Savitzky–Golay method is explained mathematically and developed through the case of constructing a 5‐point quadratic polynomial convolution kernel suitable for differentiating spectra of adequate signal to noise. The concept of differentiation of spectra where signal to noise is less than adequate is also developed. Two alternative strategies to Savitzky–Golay differentiation are presented, which fit curves to data that allow derivatives to be obtained where Savitzky–Golay would otherwise fail. These alternative methods involve constructing a parametric curve that fits data over the entire energy interval of interest. Derivatives of spectra are then obtained by differentiating these parametric curves directly. A comparison of results for different materials for which specific sp 2 ‐ vs sp 3 ‐hybridized carbon proportions are of interest is used to emphasize the importance of characterizing methods used to differentiate spectra and understanding the characteristics of instrumentation used to measure spectra. The case for using Principal Component Analysis noise reduction with C KLL spectra is made for spectra collected from a heterogeneous graphene oxide sample.
Bridgman crystal growth can be conducted in the so-called "detached" solidification regime, where the growing crystal is detached from the crucible wall. A small gap between the growing crystal and the crucible wall, of the order of 100 micrometers or less, can be maintained during the process. A meniscus is formed at the bottom of the melt between the crystal and crucible wall. Under proper conditions, growth can proceed without collapsing the meniscus. The meniscus shape plays a key role in stabilizing the process. Thermal and other process parameters can also affect the geometrical steady-state stability conditions of solidification. The dynamic stability theory of the shaped crystal growth process has been developed by Tatarchenko. It consists of finding a simplified autonomous set of differential equations for the radius, height, and possibly other process parameters. The problem then reduces to analyzing a system of first order linear differential equations for stability. Here we apply a modified version of this theory for a particular case of detached solidification. Approximate analytical formulas as well as accurate numerical values for the capillary stability coefficients are presented. They display an unexpected singularity as a function of pressure differential. A novel approach to study the thermal field effects on the crystal shape stability has been proposed. In essence, it rectifies the unphysical assumption of the model that utilizes a perturbation of the crystal radius along the axis as being instantaneous. It consists of introducing time delay effects into the mathematical description and leads, in general, to stability over a broader parameter range. We believe that this novel treatment can be advantageously implemented in stability analyses of other crystal growth techniques such as Czochralski and float zone methods.