Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “differential equations”

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 1,009 records · Page 56

From clutter to clarity: Emergent neural operators via questionnaire metrics

Real-world datasets in chemical engineering and bioengineering processes—such as those from catalytic reactors, multiphase flows, polymerization reactors, bioreactors, and clinical trials—can often be unlabeled or disorganized, rendering the training of existing supervised learning models ineffective at learning the underlying dynamics. To salvage these datasets for decision-making, we first seek to obtain clarity from the cluttered data. Here, we present a framework for developing “structural” generative models, discovering emergent equations, and constructing efficient emulators from scrambled datasets by integrating unsupervised organizational learning techniques (Questionnaires) with advanced deep learning architectures (Deep Hidden Physics Models and Deep Operator Networks). Our approach is demonstrated on two illustrative model systems: (a) a 1D advection–diffusion partial differential equation representing a winding underground pipe and (b) an ensemble of Stuart–Landau oscillators, an agent-based system of coupled ordinary differential equations. In both cases, we successfully reconstruct meaningful spatial, temporal, and parameter embeddings from scrambled data, enabling good predictions of system dynamics. As a result, we highlight the framework’s potential for broader applications, enabling data-driven system identification in fields with inherently disorganized or hidden parameter spaces.

42 ENGINEERING↗

On Properties of Adjoint Systems for Evolutionary PDEs

We investigate the geometric structure of adjoint systems associated with evolutionary partial differential equations at the fully continuous, semi-discrete, and fully discrete levels and the relations between these levels. We show that the adjoint system associated with an evolutionary partial differential equation has an infinite-dimensional Hamiltonian structure, which is useful for connecting the fully continuous, semi-discrete, and fully discrete levels. We subsequently address the question of discretize-then-optimize versus optimize-then-discrete for both semi-discretization and time integration, by characterizing the commutativity of discretize-then-optimize methods versus optimize-then-discretize methods uniquely in terms of an adjoint-variational quadratic conservation law. For Galerkin semi-discretizations and one-step time integration methods in particular, we explicitly construct these commuting methods by using structure-preserving discretization techniques.

97 MATHEMATICS AND COMPUTING↗

Finite element solution theory for three-dimensional boundary flows

A finite element algorithm is derived for the numerical solution of a three-dimensional flow field described by a system of initial-valued, elliptic boundary value partial differential equations. The familiar three-dimensional boundary layer equations belong to this description when diffusional processes in only one coordinate direction are important. The finite element algorithm transforms the original description into large order systems of ordinary differential equations written for the dependent variables discretized at node points of an arbitrarily irregular computational lattice. The generalized elliptic boundary conditions is piecewise valid for each dependent variable on boundaries that need not explicitly coincide with coordinate surfaces. Solutions for sample problems in laminar and turbulent boundary flows illustrate favorable solution accuracy, convergence, and versatility.

Baker, A. J.↗

Theoretical predictions of latitude dependencies in the solar wind

Results are presented which were obtained with the Winge-Coleman model for theoretical predictions of latitudinal dependencies in the solar wind. A first-order expansion is described which allows analysis of first-order latitudinal variations in the coronal boundary conditions and results in a second-order partial differential equation for the perturbation stream function. Latitudinal dependencies are analytically separated out in the form of Legendre polynomials and their derivative, and are reduced to the solution of radial differential equations. This analysis is shown to supply an estimate of how large the coronal variation in latitude must be to produce an 11 km/sec/deg gradient in the radial velocity of the solar wind, assuming steady-state processes.

Winge, C. R., Jr.↗

An adaptive grid algorithm for one-dimensional nonlinear equations

Richards' equation, which models the flow of liquid through unsaturated porous media, is highly nonlinear and difficult to solve. Step gradients in the field variables require the use of fine grids and small time step sizes. The numerical instabilities caused by the nonlinearities often require the use of iterative methods such as Picard or Newton interation. These difficulties result in large CPU requirements in solving Richards equation. With this in mind, adaptive and multigrid methods are investigated for use with nonlinear equations such as Richards' equation. Attention is focused on one-dimensional transient problems. To investigate the use of multigrid and adaptive grid methods, a series of problems are studied. First, a multigrid program is developed and used to solve an ordinary differential equation, demonstrating the efficiency with which low and high frequency errors are smoothed out. The multigrid algorithm and an adaptive grid algorithm is used to solve one-dimensional transient partial differential equations, such as the diffusive and convective-diffusion equations. The performance of these programs are compared to that of the Gauss-Seidel and tridiagonal methods. The adaptive and multigrid schemes outperformed the Gauss-Seidel algorithm, but were not as fast as the tridiagonal method. The adaptive grid scheme solved the problems slightly faster than the multigrid method. To solve nonlinear problems, Picard iterations are introduced into the adaptive grid and tridiagonal methods. Burgers' equation is used as a test problem for the two algorithms. Both methods obtain solutions of comparable accuracy for similar time increments. For the Burgers' equation, the adaptive grid method finds the solution approximately three times faster than the tridiagonal method. Finally, both schemes are used to solve the water content formulation of the Richards' equation. For this problem, the adaptive grid method obtains a more accurate solution in fewer work units and less computation time than required by the tridiagonal method. The performance of the adaptive grid method tends to degrade as the solution process proceeds in time, but still remains faster than the tridiagonal scheme.

Gutierrez, William E.↗

Perturbed Equations of Motion for Formation Flight Near the Sun-Earth L2 Point

NASA is planning missions to the vicinity of the Sun-Earth L(sub 2) point, some involving a distributed system of telescope spacecraft, configured in a plane about a hub. Several sets of differential equations are written for the formation flight of such telescopes relative to the hub, with varying levels of fidelity. Effects are cast as additive perturbations to the circular restricted three-body problem, expanded in terms of the system distanced, to an accuracy of 10-20 m. These include Earth's orbital eccentricity, lunar motion, solar radiation pressure, and small thrusting forces. Simulations validating the expanded differential equations are presented.

Luquette, Richard↗

Simulation of cryogenic turbopump annular seals

San Andres employed the NBS software package MIPROPS to account for density's dependence on pressure in the simulation of liquid annular seals. His example on a LH2 seal showed a significant change in the mass coefficient compared to a constant density model. San Andres, Yang, and Childs extended this analysis by including the pressure and temperature dependence of density, specific heat, viscosity, volumetric expansion, and thermal conductivity in a coupled solution of the energy, momentum, and continuity equations. Their example showed very significant changes in stiffness and inertia for a high speed (38,000 rpm), large L/D ratio (0.5) LOX seal, as compared to their constant temperature results. The current research rederived the San Andres-Yang-Childs (SYC) analysis and extended it to include not only the Moody friction model of SYC but also the Hir's friction model. The derivation begins with obtaining the local differential equations of continuity, momentum, and energy conservation in the seal. These equations are averaged across the film thickness to obtain the resulting 'bulk flow' differential equations. Shear stress and convective heat loss through the stator (seal) and rotor are related to the Moody and Hir's friction factor model. The Holman analogy is employed to relate heat conduction in or out of the fluid film's boundary layer to the friction induced shear stress.

Palazzolo, Alan B.↗

Viscous-fluid reaction to a torsionally oscillating spheroid Linearized steady-state weak-separation solution

For both the linearized prolate and oblate problems, perturbation series for all orders of eccentricity are obtained. The 'weak separation' method of solving partial-differential-equation boundary-value problems assumes that each linearly independent solution is a sum of products of single-variable functions. The method is applied to a hierarchy of partial differential equations. The weak-separation solutions are constructed by generalizing only the Bessel functions of the homogeneous solution.

Tompkins, D. R., Jr.↗

Condensation modes in magnetized cylindrical plasmas

In order to understand the formation and stability of solar prominences, attention is given to the condensation modes of the thermal instability in a cylindriical plasma whose magnetic field has both potential (longitudinal) and nonpotential (poloidal) components forming helically twisted field lines. It is noted that the field line twist has a significant effect on the stability of condensation modes which are unstable if field lines are straight, but become stable as the twist increases. The problem is treated in a fully self-consistent way, to derive a second-order ordinary differential equation for radiative MHD stability. The approximation made by neglecting inertial terms is valid for coronal loop conditions if there is no singular surface in the plasma, and a simplified differential equation is derived. Attention is given to stability for the m = zero, 1, and 2 modes numerically, for various loop parameters.

An, C.-H.↗

Stochastic time-optimal control problems

Two types of stochastic time-optimal controls in a one-dimensional setting are considered. Multidimensional problems, in the case of complete state information available and the system modeled by stochastic differential equations, are studied under the formulation of minimizing the expected transient-response time. The necessary condition of optimality is the satisfaction for the value function of a parabolic partial differential equation with boundary conditions. The sufficient condition of optimality is also provided, based on Dynkin's formula. Finally, three examples are given.

Zhang, W.↗

[Development of New Mathematical Methodology in Air Traffic Control for the Analysis of Hybrid Systems]

The aim of this research is to develop new mathematical methodology for the analysis of hybrid systems of the type involved in Air Traffic Control (ATC) problems. Two directions of investigation were initiated. The first used the methodology of nonlinear generalized functions, whose mathematical foundations were initiated by Colombeau and developed further by Oberguggenberger; it has been extended to apply to ordinary differential. Systems of the type encountered in control in joint work with the PI and M. Oberguggenberger. This involved a 'mixture' of 'continuous' and 'discrete' methodology. ATC clearly involves mixtures of two sorts of mathematical problems: (1) The 'continuous' dynamics of a standard control type described by ordinary differential equations (ODE) of the form: {dx/dt = f(x, u)} and (2) the discrete lattice dynamics involved of cellular automata. Most of the CA literature involves a discretization of a partial differential equation system of the type encountered in physics problems (e.g. fluid and gas problems). Both of these directions requires much thinking and new development of mathematical fundamentals before they may be utilized in the ATC work. Rather than consider CA as 'discretization' of PDE systems, I believe that the ATC applications will require a completely different and new mathematical methodology, a sort of discrete analogue of jet bundles and/or the sheaf-theoretic techniques to topologists. Here too, I have begun work on virtually 'virgin' mathematical ground (at least from an 'applied' point of view) which will require considerable preliminary work.

Hermann, Robert↗

A finite difference method for the solution of the transonic flow around harmonically oscillating wings

A finite difference method for the solution of the transonic flow about a harmonically oscillating wing is presented. The partial differential equation for the unsteady transonic flow was linearized by dividing the flow into separate steady and unsteady perturbation velocity potentials and by assuming small amplitudes of harmonic oscillation. The resulting linear differential equation is of mixed type, being elliptic or hyperbolic whereever the steady flow equation is elliptic or hyperbolic. Central differences were used for all derivatives except at supersonic points where backward differencing was used for the streamwise direction. Detailed formulas and procedures are described in sufficient detail for programming on high speed computers. To test the method, the problem of the oscillating flap on a NACA 64A006 airfoil was programmed. The numerical procedure was found to be stable and convergent even in regions of local supersonic flow with shocks.

Ehlers, E. F.↗

Grid spacing control with variation diminishing splines

Methods used to specify and control two and three dimensional grids on which numerical solutions of partial differential equations may be obtained were studied. While initially focusing on grid generation, the research will evolve into a consideration of the interaction of grid generation with the solution of a partial differential equation. The multisurface method of grid generation was used to continuously patch a grid onto an existing grid. In the resulting grid the elements of the Jacobian matrix must be continuous across the boundary between the original grid and the patched grid. Programs were written which accept as input the coordinates of the original grid and the desired new boundary, and then use the three-surface or four-surface version of the multisurface method to extend the original grid out to a new boundary.

Smith, P. W.↗

Computational problems on composite grids

Most currently used algorithms for the numerical solution of the partial differential equations encountered in fluid flow problems can be implemented on composite grid systems. Finite volume formulations are easier to derive on composite grids, and may in principle be derived for partial differential equations of all types. Except when using Alternating Difference Implicit-type schemes, the overlapping of grids is an alternative to the more common grid construction procedure where grid lines continue smoothly from one subregion to the next. In the solution of model problems, the correct choice of an interpolation formula has been found able to reduce errors by a factor of two.

Mastin, C. W.↗

Anticipation of the landing shock phenomenon in flight simulation

An aircraft landing may be described as a controlled crash because a runway surface is intercepted. In a simulation model the transition from aerodynamic flight to weight on wheels involves a single computational cycle during which stiff differential equations are activated; with a significant probability these initial conditions are unrealistic. This occurs because of the finite cycle time, during which large restorative forces will accompany unrealistic initial oleo compressions. This problem was recognized a few years ago at Ames Research Center during simulation studies of a supersonic transport. The mathematical model of this vehicle severely taxed computational resources, and required a large cycle time. The ground strike problem was solved by a described technique called anticipation equations. This extensively used technique has not been previously reported. The technique of anticipating a significant event is a useful tool in the general field of discrete flight simulation. For the differential equations representing a landing gear model stiffness, rate of interception and cycle time may combine to produce an unrealistic simulation of the continuum.

Mcfarland, Richard E.↗

Quadrature imposition of compatibility conditions in Chebyshev methods

Often, in solving an elliptic equation with Neumann boundary conditions, a compatibility condition has to be imposed for well-posedness. This condition involves integrals of the forcing function. When pseudospectral Chebyshev methods are used to discretize the partial differential equation, these integrals have to be approximated by an appropriate quadrature formula. The Gauss-Chebyshev (or any variant of it, like the Gauss-Lobatto) formula can not be used here since the integrals under consideration do not include the weight function. A natural candidate to be used in approximating the integrals is the Clenshaw-Curtis formula, however it is shown that this is the wrong choice and it may lead to divergence if time dependent methods are used to march the solution to steady state. The correct quadrature formula is developed for these problems. This formula takes into account the degree of the polynomials involved. It is shown that this formula leads to a well conditioned Chebyshev approximation to the differential equations and that the compatibility condition is automatically satisfied.

Gottlieb, D.↗

Half-life determination of heavy ions in a storage ring considering feeding and depleting background processes

Heavy-ion storage rings have relatively large momentum acceptance which allows for multiple ion species to circulate at the same time. This needs to be considered in radioactive decay measurements of highly charged ions, where atomic charge exchange reactions can significantly alter the intensities of parent and daughter ions. In this study, we investigate this effect using the decay curves of ion numbers in the recent 205 Tl 81+ bound-state beta decay experiment conducted using the Experimental Storage Ring at GSI Darmstadt. To understand the intricate dynamics of ion numbers, we present a set of differential equations that account for various atomic and nuclear reaction processes—bound state beta decay, atomic electron recombination and capture, and electron ionization. By incorporating appropriate boundary conditions, we develop a set of differential equations that accurately simulate the decay curves of various simultaneously stored ions in the storage ring: 205 Tl 81+ , 205 Pb 81+ , 205 Pb 82+ , 200 Hg 79+ , and 200 Hg 80+ . Through a quantitative comparison between simulations and experimental data, we provide insights into the detailed reaction mechanisms governing stored heavy ions within the storage ring. Our approach effectively models charge-changing processes, reduces the complexity of the experimental setup, and provides a simpler method for measuring the decay half-lives of highly charged ions in storage rings.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

On the stability of numerical integration routines.

Numerical integration methods for the solution of initial value problems for ordinary vector differential equations may be modelled as discrete time feedback systems. The stability criteria discovered in modern control theory are applied to these systems and criteria involving the routine, the step size and the differential equation are derived. Linear multistep, Runge-Kutta, and predictor-corrector methods are all investigated.

Glover, K.↗