Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical relativity”

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 1,099 records · Page 61

Materials processing in a centrifuge - Numerical modeling of macrogravity effects

The fluid mechanics associated with crystal growth processes on a centrifuge is investigated. A simple scaling analysis is used to examine the relative magnitudes of the forces acting on the system and good agreement is obtained with previous studies. A two-dimensional model of crystal growth on a centrifuge is proposed and calculations are undertaken to help in understanding the fundamental transport processes within the crystal growth cell. Results from three-dimensional calculations of actual centrifuge-based crystal growth systems are presented both for the thermodynamically stable and unstable configurations. The calculations show the existence of flow bifurcations in certain configurations but not in all instances. The numerical simulations also show that the centrifugal force is the dominant stabilizing force on fluid convection in the stable configuration. The stabilizing influence of the Coriolis force is found to be only secondary in nature. No significant impact of gravity gradient is found in the calculations. Simulations of unstable configurations show that the Coriolis force has a stabilizing influence on fluid motion by delaying the onset of unsteady convection. Detailed flow and thermal field characteristics are presented for all the different cases that are simulated.

Ramachandran, N.↗

Numerical algorithms for computations of feedback laws arising in control of flexible systems

Several continuous models will be examined, which describe flexible structures with boundary or point control/observation. Issues related to the computation of feedback laws are examined (particularly stabilizing feedbacks) with sensors and actuators located either on the boundary or at specific point locations of the structure. One of the main difficulties is due to the great sensitivity of the system (hyperbolic systems with unbounded control actions), with respect to perturbations caused either by uncertainty of the model or by the errors introduced in implementing numerical algorithms. Thus, special care must be taken in the choice of the appropriate numerical schemes which eventually lead to implementable finite dimensional solutions. Finite dimensional algorithms are constructed on a basis of a priority analysis of the properties of the original, continuous (infinite diversional) systems with the following criteria in mind: (1) convergence and stability of the algorithms and (2) robustness (reasonable insensitivity with respect to the unknown parameters of the systems). Examples with mixed finite element methods and spectral methods are provided.

Lasiecka, Irena↗

Inviscid, nonadiabatic flow fields over blunt, sonic corner bodies for outer planet entry conditions by a method of integral relations

An investigation has been made into the ability of a method of integral relations to calculate inviscid zero degree angle of attack, radiative heating distributions over blunt, sonic corner bodies for some representative outer planet entry conditions is investigated. Comparisons have been made with a more detailed numerical method, a time asymptotic technique, using the same equilibrium chemistry and radiation transport subroutines. An effort to produce a second order approximation (two-strip) method of integral relations code to aid in this investigation is also described and a modified two-strip routine is presented. Results indicate that the one-strip method of integral relations cannot be used to obtain accurate estimates of the radiative heating distribution because of its inability to resolve thermal gradients near the wall. The two-strip method can sometimes be used to improve these estimates; however, the two-strip method has only a small range of conditions over which it will yield significant improvement over the one-strip method.

Gnoffo, P. A.↗

A Newton/upwind method and numerical study of shock wave/boundary layer interactions

The objective of the paper is two-fold. First, an upwind/central differencing method for solving the steady Navier-Stokes equations is described. The symmetric line relation method is used to solve the resulting algebraic system to achieve high computational efficiency. The grid spacings used in the calculations are determined from the triple-deck theory, in terms of Mach and Reynolds numbers and other flow parameters. Thus the accuracy of the numerical solutions is improved by comparing them with experimental, analytical, and other computational results. Secondly, the shock wave/boundary layer interactions are studied numerically, with special attention given to the flow separation. The concept of free interaction is confirmed. Although the separated region varies with Mach and Reynolds numbers, it is found that the transverse velocity component behind the incident shock, which has not been identified heretofore, is also an important parameter. A small change of this quantity is sufficient to eliminate the flow separation entirely.

Liou, Meng-Sing↗

Optimization methods and silicon solar cell numerical models

An optimization algorithm for use with numerical silicon solar cell models was developed. By coupling an optimization algorithm with a solar cell model, it is possible to simultaneously vary design variables such as impurity concentrations, front junction depth, back junction depth, and cell thickness to maximize the predicted cell efficiency. An optimization algorithm was developed and interfaced with the Solar Cell Analysis Program in 1 Dimension (SCAP1D). SCAP1D uses finite difference methods to solve the differential equations which, along with several relations from the physics of semiconductors, describe mathematically the performance of a solar cell. A major obstacle is that the numerical methods used in SCAP1D require a significant amount of computer time, and during an optimization the model is called iteratively until the design variables converge to the values associated with the maximum efficiency. This problem was alleviated by designing an optimization code specifically for use with numerically intensive simulations, to reduce the number of times the efficiency has to be calculated to achieve convergence to the optimal solution.

Girardini, K.↗

Impact crater scaling laws

Impact craters are numerous on planetary bodies and furnish important information about the composition and past histories of those bodies. The interpretation of that information requires knowledge about the fundamental aspects of impact cratering mechanics. Since the typical conditions of impacts are at a size scale and velocity far in excess of experimental capabilities, direct simulations are precluded. Therefore, one must rely on extrapolation from experiments of relatively slow impacts of very small bodies, using physically based scaling laws, or must study the actual cases of interest using numerical code solutions of the fundamental physical laws that govern these processes. A progress report is presented on research on impact cratering scaling laws, on numerical studies that were designed to investigate those laws, and on various applications of the scaling laws developed by the author and his colleagues. These applications are briefly reviewed.

Holsapple, K. A.↗

Studies in Astronomical Time Series Analysis: VII. An Enquiry Concerning Non-Linearity, the RMS-Mean Flux Relation, and log-Normal Flux Distributions

A broad and widely used class of stationary, linear, additive time series models can have statistical properties which many authors have asserted imply that the underlying process must be non-linear, non-stationary, multiplicative, or inconsistent with shot noise. This result is demonstrated with exact and numerical evaluation of the model flux distribution function and dependence of flux standard deviation on mean flux (here and in the literature called the rms-flux relation). These models can: (1) exhibit normal, log-normal or other flux distributions; (2) show linear or slightly non-linear rms-mean flux dependencies; as well as (3) match arbitrary second order statistics of the time series data. Accordingly the above assertions cannot be made on the basis of statistical time series analysis alone. Also idscussed are ambiguities in the meaning of terms relevant to this study – linear, stationary and multiplicative – and functions that can transform observed fluxes to a normal distribution as well or better than the logarithm.

Jeffrey D Scargle↗

The numerical analysis of the rotational theory for the formation of lunar globules

The morphology of lunar globules is studied through the application of a numerical analysis of their rotation in space during cooling. It is assumed that molten rock is shot from the surface of the moon, solidifies in space above the moon and then falls back to the surface. The rotational theory studied makes the following assumptions: the volume of the molten rock does not change during cooling; the angular momentum is conserved; there are no internal motions because of the high viscosity of the molten rock, i.e., in equilibrium the globule is rotating as a rigid body; finally, the kinetic reaction of the globule to the forces is fast relative to the rate of cooling, i.e., the globule reaches equilibrium at constant energy. These assumptions are subjected to numerical analysis yielding good agreement between the actual globule shapes and the numerical results, but leaving some doubt as to the validity of the rotational theory due to the failure to establish the existence of true local minima and an incomplete understanding of the thermokentics.

Ross, J.↗

Developing Lightning Prediction Tools for the CCAFS Dual-Polarimetric Radar

NASA Marshall Space Flight Center and the University of Alabama Huntsville are collaborating with the 45th Weather Squadron (45WS) to develop improved lightning prediction capabilities for the new C-band dual-polarimetric weather radar being acquired for use by 45WS and launch weather forecasters at Cape Canaveral Air Force Station (CCAFS). In particular, these algorithms will focus on lightning onset, cessation and combined lightning-radar applications for convective winds assessment. Research using radar reflectivity (Z) data for prediction of lightning onset has been extensively discussed in the literature and subsequently applied by launch weather forecasters as it pertains to lightning nowcasting. Currently the forecasters apply a relatively straight forward but effective temperature-Z threshold algorithm for assessing the likelihood of lightning onset in a given storm. In addition, a layered VIL above the freezing level product is used as automated guidance for the onset of lightning. Only limited research and field work has been conducted on lightning cessation using Z and vertically-integrated Z for determining cessation. Though not used operationally vertically-integrated Z (basis for VIL) has recently shown promise as a tool for use in nowcasting lightning cessation. The work discussed herein leverages and expands upon these and similar reflectivity-threshold approaches via the application/addition of over two decades of polarimetric radar research focused on distinct multi-parameter radar signatures of ice/mixed-phase initiation and ice-crystal orientation in highly electrified convective clouds. Specifically, our approach is based on numerous previous studies that have observed repeatable patterns in the behavior of the vertical hydrometeor column as it relates to the temporal evolution of differential reflectivity and depolarization (manifested in either LDR or p(sub hv)), development of in-situ mixed and ice phase microphysics, electric fields, and ensuing lightning in the sub-tropical/tropical convection typical of the southeastern U.S., Maritime Continent, and southwestern Amazon. The polarimetric signatures detected in this setting provide a basis for automated 3-D detection of hydrometeor types in fuzzy logic hydrometeor identification algorithms (HID). Our working hypothesis is that improvement in lightning onset warning lead time and specificity for a given storm, relative to application of a Z-threshold algorithm, should arise as a consequence of the ability of dual-polarimetric radar to unambiguously detect and identify (through HID algorithms) the updraft elevation of rain-water cores above the freezing level and subsequent onset of drop freezing, riming, and robust mixed phase processes leading to significant charge separation and lightning. This type of algorithm, though dependent on the quality of the polarimetric data should be less susceptible to variable Z-calibration that can impact a given Z-threshold approach. To facilitate development of the algorithm while the 45WS dual-pol radar is in its current test stages and to evaluate the impact of polarimetric data quality (e.g., modified scan parameters and sampling) on the ensuing algorithms, we are using the ARMOR C-band dual-pol radar in Huntsville combined with N. Alabama LMA data and ARMOR HID algorithms [NCAR algorithm modified for application at C-band] in a testbed fashion. For lightning cessation we are revisiting the application of differential propagation phase variables for the monitoring of ice crystal alignment driven by in-cloud electric fields combined with metrics of ice water path (i.e., vertically integrated reflectivity). Importantly it should be noted that this approach is still very much a research topic and as such, we will explore operational applications that involve radar frequencies other than C-Band by using the UAH MAX X-band dual-pol radar in slow staring modes.

Petersen, W. A.↗

Doppler frequency in interplanetary radar and general relativity

The change of frequency of an interplanetary radar signal sent from the earth to another planet or to a space probe is worked out according to general relativity. The Schwarzschild spacetime is employed and its null geodesics control the motion of the signals. Exact Doppler frequency formulas are derived for one-way and two-way radar in terms of an arbitrary Schwarzschild radial coordinate. A reduction to the special relativity case is used to interpret the formulas in terms of the relative radial velocity of emitter and target. The general relativity corrections are worked out approximately for each of three possible Schwarzschild radial coordinates, and a numerical example is given. The amount of the correction is different according as one or the other of the Schwarzschild coordinates is identified with the radius vector deduced from classical celestial mechanics. The identification problem is discussed.

Mcvittie, G. C.↗

Optimal policies for identification of stochastic linear systems

The problem of designing closed-loop policies for identification of multiinput-multioutput linear discrete-time systems with random time-varying parameters is considered in this paper using a Bayesian approach. A sensitivity index gives a measure of performance for the closed-loop laws. The computation of the optimal laws is shown to be nontrivial, an exercise in stochastic control, but open-loop, affine, and open-loop feedback optimal inputs are shown to yield tractable problems. Numerical examples are given. For time-invariant systems, the criterion considered is shown to be related to the trace of the information matrix associated with the system.

Lopez-Toledo, A. A.↗

Solar generated quasi-biennial geomagnetic variation

The existence of highly correlated quasi-biennial variations in the geomagnetic field and in solar activity is demonstrated. The analysis uses a numerical filter technique applied to monthly averages of the geomagnetic horizontal component and of the Zurich relative sunspot number. Striking correlations are found between the quasi-biennial geomagnetic variations determined from several magnetic observatories located at widely different longitudes, indicating a worldwide nature of the obtained variation. The correlation coefficient between the filtered Dst index and the filtered relative sunspot number is found to be -0.79 at confidence level greater than 99% with a time-lag of 4 months, with solar activity preceding the Dst variation. The correlation between the unfiltered data of Dst and of the sunspot number is also high with a similar time-lag. Such a timelag has not been discussed in the literature, and a further study is required to establish the mode of sun-earth relationship that gives this time delay.

Sugiura, M.↗

A rapid solver for hyperbolic systems of equations

A numerical method is proposed for solving the time-dependent compressible Navier-Stokes equations in two dimensions. The equations are time-split into a hyperbolic part and a parabolic part. This paper describes an explicit numerical method for solving the hyperbolic (inviscid) part. The hyperbolic operator is explicit, conservative, uses characteristic relations to predict convection and pressure fields, and is stable under a specified condition.

Maccormack, R. W.↗

Application of LANDSAT images to wetland study and land use classification in west Tennessee, part 1

The author has identified the following significant results. densitometric analysis was performed on LANDSAT data to permit numerical classification of objects observed in the imagery on the basis of measurements of optical density. Relative light transmission measurements were taken on four types of scene elements in each of three LANDSAT black and white bands in order to determine which classification could be distinguished. The analysis of band 6 determined forest and agricultural classifications, but not the urban and wetlands. Both bands 4 and 5 showed a significant difference existed between the confirmed classification of wetlands-agriculture, and urban areas. Therefore, the combination of band 6 with either 4 or 5 would permit the separation of the urban from the wetland classification. To enhance the urban and wetland boundaries, the LANDSAT black and white bands were combined in a multispectral additive color viewer. Several combinations of filters and light intensities were used to obtain maximum discrimination between points of interest. The best results for enhancing wetland boundaries and urban areas were achieved by using a color composite (a blue, green, and red filter on bands 4, 5 and 6 respectively).

Shahrokhi, F.↗

Solar-generated quasi-biennial geomagnetic variation

The existence of highly correlated quasi-biennial variations in the geomagnetic field and in solar activity is demonstrated. The analysis uses a numerical filter technique applied to monthly averages of the geomagnetic horizontal component and of the Zurich relative sunspot number. Striking correlations are found between the quasi-biennial geomagnetic variations determined from several magnetic observatories located at widely different longitudes, indicating the worldwide nature of the obtained variation. The correlation coefficient between the filtered Dst index and the filtered relative sunspot number is found to be -0.79 at confidence level greater than 99% with a time lag of 4 months, solar activity preceding the Dst variation. The correlation between the unfiltered data of Dst and those of the sunspot number is also high with a similar time lag. Such a time lag has not been discussed in the literature, and a further study is required to establish the mode of sun-earth relationship that gives this time delay.

Sugiura, M.↗

The influence of a local wall deformation on the development of natural instabilities in a laminar boundary layer

The natural instabilities which propagate in the laminar boundary layer of a flat plate composed of intermittent wave trains are described. A spectral analysis determines the frequency range and gives a frequency and the harmonic 2 only if there is a wall deformation. This analysis provides the amplitude modulation spectrum of the instabilities. Plots of the evolution of power spectral density are compared with the numerical results obtained from the resolve of the Orr-Sommerfeld equation, while the harmonic is related to a micro-recirculating flow near the wall deformation.

Burnel, S.↗

Strain energy release rate determination of stress intensity factors by finite element methods

The stiffness derivative finite element technique is used to determine the Mode I stress intensity factors for three-crack configurations. The geometries examined include the double edge notch, single edge notch, and the center crack. The results indicate that when the specified guidelines of the Stiffness Derivative Method are used, a high degree of accuracy can be achieved with an optimized, relatively coarse finite element mesh composed of standard, four-node, plane strain, quadrilateral elements. The numerically generated solutions, when compared with analytical ones, yield results within 0.001 percent of each other for the double edge crack, 0.858 percent for the single edge crack, and 2.021 percent for the center crack.

Walsh, R. M., Jr.↗

Response of hot element wall shear stress gages in laminar oscillating flows

An experimental investigation of the time-dependent response of hot element wall shear stress gages in unsteady periodic air flows is reported. The study has focused on wall shear stress in laminar oscillating flows produced on a flat plate by a free stream velocity composed of a mean component and a superposed sinusoidal variation. Two types of hot element gages, platinum film and flush wire, were tested for values of reduced frequency ranging from 0.14 to 2.36. Values of the phase angle of the wall shear stress variation relative to the free stream velocity, as indicated by the hot element gages, are compared with numerical prediction. The comparisons show that the gages indicate a wall shear stress variation that lags the true variation, and that the gages will also not indicate the correct wall shear stress variation in periodic turbulent flows.

Cook, W. J.↗