Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “EXACT SOLUTIONS”

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 307 records · Page 17

The method of Ritz applied to the equation of Hamilton

Without any reference to the theory of differential equations, the initial value problem of the nonlinear, nonconservative double pendulum system is solved by the application of the method of Ritz to the equation of Hamilton. Also shown is an example of the reduction of the traditional eigenvalue problem of linear, homogeneous, differential equations of motion to the solution of a set of nonhomogeneous algebraic equations. No theory of differential equations is used. Solution of the time-space path of the linear oscillator is demonstrated and compared to the exact solution.

Bailey, C. D.↗

Buckling and vibration of periodic lattice structures

Lattice booms and platforms composed of flexible members or large diameter rings which may be stiffened by cables in order to support membrane-like antennas or reflector surfaces are the main components of some large space structures. The nature of these structures, repetitive geometry with few different members, makes possible relatively simple solutions for buckling and vibration of a certain class of these structures. Each member is represented by a stiffness matrix derived from the exact solution of the beam column equation. This transcendental matrix gives the current member stiffness at any end load or frequency. Using conventional finite element techniques, equilibrium equations can be written involving displacements and rotations of a typical node and its neighbors. The assumptions of a simple trigonometric mode shape is found to satisfy these equations exactly; thus the entire problem is governed by just one 6 x 6 matrix equation involving the amplitude of the displacement and rotation mode shapes. The boundary conditions implied by this solution are simple supported ends for the column type configurations.

Anderson, M. S.↗

A theoretical investigation of the sound radiation fields associated with a Bellmouth inlet

Analytical results are obtained by numerical integration of a cylindrically symmetric integral representation of the external solutions of the Helmholtz equation. The accuracy of this method is checked by comparisons of computed results with 'exact' solutions generated by the point source method. In all cases, the average error for the amplitude and phase of the points calculated in the field is found to be less than ten percent. Theoretical studies which model experiments run for the NASA Langely Bellmouth inlet configuration are presented and comparisons are made with the experimental results. In all cases very good agreement is obtained between the experimental and theoretically calculated values.

Meyer, W. L.↗

Static Load Distribution in Ball Bearings

A numerical procedure for computing the internal loading distribution in statically loaded, single-row, angular-contact ball bearings when subjected to a known combined radial and thrust load is presented. The combined radial and thrust load must be applied in order to avoid tilting between inner and outer rings. The numerical procedure requires the iterative solution of Z + 2 simultaneous nonlinear equations - where Z is the number of the balls - to yield an exact solution for axial and radial deflections, and contact angles. Numerical results for a 218 angular-contact ball bearing have been compared with those from the literature and show significant differences in the magnitudes of the ball loads, contact angles, and the extent of the loading zone.

Ricci, Mario↗

Analytical solutions of the Arrhenius-Semenov problem for constant volume burn

Analytical solutions to the Semenov thermal ignition problem for constant volume burn governed by Arrhenius reaction kinetics are derived. Specifically, an approximate analytical solution technique for the Arrhenius-Semenov differential equation is derived for reaction orders n ϵ R> 0 and exact solutions are also constructed for reaction orders n ϵ N : n ≤ 3. The approximation technique relies on expansion of the respective nondominant terms in the differential equation at the lower and upper bounds of the reaction progress variable in order to create a pair of integrable series. The two integrated series are then connected to create a single continuous analytical solution. Excellent agreement is observed between the analytical approximation and solutions obtained numerically. The presented approximation constitutes a simple and robust strategy for solving the Arrhenius-Semenov problem analytically.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Reduced and simplified chemical kinetics for air dissociation using Computational Singular Perturbation

The Computational Singular Perturbation CSP methods is employed (1) in the modeling of a homogeneous isothermal reacting system and (2) in the numerical simulation of the chemical reactions in a hypersonic flowfield. Reduced and simplified mechanisms are constructed. The solutions obtained on the basis of these approximate mechanisms are shown to be in very good agreement with the exact solution based on the full mechanism. Physically meaningful approximations are derived. It is demonstrated that the deduction of these approximations from CSP is independent of the complexity of the problem and requires no intuition or experience in chemical kinetics.

Goussis, D. A.↗

Model of brittle matrix composite toughening based on discrete fiber reinforcement

An analytical approach for the analysis of the effectiveness of fiber reinforcement in brittle matrix composites is presented. The analytical method allows consideration of discrete fiber distribution and examination of the development of crack growth parameters on microscale. The problem associated with the bridging zone development is addressed here; therefore, the bridging zone is considered to be smaller than the main preexisting crack, and the small scale approach is used. The mechanics of the reinforcement is accurately accounted for in the process zone of a growing crack. Closed form solutions characterizing the initial failure process are presented for linear and nonlinear force - fiber pullout displacement relationships. The implicit exact solution for the extended bridging zone is presented as well.

Rubinstein, Asher A.↗

Integrating Multi-Source Data for Bi-Level Traffic Simulator Calibration: A Literature Review and Highway Case Study

Traffic simulation serves as a powerful tool for pre-evaluating policies and technologies. In this context, simulation-based Dynamic traffic assignment (DTA) models are capable of capturing traffic dynamics. They are well-known as critical tools in controlling and predicting traffic situations. The reliability of simulation results heavily depends on the calibration process. Most studies in the literature formulate and calibrate simulators based on a single source of collected data or multiple data sets with the same spatiotemporal characteristics. However, in practice, traffic data is collected by various tools with usually different spatial and temporal resolutions. This study introduces a novel approach to taking into account diverse input data from a variety of sources. An iterative bi-level solution is proposed. to equally treat traffic flow and speed data. The upper level solves flow calibration with the exact solution method, and the lower level calibrates the speed with the simultaneous perturbation stochastic approximation (SPSA) algorithm. Subsequently, the effectiveness of the proposed model is investigated using data from a six-mile section of Nashville's I-24 highway in Tennessee. The results demonstrate that our proposed model creates an effective feedback loop between the optimizer and the simulator for calibrating flow and speed to reduce the error between simulated and real data.

42 ENGINEERING↗

Modal propagation angles in a cylindrical duct with flow and their relation to sound radiation

The angles of propagation for the wave front making up a duct mode are presented with the Mach number in the duct. Approximate equations are derived to provide simple utilitarian expressions. These expressions are valid only near the outer wall which is the most important region since the bulk of the acoustic intensity is located there. Exact solutions using Hankel functions are given in an appendix. These data corroborate the approximate solution accuracy near the outer wall. The axial propagation angle is used to infer information about the far-field radiation pattern. The resultant axial angle of propagation in the duct is shown to agree exactly with the peak of the principal lobe to far-field radiation obtained from formal radiation calculations when the Mach number is uniform everywhere. The obtained solution is extended to cover the case of different Mach numbers inside and outside the duct for which exact calculations have not been available for engine inlet configurations.

Rice, E. J.↗

Pressure falloff testing to characterize CO 2 plume and dry-out zone during CO 2 injection in saline aquifers

The fluid mobilities and the lateral extent of the CO 2 plume and the dry-out zone are among important unknowns required to effectively manage CO 2 storage in saline aquifers projects. These unknowns can be derived from pressure falloff testing which has been widely used to characterize the subsurface. This paper provides a three-region analytical solution for falloff pressure during the infinite-acting period of CO 2 injection in saline aquifers. The solution is derived using Laplace transformation technique of the governing diffusivity equations and the accompanying initial and boundary conditions. The exact solution is validated against numerical results obtained from a black-oil reservoir simulator for two cases. The solution is approximated into simple forms in real-time domain considering different time periods. Then, a graphical semi-log analysis technique based on the approximate solutions is provided as an interpretation tool of falloff tests where the fluid mobilities and the extents of dry-out zone and CO 2 plume can be inferred. The interpretation technique is applied to synthetic falloff test data to investigate its potential and reliability.

02 PETROLEUM↗

Micromechanical analysis of the failure process in ceramic matrix composites

An analysis of the effectiveness of fiber reinforcement in brittle matrix composites is presented. The analytical method allows consideration of discrete fiber distribution and examination of the development of crack growth parameters on the microscale. The problem associated with bridging zone development is addressed here; therefore, the bridging zone is considered to be smaller than the main preexisting crack, and the small scale approach is used. The mechanics of the reinforcement is accurately accounted for in the process zone of a growing crack. Closed form solutions characterizing the initial failure process are presented for linear and nonlinear force-fiber pullout displacement relationships. The implicit exact solution for the extended bridging zone is presented as well.

Rubinstein, Asher A.↗

Code verification for practically singular equations

We report the method-of-moments implementation of the electric-field integral equation (EFIE) yields many code-verification challenges due to the various sources of numerical error and their possible interactions. Matters are further complicated by singular integrals, which arise from the presence of a Green's function. To address these singular integrals, an approach is presented in wherein both the solution and Green's function are manufactured. Because the arising equations are poorly conditioned, they are reformulated as a set of constraints for an optimization problem that selects the solution closest to the manufactured solution. In this paper, we demonstrate how, for such practically singular systems of equations, computing the truncation error by inserting the exact solution into the discretized equations cannot detect certain orders of coding errors. On the other hand, the discretization error from the optimal solution is a more sensitive metric that can detect orders less than those of the expected convergence rate.

97 MATHEMATICS AND COMPUTING↗

Three-dimensional unsteady lifting surface theory in the subsonic range

The methods of the unsteady lifting surface theory are surveyed. Linearized Euler's equations are simplified by means of a Galileo-Lorentz transformation and a Laplace transformation so that the time and the compressibility of the fluid are limited to two constants. The solutions to this simplified problem are represented as integrals with a differential nucleus; these results in tolerance conditions, for which any exact solution must suffice. It is shown that none of the existing three-dimensional lifting surface theories in subsonic range satisfy these conditions. An oscillating elliptic lifting surface which satisfies the tolerance conditions is calculated through the use of Lame's functions. Numerical examples are calculated for the borderline cases of infinitely stretched elliptic lifting surfaces and of circular lifting surfaces. Out of the harmonic solutions any such temporal changes of the down current are calculated through the use of an inverse Laplace transformation.

Kuessner, H. G.↗

Wave interactions in a three-dimensional attachment line boundary layer

The 3-D boundary layer on a swept wing can support different types of hydrodynamic instability. Attention is focused on the so-called spanwise contamination problem, which occurs when the attachment line boundary layer on the leading edge becomes unstable to Tollmien-Schlichting waves. In order to gain insight into the interactions important in that problem, a simplified basic state is considered. This simplified flow corresponds to the swept attachment line boundary layer on an infinite flat plate. The basic flow here is an exact solution of the Navier-Stokes equations and its stability to 2-D waves propagating along the attachment can be considered exactly at finite Reynolds number. This has been done in the linear and weakly nonlinear regimes. The corresponding problem is studied for oblique waves and their interaction with 2-D waves is investigated. In fact, oblique modes cannot be described exactly at finite Reynolds number so it is necessary to make a high Reynolds number approximation and use triple deck theory. It is shown that there are two types of oblique wave which, if excited, cause the destabilization of the 2-D mode and the breakdown of the disturbed flow at a finite distance from the leading edge. First, a low frequency mode related to the viscous stationary crossflow mode is a possible cause of breakdown. Second, a class of oblique wave with frequency comparable with that of the 2-D mode is another cause of breakdown. It is shown that the relative importance of the modes depends on the distance from the attachment line.

Hall, Philip↗

Numerical simulations of high-speed flows about waveriders with sharp leading edges

A procedure is developed for the numerical simulation of stagnation-free inviscid supersonic and hypersonic flows about waveriders with sharp leading edges. The numerical approach involves the development of a specialized grid generator (named HYGRID), an algebraic solution-adaptive grid scheme, and a modified flow solving method. A comparison of the results obtained for several waverider geometries with exact solutions, other numerical solutions, and experimental results demonstrated the ability of the new procedure to produce stagnation-free Euler solutions about sharp-edged configurations and to describe the physics of the flow in these regions.

Jones, Kevin D.↗

Analytic solution of the box diffusion model for a global ocean

The paper presents an analytic solution for the temperature of the mixed layer for a box diffusion model of a global ocean. It is shown that the exact solution approaches the appropriate limits of a mixed layer ocean model and of a diffusive ocean without a mixed layer on top. It is also shown that, over a large range of parameters of the box diffusion model, the difference between the temperature calculated using a box diffusion model and the temperature calculated using a purely diffusive ocean is buried in the noise.

Lebedeff, Sergej A.↗

A minimally invasive, efficient method for propagation of full-field uncertainty in solid dynamics

In this work, we present a minimally invasive method for forward propagation of material property uncertainty to full-field quantities of interest in solid dynamics. Full-field uncertainty quantification enables the design of complex systems where quantities of interest, such as failure points, are not known a priori. The method, motivated by the well-known probability density function (PDF) propagation method of turbulence modeling, uses an ensemble of solutions to provide the joint PDF of desired quantities at every point in the domain. A small subset of the ensemble is computed exactly, and the remainder of the samples are computed with approximation of the evolution equations based on those exact solutions. Although the proposed method has commonalities with traditional interpolatory stochastic collocation methods applied directly to quantities of interest, it is distinct and exploits the parameter dependence and smoothness of the driving term of the evolution equations. The implementation is model independent, storage and communication efficient, and straightforward. We demonstrate its efficiency, accuracy, scaling with dimension of the parameter space, and convergence in distribution with two problems: a quasi-one-dimensional bar impact, and a two material notched plate impact. For the bar impact problem, we provide an analytical solution to PDF of the solution fields for method validation. With the notched plate problem, we also demonstrate good parallel efficiency and scaling of the method.

42 ENGINEERING↗