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 793 records · Page 44

Stokes problems for moving half-planes

New exact solutions of the Navier-Stokes equations are obtained for the unbounded and bounded oscillatory and impulsive tangential edgewise motion of touching half-infinite plates in their own plane. In contrast to Stokes classical solutions for the harmonic and impulsive motion of an infinite plane wall, where the solutions are separable or have a simple similarity form, the present solutions have a two-dimensional structure in the near region of the contact between the half-infinite plates. Nevertheless, it is possible to obtain relatively simple closed-form solutions for the flow field in each case by defining new variables which greatly simplify the r- and theta- dependence of the solutions in the vicinity of the contact region. These solutions for flow in a half-infinite space are then extended to bounded flows in a channel using an image superposition technique. The impulsive motion has application to the motion near geophysical faults, whereas the oscillatory motion has arisen in the design of a novel oscillating half-plate flow chamber for examining the effect of fluid shear stress on cultured cell monolayers.

NASA Discipline Musculoskeletal↗

A Study of Supersonic Surface Sources: The Ffowcs Williams-Hawkings Equation and the Kirchhoff Formula

In this paper we address the mathematical problem of noise generation from high speed moving surfaces. The problem we are solving is the linear wave equation with sources on a moving surface. The Ffowcs Williams-Hawkings (FW-H) equation as well as the govern- ing equation for deriving the Kirchhoff formula for moving surfaces are both this type of partial differential equation. We give a new exact solution of this problem here in closed form which is valid for subsonic and supersonic motion of the surface but it is particularly suitable for supersonically moving surfaces. This new solution is the simplest of all high speed formulations of Langley and is denoted formulation 4 following the tradition of numbering of our major results for the prediction of the noise of rotating blades. We show that for a smooth surface moving at supersonic speed, our solution has only removable singularities. Thus it can be used for numerical work.

Farassat, F.↗

Implicit/Multigrid Algorithms for Incompressible Turbulent Flows on Unstructured Grids

An implicit code for computing inviscid and viscous incompressible flows on unstructured grids is described. The foundation of the code is a backward Euler time discretization for which the linear system is approximately solved at each time step with either a point implicit method or a preconditioned Generalized Minimal Residual (GMRES) technique. For the GMRES calculations, several techniques are investigated for forming the matrix-vector product. Convergence acceleration is achieved through a multigrid scheme that uses non-nested coarse grids that are generated using a technique described in the present paper. Convergence characteristics are investigated and results are compared with an exact solution for the inviscid flow over a four-element airfoil. Viscous results, which are compared with experimental data, include the turbulent flow over a NACA 4412 airfoil, a three-element airfoil for which Mach number effects are investigated, and three-dimensional flow over a wing with a partial-span flap.

Anderson, W. Kyle↗

Three-Dimensional Postbuckling Analysis of Curved Beams

Presented here is a method of solving highly flexible curved beam undergoing huge static or quasi-static deformations. A geometrically exact beam theory based on the use of Jaumann stresses and strains and exact coordinate transformation is presented in terms of 17 first-order ordinary differential equations, and a multiple shooting method is used to solve the corresponding nonlinear two-point boundary value problems. The geometrically exact beam theory accounts far large rotations, large displacements, initial curvatures, extensionality, and transverse shear strains. Four examples are used to demonstrate this method, including a rotating clamped-free beam under the influence of gravity and centrifugal forces, an L-frame subjected to an in-plane tip load, a circular arch subjected to a concentrated load, and a clamped-hinged helical spring subjected to an axial displacement. Results show that the combination of the multiple shooting method and the geometrically exact beam theory works very well. Moreover, the obtained numerically exact solutions can be used to verify the accuracy of nonlinear finite element codes for nonlinear analysis of complex structures.

Pai, P. Frank↗

Hybrid Particle-Element Simulation of Impact on Composite Orbital Debris Shields

This report describes the development of new numerical methods and new constitutive models for the simulation of hypervelocity impact effects on spacecraft. The research has included parallel implementation of the numerical methods and material models developed under the project. Validation work has included both one dimensional simulations, for comparison with exact solutions, and three dimensional simulations of published hypervelocity impact experiments. The validated formulations have been applied to simulate impact effects in a velocity and kinetic energy regime outside the capabilities of current experimental methods. The research results presented here allow for the expanded use of numerical simulation, as a complement to experimental work, in future design of spacecraft for hypervelocity impact effects.

Fahrenthold, Eric P.↗

A Kernel-Free Particle-Finite Element Method for Hypervelocity Impact Simulation

An improved hybrid particle-finite element method has been developed for the simulation of hypervelocity impact problems. Unlike alternative methods, the revised formulation computes the density without reference to any kernel or interpolation functions, for either the density or the rate of dilatation. This simplifies the state space model and leads to a significant reduction in computational cost. The improved method introduces internal energy variables as generalized coordinates in a new formulation of the thermomechanical Lagrange equations. Example problems show good agreement with exact solutions in one dimension and good agreement with experimental data in a three dimensional simulation.

Park, Young-Keun↗

Model of Atmospheric Links on Optical Communications from High Altitude

Optical communication links have the potential to solve many of the problems of current radio and microwave links to satellites and high-altitude aircraft. The higher frequency involved in optical systems allows for significantly greater signal bandwidth, and thus information transfer rate, in excess of 10 Gbps, and the highly directional nature of laser-based signals eliminates the need for frequency-division multiplexing seen in radio and microwave links today. The atmosphere, however, distorts an optical signal differently than a microwave signal. While the ionosphere is one of the most significant sources of noise and distortion in a microwave or radio signal, the lower atmosphere affects an optical signal more significantly. Refractive index fluctuations, primarily caused by changes in atmospheric temperature and density, distort the incoming signal in both deterministic and nondeterministic ways. Additionally, suspended particles, such as those in haze or rain, further corrupt the transmitted signal. To model many of the atmospheric effects on the propagating beam, we use simulations based on the beam-propagation method. This method, developed both for simulation of signals in waveguides and propagation in atmospheric turbulence, separates the propagation into a diffraction and refraction problem. The diffraction step is an exact solution, within the limits of numerical precision, to the problem of propagation in free space, and the refraction step models the refractive index variances over a segment of the propagation path. By applying refraction for a segment of the propagation path, then diffracting over that same segment, this method forms a good approximation to true propagation through the atmospheric medium. Iterating over small segments of the total propagation path gives a good approximation to the problem of propagation over the entire path. Parameters in this model, such as initial beam profile and atmospheric constants, are easily modified in a simulation such as this, which allows for the rapid analysis of different propagation scenarios. Therefore, this method allows the development of a near-optimal system design for a wide range of situations, typical of what would be seen in different atmospheric conditions over a receiving ground station. A simulation framework based upon this model was developed in FORTRAN, and for moderate grid sizes and propagation distances these simulations are computable in reasonable time on a standard workstation. This presentation will discuss results thus far.

Subich, Christopher↗

Nonlinear Evolution of Alfvenic Wave Packets

Alfven waves are a ubiquitous feature of the solar wind. One approach to studying the evolution of such waves has been to study exact solutions to approximate evolution equations. Here we compare soliton solutions of the Derivative Nonlinear Schrodinger evolution equation (DNLS) to solutions of the compressible MHD equations.

nonlinear evolution alfvenic wave packets solar wi↗

The SMM Model as a Boundary Value Problem Using the Discrete Diffusion Equation

A generalized single step stepwise mutation model (SMM) is developed that takes into account an arbitrary initial state to a certain partial difference equation. This is solved in both the approximate continuum limit and the more exact discrete form. A time evolution model is developed for Y DNA or mtDNA that takes into account the reflective boundary modeling minimum microsatellite length and the original difference equation. A comparison is made between the more widely known continuum Gaussian model and a discrete model, which is based on modified Bessel functions of the first kind. A correction is made to the SMM model for the probability that two individuals are related that takes into account a reflecting boundary modeling minimum microsatellite length. This method is generalized to take into account the general n-step model and exact solutions are found. A new model is proposed for the step distribution.

Campbell, Joel↗

Optical Characterization of Metallic Aerosols

Airborne metallic particulates from industry and urban sources are highly conducting aerosols. The characterization of these pollutant particles is important for environment monitoring and protection. Because these metallic particulates are highly reflective, their effect on local weather or regional radiation budget may also need to be studied. In this work, light scattering characteristics of these metallic aerosols are studied using exact solutions on perfectly conducting spherical and cylindrical particles. It is found that for perfectly conducting spheres and cylinders, when scattering angle is larger than approx. 90 deg. the linear polarization degree of the scattered light is very close to zero. This light scattering characteristics of perfectly conducting particles is significantly different from that of other aerosols. When these perfectly conducting particles are immersed in an absorbing medium, this light scattering characteristics does not show significant change. Therefore, measuring the linear polarization of scattered lights at backward scattering angles can detect and distinguish metallic particulates from other aerosols. This result provides a great potential of metallic aerosol detection and monitoring for environmental protection.

Sun, Wenbo↗

Magnetofluid Turbulence in the Solar Wind

The solar wind shows striking characteristics that suggest that it is a turbulent magnetofluid, but the picture is not altogether simple. From the earliest observations, a strong correlation between magnetic fluctuations and plasma velocity fluctuations was noted. The high corrections suggest that the fluctuations are Alfven waves. In addition, the power spectrum of the magnetic fluctuation showed evidence of an inertial range that resembled that seen in fully-developed fluid turbulence. Alfven waves, however, are exact solutions of the equations of incompressible magnetohydrodynamics. Thus, there was a puzzle: how can a magnetofluid consisting of Alfven waves be turbulent? The answer lay in the role of velocity shears in the solar wind that could drive turbulent evolution. Puzzles remain: for example, the power spectrum of the velocity fluctuations is less steep than the slope of the magnetic fluctuations, nor do we understand even now why the solar wind appears to be nearly incompressible with a -5/3 power-spectral index.

Goldstein, Melvyn L.↗

Cosmological Inflation: A Personal Perspective

We present a brief review of Cosmological Inflation from the personal perspective of the speaker who almost 30 years ago proposed a way of resolving the problem of Cosmological Horizon by employing certain notions and developments from the field of High Energy Physics. Along with a brief introduction of the Horizon and Flatness problems of standard cosmology, this lecture concentrates on personal reminiscing of the notions and ideas that prevailed and influenced the author's thinking at the time. The lecture then touches upon some more recent developments related to the subject including exact solutions to conformal gravity that provide a first principles emergence of a characteristic acceleration in the universe and concludes with some personal views concerning the direction that the cosmology field has taken in the past couple of decades and certain speculations some notions that may indicate future directions of research.

Kazanas, Demos↗

Asymptotic Energies and QED Shifts for Rydberg States of Helium

This paper reviews progress that has been made in obtaining essentially exact solutions to the nonrelativistic three-body problem for helium by a combination of variational and asymptotic expansion methods. The calculation of relativistic and quantum electrodynamic corrections by perturbation theory is discussed, and in particular, methods for the accurate calculation of the Bethe logarithm part of the electron self energy are presented. As an example, the results are applied to the calculation of isotope shifts for the short-lived 'halo' nucleus He-6 relative to He-4 in order to determine the nuclear charge radius of He-6 from high precision spectroscopic measurements carried out at the Argonne National Laboratory. The results demonstrate that the high precision that is now available from atomic theory is creating new opportunities to create novel measurement tools, and helium, along with hydrogen, can be regarded as a fundamental atomic system whose spectrum is well understood for all practical purposes.

Drake, G.W.F.↗

Turbulence and Global Properties of the Solar Wind

The solar wind shows striking characteristics that suggest that it is a turbulent magnetofluid, but the picture is not altogether simple. From the earliest observations, a strong correlation between magnetic fluctuations and plasma velocity fluctuations as noted. The high corrections suggest that the fluctuations are Alven waves. In addition, the power spectrum of the magnetic fluctuation showed evidence of an inertial range that resembled that seen in fully-developed fluid turbulence. Alfven waves, however, are exact solutions of the equations of an incompressible magnetohydrodynamics. Thus, there was a puzzle: how can a magnetofluid consisting of Alfven waves be turbulent? The answer lay in the role of velocity shears in the solar wind that could drive turbulent evolution. Puzzles remain: for example, the power spectrum of the velocity fluctuations is less steep than the slope of the magnetic fluctuations. The plasma in the magnetic tail of Earth's magnetosphere also shows aspects of turbulence, as does the plasma in the dayside magnetosphere near the poles--the so-called dayside cusps.

Goldstein, Melvyn L.↗

A Refined Zigzag Beam Theory for Composite and Sandwich Beams

A new refined theory for laminated composite and sandwich beams that contains the kinematics of the Timoshenko Beam Theory as a proper baseline subset is presented. This variationally consistent theory is derived from the virtual work principle and employs a novel piecewise linear zigzag function that provides a more realistic representation of the deformation states of transverse-shear flexible beams than other similar theories. This new zigzag function is unique in that it vanishes at the top and bottom bounding surfaces of a beam. The formulation does not enforce continuity of the transverse shear stress across the beam s cross-section, yet is robust. Two major shortcomings that are inherent in the previous zigzag theories, shear-force inconsistency and difficulties in simulating clamped boundary conditions, and that have greatly limited the utility of these previous theories are discussed in detail. An approach that has successfully resolved these shortcomings is presented herein. Exact solutions for simply supported and cantilevered beams subjected to static loads are derived and the improved modelling capability of the new zigzag beam theory is demonstrated. In particular, extensive results for thick beams with highly heterogeneous material lay-ups are discussed and compared with corresponding results obtained from elasticity solutions, two other zigzag theories, and high-fidelity finite element analyses. Comparisons with the baseline Timoshenko Beam Theory are also presented. The comparisons clearly show the improved accuracy of the new, refined zigzag theory presented herein over similar existing theories. This new theory can be readily extended to plate and shell structures, and should be useful for obtaining relatively low-cost, accurate estimates of structural response needed to design an important class of high-performance aerospace structures.

Tessler, Alexander↗

Simple Formulas and Results for Buckling-Resistance and Stiffness Design of Compression-Loaded Laminated-Composite Cylinders

Simple formulas for the buckling stress of homogeneous, specially orthotropic, laminated-composite cylinders are presented. The formulas are obtained by using nondimensional parameters and equations that facilitate general validation, and are validated against the exact solution for a wide range of cylinder geometries and laminate constructions. Results are presented that establish the ranges of the nondimensional parameters and coefficients used. General results, given in terms of the nondimensional parameters, are presented that encompass a wide range of geometries and laminate constructions. These general results also illustrate a wide spectrum of behavioral trends. Design-oriented results are also presented that provide a simple, clear indication of laminate composition on critical stress, critical strain, and axial stiffness. An example is provided to demonstrate the application of these results to thin-walled column designs.

Nemeth, Michael P.↗

Development and Verification of the Charring, Ablating Thermal Protection Implicit System Simulator

The development and verification of the Charring Ablating Thermal Protection Implicit System Solver (CATPISS) is presented. This work concentrates on the derivation and verification of the stationary grid terms in the equations that govern three-dimensional heat and mass transfer for charring thermal protection systems including pyrolysis gas flow through the porous char layer. The governing equations are discretized according to the Galerkin finite element method (FEM) with first and second order fully implicit time integrators. The governing equations are fully coupled and are solved in parallel via Newton s method, while the linear system is solved via the Generalized Minimum Residual method (GMRES). Verification results from exact solutions and Method of Manufactured Solutions (MMS) are presented to show spatial and temporal orders of accuracy as well as nonlinear convergence rates.

Amar, Adam J.↗