Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Variational 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.

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

At least 811 records · Page 45

Magnetospheric electric fields deduced from drifting whistler paths.

Results of a study in which the amplitude of the E-W component E sub w of the convection electric field in the nightside magnetosphere has been inferred from the observed cross-L motions of whistler ducts within the plasmasphere, and several ducts distributed over 1 to 2 earth radii in L space and over plus or minus deg 15 around the longitude of the Eights, Antarctica, whistler station have been tracked simultaneously. The method appears capable of resolving fluctuations in E sub w with period T equal to approximately 15 min and rms amplitude as low as 0.05 mV/m. For variations with T greater than 1 hour the method has a sensitivity of the order of 0.01 mV/m. Three case studies are presented, two of which illustrate convection activity associated with relatively isolated substorms. In these two cases E sub w reversed from westward to eastward for a period following the decay of substorm bay activity. In the third case the substorm bay activity was prolonged, and E sub w remained westward and at enhanced levels until local dawn. Evidence was found that, at least in a limited longitudinal sector, perturbing substorm E sub w fields can penetrate deep within the plasmasphere. In two of the case studies comparisons of E sub w and the interplanetary magnetic-field theta component show evidence of a possible relation based on brief (less than or equal to 1 hour) southward excursions but not on long preceding southward events.

Carpenter, D. L.↗

Transition Measurements on the SWiFT Model in the National Transonic Facility

Boundary layer transition locations on the Swept Wing Flow Test (SWiFT) Hybrid Wing Body model were measured in the National Transonic Facility (NTF) at the NASA Langley Research Center using several different methods. These methods include sublimating chemical flow visualization for a limited number of low Reynolds number conditions, as well as mean static pressure and dynamic pressure measurements. Results are presented for a range of conditions in the NTF, demonstrating Mach and Reynolds number effects. The transition front locations obtained from the different techniques agree well. A variation of the static pressure transition detection method is described in which continuous angle of attack polars are performed. This approach enables detection of the transition front at more challenging (high Reynolds number) conditions than is possible with sparser data from pitch-pause angle-of-attack polars. The trip dot effectiveness is also verified using all of the available diagnostics. Additionally, low Reynolds number results are compared with those acquired on the same model at the Aircraft Research Association (ARA) facility, and these results show good agreement overall, suggesting similar freestream flow quality conditions of the two facilities.

National Transonic Facility↗

Transition Measurements on the SWiFT Model in the National Transonic Facility

Boundary layer transition locations on the Swept Wing Flow Test (SWiFT) Hybrid Wing Body model were measured in the National Transonic Facility (NTF) at the NASA Langley Research Center using several different methods. These methods include sublimating chemical flow visualization for a limited number of low Reynolds number conditions, as well as mean static pressure and dynamic pressure measurements. Results are presented for a range of conditions in the NTF, demonstrating Mach and Reynolds number effects. The transition front locations obtained from the different techniques agree well. A variation of the static pressure transition detection method is described in which continuous angle of attack polars are performed. This approach enables detection of the transition front at more challenging (high Reynolds number) conditions than is possible with sparser data from pitch-pause angle-of-attack polars. The trip dot effectiveness is also verified using all of the available diagnostics. Additionally, low Reynolds number results are compared with those acquired on the same model at the Aircraft Research Association (ARA) facility, and these results show good agreement overall, suggesting similar freestream flow quality conditions of the two facilities.

National Transonic Facility↗

Evaluation of an improved finite-element thermal stress calculation technique

A procedure for generating accurate thermal stresses with coarse finite element grids (Ojalvo's method) is described. The procedure is based on the observation that for linear thermoelastic problems, the thermal stresses may be envisioned as being composed of two contributions; the first due to the strains in the structure which depend on the integral of the temperature distribution over the finite element and the second due to the local variation of the temperature in the element. The first contribution can be accurately predicted with a coarse finite-element mesh. The resulting strain distribution can then be combined via the constitutive relations with detailed temperatures from a separate thermal analysis. The result is accurate thermal stresses from coarse finite element structural models even where the temperature distributions have sharp variations. The range of applicability of the method for various classes of thermostructural problems such as in-plane or bending type problems and the effect of the nature of the temperature distribution and edge constraints are addressed. Ojalvo's method is used in conjunction with the SPAR finite element program. Results are obtained for rods, membranes, a box beam and a stiffened panel.

Camarda, C. J.↗

Learning process mapping heuristics under stochastic sampling overheads

A statistical method was developed previously for improving process mapping heuristics. The method systematically explores the space of possible heuristics under a specified time constraint. Its goal is to get the best possible heuristics while trading between the solution quality of the process mapping heuristics and their execution time. The statistical selection method is extended to take into consideration the variations in the amount of time used to evaluate heuristics on a problem instance. The improvement in performance is presented using the more realistic assumption along with some methods that alleviate the additional complexity.

Ieumwananonthachai, Arthur↗

Clutter filter design considerations for Airborne Doppler radar detection of windshear

The problem of clutter rejection when processing down-looking Doppler radar returns from a low altitude airborne platform is a paramount problem. With radar as a remote sensor for detecting and predicting windshear in the vicinity of an urban airport, dynamic range requirements can exceed 50 dB because of high clutter to signal ratios. This presentation describes signal processing considerations in the presence of distributed and/or discrete clutter interference. Previous analyses have considered conventional range cell processing of radar returns from a rigidly mounted radar platform using either the Fourier or the pulse-pair method to estimate average windspeed and windspeed variation within a cell. Clutter rejection has been based largely upon analyzing a particular environment in the vicinity of the radar and employing a variety of techniques to reduce interference effects including notch filtering, Fourier domain line editing, and use of clutter maps. For the airborne environment the clutter characteristics may be somewhat different. Conventional clutter rejection methods may have to be changed and new methods will probably be required to provide useful signal to noise ratios. Various considerations are described. A major thrust has been to evaluate the effect of clutter rejection filtering upon the ability to derive useful information from the post-filter radar data. This analysis software is briefly described. Finally, some ideas for future analysis are considered including the use of adaptive filtering for clutter rejection and the estimation of windspeed spatial gradient directly from radar returns as a means of reducing the effects of clutter on the determination of a windshear hazard.

Baxa, Ernest G., Jr.↗

Investigation of mounting discrete chip components for hybrid microelectronic applications.

In the assembly of hybrid microcircuits, passive and active discrete chip components are attached by epoxy cements and by solders. This study was initiated to determine the integrity of such attachments. Alumina substrates were subjected to a group of thermal and mechanical test environments defined by test methods in MIL-STD-883. Lot samples of substrates with assembly variations were subjected to separately applied testing and to serial testing. These test methods were applied sequentially to the same lot of substrates. Epoxy cements were generally satisfactory with respect to separate tests, although sequential tests severely degraded epoxy bond strengths of all types. Bond strengths to unglazed substrates were slightly superior to those of the same components on glazed types. Solder joints were vulnerable to catastrophic failure in several environments, especially temperature cycling and thermal shock. Vibration and impact shock also produced solder joint failures.

Caruso, S. V.↗

Development of computational methods for heavy lift launch vehicles

The research effort has been focused on the development of an advanced flow solver for complex viscous turbulent flows with shock waves. The three-dimensional Euler and full/thin-layer Reynolds-averaged Navier-Stokes equations for compressible flows are solved on structured hexahedral grids. The Baldwin-Lomax algebraic turbulence model is used for closure. The space discretization is based on a cell-centered finite-volume method augmented by a variety of numerical dissipation models with optional total variation diminishing limiters. The governing equations are integrated in time by an implicit method based on lower-upper factorization and symmetric Gauss-Seidel relaxation. The algorithm is vectorized on diagonal planes of sweep using two-dimensional indices in three dimensions. A new computer program named CENS3D has been developed for viscous turbulent flows with discontinuities. Details of the code are described in Appendix A and Appendix B. With the developments of the numerical algorithm and dissipation model, the simulation of three-dimensional viscous compressible flows has become more efficient and accurate. The results of the research are expected to yield a direct impact on the design process of future liquid fueled launch systems.

Yoon, Seokkwan↗

Development of advanced Navier-Stokes solver

The objective of research was to develop and validate new computational algorithms for solving the steady and unsteady Euler and Navier-Stokes equations. The end-products are new three-dimensional Euler and Navier-Stokes codes that are faster, more reliable, more accurate, and easier to use. The three-dimensional Euler and full/thin-layer Reynolds-averaged Navier-Stokes equations for compressible/incompressible flows are solved on structured hexahedral grids. The Baldwin-Lomax algebraic turbulence model is used for closure. The space discretization is based on a cell-centered finite-volume method augmented by a variety of numerical dissipation models with optional total variation diminishing limiters. The governing equations are integrated in time by an implicit method based on lower-upper factorization and symmetric Gauss-Seidel relaxation. The algorithm is vectorized on diagonal planes of sweep using two-dimensional indices in three dimensions. Convergence rates and the robustness of the codes are enhanced by the use of an implicit full approximation storage multigrid method.

Yoon, Seokkwan↗

Useful extremum principle for the variational calculation of matrix elements

Variational principles are considered for the approximate evaluation of the diagonal matrix elements of an arbitrary known linear Hermitian operator. A method is derived that is immediately applicable to the variational determination of both the off-diagonal and diagonal matrix elements of normal and modified Green's functions.

Gerjuoy, E.↗

Variational approach to probabilistic finite elements

Probabilistic finite element method (PFEM), synthesizing the power of finite element methods with second-moment techniques, are formulated for various classes of problems in structural and solid mechanics. Time-invariant random materials, geometric properties, and loads are incorporated in terms of their fundamental statistics viz. second-moments. Analogous to the discretization of the displacement field in finite element methods, the random fields are also discretized. Preserving the conceptual simplicity, the response moments are calculated with minimal computations. By incorporating certain computational techniques, these methods are shown to be capable of handling large systems with many sources of uncertainties. By construction, these methods are applicable when the scale of randomness is not very large and when the probabilistic density functions have decaying tails. The accuracy and efficiency of these methods, along with their limitations, are demonstrated by various applications. Results obtained are compared with those of Monte Carlo simulation and it is shown that good accuracy can be obtained for both linear and nonlinear problems. The methods are amenable to implementation in deterministic FEM based computer codes.

Belytschko, T.↗

Variational approach to probabilistic finite elements

Probabilistic finite element methods (PFEM), synthesizing the power of finite element methods with second-moment techniques, are formulated for various classes of problems in structural and solid mechanics. Time-invariant random materials, geometric properties and loads are incorporated in terms of their fundamental statistics viz. second-moments. Analogous to the discretization of the displacement field in finite element methods, the random fields are also discretized. Preserving the conceptual simplicity, the response moments are calculated with minimal computations. By incorporating certain computational techniques, these methods are shown to be capable of handling large systems with many sources of uncertainties. By construction, these methods are applicable when the scale of randomness is not very large and when the probabilistic density functions have decaying tails. The accuracy and efficiency of these methods, along with their limitations, are demonstrated by various applications. Results obtained are compared with those of Monte Carlo simulation and it is shown that good accuracy can be obtained for both linear and nonlinear problems. The methods are amenable to implementation in deterministic FEM based computer codes.

Belytschko, T.↗

Small-scale experiments in STOVL ground effects

A series of tests was completed in which suckdown and fountain forces and pressures were measured on circular plates and twin-tandem-jet generic STOVL (short takeoff and vertical landing) configurations. The tests were conducted using a small-scale hover rig, for jet pressure ratios up to 6 and jet temperatures up to 700 F. The measured suckdown force on a circular plate with a central jet was greater than that found with a commonly used empirical prediction method. The present data showed better agreement with other sets of data. The tests of the generic STOVL configurations were conducted to provide force and pressure data with a parametric variation of parameters so that an empirical prediction method could be developed. The effects of jet pressure ratio and temperature were found to be small. Lift improvement devices were shown to substantially reduce the net suckdown forces.

Corsiglia, Victor R.↗

Small-scale experiments in STOVL ground effects

A series of tests has been completed in which suckdown and fountain forces and pressures were measured on circular plates and twin-tandem-jet generic STOVL (short takeoff and vertical landing) configurations. The tests were conducted using a small-scale hover rig, for jet pressure ratios up to 6 and jet temperatures up to 700 F. The measured suckdown force on a circular plate with a central jet was greater than that found with a commonly used empirical prediction method. The present data showed better agreement with other sets of data. The tests of the generic STOVL configurations were conducted to provide force and pressure data with a parametric variation of parameters so that an empirical prediction method cold be developed. The effects of jet pressure ratio and temperature were found to be small. Lift improvement devices were shown to substantially reduce the net suckdown forces.

Corsiglia, Victor R.↗

Geophysical interpretation of the magnetic anomalies of the Earth derived from MAGSAT data

The ambiguities about geophysical implications based on the correlation of scalar magnetic anomalies and geological features were investigated. A method was developed to convert scalar magnetic anomalies into a map of the lateral variations of magnetic susceptibility of the lithosphere. This map is directly correlated with the causative sources. The method is based on spherical harmonic analysis of lateral variations seen on the scalar magnetic anomaly map and those of the lithospheric magnetic susceptibility. The harmonic coefficients are related through the fundamental causality relationship governing a magnetized body and its associated scalar magnetic anomaly. The main features of the resulting magnetic susceptibility anomalies are outlined.

Arkani-Hamed, J.↗

(abstract) Optimal Low Thrust Trajectories Using Differential Inclusion Concepts

Low thrust propulsion systems typically have their greatest benefit for high energy missions or missions with large post-launch maneuver requirements. Missions which have been examined include main belt asteroid rendezvous, comet rendezvous, outer planet and Mercury orbiters, Pluto flyby, and solar probe missions. Low thrust mission design software used to determine these trajectories is based on two distinct formulations of the optimal control problem: the indirect and direct methods. The traditional approach (indirect) is to use the calculus of variations to obtain first order necessary conditions on the states and costates. In contrast, direct methods are conceptually different in that no explicit integration takes place. A direct method based on differential inclusion concepts has been developed and successfully used to compute low thrust trajectories. This new approach removes explicit control dependence from the problem thereby reducing the dimension of the parameter space for the nonlinear programming problem. Also when compared to other direct methods, fewer nonlinear constraints are required to represent the dynamics of the problem.

low↗

Solar tri-diurnal variation of cosmic rays in a wide range of rigidity

Solar tri-diurnal variations of cosmic rays have been analyzed in a wide range of rigidity, using data from neutron monitors, and the surface and underground muon telescopes for the period 1978-1983. The rigidity spectrum of the anisotropy in space is assumed to be of power-exponential type as (P/gamma P sub o) to the gamma exp (gamma-P/P sub o). By means of the best-fit method between the observed and the expected variations, it is obtained that the spectrum has a peak at P (=gamma P sub o) approx = 90 GV, where gamma=approx 3.0 and P sub o approx. 30 GV. The phase in space of the tri-diurnal variation is also obtained as 7.0 hr (15 hr and 23 hr LT), which is quite different from that of approx. 1 hr. arising from the axisymmetric distribution of cosmic rays with respect to the IMF.

Mori, S.↗

Buckling of cylindrical shells with random imperfections

Buckling of long cylindrical shells with random imperfections subjected to axial loads is treated using the method of truncated hierarchy. A system of homogeneous variational equations is set up to examine the existence of bifurcation in the neighborhood of the equilibrium state. These equations are linear and involve stationary random coefficients. The truncated-hierarchy method is applied, and characteristic equations are obtained. Various exponential cosine correlation functions associated with asymmetric imperfections are examined numerically. Qualitatively, the results obtained are as anticipated.

Fersht, R. S.↗