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 937 records · Page 52

A solution procedure for behavior of thick plates on a nonlinear foundation and postbuckling behavior of long plates

Approximate solutions for three nonlinear orthotropic plate problems are presented: (1) a thick plate attached to a pad having nonlinear material properties which, in turn, is attached to a substructure which is then deformed; (2) a long plate loaded in inplane longitudinal compression beyond its buckling load; and (3) a long plate loaded in inplane shear beyond its buckling load. For all three problems, the two dimensional plate equations are reduced to one dimensional equations in the y-direction by using a one dimensional trigonometric approximation in the x-direction. Each problem uses different trigonometric terms. Solutions are obtained using an existing algorithm for simultaneous, first order, nonlinear, ordinary differential equations subject to two point boundary conditions. Ordinary differential equations are derived to determine the variable coefficients of the trigonometric terms.

Stein, M.↗

Design of Flight Vehicle Management Systems

As the operation of large systems becomes ever more dependent on extensive automation, the need for an effective solution to the problem of design and validation of the underlying software becomes more critical. Large systems possess much detailed structure, typically hierarchical, and they are hybrid. Information processing at the top of the hierarchy is by means of formal logic and sentences; on the bottom it is by means of simple scalar differential equations and functions of time; and in the middle it is by an interacting mix of nonlinear multi-axis differential equations and automata, and functions of time and discrete events. The lecture will address the overall problem as it relates to flight vehicle management, describe the middle level, and offer a design approach that is based on Differential Geometry and Discrete Event Dynamic Systems Theory.

Meyer, George↗

Nonlinear Control and Discrete Event Systems

As the operation of large systems becomes ever more dependent on extensive automation, the need for an effective solution to the problem of design and validation of the underlying software becomes more critical. Large systems possesses much detailed structure, typically hierarchical, and they are hybrid. Information processing at the top of the hierarchy is by means of formal logic and sentences; on the bottom it is by means of simple scalar differential equations and functions of time; and in the middle it is by an interacting mix of nonlinear multi-axis differential equations and automata, and functions of time and discrete events. The lecture will address the overall problem as it relates to flight vehicle management, describe the middle level, and offer a design approach that is based on Differential Geometry and Discrete Event Dynamic Systems Theory.

Meyer, George↗

Validation and parameterization of a novel physics-constrained neural dynamics model applied to turbulent fluid flow

We report, in fluid physics, data-driven models to enhance or accelerate time to solution are becoming increasingly popular for many application domains, such as alternatives to turbulence closures, system surrogates, or for new physics discovery. In the context of reduced order models of high-dimensional time-dependent fluid systems, machine learning methods grant the benefit of automated learning from data, but the burden of a model lies on its reduced-order representation of both the fluid state and physical dynamics. In this work, we build a physics-constrained, data-driven reduced order model for Navier–Stokes equations to approximate spatiotemporal fluid dynamics in the canonical case of isotropic turbulence in a triply periodic box. The model design choices mimic numerical and physical constraints by, for example, implicitly enforcing the incompressibility constraint and utilizing continuous neural ordinary differential equations for tracking the evolution of the governing differential equation. We demonstrate this technique on a three-dimensional, moderate Reynolds number turbulent fluid flow. In assessing the statistical quality and characteristics of the machine-learned model through rigorous diagnostic tests, we find that our model is capable of reconstructing the dynamics of the flow over large integral timescales, favoring accuracy at the larger length scales. More significantly, comprehensive diagnostics suggest that physically interpretable model parameters, corresponding to the representations of the fluid state and dynamics, have attributable and quantifiable impact on the quality of the model predictions and computational complexity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Kinematic Flow and the Emergence of Time

Perhaps the most basic question we can ask about cosmological correlations is how their strength changes as we smoothly vary kinematic parameters. The answer is encoded in differential equations that govern this evolution in kinematic space. In this Letter, we introduce a new perspective on these differential equations. We show that, in the simplified setting of conformally coupled scalars in power-law Friedmann-Robertson-Walker spacetimes, the equations for arbitrary tree-level processes can be obtained from a small number of simple combinatorial rules. While this “kinematic flow” is defined purely in terms of boundary data, it reflects the physics of bulk time evolution. The unexpected regularity of the equations suggests the existence of an autonomously defined mathematical structure from which cosmological correlations and the time evolution of the associated spacetime emerge.

79 ASTRONOMY AND ASTROPHYSICS↗

A method for phenomenological analysis of ecological data.

The experimental meaning of the phenomenological differential equations for a competing population is reviewed. It is concluded that it is virtually impossible to construct the differential equations precise enough for studying stability. We consider instead a method of phenomenological analysis which can be applied to a set of population curves. We suggest an ecological index calculated from the population curves, which indicates a group property of the entire system. As a function of time, the index is presumably insensitive to Volterra type fluctuations. A marked increase of the index's value however indicates a marked change of the environmental conditions. It is not easy to deduce the group property from the population curves alone, because a change in population is in general due to the superposition of external disturbances and Volterra fluctuations.

Huang, H.-W.↗

Spline approximations for nonlinear hereditary control systems

A sline-based approximation scheme is discussed for optimal control problems governed by nonlinear nonautonomous delay differential equations. The approximating framework reduces the original control problem to a sequence of optimization problems governed by ordinary differential equations. Convergence proofs, which appeal directly to dissipative-type estimates for the underlying nonlinear operator, are given and numerical findings are summarized.

Daniel, P. L.↗

A new approach to optimization for aerodynamic applications

A new approach to solving optimization problems that involve nonlinear partial differential equations is presented. The approach eliminates the need for an inner-outer iterative procedure, solving the partial differential equation only once, thereby reducing the cost of computation to an extent which would allow its use as a practical tool in optimization problems. The approach is tested on a single design parameter problem through the use of a specially developed scheme. The results indicate comparable convergence properties for the present iterative process and the standard iterative scheme. The presented ideas are also applicable to multidesign parameter problems.

Rizk, M. H.↗

Switch-on shocks

It is argued that parallel propagating Alfven waves steepen to second order in the wave amplitude, rather than first order, as is the case for oblique MHD waves. An analytical solution for the structure of weak switch-on shocks is derived. The Rankine-Hugoniot relations for finite amplitude switch-on shocks are obtained in a compact form. The differential equation for resistive switch-on shocks is solved, and the resistive critical Mach number is obtained. The differential equation for resistive switch-on shocks with finite ion inertia is identical to the one obtained in the pure resistive case. It is argued that efficient electron heat conduction can significantly extend the range of upstream parameters for which switch-on shocks are expected.

Kennel, C. F.↗

Almost periodic solutions to difference equations

The theory of Massera and Schaeffer relating the existence of unique almost periodic solutions of an inhomogeneous linear equation to an exponential dichotomy for the homogeneous equation was completely extended to discretizations by a strongly stable difference scheme. In addition it is shown that the almost periodic sequence solution will converge to the differential equation solution. The preceding theory was applied to a class of exponentially stable partial differential equations to which one can apply the Hille-Yoshida theorem. It is possible to prove the existence of unique almost periodic solutions of the inhomogeneous equation (which can be approximated by almost periodic sequences) which are the solutions to appropriate discretizations. Two methods of discretizations are discussed: the strongly stable scheme and the Lax-Wendroff scheme.

Bayliss, A.↗

Some MACSYMA program for solving difference equations

A set of MACSYMA programs are described for finding closed form solutions to linear recurrence relations in equations having either constant or variable coefficients. In the homogenous case, a polymonial equation is obtained and the solution to the recurrence relation can be written as a linear combination of the roots of the polynomial. Exponential generating functions are used to solve variable coefficient relations. Taking successive derivatives and using the recurrence relation, an ordinary differential equation is obtained. Expanding the solution to the differential equation in a Taylor series, shows that the nth term of the series is the solution to the recurrence relation. For second order recurrences, a check is made for those that can be solved in terms of Bessel functions.

Ivie, J.↗

An analysis of turbulent diffusion flame in axisymmetric jet

The kinetic theory of turbulent flow was employed to study the mixing limited combustion of hydrogen in axisymmetric jets. The integro-differential equations in two spatial and three velocity coordinates describing the combustion were reduced to a set of hyperbolic partial differential equations in the two spatial coordinates by a binodal approximation. The MacCormick's finite difference method was then employed for solution. The flame length was longer than that predicted by the flame-sheet analysis, and was found to be in general agreement with a recent experimental result. Increase of the turbulence energy and scale resulted in an enhancement of the combustion rate and, hence, in a shorter flame length. Details of the numerical method as well as of the physical findings are discussed.

Chung, P. M.↗

Theory of biaxial graded-index optical fiber

The problem of wave propagation in a biaxial graded-index fiber with circular symmetry is considered. The problem is formulated in terms of four first-order differential equations for the tangential components of the electric and magnetic fields. A general solution method for solving systems of differential equations is presented. This solution method is then used to solve the system of equations for a particular example of a biaxial graded-index fiber. Numerical results for the propagation constant in the fiber are also given.

Kawalko, Stephen F.↗

A method for the analysis of biaxial graded-index optical fibers

The problem of wave propagation is a biaxial graded-index fiber with circular symmetry is considered. The problem is formulated in terms of four first-order differential equations for the tangential components of the electric and magnetic fields. A general solution method for solving systems of differential equations is presented. This solution method is then used to solve the system of equations for a particular example of a biaxial graded-index fiber. Numerical results for the propagation constant in the fiber are given.

Kawalko, Stephen F.↗

Inclusion of transverse shear deformation in the exact buckling and vibration analysis of composite plate assemblies

The problem considered is the development of the necessary plate stiffnesses for use in the general purpose program VICONOPT for buckling and vibration of composite plate assemblies. The required stiffnesses include the effects of transverse shear deformation and are for sinusoidal response along the plate length as required in VICONOPT. The method is based on the exact solution of the plate differential equations for a composite laminate having fully populated A, B, and D stiffness matrices which leads to an ordinary differential equation of tenth order.

Anderson, Melvin S.↗

Lagrangian density for collisional plasma

For the purpose of deriving appropriate Lagrangians for plasma equations that include effects of energy loss, the paper examines the inverse problem of the calculus of variations for systems of first- and second-order quasi-linear partial differential equations. This results in convenient forms of the sufficient conditions under which the given differential equations are Euler-Lagrange equations of a Lagrangian. These conditions are then applied to determine the necessary transformation that converts equations, apparently not already in it, into Euler-Lagrange form. The appropriate Lagrangian for a warm collisional plasma is obtained, and the Lagrangian is derived for a resistive transmission line.

Peng, Y.-K. M.↗

Stability of Detached Solidification

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.

Mazuruk, K.↗