SEARCH · Engineering Papers

Results for “numerical integration”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 415 records · Page 23

All sources NASA NTRS DOE OSTI

NASA NTRS1974

A procedure for verifying the structural integrity of an existing pressurized wind tunnel

The present work describes the approach employed in verifying the structural integrity at maximum operating pressure of a low turbulence pressurized wind tunnel built in 1939-1940. The approach utilized existing ASME Code guidelines and fracture mechanics techniques in combination with the latest nondestructive examination methods to validate the desired life of the facility while minimizing reconstruction or major modifications. Detailed stress analyses using a finite element computer program and a numerical integration program were performed on critical areas of the tunnel. Critical vessel areas were evaluated for remaining fatigue life. Fracture mechanics analyses of the shell were performed at various regions, and critical flaw sizes at each location were determined. Analysis of the test data permitted the remaining useful life of the tunnel to be estimated at 8 to 10 years.

Taylor, J. T.↗

DOE OSTI2021

Numerical characterization of quasi-steady thermal load for thin crystal at cryogenic temperature with nondiffusive heat transfer

With dramatically improved brightness and repetition rate, the thermal load for crystal optics in x-ray free-electron laser applications has also significantly increased. To mitigate the thermal load, one effective method is cryogenic cooling. However, the emerging nondiffusive heat transfer phenomenon at cryogenic temperature may cause design failure if overlooked. To evaluate the optical performance of thin crystal optics under thermal load at cryogenic temperature, an integrated numerical tool is presented and applied to characterize the thermal load on thin crystals with nondiffusive phenomena accounted. Significant thermally induced distortion of the rocking curve is observed from numerical simulation, leading to potential seed power reduction in hard x-ray self-seeding application. Cryogenic cooling is confirmed necessary by simulation to handle the thermal load at high repetition rate operation.

42 ENGINEERING↗

NASA NTRS1976

Algorithms for computing the geopotential using a simple density layer

Several algorithms have been developed for computing the potential and attraction of a simple density layer. These are numerical cubature, Taylor series, and a mixed analytic and numerical integration using a singularity-matching technique. A computer program has been written to combine these techniques for computing the disturbing acceleration on an artificial earth satellite. A total of 1640 equal-area, constant surface density blocks on an oblate spheroid are used. The singularity-matching algorithm is used in the subsatellite region, Taylor series in the surrounding zone, and numerical cubature on the rest of the earth.

Morrison, F.↗

NASA NTRS1982

A study of viscous flow in stator and rotor passages

A numerical technique to predict the viscous behavior in the three-dimensional domain of rotating and stationary blade passages is presented. The analysis is based upon the numerical integration of the incompressible Navier-Stokes equations. An approach is taken to combine the three momentum equations and the continuity equation into two sets of vorticity transport-stream function equations. These equations are expressed along families of arbitrarily defined orthogonally intersecting surfaces placed within the blade passage. Unlike the stream surfaces, the arbitrary surfaces allow fluid to flow across. The numerical code employs a nonorthogonal body fitted coordinate system on each of the intersecting surfaces. The solution process involves the iterative combination of the solutions on the two families of intersecting surfaces for a complete solution. The results of flow in thin blade stationary and rotating polar ducts is presented.

Sheoran, Y.↗

NASA NTRS1988

Dynamic analysis of nonlinear rotor-housing systems

Nonlinear analysis methods are developed which will enable the reliable prediction of the dynamic behavior of the space shuttle main engine (SSME) turbopumps in the presence of bearing clearances and other local nonlinearities. A computationally efficient convolution method, based on discretized Duhamel and transition matrix integral formulations, is developed for the transient analysis. In the formulation, the coupling forces due to the nonlinearities are treated as external forces acting on the coupled subsystems. Iteration is utilized to determine their magnitudes at each time increment. The method is applied to a nonlinear generic model of the high pressure oxygen turbopump (HPOTP). As compared to the fourth order Runge-Kutta numerical integration methods, the convolution approach proved to be more accurate and more highly efficient. For determining the nonlinear, steady-state periodic responses, an incremental harmonic balance method was also developed. The method was successfully used to determine dominantly harmonic and subharmonic responses fo the HPOTP generic model with bearing clearances. A reduction method similar to the impedance formulation utilized with linear systems is used to reduce the housing-rotor models to their coordinates at the bearing clearances. Recommendations are included for further development of the method, for extending the analysis to aperiodic and chaotic regimes and for conducting critical parameteric studies of the nonlinear response of the current SSME turbopumps.

Noah, Sherif T.↗

NASA NTRS2010

An Improved Neutron Transport Algorithm for HZETRN

Long term human presence in space requires the inclusion of radiation constraints in mission planning and the design of shielding materials, structures, and vehicles. In this paper, the numerical error associated with energy discretization in HZETRN is addressed. An inadequate numerical integration scheme in the transport algorithm is shown to produce large errors in the low energy portion of the neutron and light ion fluence spectra. It is further shown that the errors result from the narrow energy domain of the neutron elastic cross section spectral distributions, and that an extremely fine energy grid is required to resolve the problem under the current formulation. Two numerical methods are developed to provide adequate resolution in the energy domain and more accurately resolve the neutron elastic interactions. Convergence testing is completed by running the code for various environments and shielding materials with various energy grids to ensure stability of the newly implemented method.

Slaba, Tony C.↗

DOE OSTI2021

Ab Initio Direct Dynamics

The reactivity and dynamics of molecular systems can be explored computationally by classical trajectory calculations. The traditional approach involves fitting a functional form of a potential energy surface (PES) to the energies from a large number of electronic structure calculations and then integrating numerous trajectories on this fitted PES to model the molecular dynamics. The ever-decreasing cost of computing and continuing advances in computational chemistry software have made it possible to use electronic structure calculations directly in molecular dynamics simulations without first having to construct a fitted PES. In this “on-the-fly” approach, every time the energy and its derivatives are needed for the integration of the equations of motion, they are obtained directly from quantum chemical calculations. This approach started to become practical in the mid-1990s as a result of increased availability of inexpensive computer resources and improved computational chemistry software. The application of direct dynamics calculations has grown rapidly over the last 25 years and would require a lengthy review article. The present Account is limited to some of our contributions to methods development and various applications. To improve the efficiency of direct dynamics calculations, we developed a Hessian-based predictor-corrector algorithm for integrating classical trajectories. Hessian updating made this even more efficient. Furthermore, this approach was also used to improve algorithms for following the steepest descent reaction paths. For larger molecular systems, we developed an extended Lagrangian approach in which the electronic structure is propagated along with the molecular structure. Strong field chemistry is a rapidly growing area, and to improve the accuracy of molecular dynamics in intense laser fields, we included the time-varying electric field in a novel predictor-corrector trajectory integration algorithm. Since intense laser fields can excite and ionize molecules, we extended our studies to include electron dynamics. Specifically, we developed code for time-dependent configuration interaction electron dynamics to simulate strong field ionization by intense laser pulses. Our initial application of ab initio direct dynamics in 1994 was to CH 2 O → H 2 + CO; the calculated vibrational distributions in the products were in very good agreement with experiment. In the intervening years, we have used direct dynamics to explore energy partitioning in various dissociation reactions, unimolecular dissociations yielding three fragments, reactions with branching after the transition state, nonstatistical dynamics of chemically activated molecules, dynamics of molecular fragmentation by intense infrared laser pulses, selective activation of specific dissociation channels by aligned intense infrared laser fields, angular dependence of strong field ionization, and simulation of sequential double ionization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

A scalable matrix-free spectral element approach for unsteady PDE constrained optimization using PETSc/TAO

In this work, we provide a new approach for the efficient matrix-free application of the transpose of the Jacobian for the spectral element method for the adjoint-based solution of partial differential equation (PDE) constrained optimization. This results in optimizations of nonlinear PDEs using explicit integrators where the integration of the adjoint problem is not more expensive than the forward simulation. Solving PDE constrained optimization problems entails combining expertise from multiple areas, including simulation, computation of derivatives, and optimization. The Portable, Extensible Toolkit for Scientific computation (PETSc) together with its companion package, the Toolkit for Advanced Optimization (TAO), is an integrated numerical software library that contains an algorithmic/software stack for solving linear systems, nonlinear systems, ordinary differential equations, differential algebraic equations, and large-scale optimization problems and, as such, is an ideal tool for performing PDE-constrained optimization. This paper describes an efficient approach in which the software stack provided by PETSc/TAO can be used for large-scale nonlinear time-dependent problems. Time integration can involve a range of high-order methods, both implicit and explicit. The PDE-constrained optimization algorithm used is gradient-based and seamlessly integrated with the simulation of the physical problem.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS2018

Machine Learning Application to Atmospheric Chemistry Modeling

Atmospheric chemistry is a high-dimensionality, large-data problem and thus may be suited to machine-learning algorithms. We show here the potential of a random forest regression algorithm to replace the gas-phase chemistry solver in the GEOS-Chem chemistry model. In this proof-of-concept study, we used one month of model output to train random forest regression models to predict the concentrations of each long-lived chemical species after integration based upon the physical and chemical conditions before the chemical integration. The choice of prediction type has a strong impact on the skill of the regression model. We find best results from predicting the change in concentration for very long-lived species and the absolute concentration for shorter lived species. The skill of the machine learning algorithm is further improved by using a family approach for NO and NO2 rather than treating them independently.By replacing the numerical integrator with the random forest algorithm and running this model for one month, we find that the model is able to reproduce many of the features of the reference chemistry simulation. Replacing the integration methodology with a machine learning algorithm has the potential to be substantially faster. There are a wide range of applications for such an approach, e.g. to generate boundary conditions, for use in air quality forecasts or chemical data assimilation systems, etc.

Keller, Christoph A.↗

NASA NTRS2015

Efficient Trajectory Propagation for Orbit Determination Problems

Regularized formulations of orbital motion apply a series of techniques to improve the numerical integration of the orbit. Despite their advantages and potential applications little attention has been paid to the propagation of the partial derivatives of the corresponding set of elements or coordinates, required in many orbit-determination scenarios and optimization problems. This paper fills this gap by presenting the general procedure for integrating the state-transition matrix of the system together with the nominal trajectory using regularized formulations and different sets of elements. The main difficulty comes from introducing an independent variable different from time, because the solution needs to be synchronized. The correction of the time delay is treated from a generic perspective not focused on any particular formulation. The synchronization using time-elements is also discussed. Numerical examples include strongly-perturbed orbits in the Pluto system, motivated by the recent flyby of the New Horizons spacecraft, together with a geocentric flyby of the NEAR spacecraft.

numerical methods↗

NASA NTRS1978

On the attitude motion of an orbiting rigid body under the influence of gravity gradient torque

An investigation is made of the rotational motion of a rigid body orbiting the earth, under the influence of the geogravity-gradient torque. The attitude dynamics are formulated as perturbations from a nominal case. The perturbed equations are solved asymptotically, by the multiple scales technique. The independent variable, time, is extended into a space of higher dimension by means of new scales, fast, slow, etc. Integration of the equations is carried out separately in the new variables. The rapid and slow aspects of the attitude dynamics are systematically separated, resulting in a more efficient computer implementation and enhanced physical insight. The theory is applied to predict the attitude dynamics of an asymmetric rigid body satellite. A comparison is made of the maximum errors as the step size increases as predicted by the multiple scales solution and direct numerical integration. An improvement in computational speed of an order of magnitude is demonstrated.

Tao, Y. C.↗

NASA NTRS2011

G-DYN Multibody Dynamics Engine

G-DYN is a multi-body dynamic simulation software engine that automatically assembles and integrates equations of motion for arbitrarily connected multibody dynamic systems. The algorithm behind G-DYN is based on a primal-dual formulation of the dynamics that captures the position and velocity vectors (primal variables) of each body and the interaction forces (dual variables) between bodies, which are particularly useful for control and estimation analysis and synthesis. It also takes full advantage of the spare matrix structure resulting from the system dynamics to numerically integrate the equations of motion efficiently. Furthermore, the dynamic model for each body can easily be replaced without re-deriving the overall equations of motion, and the assembly of the equations of motion is done automatically. G-DYN proved an essential software tool in the simulation of spacecraft systems used for small celestial body surface sampling, specifically in simulating touch-and-go (TAG) maneuvers of a robotic sampling system from a comet and asteroid. It is used extensively in validating mission concepts for small body sample return, such as Comet Odyssey and Galahad New Frontiers proposals.

Acikmese, Behcet↗

NASA NTRSDate not supplied

Numerical Modeling of No Vent Filling of a Cryogenic Tank

This paper presents a multi-node finite volume model of the chilling and filling of a cryogenic tank using the Generalized Fluid System Simulation Program (GFSSP), a general purpose flow network code. Vented Chill and No Vent Fill (VCNVF) tests were conducted at Marshall Space Flight Center where a flight-like tank was filled with liquid nitrogen from a supply tank. In the VCNVF tests, the tank was partially chilled while the vent valve was open. After the partial chilling, the vent valve was closed and the tank was filled without any venting. An integrated numerical model of the test set up was developed. The model included the transfer line from the supply tank, the target tank with the spray nozzle and solid walls, and the discharge line with the vent valve. The tank was discretized into multiple fluid nodes and branches to represent the ullage and liquid nitrogen and multiple solid nodes to represent the tank wall and structure. The heat transfer between solid to fluid was calculated from pool boiling correlations which include film, transition, and nucleate boiling, as well as natural convection during pre- and post-boiling. The model also accounts for the condensation of vapor in the tank when it comes in contact with the liquid spray. The predicted pressure, resident mass, wall and ullage temperature in the tank were compared with the test data.

Alok Majumdar↗

NASA NTRSDate not supplied

Numerical Modeling of No Vent Filling of a Cryogenic Tank

Alok Majumdar↗

DOE OSTI2024

FLOWERS AEP: An Analytical Model for Wind Farm Layout Optimization

Annual energy production (AEP) is commonly used in objective functions for wind farm layout optimization. AEP is proportional to wind farm power production integrated over an annual distribution of free-stream wind conditions. Physics-based estimates of wind farm power production typically rely on low-fidelity engineering wake models that approximate the steady-state wind farm flow field. AEP estimates are then obtained by performing independent simulations for discrete wind conditions and using rectangular quadrature to account for each condition's expected frequency of occurrence. Depending on the number of simulated discrete wind conditions, this numerical integral could be hampered by poor accuracy or high computational costs. The FLOWERS AEP model instead poses an analytical integral of the engineering wake model over the variable wind conditions, yielding a closed-form, analytical function for wind farm AEP. This paper derives the analytical functions for FLOWERS AEP and its derivatives with respect to turbine position, which are useful for gradient-based wind farm layout optimization, in nondimensional form. We then analyze the benefits of the FLOWERS AEP model over conventional reference models, focusing on its low cost, adequate wake loss predictions, and smooth design space. Although the FLOWERS approach is found to predict the exact value of AEP with some error relative to the reference model (within 14% on average), it dramatically reduces computation time by an order of magnitude, produces a qualitatively similar design space at relatively low resolution, and yields comparable optimal layouts. This significant speed improvement is critical in layout optimization applications, where determining an optimal layout in an efficient manner is more important than precise AEP prediction.

17 WIND ENERGY↗

NASA NTRS1996

Kuiper Belt Dust Grains as a Source of Interplanetary Dust Particles

The recent discovery of the so-called Kuiper belt objects has prompted the idea that these objects produce dust grains that may contribute significantly to the interplanetary dust population. In this paper, the orbital evolution of dust grains, of diameters 1 to 9 microns, that originate in the region of the Kuiper belt is studied by means of direct numerical integration. Gravitational forces of the Sun and planets, solar radiation pressure, as well as Poynting-Robertson drag and solar wind drag are included. The interactions between charged dust grains and solar magnetic field are not considered in the model. Because of the effects of drag forces, small dust grains will spiral toward the Sun once they are released from their large parent bodies. This motion leads dust grains to pass by planets as well as encounter numerous mean motion resonances associated with planets. Our results show that about 80% of the Kuiper belt grains are ejected from the Solar System by the giant planets, while the remaining 20% of the grains evolve all the way to the Sun. Surprisingly, the latter dust grains have small orbital eccentricities and inclinations when they cross the orbit of the Earth. This makes them behave more like asteroidal than cometary-type dust particles. This also enhances their chances of being captured by the Earth and makes them a possible source of the collected interplanetary dust particles; in particular, they represent a possible source that brings primitive/organic materials from the outer Solar System to the Earth. When collisions with interstellar dust grains are considered, however, Kuiper belt dust grains around 9 microns appear likely to be collisionally shattered before they can evolve toward the inner part of the Solar System. The collision destruction can be applied to Kuiper belt grains up to about 50 microns. Therefore, Kuiper belt dust grains within this range may not be a significant part of the interplanetary dust complex in the inner Solar System.

Liou, Jer-Chyi↗

NASA NTRS1971

Techniques of orbital decay and long-term ephemeris prediction for satellites in earth orbit

In the special perturbation method, Cowell and variation-of-parameters formulations of the motion equations are implemented and numerically integrated. Variations in the orbital elements due to drag are computed using the 1970 Jacchia atmospheric density model, which includes the effects of semiannual variations, diurnal bulge, solar activity, and geomagnetic activity. In the general perturbation method, two-variable asymptotic series and automated manipulation capabilities are used to obtain analytical solutions to the variation-of-parameters equations. Solutions are obtained considering the effect of oblateness only and the combined effects of oblateness and drag. These solutions are then numerically evaluated by means of a FORTRAN program in which an updating scheme is used to maintain accurate epoch values of the elements. The atmospheric density function is approximated by a Fourier series in true anomaly, and the 1970 Jacchia model is used to periodically update the Fourier coefficients. The accuracy of both methods is demonstrated by comparing computed orbital elements to actual elements over time spans of up to 8 days for the special perturbation method and up to 356 days for the general perturbation method.

Barry, B. F.↗

NASA NTRS2011

Designing the Space Shuttle Propulsion System

The major elements of the Space Shuttle Main Propulsion System include two reusable solid rocket motors integrated into recoverable solid rocket boosters, an expendable external fuel and oxidizer tank, and three reusable Space Shuttle Main Engines. Both the solid rocket motors and space shuttle main engines ignite prior to liftoff, with the solid rocket boosters separating about two minutes into flight. The external tank separates after main engine shutdown and is safely expended in the ocean. The SSME's, integrated into the Space Shuttle Orbiter aft structure, are reused after post landing inspections. Both the solid rocket motors and the space shuttle main engine throttle during early ascent flight to limit aerodynamic loads on the structure. The configuration is called a stage and a half as all the propulsion elements are active during the boost phase, and the SSME's continue operation to achieve orbital velocity approximately eight and a half minutes after liftoff. Design and performance challenges were numerous, beginning with development work in the 1970 s. The solid rocket motors were large, and this technology had never been used for human space flight. The SSME s were both reusable and very high performance staged combustion cycle engines, also unique to the Space Shuttle. The multi body side mount configuration was unique and posed numerous integration and interface challenges across the elements. Operation of the system was complex and time consuming. This paper discusses a number of the system level technical challenges including development and operations.

Owen, James↗