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 901 records · Page 50

Similar solutions for the compressible laminar boundary layer with heat transfer and pressure gradient

Stewartson's transformation is applied to the laminar compressible boundary-layer equations and the requirement of similarity is introduced, resulting in a set of ordinary nonlinear differential equations previously quoted by Stewartson, but unsolved. The requirements of the system are Prandtl number of 1.0, linear viscosity-temperature relation across the boundary layer, an isothermal surface, and the particular distributions of free-stream velocity consistent with similar solutions. This system admits axial pressure gradients of arbitrary magnitude, heat flux normal to the surface, and arbitrary Mach numbers. The system of differential equations is transformed to integral system, with the velocity ratio as the independent variable. For this system, solutions are found by digital computation for pressure gradients varying from that causing separation to the infinitely favorable gradient and for wall temperatures from absolute zero to twice the free-stream stagnation temperature. Some solutions for separated flows are also presented.

Cohen, Clarence B↗

PolyODENet: Deriving mass-action rate equations from incomplete transient kinetics data

Kinetics of a reaction network that follows mass-action rate laws can be described with a system of ordinary differential equations (ODEs) with polynomial right-hand side. However, it is challenging to derive such kinetic differential equations from transient kinetic data without knowing the reaction network, especially when the data are incomplete due to experimental limitations. We introduce a program, PolyODENet, toward this goal. Based on the machine-learning method Neural ODE, PolyODENet defines a generative model and predicts concentrations at arbitrary time. As such, it is possible to include unmeasurable intermediate species in the kinetic equations. Importantly, we have implemented various measures to apply physical constraints and chemical knowledge in the training to regularize the solution space. Here, using simple catalytic reaction models, we demonstrate that PolyODENet can predict reaction profiles of unknown species and doing so even reveal hidden parts of reaction mechanisms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Application of the Green's function method for 2- and 3-dimensional steady transonic flows

A Time-Domain Green's function method for the nonlinear time-dependent three-dimensional aerodynamic potential equation is presented. The Green's theorem is being used to transform the partial differential equation into an integro-differential-delay equation. Finite-element and finite-difference methods are employed for the spatial and time discretizations to approximate the integral equation by a system of differential-delay equations. Solution may be obtained by solving for this nonlinear simultaneous system of equations in time. This paper discusses the application of the method to the Transonic Small Disturbance Equation and numerical results for lifting and nonlifting airfoils and wings in steady flows are presented.

Tseng, K.↗

Singular perturbation equations for flexible satellites

The dynamic model of a flexible satellite with n generalized structural position coordinates requires the solution of a set of N coupled nonlinear ordinary or partial differential equations. For single composite body satellites, N is equal to (n + 3). For dual-spin systems, N is equal to (n + 9). These equations usually involve time derivatives up to the second order. For large values of n, linearization of the system has so far been the only practicable way of solution. The present study shows a method of avoiding this linearization by reducing the N equations to three or nine nonlinear, coupled, first-order ordinary differential equations involving only the angular velocities of the composite bodies. The solutions for these angular velocities lead to linear equations in the n generalized structural position coordinates, which can then be solved by known methods.

Huang, T. C.↗

Variational objective analyses for cyclone studies

The basic analysis equations, i.e., the two horizontal momentum equations, the hydrostatic equation, and the integrated continuity equation were derived for the nonlinear vertical coordinate, nondimensionalized, and expressed in finite differences on a staggered grid. Special care was taken to transform the hydrostatic equation and the pressure gradient terms of the horizontal momentum equations to nearly eliminate truncation error over steeply sloping terrain. This formulation also eliminated explicit reference to orographically induced variations in the thermodynamic variables so that the variational adjustments are on the scale of the meteorological perturbations. The analysis equations were subjected to the Euler-Lagrange operations as expressed for finite differences and an additional set of five partial differential equations was derived, bringing to nine the number of equations in Model I. Higher order terms, terms containing observed quantities, and terms containing none of the variables to be adjusted were grouped into forcing functions and the equations were solved for the zero order terms. Zero order variables were eliminated between these equations and there resulted two diagnostic equations which take the form of general linear second order partial differential equations with nonconstant coefficients.

Achtemeier, G. L.↗

The extended auxiliary equation mapping method to determine novel exact solitary wave solutions of the nonlinear fractional PDEs

Abstract In this paper, some new nonlinear fractional partial differential equations (PDEs) have been considered.Three models are including the space-time fractional-order Boussinesq equation, space-time (2 + 1)-dimensional breaking soliton equations, and space-time fractional-order SRLW equation describe the behavior of these equations in the diverse applications. Meanwhile, the fractional derivatives in the sense of β -derivative are defined. Some fractional PDEs will convert to the considered ordinary differential equations by the help of transformation of β -derivative. These equations are analyzed utilizing an integration scheme, namely, the extended auxiliary equation mapping method. The different kinds of traveling wave solutions, solitary, topological, dark soliton, periodic, kink, and rational, fall out as a by-product of this scheme. Finally, the existence of the solutions for the constraint conditions is also shown. The outcome indicates that some fractional PDEs are used as a growing finding in the engineering sciences, mathematical physics, and so forth.

Engineering↗

On the integrability cases of the equation of motion for a satellite in an axially symmetric gravitational field.

Investigation of two cases of integrability of a second-order differential equation describing the projection of an axisymmetric satellite orbit on to a plane perpendicular to the rotation axis. It is demonstrated that for these two cases the integration can be carried out either by quadratures or reduced to a first-order differential equation. Analytical and physical properties are expressed, and it is shown that the equation can be derived from the classical plane eikonal equation of geometric optics.

Ghaffari, A.↗

Quadrature imposition of compatibility conditions in Chebyshev methods

This paper handles a problem encountered in the solution by way of Chebyshev spectral collocation discretization of the Poisson equation with homogeneous Neumann boundary conditions. When pseudospectral Chebyshev methodology is utilized to discretize the partial differential equation these integrals have to be approximated by an appropriate quadrature formula. The correct quadrature formula is developed for these problems, taking 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.↗

Modeling of outgassing and matrix decomposition in carbon-phenolic composites

Work done in the period Jan. - June 1994 is summarized. Two threads of research have been followed. First, the thermodynamics approach was used to model the chemical and mechanical responses of composites exposed to high temperatures. The thermodynamics approach lends itself easily to the usage of variational principles. This thermodynamic-variational approach has been applied to the transpiration cooling problem. The second thread is the development of a better algorithm to solve the governing equations resulting from the modeling. Explicit finite difference method is explored for solving the governing nonlinear, partial differential equations. The method allows detailed material models to be included and solution on massively parallel supercomputers. To demonstrate the feasibility of the explicit scheme in solving nonlinear partial differential equations, a transpiration cooling problem was solved. Some interesting transient behaviors were captured such as stress waves and small spatial oscillations of transient pressure distribution.

Mcmanus, Hugh L.↗

Latent Neural ODE for Integrating Multi-Timescale Measurements in Smart Distribution Grids

Under a smart grid paradigm, there has been an increase in sensor installations to enhance situational awareness. The measurements from these sensors can be leveraged for real-time monitoring, control, and protection. However, these measurements are typically irregularly sampled. These measure-ments may also be intermittent due to communication bandwidth limitations. To tackle this problem, this paper proposes a novel latent neural ordinary differential equations (LODE) approach to aggregate the unevenly sampled multivariate time-series measurements. The proposed approach is flexible in performing both imputations and predictions while being computationally efficient. Simulation results on IEEE 37 bus test systems illustrate the efficiency of the proposed approach.

multi time-scale measurements↗

Leveraging Multitime Hamilton–Jacobi PDEs for Certain Scientific Machine Learning Problems

Hamilton-Jacobi partial differential equations (HJ PDEs) have deep connections with a wide range of fields, including optimal control, differential games, and imaging sciences. By considering the time variable to be a higher dimensional quantity, HJ PDEs can be extended to the multi-time case. In this paper, we establish a novel theoretical connection between specific optimization problems arising in machine learning and the multi-time Hopf formula, which corresponds to a representation of the solution to certain multi-time HJ PDEs. Through this connection, we increase the interpretability of the training process of certain machine learning applications by showing that when we solve these learning problems, we also solve a multi-time HJ PDE and, by extension, its corresponding optimal control problem. As a first exploration of this connection, we develop the relation between the regularized linear regression problem and the Linear Quadratic Regulator (LQR). We then leverage our theoretical connection to adapt standard LQR solvers (namely, those based on the Riccati ordinary differential equations) to design new training approaches for machine learning. Lastly, we provide some numerical examples that demonstrate the versatility and possible computational advantages of our Riccati-based approach in the context of continual learning, post-training calibration, transfer learning, and sparse dynamics identification.

97 MATHEMATICS AND COMPUTING↗

Simulation of a steady-state integrated human thermal system.

The mathematical model of an integrated human thermal system is formulated. The system consists of an external thermal regulation device on the human body. The purpose of the device (a network of cooling tubes held in contact with the surface of the skin) is to maintain the human body in a state of thermoneutrality. The device is controlled by varying the inlet coolant temperature and coolant mass flow rate. The differential equations of the model are approximated by a set of algebraic equations which result from the application of the explicit forward finite difference method to the differential equations. The integrated human thermal system is simulated for a variety of combinations of the inlet coolant temperature, coolant mass flow rate, and metabolic rates. Two specific cases are considered: (1) the external thermal regulation device is placed only on the head and (2) the devices are placed on the head and the torso. The results of the simulation indicate that when the human body is exposed to hot environment, thermoneutrality can be attained by localized cooling if the operating variables of the external regulation device(s) are properly controlled.

Hsu, F. T.↗

Simulation of a steady-state integrated human thermal system.

The mathematical model of an integrated human thermal system is formulated. The system consists of an external thermal regulation device on the human body. The purpose of the device (a network of cooling tubes held in contact with the surface of the skin) is to maintain the human body in a state of thermoneutrality. The device is controlled by varying the inlet coolant temperature and coolant mass flow rate. The differential equations of the model are approximated by a set of algebraic equations which result from the application of the explicit forward finite difference method to the differential equations. The integrated human thermal system is simulated for a variety of combinations of the inlet coolant temperature, coolant mass flow rate, and metabolic rates.

Hsu, F. T.↗

Finite element or Galerkin type semidiscrete schemes

A finite element of Galerkin type semidiscrete method is proposed for numerical solution of a linear hyperbolic partial differential equation. The question of stability is reduced to the stability of a system of ordinary differential equations for which Dahlquist theory applied. Results of separating the part of numerical solution which causes the spurious oscillation near shock-like response of semidiscrete scheme to a step function initial condition are presented. In general all methods produce such oscillatory overshoots on either side of shocks. This overshoot pathology, which displays a behavior similar to Gibb's phenomena of Fourier series, is explained on the basis of dispersion of separated Fourier components which relies on linearized theory to be satisfactory. Expository results represented.

Durgun, K.↗

A new approach for designing self-organizing systems and application to adaptive control

There is tremendous interest in the design of intelligent machines capable of autonomous learning and skillful performance under complex environments. A major task in designing such systems is to make the system plastic and adaptive when presented with new and useful information and stable in response to irrelevant events. A great body of knowledge, based on neuro-physiological concepts, has evolved as a possible solution to this problem. Adaptive resonance theory (ART) is a classical example under this category. The system dynamics of an ART network is described by a set of differential equations with nonlinear functions. An approach for designing self-organizing networks characterized by nonlinear differential equations is proposed.

Ramamoorthy, P. A.↗

Implementing abstract multigrid or multilevel methods

Multigrid methods can be formulated as an algorithm for an abstract problem that is independent of the partial differential equation, domain, and discretization method. In such an abstract setting, problems not arising from partial differential equations can be treated. A general theory exists for linear problems. The general theory was motivated by a series of abstract solvers (Madpack). The latest version was motivated by the theory. Madpack now allows for a wide variety of iterative and direct solvers, preconditioners, and interpolation and projection schemes, including user callback ones. It allows for sparse, dense, and stencil matrices. Mildly nonlinear problems can be handled. Also, there is a fast, multigrid Poisson solver (two and three dimensions). The type of solvers and design decisions (including language, data structures, external library support, and callbacks) are discussed. Based on the author's experiences with two versions of Madpack, a better approach is proposed. This is based on a mixed language formulation (C and FORTRAN + preprocessor). Reasons for not using FORTRAN, C, or C++ (individually) are given. Implementing the proposed strategy is not difficult.

Douglas, Craig C.↗

Numerical Stability In Hyperbolic Boundary-Value Problems

Technical memorandum discusses stability of numerical solutions involving semidiscrete approximations to hyperbolic partial differential equations in initial-and-boundary-value problems. Topic of practical significance because hyperbolic partial differential equations arise in mathematical modeling of waves and blasts. Solutions often needed over restricted regions of unbounded spaces. Outer boundaries artificial, introduced only to limit domains of numerical solutions. Conditions at such artificial boundaries cause numerical instabilities that degrade computed solutions.

Warming, Robert F.↗