SEARCH · Engineering Papers

Results for “CALCULATION METHODS”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,135 records · Page 63

All sources NASA NTRS DOE OSTI

NASA NTRS1990

Prediction of vortex shedding from forebodies with chines

An engineering prediction method to calculate vortex shedding from noncircular forebodies with sharp chine edges in subsonic flow at large incidence angles is presented. The forebody is represented by two- and three-dimensional singularities, and the lee side vortex wake is modeled by discrete vortices in crossflow planes along the body. The computational procedure is described, and comparisons of measured and predicted surface pressure distributions and predicted flow field vectors are presented to illustrate the method.

Mendenhall, Michael R.↗

NASA NTRS1991

Nonequilibrium turbulence modeling effects on transonic vortical flows about delta wings

The Johnson-King turbulence model that is a viable method for calculating two dimensional transonic separated flows was extended into three dimensions. The implementation was done for Navier Stokes flow solvers written in general curvilinear coordinates. The present approach used in turbulence modeling is based on streamwise integration of an ordinary differential equation (o.d.e.) that governs the maximum Reynolds shear stress behavior. Streamwise integration of the o.d.e. approach was found to offer great mathematical simplicity and economy for three dimensional Navier Stokes methods. Thus, the new method is quick, simple, and very cheap. The new method was first checked against the data of a well known transonic axisymmetric bump experiment, and a good agreement was obtained. Later, the new method was used to compute the flow around a low aspect ratio wing in a transonic wind tunnel. Finally it was employed to study the nonequilibrium turbulence effects on the transonic vortical flows about a 65 deg sweep round leading edge delta wing.

Unver Kaynak↗

NASA NTRS1999

Optical model methods of predicting nuclide production cross sections from heavy ion fragmentation

Quantum mechanical optical potential methods for calculating inclusive isotope and element production cross sections from the fragmenting of heavy nuclei by intermediate- and high-energy protons and heavy ions are presented based upon a modified abrasion-ablation-FSI (frictional spectator interaction) collision model. The abrasion stage is treated as a quantum mechanical knockout process that leaves the residual prefragment in an excited state. Prefragment excitation energies are estimated using a combined liquid drop and FSI method. In ablation the prefragment deexcites by particle and photon emission to produce the final fragment. Contributions from electromagnetic dissociation to single nucleon removal cross sections are incorporated using a Weiszacker-Williams theory that includes electric dipole and electric quadrupole interactions. Estimates of elemental and isotopic production cross sections are in good agreement with published cross section measurements for a variety of projectile-target-beam energy combinations.

NASA Discipline Radiation Health↗

NASA NTRS1985

Interlaminar G(IIc) evaluation of toughened-resin matrix composites using the end-notched flexure test

The usefulness of the end-notched flexure (ENF) test for measuring G(IIc) for toughened composites was evaluated. Calculations were made using the measured flexural compliance from the ENF test and the measured axial modulus. Tests were conducted on three different composite materials, using laminates with two different span lengths and insert sizes. The error associated with neglecting nonlinear terms in the displacement-curvature relationship was analyzed to demonstrate the influence of specimen geometry on G(IIc) measurement. It is concluded that G(IIc) measurement may vary with the method of calculation, span length, or insert thickness, but these effects are minor compared to the large differences in measurements obtained for different materials. Hence, the ENF test is a valid one for screening materials for improvements in interlaminar shear fracture toughness.

Murri, G. B.↗

NASA NTRS1996

Stress and Damage in Polymer Matrix Composite Materials Due to Material Degradation at High Temperatures

This report describes analytical methods for calculating stresses and damage caused by degradation of the matrix constituent in polymer matrix composite materials. Laminate geometry, material properties, and matrix degradation states are specified as functions of position and time. Matrix shrinkage and property changes are modeled as functions of the degradation states. The model is incorporated into an existing composite mechanics computer code. Stresses, strains, and deformations at the laminate, ply, and micro levels are calculated, and from these calculations it is determined if there is failure of any kind. The rationale for the model (based on published experimental work) is presented, its integration into the laminate analysis code is outlined, and example results are given, with comparisons to existing material and structural data. The mechanisms behind the changes in properties and in surface cracking during long-term aging of polyimide matrix composites are clarified. High-temperature-material test methods are also evaluated.

McManus, Hugh L.↗

NASA NTRS1989

Power flows and Mechanical Intensities in structural finite element analysis

The identification of power flow paths in dynamically loaded structures is an important, but currently unavailable, capability for the finite element analyst. For this reason, methods for calculating power flows and mechanical intensities in finite element models are developed here. Formulations for calculating input and output powers, power flows, mechanical intensities, and power dissipations for beam, plate, and solid element types are derived. NASTRAN is used to calculate the required velocity, force, and stress results of an analysis, which a post-processor then uses to calculate power flow quantities. The SDRC I-deas Supertab module is used to view the final results. Test models include a simple truss and a beam-stiffened cantilever plate. Both test cases showed reasonable power flow fields over low to medium frequencies, with accurate power balances. Future work will include testing with more complex models, developing an interactive graphics program to view easily and efficiently the analysis results, applying shape optimization methods to the problem with power flow variables as design constraints, and adding the power flow capability to NASTRAN.

Hambric, Stephen A.↗

NASA NTRS1994

Improving Diffusing S-duct Performance by Secondary Flow Control

The objective of this research was to study ways to reduce inlet flow distortion (i.e., total pressure nonuniformity) and improve total pressure recovery in a diffusing S-duct. This was accomplished by controlling the development of secondary flows within the duct through the use of tapered-fin type vortex generators. Reported are results for the bare duct and seven different configurations of vortex generators. Data presented for each configuration include surface static pressure, surface flow visualization, and exit plane total pressure and transverse velocity. The performance of each configuration was assessed by calculating total pressure recovery and inlet distortion descriptors from the data and comparing them to the values for the bare duct. The best configuration tested reduced distortion (as measured by the DC(45) and DC(90) descriptors) by more than 50 percent while improving total pressure recovery by 0.5 percent. These results should provide valuable guidance in designing vortex generator installations in ducts and for assessing the accuracy of computational fluid dynamics (CFD) methods to calculate duct flows with installed vortex generators.

Reichert, Bruce A.↗

NASA NTRS1982

An efficient full potential implicit method based on characteristics for analysis of supersonic flows

A nonlinear aerodynamic prediction technique based on the full potential equation in conservation form has been developed for the treatment of supersonic flows. The method uses the theory of characteristic signal propagation to accurately simulate the flow structure, which includes shock waves and mixed elliptic-hyperbolic crossflow. An implicit approximate factorization scheme is employed to solve the finite-differenced equation. The necessary body-fitted grid system in every marching plane is generated numerically, using an elliptic grid solver. Results are shown for conical and nonconical wing-body combinations and compared with experimental data and Euler calculations. The method demonstrates an enormous savings in execution time and memory requirements over Euler methods.

Shankar, V.↗

NASA NTRS1988

A diagonally inverted LU implicit multigrid scheme for the 3-D Navier-Stokes equations and a two equation model of turbulence

An LU implicit multigrid algorithm is developed to calculate 3-D compressible viscous flows. This scheme solves the full 3-D Reynolds-Averaged Navier-Stokes equation with a two-equation kappa-epsilon model of turbulence. The flow equations are integrated by an efficient, diagonally inverted, LU implicit multigrid scheme while the kappa-epsilon equations are solved, uncoupled from the flow equations, by a block LU implicit algorithm. The flow equations are solved within the framework of the multigrid method using a four-grid level W-cycle, while the kappa-epsilon equations are iterated only on the finest grid. This treatment of the Reynolds-Averaged Navier-Stokes equations proves to be an efficient method for calculating 3-D compressible viscous flows.

Yokota, Jeffrey W.↗

NASA NTRS1989

A diagonally inverted LU implicit multigrid scheme for the 3-D Navier-Stokes equations and a two equation model of turbulence

Yokota, Jeffrey W.↗

NASA NTRS1995

A well-posed numerical method to track isolated conformal map singularities in Hele-Shaw flow

We present a new numerical method for calculating an evolving 2D Hele-Shaw interface when surface tension effects are neglected. In the case where the flow is directed from the less viscous fluid into the more viscous fluid, the motion of the interface is ill-posed; small deviations in the initial condition will produce significant changes in the ensuing motion. This situation is disastrous for numerical computation, as small round-off errors can quickly lead to large inaccuracies in the computed solution. Our method of computation is most easily formulated using a conformal map from the fluid domain into a unit disk. The method relies on analytically continuing the initial data and equations of motion into the region exterior to the disk, where the evolution problem becomes well-posed. The equations are then numerically solved in the extended domain. The presence of singularities in the conformal map outside of the disk introduces specific structures along the fluid interface. Our method can explicitly track the location of isolated pole and branch point singularities, allowing us to draw connections between the development of interfacial patterns and the motion of singularities as they approach the unit disk. In particular, we are able to relate physical features such as finger shape, side-branch formation, and competition between fingers to the nature and location of the singularities. The usefulness of this method in studying the formation of topological singularities (self-intersections of the interface) is also pointed out.

Baker, Gregory↗

NASA NTRS1996

Artificial Boundary Conditions for Computation of Oscillating External Flows

In this paper, we propose a new technique for the numerical treatment of external flow problems with oscillatory behavior of the solution in time. Specifically, we consider the case of unbounded compressible viscous plane flow past a finite body (airfoil). Oscillations of the flow in time may be caused by the time-periodic injection of fluid into the boundary layer, which in accordance with experimental data, may essentially increase the performance of the airfoil. To conduct the actual computations, we have to somehow restrict the original unbounded domain, that is, to introduce an artificial (external) boundary and to further consider only a finite computational domain. Consequently, we will need to formulate some artificial boundary conditions (ABC's) at the introduced external boundary. The ABC's we are aiming to obtain must meet a fundamental requirement. One should be able to uniquely complement the solution calculated inside the finite computational domain to its infinite exterior so that the original problem is solved within the desired accuracy. Our construction of such ABC's for oscillating flows is based on an essential assumption: the Navier-Stokes equations can be linearized in the far field against the free-stream back- ground. To actually compute the ABC's, we represent the far-field solution as a Fourier series in time and then apply the Difference Potentials Method (DPM) of V. S. Ryaben'kii. This paper contains a general theoretical description of the algorithm for setting the DPM-based ABC's for time-periodic external flows. Based on our experience in implementing analogous ABC's for steady-state problems (a simpler case), we expect that these boundary conditions will become an effective tool for constructing robust numerical methods to calculate oscillatory flows.

Tsynkov, S. V.↗

NASA NTRS2017

Determination of the Meteor Limiting Magnitude

We present our method to calculate the meteor limiting magnitude. The limiting meteor magnitude defines the faintest magnitude at which all meteors are still detected by a given system. An accurate measurement of the limiting magnitude is important in order to calculate the meteoroid flux from a meteor shower or sporadic source. Since meteor brightness is linked to meteor mass, the limiting magnitude is needed to calculate the limiting mass of the meteor flux measurement. The mass distribution of meteoroids is thought to follow a power law, thus being slightly off in the limiting magnitude can have a significant effect on the measured flux. Sky conditions can change on fairly short timescales; therefore one must monitor the meteor limiting magnitude at regular intervals throughout the night, rather than just measuring it once. We use the stellar limiting magnitude as a proxy of the meteor limiting magnitude. Our method for determining the stellar limiting magnitude and how we transform it into the meteor limiting magnitude is presented. These methods are currently applied to NASA's wide-field meteor camera network to determine nightly fluxes, but are applicable to other camera networks.

A Kingery↗

NASA NTRS1966

On the calculation of nonlinear vibrations of flexible plates and shallow shells by the small parameter method

Small parameter method used for calculating periodic vibrations of flexible shallow shells and plates

PLATE THEORY↗

NASA NTRS1973

Static penetration resistance of soils

Model test results were used to define the failure mechanism associated with the static penetration resistance of cohesionless and low-cohesion soils. Knowledge of this mechanism has permitted the development of a new analytical method for calculating the ultimate penetration resistance which explicitly accounts for penetrometer base apex angle and roughness, soil friction angle, and the ratio of penetration depth to base width. Curves relating the bearing capacity factors to the soil friction angle are presented for failure in general shear. Strength parameters and penetrometer interaction properties of a fine sand were determined and used as the basis for prediction of the penetration resistance encountered by wedge, cone, and flat-ended penetrometers of different surface roughness using the proposed analytical method. Because of the close agreement between predicted values and values measured in laboratory tests, it appears possible to deduce in-situ soil strength parameters and their variation with depth from the results of static penetration tests.

Durgunoglu, H. T.↗

NASA NTRS2001

Computationally efficient control allocation

A computationally efficient method for calculating near-optimal solutions to the three-objective, linear control allocation problem is disclosed. The control allocation problem is that of distributing the effort of redundant control effectors to achieve some desired set of objectives. The problem is deemed linear if control effectiveness is affine with respect to the individual control effectors. The optimal solution is that which exploits the collective maximum capability of the effectors within their individual physical limits. Computational efficiency is measured by the number of floating-point operations required for solution. The method presented returned optimal solutions in more than 90% of the cases examined; non-optimal solutions returned by the method were typically much less than 1% different from optimal and the errors tended to become smaller than 0.01% as the number of controls was increased. The magnitude of the errors returned by the present method was much smaller than those that resulted from either pseudo inverse or cascaded generalized inverse solutions. The computational complexity of the method presented varied linearly with increasing numbers of controls; the number of required floating point operations increased from 5.5 i, to seven times faster than did the minimum-norm solution (the pseudoinverse), and at about the same rate as did the cascaded generalized inverse solution. The computational requirements of the method presented were much better than that of previously described facet-searching methods which increase in proportion to the square of the number of controls.

Durham, Wayne↗

NASA NTRS1984

Costs of measuring leaf area index of corn

The magnitude of plant-to-plant variability of leaf area of corn plants selected from uniform plots was examined and four representative methods for measuring leaf area index (LAI) were evaluated. The number of plants required and the relative costs for each sampling method were calculated to detect 10, 20, and 50% differences in LAI using 0.05 and 0.01 tests of significance and a 90% probability of success (beta = 0.1). The natural variability of leaf area per corn plant was nearly 10%. Additional variability or experimental error may be introduced by the measurement technique employed and by nonuniformity within the plot. Direct measurement of leaf area with an electronic area meter had the lowest CV, required that the fewest plants be sampled, but required approximately the same amount of time as the leaf area/weight ratio method to detect comparable differences. Indirect methods based on measurements of length and width of leaves required more plants but less total time than the direct method. Unless the coefficients for converting length and width to area are verified frequently, the indirect methods may be biased. When true differences in LAI among treatments exceed 50% of mean, all four methods are equal. The method of choice depends on the resources available, the differences to be detected, and what additional information, such as leaf weight or stalk weight, is also desired.

Daughtry, C. S. T.↗

NASA NTRS1977

A user's guide for V174, a program using a finite difference method to analyze transonic flow over oscillating wings

The design and usage of a pilot program using a finite difference method for calculating the pressure distributions over harmonically oscillating wings in transonic flow are discussed. The procedure used is based on separating the velocity potential into steady and unsteady parts and linearizing the resulting unsteady differential equation for small disturbances. The steady velocity potential which must be obtained from some other program, is required for input. The unsteady differential equation is linear, complex in form with spatially varying coefficients. Because sinusoidal motion is assumed, time is not a variable. The numerical solution is obtained through a finite difference formulation and a line relaxation solution method.

Butler, T. D.↗