Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Analytic gradients”

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 505 records · Page 28

Reactive transport modeling of organic carbon degradation in marine methane hydrate systems

Abstract Natural methane hydrate has often been observed in sand layers that contain no particulate organic carbon (POC), but are surrounded by organic-rich, fine-grained marine muds. In this paper, we develop a reactive transport model (RTM) of a microbially-mediated set of POC degradation reactions, including hydrolysis of POC driven by extracellular enzymes, fermentation of the resulting high-molecular weight dissolved organic carbon (HMW-DOC), and methanogenesis that consumes low-molecular weight dissolved organic carbon (LMW-DOC). These processes are mediated by two groups of microbes, fermenters and methanogens that are heterogeneously distributed in different lithologies, with the largest numbers of microbes in the large pores of coarse-grained layers. We find that the RTM can reproduce methane hydrate occurrences observed in two different geological environments, at Walker Ridge Site 313-H (Gulf of Mexico) and IODP Site U1325 (Cascadia Margin). We also find that microbes can degrade POC even if they are physically separated, as extracellular enzymes and DOC can diffuse away from where they are produced by microbes. Microbial activity is highest at relatively early times after burial at shallow depths and near lithological boundaries, where concentration gradients transport solutes to intervals that contain the most microbes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Stability of beta-induced Alfvén eigenmodes (BAE) in DIII-D

Although the stability of ellipticity, toroidal and reversed-shear Alfvén eigenmodes (EAE, TAE, RSAE) are relatively well understood, less is known about the stability of lower-frequency modes such as the beta-induced Alfvén eigenmode (BAE) but, because they are often unstable in present devices and are implicated in fast-ion transport, understanding their stability is vital. BAE stability is studied in primarily weak or reversed shear DIII-D plasmas with sub-Alfvénic deuterium beams. Modes are classified based on electron cyclotron emission, beam emission spectroscopy, magnetics, and interferometer data. The study is limited to the initial two seconds of the discharge, where the evolving q profile provides an effective scan of the dependence of stability upon q. In a dedicated experiment, BAEs are unstable at times in the discharge when the minimum of the safety factor q min is close to a rational number. The observed mode frequencies are usually close to analytic estimates of the BAE accumulation point and the eigenfunction peaks in the vicinity of q min . Unstable BAEs usually occur in bursts that chirp rapidly in frequency. To isolate the importance of thermal and beam gradients in driving the modes, the beam and electron cyclotron heating power is altered for 50–100 ms durations in reproducible discharges. As expected from the resonance condition, BAEs depend sensitively on the beam power and injection geometry. Modes only persist for ~25 ms because the anisotropic beam population only interacts strongly with the modes over a relatively narrow range of q. Lastly, a database of over 1000 beam-heated discharges shows that BAEs are more likely to be unstable when the poloidal beta exceeds 0.5.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Turbulence-level dependence of cosmic ray parallel diffusion

ABSTRACT Understanding the transport of energetic cosmic rays belongs to the most challenging topics in astrophysics. Diffusion due to scattering by electromagnetic fluctuations is a key process in cosmic ray transport. The transition from a ballistic to a diffusive-propagation regime is presented in direct numerical calculations of diffusion coefficients for homogeneous magnetic field lines subject to turbulent perturbations. Simulation results are compared with theoretical derivations of the parallel diffusion coefficient’s dependences on the energy and the fluctuation amplitudes in the limit of weak turbulence. The present study shows that the widely used extrapolation of the energy scaling for the parallel diffusion coefficient to high turbulence levels predicted by quasi-linear theory does not provide a universally accurate description in the resonant-scattering regime. It is highlighted here that the numerically calculated diffusion coefficients can be polluted for low energies due to missing resonant interaction possibilities of the particles with the turbulence. Five reduced-rigidity regimes are established, which are separated by analytical boundaries derived in this work. Consequently, a proper description of cosmic ray propagation can only be achieved by using a turbulence-level-dependent diffusion coefficient and can contribute to solving the Galactic cosmic ray gradient problem.

79 ASTRONOMY AND ASTROPHYSICS↗

Improved method for calculating transonic velocities on blade-to-blade stream surfaces of a turbomachine

A method was developed to improve the accuracy of an existing computer program used to calculate transonic velocities on a blade-to-blade surface of a turbomachine. The method eliminates problems encountered in obtaining solutions with the velocity gradient equation when large gradients in velocity occur through the blade row. With the improved method, results indicate that the transonic solution can be obtained by scaling the velocities obtained at the reduced mass flow rate where all velocities are subsonic thereby eliminating the need for a solution of the velocity gradient equation. Solutions obtained with the scaling method on a two dimensional compressor cascade and an axial turbine stator show good agreement with experimental data. The results obtained for the stationary blade rows and comparison of analytical results obtained with and without the present method suggest that the method will yield an improved solution for centrifugal compressor impellers.

Wood, J. R.↗

Measurement of oil film thickness for application to elastomeric Stirling engine rod seals

The rod seal in the Stirling engine has the function of separating high pressure gas from low or ambient pressure oil. An experimental apparatus was designed to measure the oil film thickness distribution for an elastomeric seal in a reciprocating application. Tests were conducted on commercial elastomeric seals having a 76 mm rod and a 3.8 mm axial width. Test conditions included 70 and 90 seal durometers, a sliding velocity of 0.8 m/sec, and a zero pressure gradient across the seal. An acrylic cylinder and a typical synthetic base automotive lubricant were used. The experimental results showed that the effect of seal hardness on the oil film thickness is considerable. A comparison between analytical and experimental oil film profiles for an elastomeric seal during relatively high speed reciprocating motion showed an overall qualitative agreement.

Krauter, A. I.↗

A current generation by Compton scattering in a relativistic plasma with velocity shear and temperature gradient

Current generation by Thomson scattering in a non-relativistic plasma with the velocity shear and the temperature gradient (Hinata and Daneshvar, 1983) is extended to a relativistic plasma by replacing Thomson cross section by the Klein-Nishina formula. Because of the energy dependence of the cross-section, a numerical rather than analytic result is presented. The present calculation may be applied to a supernova implosion where the temperature may reach several MeV and a strong differential rotation is expected. It may also find applications in the early universe, and laser-pellet interaction.

Hinata, S.↗

Plasma electron collection through biased slits in a dielectric

A large number of experimental and analytical efforts have been directed toward understanding the plasma sheath growth and discharge phenomena which lead to high-voltage solar array/space plasma interactions. An important question which has not been addressed is how the voltage gradient in the plasma sheath near the surface of such an array may affect these interactions. The purpose of the experimental study described in this paper is to examine the merging of the sheaths around biased slits in a dielectric and how this affects the collection of electrons through these slits. The data, which are obtained by emissive probes and direct measurement of the current collected through the slits, indicate that when the sheaths merge the current collection by the slits is significantly altered with the most positive slit collecting more electrons than it otherwise would. Therefore, the effect of a voltage gradient in the sheath around a solar array should be considered when evaluating solar array performance.

Carruth, M. R., Jr.↗

Computer-aided design analysis of 57-mm, angular-contact, cryogenic turbopump bearings

The Space Shuttle main engine high-pressure oxygen turbopumps have not experienced the sevice life required of them. This insufficiency has been due in part to the shortened life of the bearings. To improve the life of the existing turbopump bearings, an effort is under way to investigate bearing modifications that could be retrofitted into the present bearing cavity. Several bearing parameters were optimized using the computer program SHABERTH, which performs a thermomechanical simulation of a load support system. The computer analysis showed that improved bearing performance is feasible if low friction coefficients can be attained. Bearing geometries were optimized considering heat generation, equilibrium temperatures, and relative life. Thermal gradients through the bearings were found to be lower with liquid lubrication than with solid film lubrication, and a liquid oxygen coolant flowrate of approximately 4.0 kg/s was found to be optimal. This paper describes the analytical modeling used to determine these feasible modifications to improve bearing performance.

Armstrong, Elizabeth S.↗

Nonlinear dynamic analysis of quasi-symmetric anisotropic structures

An efficient computational method for the nonlinear dynamic analysis of quasi-symmetric anisotropic structures is proposed. The application of mixed models simplifies the analytical development and improves the accuracy of the response predictions, and operator splitting allows the reduction of the analysis model of the quasi-symmetric structure to that of the corresponding symmetric structure. The preconditoned conjugate gradient provides a stable and effective technique for generating the unsymmetric response of the structure as the sum of a symmetrized response plus correction modes. The effectiveness of the strategy is demonstrated with the example of a laminated anisotropic shallow shell of quadrilateral planform subjected to uniform normal loading.

Noor, Ahmed K.↗

Analytical estimates of radial segregation in Bridgman growth from low-level steady and periodic accelerations

Approximate analytical solutions are developed that are mathematically tractable but provide the elements of CFD models for estimating convective flows subjected to low-level accelerations. The estimates are applicable to dilute systems that have horizontal temperature gradients in the vertical Bridgman configuration. The Navier-Stokes momentum equation is solved in 1D and subsequently in 2D by the first-order perturbation method. The analysis used is the case of plane-front solidification of a dilute system in the vertical thermally stable configuration. The model is found to predict the degree of radial segregation within a factor of two for a range of material and processing parameters. The method provides detailed information on the effects of processing on solute distribution in the grown crystal which are of interest in the orbital experiments emphasizing the control of radial segregation.

Naumann, Robert J.↗

A new approach to mixed H2/H infinity controller synthesis using gradient-based parameter optimization methods

In the past few years, the mixed H(sub 2)/H-infinity control problem has been the object of much research interest since it allows the incorporation of robust stability into the LQG framework. The general mixed H(sub 2)/H-infinity design problem has yet to be solved analytically. Numerous schemes have considered upper bounds for the H(sub 2)-performance criterion and/or imposed restrictive constraints on the class of systems under investigation. Furthermore, many modern control applications rely on dynamic models obtained from finite-element analysis and thus involve high-order plant models. Hence the capability to design low-order (fixed-order) controllers is of great importance. In this research a new design method was developed that optimizes the exact H(sub 2)-norm of a certain subsystem subject to robust stability in terms of H-infinity constraints and a minimal number of system assumptions. The derived algorithm is based on a differentiable scalar time-domain penalty function to represent the H-infinity constraints in the overall optimization. The scheme is capable of handling multiple plant conditions and hence multiple performance criteria and H-infinity constraints and incorporates additional constraints such as fixed-order and/or fixed structure controllers. The defined penalty function is applicable to any constraint that is expressible in form of a real symmetric matrix-inequity.

Ly, Uy-Loi↗

Personalized Tucker Decomposition: Modeling Commonality and Peculiarity on Tensor Data

In this paper, we propose a personalized Tucker decomposition (perTucker) to address the limitations of traditional tensor decomposition methods in capturing heterogeneity across different datasets. perTucker decomposes tensor data into shared global components and personalized local components. We introduce an order orthogonality assumption and develop a proximal gradient regularized block coordinate descent algorithm guaranteed to converge to a stationary point. The unique and common representations learned by perTucker reveal intrinsic statistical patterns in data and provide valuable information for a wide range of downstream analytics, including anomaly detection, source classification, and clustering. We demonstrate perTucker’s effectiveness through a simulation study and two case studies on solar flare detection and tonnage signal classification.

14 SOLAR ENERGY↗

Computation of Direct Sensitivities of Spatial Multibody Systems With Joint Friction

Abstract Friction exists in most mechanical systems and may have a major influence on the dynamic performance of the system. The incorporation of friction in dynamic systems has been a subject of active research for several years owing to its high nonlinearity and its dependence on several parameters. Consequently, optimization of dynamic systems with friction becomes a challenging task. Gradient-based optimization of dynamical systems is a prominent technique for optimal design and requires the computation of model sensitivities with respect to the design parameters. The novel contribution of this paper is the derivation of the analytical methodology for the computation of direct sensitivities for smooth multibody systems with joint friction using the Lagrangian index-1 formulation. System dynamics have been computed using two different friction models; the Brown and McPhee, and the Gonthier et al. model. The methodology proposed to obtain model sensitivities has also been validated using the complex finite difference method. A case study has been conducted on a spatial multibody system to observe the effect of friction on the dynamics and model sensitivities, compare sensitivities with respect to different parameters and demonstrate the numerical and validation aspects. Since design parameters can have very different magnitudes and units, the sensitivities have been scaled with the parameters for comparison. Finally, a discussion has been presented on the interpretation of the case study results. Due to the incorporation of joint friction, ‘jumps’ or discontinuities are observed in the model sensitivities akin to those observed for hybrid dynamical systems.

Engineering↗

Application of Dorodnitsyn's Integral Method to Nonequilibrium Flows over Pointed Bodies

Dorodnitsyn's integral method is used to obtain an approximate solution to the supersonic nonequilibrium flow over pointed bodies with attached shock waves. The partial differential equations governing the flow are converted to an approximate set of ordinary equations, which are solved by numerical integration starting at the body tip. Detailed analytical and numerical results for the first approximation are presented, considering the vibrational relaxation of a diatomic gas over a wedge or cone. It is shown that the first approximation yields: (1) The exact flow-variable gradients at the wedge tip (2) Expressions for the flow-variable gradients at the cone tip which are in agreement with extrapolations of characteristics calculations (3) A good approximate algebraic solution for frozen or equilibrium conical flow (4) An approximate expression for the nonequilibrium-flow stream function which affords a means of obtaining variations across the shock layer of the temperature and vibrational energy. Numerical results for both the wedge and cone compare favorably with identical cases computed by the method of characteristics.

NONEQUILIBRIUM FLOW↗

Orbiting chains and rings

Two very different, highly flexible, space structures proposed during the last decade which involve analytical solution of certain partial differential equations are discussed. The first structure is an array, or hanging chain of aluminum beads which would serve as a convenient communicator if, under the influence of the Earth's gravity gradient, it assumes a local vertical orientation as it circles the Earth. Two passive schemes have been proposed for damping the rigid modes: twist the wire at the ends to provide non-zero moment of inertia about the vertical, thereby inducing relative motion of the two tips during rigid pitch or roll, and thus exercising a damper; and introduce weak lossy springs between the end sections and the main section; thereby providing linear coupling between the springs and all the in-plane (pitch) modes. The second structure is a complete ring of satellites cabled together at synchronous altitude. By a slight increase in altitude, the cable is in tension. This configuration is, however, unstable, and an active feedback control scheme is required to stabilize it. A possible scheme involves local cable length adjustment based on measurement of local altitude and shape variations and their rates.

Breakwell, J. V.↗

Analytical models of helical wind-type astrophysical flows

Three classes of analytic solutions of the basic hydrodynamic equations assumed to govern plasma flow in the atmosphere of a rotating astrophysical object are presented. Attention is focused on the balance of the inertial and gravitational forces with suitable pressure gradients in the presence of rotation. The solutions are written in terms of dimensionless physical parameters in order to facilitate their direct application to specific astrophysical flows, such as winds from massive premain sequence objects and T Tauri stars, and bipolar flows in young stellar objects.

Tsinganos, K.↗

Boost-phase discrimination research activities

Theoretical research in two areas was performed. The aerothermodynamics research focused on the hard-body and rocket plume flows. Analytical real gas models to describe finite rate chemistry were developed and incorporated into the three-dimensional flow codes. New numerical algorithms capable of treating multi-species reacting gas equations and treating flows with large gradients were also developed. The computational chemistry research focused on the determination of spectral radiative intensity factors, transport properties and reaction rates. Ab initio solutions to the Schrodinger equation provided potential energy curves transition moments (radiative probabilities and strengths) and potential energy surfaces. These surfaces were then coupled with classical particle reactive trajectories to compute reaction cross-sections and rates.

Cooper, David M.↗

Using adaptive grid in modeling rocket nozzle flow

The mechanical behavior of a rocket motor internal flow field results in a system of nonlinear partial differential equations which cannot be solved analytically. However, this system of equations called the Navier-Stokes equations can be solved numerically. The accuracy and the convergence of the solution of the system of equations will depend largely on how precisely the sharp gradients in the domain of interest can be resolved. With the advances in computer technology, more sophisticated algorithms are available to improve the accuracy and convergence of the solutions. An adaptive grid generation is one of the schemes which can be incorporated into the algorithm to enhance the capability of numerical modeling. It is equivalent to putting intelligence into the algorithm to optimize the use of computer memory. With this scheme, the finite difference domain of the flow field called the grid does neither have to be very fine nor strategically placed at the location of sharp gradients. The grid is self adapting as the solution evolves. This scheme significantly improves the methodology of solving flow problems in rocket nozzles by taking the refinement part of grid generation out of the hands of computational fluid dynamics (CFD) specialists and place it into the computer algorithm itself.

Chow, Alan S.↗