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,081 records · Page 60

Response of a Rotating Propeller to Aerodynamic Excitation

The flexural vibration of a rotating propeller blade with clamped shank is analyzed with the object of presenting, in matrix form, equations for the elastic bending moments in forced vibration resulting from aerodynamic forces applied at a fixed multiple of rotational speed. Matrix equations are also derived which define the critical speeds end mode shapes for any excitation order and the relation between critical speed and blade angle. Reference is given to standard works on the numerical solution of matrix equations of the forms derived. The use of a segmented blade as an approximation to a continuous blade provides a simple means for obtaining the matrix solution from the integral equation of equilibrium, so that, in the numerical application of the method presented, the several matrix arrays of the basic physical characteristics of the propeller blade are of simple form, end their simplicity is preserved until, with the solution in sight, numerical manipulations well-known in matrix algebra yield the desired critical speeds and mode shapes frame which the vibration at any operating condition may be synthesized. A close correspondence between the familiar Stodola method and the matrix method is pointed out, indicating that any features of novelty are characteristic not of the analytical procedure but only of the abbreviation, condensation, and efficient organization of the numerical procedure made possible by the use of classical matrix theory.

Arnoldi, Walter E.↗

The Nature of Global Large-scale Sea Level Variability in Relation to Atmospheric Forcing: A Modeling Study

The relation between large-scale sea level variability and ocean circulation is studied using a numerical model. A global primitive equaiton model of the ocean is forced by daily winds and climatological heat fluxes corresponding to the period from January 1992 to February 1996. The physical nature of the temporal variability from periods of days to a year, are examined based on spectral analyses of model results and comparisons with satellite altimetry and tide gauge measurements.

ocean circulation circulation models ocean models ↗

Practical Considerations for Optimizing Position Sensitivity in Arrays of Position-sensitive TES's

We are developing Position-Sensitive Transitions-Edge Sensors (PoST's) for future X-ray astronomy missions such as NASA's Constellation-X. The PoST consists of one or more Transitions Edge Sensors (TES's) thermally connected to a large X-ray absorber, which through heat diffusion, gives rise to position dependence. The development of PoST's is motivated by the desire to achieve the largest the focal-plan coverage with the fewest number of readout channels. In order to develop a practical array, consisting of an inner pixellated core with an outer array of large absorber PoST's, we must be able to simultaneously read out all (-1800) channels in the array. This is achievable using time division multiplexing (TDM), but does set stringent slew rate requirements on the array. Typically, we must damp the pulses to reduce the slew rate of the input signal to the TDM. This is achieved by applying a low-pass analog filter with large inductance to the signal. This attenuates the high frequency components of the signal, essential for position discrimination in PoST's, relative to the white noise of the readout chain and degrades the position sensitivity. Using numerically simulated data, we investigate the position sensing ability of typical PoST designs under such high inductance conditions. We investigate signal-processing techniques for optimal determination of the event position and discuss the practical considerations for real-time implementation.

Smith, Stephen J.↗

Higher-Order Finite Elements for Computing Thermal Radiation

Two variants of the finite-element method have been developed for use in computational simulations of radiative transfers of heat among diffuse gray surfaces. Both variants involve the use of higher-order finite elements, across which temperatures and radiative quantities are assumed to vary according to certain approximations. In this and other applications, higher-order finite elements are used to increase (relative to classical finite elements, which are assumed to be isothermal) the accuracies of final numerical results without having to refine computational meshes excessively and thereby incur excessive computation times. One of the variants is termed the radiation sub-element (RSE) method, which, itself, is subject to a number of variations. This is the simplest and most straightforward approach to representation of spatially variable surface radiation. Any computer code that, heretofore, could model surface-to-surface radiation can incorporate the RSE method without major modifications. In the basic form of the RSE method, each finite element selected for use in computing radiative heat transfer is considered to be a parent element and is divided into sub-elements for the purpose of solving the surface-to-surface radiation-exchange problem. The sub-elements are then treated as classical finite elements; that is, they are assumed to be isothermal, and their view factors and absorbed heat fluxes are calculated accordingly. The heat fluxes absorbed by the sub-elements are then transferred back to the parent element to obtain a radiative heat flux that varies spatially across the parent element. Variants of the RSE method involve the use of polynomials to interpolate and/or extrapolate to approximate spatial variations of physical quantities. The other variant of the finite-element method is termed the integration method (IM). Unlike in the RSE methods, the parent finite elements are not subdivided into smaller elements, and neither isothermality nor other unrealistic physical conditions are assumed. Instead, the equations of radiative heat transfer are integrated numerically over the parent finite elements by use of a computationally efficient Gaussian integration scheme.

Gould, Dana C.↗

Operational Dust Prediction

Over the last few years, numerical prediction of dust aerosol concentration has become prominent at several research and operational weather centres due to growing interest from diverse stakeholders, such as solar energy plant managers, health professionals, aviation and military authorities and policymakers. Dust prediction in numerical weather prediction-type models faces a number of challenges owing to the complexity of the system. At the centre of the problem is the vast range of scales required to fully account for all of the physical processes related to dust. Another limiting factor is the paucity of suitable dust observations available for model, evaluation and assimilation. This chapter discusses in detail numerical prediction of dust with examples from systems that are currently providing dust forecasts in near real-time or are part of international efforts to establish daily provision of dust forecasts based on multi-model ensembles. The various models are introduced and described along with an overview on the importance of dust prediction activities and a historical perspective. Assimilation and evaluation aspects in dust prediction are also discussed.

assimilation↗

Bayesian Estimation of Earth’s Undiscovered Mineralogical Diversity Using Noninformative Priors

Recently, statistical distributions have been explored to provide estimates of the mineralogical diversity of Earth, and Earth-like planets. In this paper, a Bayesian approach is introduced to estimate Earth’s undiscovered mineralogical diversity. Samples are generated from a posterior distribution of the model parameters using Markov chain Monte Carlo simulations such that estimates and inference are directly obtained. It was previously shown that the mineral species frequency distribution conforms to a generalized inverse Gauss–Poisson (GIGP) large number of rare events model. Even though the model fit was good, the population size estimate obtained by using this model was found to be unreasonably low by mineralogists. In this paper, several zero-truncated, mixed Poisson distributions are fitted and compared, where the Poisson-lognormal distribution is found to provide the best fit. Subsequently, the population size estimates obtained by Bayesian methods are compared to the empirical Bayes estimates. Species accumulation curves are constructed and employed to estimate the population size as a function of sampling size. Finally, the relative abundances, and hence the occurrence probabilities of species in a random sample, are calculated numerically for all mineral species in Earth’s crust using the Poisson-lognormal distribution. These calculations are connected and compared to the calculations obtained in a previous paper using the GIGP model for which mineralogical criteria of an Earth-like planet were given.

Bayesian statistics↗

Understanding Relationships Between Satellite, Model, and Ground-Based Surface Temperature Characterizations From Overcast to Clear Conditions in Support of Satellite Remote Sensing of Clouds and Radiation

Accurate and consistent global estimates of cloud coverage and their properties are fundamental to long-term Earth radiation budget (ERB) monitoring efforts like the Clouds and the Earth’s Radiant Energy System (CERES) project. Cloud detection algorithms often apply thresholding approaches to identify where clouds occur by comparing satellite-measured radiances with those that are expected under cloud-free conditions. In addition, once a cloud is detected, the derivation of cloud optical and microphysical properties also requires knowledge of the background radiances below the cloud. In the infrared, knowledge of the surface emissivity and the expected skin temperature under both cloudy and cloud-free conditions is needed. These traits are generally well known over the oceans. Over land, however, comparisons between satellite-derived land surface temperature (LST) with that characterized in numerical weather analyses reveal large differences in many parts of the world, often exceeding 5 K, which can lead to significant satellite cloud detection and cloud property retrieval errors. Furthermore, clouds have a dramatic influence on the LST, and therefore characterization of that model parameter also depends on the capability of the model to accurately resolve clouds. Thus, the LST characterized in models is, at times, a poor approximation for what would otherwise be observed, thereby impeding accurate satellite cloud retrievals. As a result, we seek to develop a more robust method for estimating the LST required for satellite cloud characterizations. This effort is accomplished through a combination of surface emission/air temperature relationship studies in all-sky conditions using ground measurement stations, along with deep neural network (DNN) estimates of expected LST under overcast and cloud-free conditions. We demonstrate that substituting DNN-predicted LST for that generated by numerical models can mitigate model-inherent diurnal dependencies and reduce overall bias and uncertainty relative to satellite/ground observations by 0.5–4 K and 0.5–2 K, respectively. It is expected that this work will lead to improved satellite cloud retrievals that enhance ERB monitoring efforts.

B Scarino↗

Anomalous and Ungrouped Carbonaceous Chondrites in the US Antarctic Meteorite Collection and their Potential Relevance to Ryugu and Bennu

With two different carbonaceous asteroid sample return missions in full swing, attention has focused on what connections can be made between the asteroid samples and the wide range of carbonaceous chondrite meteorites in worldwide collections. The US Antarctic meteorite collection contains nearly 1000 carbonaceous chondrites of various types including many in well-established groups as well as ungrouped, unusual or anomalous groups [1]. Some of the latter have been included in, or are possibly related to, recently proposed new carbonaceous chondrite classifications – CA and CY chondrites [2,3]. In addition to these, there are numerous ungrouped samples that have properties intermediate between established groups (like CM and CO; [4]), distinct from any other groups [5], or have anomalous properties that might be attributable to parent body processes such as heating, fluid interaction, or impacts [6]. Some CM anomalous or ungrouped chondrites share spectral features with Ryugu, which has a small hydration peak arguably due to hydrated minerals left after either impact heating or shock in carbonaceous chondrites [7]. PCA 91008, PCA 02012, GRO 95566, and LEW 85311 are all CMs that have experienced heating or metamorphism that may be due to impacts, solar radiation, or radiogenic decay [6]. These samples all have low H contents, C/H (bulk), and low 17O [6,8] and may hold clues to understanding the mineralogy of Ryugu, or interpretations of its spectral properties. WIS 91600, on the other hand, appears to be related to several other highly altered CM [6], and shares properties with the newly proposed CY chondrites [6]. An understanding of this grouplet will also aid in the interpretation of Bennu samples which have strong hydration features. In addition, ungrouped carbonaceous chondrites may provide valuable insights into the aqueous alteration potentially recorded in Bennu and Ryugu samples; such as the relatively moderate aqueous alteration recognized and dated at 4-5 myr after CAIs in MAC 88107 (C2-ungrouped) by [9] to the extensive aqueous alteration apparent in MIL 090292 (C1-ungrouped) [10]. Finally, LON 94101/94102 is a brecciated CM chondrite with numerous lithologies. Its appearance is similar to some of the brecciated lithologies visible at the surface of Bennu and Ryugu [11]. Although CM chondrites with multiple lithologies are common (e.g., [12]), the lithologies are often difficult to resolve at the hand specimen scale and only after some detailed e-beam characterization are the subtle lithologic differences evident (e.g., [13]). LON 94102 contains visibly distinct clasts at the hand specimen scale, relatively rare for CM breccias. The largest CM chondrites by mass from the U.S. Antarctic meteorite collection may provide insights into the extent of brecciation and heterogeneity within more typical CM chondrite-like source asteroids. For example, recent curation CT scans [14] show possible clasts within in a subsplit of ALH 83100 (CM1/2) which is the largest CM chondrite or CM chondrite pairing group in the U.S. Antarctic collection with an original mass of 3.019 kg. Similar work on additional subsplits or meteorites may aid our understanding of heterogeneity within CM chondrites from Antarctica. Initial classification of CM chondrites in U.S. Antarctic meteorite collection includes preliminary pairing when petrographically similar CM chondrites that have been previously found in the same field area. The two largest CM chondrite preliminary pairing groups by mass are the ALH 83102 (CM2) pairing group with an original mass of 2.554 kg and the EET 96005 (CM2) pairing group with an original mass of 1.125 kg. However, preliminary pairing groups assigned at classification—intrinsically— do not include stones that are petrographically distinct. Rigorous pairing group studies of CM chondrites from specific field sites are needed to investigate if there is unrecognized heterogeneity in CM chondrites from Antarctica that may represent common asteroid impactors (pre-atmospheric entry asteroids/meteoroids). Detailed pairing studies have the potential to recognize initial stones with multiple lithologies and investigate if there are stones that only sample one of those respective lithologies in the collection. These groups of heated CMs, extensively hydrated CMs, intermediate between CM and CO chondrites, and brecciated samples are all potentially relevant to Bennu and Ryugu samples, where heating, hydration, and brecciation have all come into play. These bodies might also be comprised of material intermediate to CM and CO chondrites, or at least distinct from the well-established CC groups. These small and unusual groups of carbonaceous chondrites may help to unlock new information about early solar system processes and aid in the understanding of the evolution of these carbonaceous asteroids.

meteorites↗

Estimation of agronomic variables using spectral signatures

Techniques for the determination of leaf area index or leaf angle distribution from remote-sensing canopy-reflectance (CR) measurements are developed on the basis of empirical models relating CR to parameters such as soil and vegetation spectral properties, solar flux, and viewing angle. A general procedure for inverting CR models is presented and applied to the models of Suits (1972), Verhoef and Bunnik (1981), and Norman (1979) in the IR range. Numerical results for a soybean canopy are compared in a table, and the error sensitivity of the inverted models is shown to be relatively high, requiring the use of ancillary data such as soil reflectance, leaf reflectance, and leaf transmittance.

Goel, N. S.↗

High temperature electronic excitation and ionization rates in gases

The relaxation times for electronic excitation due to electron bombardment of atoms was found to be quite short, so that electron kinetic temperature (T sub e) and the electron excitation temperature (T asterisk) should equilibrate quickly whenever electrons are present. However, once equilibrium has been achieved, further energy to the excited electronic states and to the kinetic energy of free electrons must be fed in by collisions with heavy particles that cause vibrational and electronic state transitions. The rate coefficients for excitation of electronic states produced by heavy particle collision have not been well known. However, a relatively simple semi-classical theory has been developed here which is analytic up to the final integration over a Boltzmann distribution of collision energies; this integral can then be evaluated numerically by quadrature. Once the rate coefficients have been determined, the relaxation of electronic excitation energy can be evaluated and compared with the relaxation rates of vibrational excitation. Then the relative importance of these two factors, electronic excitation and vibrational excitation by heavy particle collision, on the transfer of energy to free electron motion, can be assessed.

Hansen, Frederick↗

Counterrotating prop-fan simulations which feature a relative-motion multiblock grid decomposition enabling arbitrary time-steps

Improvements are presented of a computer algorithm developed for the time-accurate flow analysis of rotating machines. The flow model is a finite volume method utilizing a high-resolution approximate Riemann solver for interface flux definitions. The numerical scheme is a block LU implicit iterative-refinement method which possesses apparent unconditional stability. Multiblock composite gridding is used to orderly partition the field into a specified arrangement of blocks exhibiting varying degrees of similarity. Block-block relative motion is achieved using local grid distortion to reduce grid skewness and accommodate arbitrary time step selection. A general high-order numerical scheme is applied to satisfy the geometric conservation law. An even-blade-count counterrotating unducted fan configuration is chosen for a computational study comparing solutions resulting from altering parameters such as time step size and iteration count. The solutions are compared with measured data.

Janus, J. Mark↗

Verification and transfer of thermal pollution model. Volume 5: Verification of 2-dimensional numerical model

The six-volume report: describes the theory of a three dimensional (3-D) mathematical thermal discharge model and a related one dimensional (1-D) model, includes model verification at two sites, and provides a separate user's manual for each model. The 3-D model has two forms: free surface and rigid lid. The former, verified at Anclote Anchorate (FL), allows a free air/water interface and is suited for significant surface wave heights compared to mean water depth; e.g., estuaries and coastal regions. The latter, verified at Lake Keowee (SC), is suited for small surface wave heights compared to depth (e.g., natural or man-made inland lakes) because surface elevation has been removed as a parameter. These models allow computation of time dependent velocity and temperature fields for given initial conditions and time-varying boundary conditions.

Lee, S. S.↗

Verification and transfer of thermal pollution model. Volume 6: User's manual for 1-dimensional numerical model

The six-volume report: describes the theory of a three dimensional (3-D) mathematical thermal discharge model and a related one dimensional (1-D) model, includes model verification at two sites, and provides a separate user's manual for each model. The 3-D model has two forms: free surface and rigid lid. The former, verified at Anclote Anchorage (FL), allows a free air/water interface and is suited for significant surface wave heights compared to mean water depth; e.g., estuaries and coastal regions. The latter, verified at Lake Keowee (SC), is suited for small surface wave heights compared to depth (e.g., natural or man-made inland lakes) because surface elevation has been removed as a parameter.

Lee, S. S.↗

Numerical Experimentation with Maximum Likelihood Identification in Static Distributed Systems

Many important issues in the control of large space structures are intimately related to the fundamental problem of parameter identification. One might also ask how well this identification process can be carried out in the presence of noisy data since no sensor system is perfect. With these considerations in mind the algorithms herein are designed to treat both the case of uncertainties in the modeling and uncertainties in the data. The analytical aspects of maximum likelihood identification are considered in some detail in another paper. The questions relevant to the implementation of these schemes are dealt with, particularly as they apply to models of large space structures. The emphasis is on the influence of the infinite dimensional character of the problem on finite dimensional implementations of the algorithms. Those areas of current and future analysis are highlighted which indicate the interplay between error analysis and possible truncations of the state and parameter spaces.

Scheid, R. E., Jr.↗

Numerical and experimental evaluations of the flow past nested chevrons

An effort is made to contribute to the development of CFD by relating the successful use of vortex dynamics in the computation of the pressure drop past a planar array of chevron-shaped obstructions. An ensemble of results was used to compute the loss coefficient k, stimulating an experimental program for the assessment of the measured loss coefficient for the same geometry. The most provocative result of this study has been the representation of kinetic energy production in terms of vorticity source terms.

Foss, J. F.↗

System related testing and analysis of FRECOPA

Results from the French Cooperative Payload (FRECOPA) system analysis are presented. It was one of the numerous experiments which were flown on the Long Duration Exposure Facility (LDEF) satellite. In our flight configuration (LEO orbit, trailing edge), the environment was a better vacuum than the leading edge, with many thermal cycles (32000) and a large amount of UV radiation (11100 equivalent sun hours). Also, the satellite was mainly bombarded by micro-particles. It saw a low atomic flux and minor doses of protons and electrons.

Durin, Christian↗

A High-Resolution Capability for Large-Eddy Simulation of Jet Flows

A large-eddy simulation (LES) code that utilizes high-resolution numerical schemes is described and applied to a compressible jet flow. The code is written in a general manner such that the accuracy/resolution of the simulation can be selected by the user. Time discretization is performed using a family of low-dispersion Runge-Kutta schemes, selectable from first- to fourth-order. Spatial discretization is performed using central differencing schemes. Both standard schemes, second- to twelfth-order (3 to 13 point stencils) and Dispersion Relation Preserving schemes from 7 to 13 point stencils are available. The code is written in Fortran 90 and uses hybrid MPI/OpenMP parallelization. The code is applied to the simulation of a Mach 0.9 jet flow. Four-stage third-order Runge-Kutta time stepping and the 13 point DRP spatial discretization scheme of Bogey and Bailly are used. The high resolution numerics used allows for the use of relatively sparse grids. Three levels of grid resolution are examined, 3.5, 6.5, and 9.2 million points. Mean flow, first-order turbulent statistics and turbulent spectra are reported. Good agreement with experimental data for mean flow and first-order turbulent statistics is shown.

DeBonis, James R.↗

Supercomputing Aspects for Simulating Incompressible Flow

The primary objective of this research is to support the design of liquid rocket systems for the Advanced Space Transportation System. Since the space launch systems in the near future are likely to rely on liquid rocket engines, increasing the efficiency and reliability of the engine components is an important task. One of the major problems in the liquid rocket engine is to understand fluid dynamics of fuel and oxidizer flows from the fuel tank to plume. Understanding the flow through the entire turbo-pump geometry through numerical simulation will be of significant value toward design. One of the milestones of this effort is to develop, apply and demonstrate the capability and accuracy of 3D CFD methods as efficient design analysis tools on high performance computer platforms. The development of the Message Passage Interface (MPI) and Multi Level Parallel (MLP) versions of the INS3D code is currently underway. The serial version of INS3D code is a multidimensional incompressible Navier-Stokes solver based on overset grid technology, INS3D-MPI is based on the explicit massage-passing interface across processors and is primarily suited for distributed memory systems. INS3D-MLP is based on multi-level parallel method and is suitable for distributed-shared memory systems. For the entire turbo-pump simulations, moving boundary capability and efficient time-accurate integration methods are built in the flow solver, To handle the geometric complexity and moving boundary problems, an overset grid scheme is incorporated with the solver so that new connectivity data will be obtained at each time step. The Chimera overlapped grid scheme allows subdomains move relative to each other, and provides a great flexibility when the boundary movement creates large displacements. Two numerical procedures, one based on artificial compressibility method and the other pressure projection method, are outlined for obtaining time-accurate solutions of the incompressible Navier-Stokes equations. The performance of the two methods is compared by obtaining unsteady solutions for the evolution of twin vortices behind a flat plate. Calculated results are compared with experimental and other numerical results. For an unsteady flow, which requires small physical time step, the pressure projection method was found to be computationally efficient since it does not require any subiteration procedure. It was observed that the artificial compressibility method requires a fast convergence scheme at each physical time step in order to satisfy the incompressibility condition. This was obtained by using a GMRES-ILU(0) solver in present computations. When a line-relaxation scheme was used, the time accuracy was degraded and time-accurate computations became very expensive.

Kwak, Dochan↗