Search NASA⌕ Search

SEARCH · Search NASA

Results for “Casting Modeling Simulation”

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 73 records · Page 4

Bridging Current Indirect and Future Direct Estimates of Convective Mass Flux Toward Quantifying the Role of Convective Cores on Anvil Area Time Tendencies

Anvil cloud spatial extent is strongly driven by convective mass flux. However, the lack of convective mass flux observations precludes global analyses of the functional relationship between anvil cloud area growth tendencies and convective areal coverage and vertical mass flux. Partly informed by analyses of a hierarchy of WRF mesoscale convective system (MCS) simulations, in this presentation, we discuss new work on estimating convective mass flux by way of an approach that casts the problem of quantifying convective mass flux into terms related to diabatic heating vertical structure, convective area and tropospheric stability. We demonstrate the usefulness of this approach at the MCS spatial scale via development of a simple, robust model that analytically connects system-scale convective mass flux and anvil cloud area growth rates. We show that the model predicts anvil area size changes in agreement with those observed in a global convective system tracking database (TOOCAN). An additional advantage of this approach for estimating convective mass flux is that it bridges future satellite mission plans for directly observing convective mass flux (e.g., AOS) with longer-record observations made by TRMM, GPM and satellite sounder mission data (the latter used for estimating stability) that, when combined, can be used for inferring convective mass fluxes at the MCS spatial scale spanning slowly changing, diverse environments. An evaluation and bridging of both direct and indirect approaches supports convective process analyses varying over longer time periods in this “era of convection.”

convective mass flux↗

Bridging current indirect and future direct estimates of convective mass flux toward quantifying the role of convective cores on anvil area time tendencies

Anvil cloud spatial extent is strongly driven by convective mass flux. However, the lack of convective mass flux observations precludes global analyses of the functional relationship between anvil cloud area growth tendencies and convective areal coverage and vertical mass flux. Partly informed by analyses of a hierarchy of WRF mesoscale convective system (MCS) simulations, in this presentation, we discuss new work on estimating convective mass flux by way of an approach that casts the problem of quantifying convective mass flux into terms related to diabatic heating vertical structure, convective area and tropospheric stability. We demonstrate the usefulness of this approach at the MCS spatial scale via development of a simple, robust model that analytically connects system-scale convective mass flux and anvil cloud area growth rates. We show that the model predicts anvil area size changes in agreement with those observed in a global convective system tracking database (TOOCAN). An additional advantage of this approach for estimating convective mass flux is that it bridges future satellite mission plans for directly observing convective mass flux (e.g., AOS) with longer-record observations made by TRMM, GPM and satellite sounder mission data (the latter used for estimating stability) that, when combined, can be used for inferring convective mass fluxes at the MCS spatial scale spanning slowly changing, diverse environments. An evaluation and bridging of both direct and indirect approaches supports convective process analyses varying over longer time periods in this “era of convection.”

convective mass flux↗

Plasma Wake Simulations and Object Charging in a Shadowed Lunar Crater During a Solar Storm

Within a permanently shadowed lunar crater the horizontal flow of solar wind is obstructed by upstream topography, forming a plasma wake that electrostatically diverts ions toward the crater floor and generates a surface potential that can reach kilovolts. In the present work kinetic plasma simulations are employed to investigate the morphology of a lunar crater wake during passage of a solar storm. Results are cast in terms of leading dimensionless ratios including the ion Mach number, ratio of crater depth to plasma Debye length, peak secondary electron yield, and electron temperature vs. electron impact energy at peak secondary yield. This small set of ratios allows generalization to a much wider range of scenarios. The kinetic simulation results are fed forward into an equivalent-circuit model of a roving astronaut. In very low-plasma-current environments triboelectric charging of the astronaut suit becomes effectively perpetual, representing a critical engineering concern for roving within shadowed lunar regions. Finally, simulated ion fluxes are used to explore sputtering and implantation processes within an idealized crater. It is suggested that the physics of plasma mini-wakes formed in the vicinity of permanently shadowed topography may play a critical role in modulating the enigmatic spatial distribution of volatiles at the lunar poles.

Zimmerman, Michael I.↗

Constitutive modeling for isotropic materials

The first year of progress on a NASA-Lewis contract with the General Electric Co is documented. The purpose of this contract (NAS3-23927) is to develop and evaluate unified constitutive equations for applications to hot-path components of aircraft gas turbine engines such as high pressure turbine blades and vanes. To accomplish this goal, uniaxial, notched, and multiaxial specimens made of conventionally cast Rene 80 are being tested under conditions that simulate engine operating conditions. To reduce the raw data, automated data reduction techniques are being developed that produce computer files containing the information needed to analyze proposed constitutive theories. Described are the analytical methods being developed to determine the parameters for these nonlinear unified theories by using the reduced data files. In another activity, a dedicated finite-element computer code is being developed to use unified theories in the structural analysis of hot-section components. This code was extensively verified for one such theory by successfully predicting the strain histories measured experimentally at the notch root of complex specimens taken from complex laboratory specimens.

Ramaswamy, V. G.↗

Probabilistic analysis for fatigue strength degradation of materials

This report presents the results of the first year of a research program conducted for NASA-LeRC by the University of Texas at San Antonio. The research included development of methodology that provides a probabilistic treatment of lifetime prediction of structural components of aerospace propulsion systems subjected to fatigue. Material strength degradation models, based on primitive variables, include both a fatigue strength reduction model and a fatigue crack growth model. Linear elastic fracture mechanics is utilized in the latter model. Probabilistic analysis is based on simulation, and both maximum entropy and maximum penalized likelihood methods are used for the generation of probability density functions. The resulting constitutive relationships are included in several computer programs, RANDOM2, RANDOM3, and RANDOM4. These programs determine the random lifetime of an engine component, in mechanical load cycles, to reach a critical fatigue strength or crack size. The material considered was a cast nickel base superalloy, one typical of those used in the Space Shuttle Main Engine.

Royce, Lola↗

Shoemaker-Levy 9 and the tidal disruption of comets

The break-up of Periodic Comet Shoemaker-Levy 9 into multiple pieces following its grazing encounter with Jupiter in July 1992 can be used to study tidally-induced fracture in comets. This spectacular event allows us not only to set limits on the size, strength and density of Shoemaker-Levy 9 itself, but provides invaluable guidance to numerical modeling of such encounters. In an extensive treatment of tidal breakup which assumed self-gravitating, homogeneous, perfectly elastic bodies, Dobrovolskis derived simple analytical expressions for the tidally-induced surface and central stresses. Both can be cast in such a way that Poisson's ratio is the only material dependent constant entering these expressions. Whether both surface and central failure must be initiated as a criterion for breakup, or either one of them is sufficient, remains a subject of disagreement. To resolve this debate, we model the details of cometary breakup using a three-dimensional Smooth Particle Hydrodynamics (SPH) code modified to simulate fracture in small solid objects. At lower stresses associated with brittle failure, we use a rate-dependent strength based on the nucleation of incipient flaws whose number density is given by a Weibull distribution.

Benz, W.↗

Current-driven plasma acceleration versus current-driven energy dissipation. III - Anomalous transport

In the present paper the linear stability description and weak turbulence theory are used to develop a second order description of wave-particle transport and anomalous dissipation. The goal is to arrive at anomalous transport coefficients that can be readily included in fluid flow codes. In particular, expressions are derived for the heating rates of ions and electrons by the unstable waves and for the electron-wave momentum exchange rate that controls the anomalous resistivity effect. Comparative calculations were undertaken assuming four different saturation models: ion trapping, electron trapping, ion resonance broadening, and thermodynamic bound. A foremost finding is the importance of the role of electron Hall parameter in scaling the level of anomalous dissipation for the parameter range of the MPD thruster plasma. Polynomial expressions of the relevant transport coefficients cast solely in terms of macroscopic parameters are also obtained for inclusion in plasma fluid codes for the self-consistent numerical simulation of real thruster flows including microturbulent effects.

Choueiri, Edgar Y.↗

Expendable launch vehicle studies

Analytical support studies of expendable launch vehicles concentrate on the stability of the dynamics during launch especially during or near the region of maximum dynamic pressure. The in-plane dynamic equations of a generic launch vehicle with multiple flexible bending and fuel sloshing modes are developed and linearized. The information from LeRC about the grids, masses, and modes is incorporated into the model. The eigenvalues of the plant are analyzed for several modeling factors: utilizing diagonal mass matrix, uniform beam assumption, inclusion of aerodynamics, and the interaction between the aerodynamics and the flexible bending motion. Preliminary PID, LQR, and LQG control designs with sensor and actuator dynamics for this system and simulations are also conducted. The initial analysis for comparison of PD (proportional-derivative) and full state feedback LQR Linear quadratic regulator) shows that the split weighted LQR controller has better performance than that of the PD. In order to meet both the performance and robustness requirements, the H(sub infinity) robust controller for the expendable launch vehicle is developed. The simulation indicates that both the performance and robustness of the H(sub infinity) controller are better than that for the PID and LQG controllers. The modelling and analysis support studies team has continued development of methodology, using eigensensitivity analysis, to solve three classes of discrete eigenvalue equations. In the first class, the matrix elements are non-linear functions of the eigenvector. All non-linear periodic motion can be cast in this form. Here the eigenvector is comprised of the coefficients of complete basis functions spanning the response space and the eigenvalue is the frequency. The second class of eigenvalue problems studied is the quadratic eigenvalue problem. Solutions for linear viscously damped structures or viscoelastic structures can be reduced to this form. Particular attention is paid to Maxwell and Kelvin models. The third class of problems consists of linear eigenvalue problems in which the elements of the mass and stiffness matrices are stochastic. dynamic structural response for which the parameters are given by probabilistic distribution functions, rather than deterministic values, can be cast in this form. Solutions for several problems in each class will be presented.

Bainum, Peter M.↗

Heat transfer distributions around nominal ice accretion shapes formed on a cylinder in the NASA Lewis icing research tunnel

Local heat transfer coefficients were obtained on irregular cylindrical shapes which typify the accretion of ice on circular cylinders in cross flow. The ice shapes were grown on a 5.1 cm (2.0 in.) diameter cylinder in the NASA Lewis Icing Research Tunnel. The shapes were 2, 5, and 15 min accumulations of glaze ice and 15 min accumulation of rime ice. Heat transfer coefficients were also measured around the cylinder with no ice accretion. These icing shapes were averaged axially to obtain a nominal shape of constant cross section for the heat transfer tests. Heat transfer coefficients around the perimeter of each shape were measured with electrically heated copper strips embedded in the surface of the model which was cast from polyurethane foam. Each strip contained a thermocouple to measure the local surface temperature. The models were run in a 15.2 x 68.6 cm (6 x 27 in.) wind tunnel at several velocities. Background turbulence in the wind tunnel was less than 0.5 percent. The models were also run with a turbulence producing grid which gave about 3.5 percent turbulence at the model location with the model removed. The effect of roughness was also simulated with sand grains glued to the surface. Results are presented as Nusselt number versus angle from the stagnation line for the smooth and rough models for both high and low levels of free stream turblence. Roughness of the surface in the region prior to flow separation plays a major role in determining the heat transfer distribution.

Vanfossen, G. J.↗

A three-dimensional incompressible flow simulation method and its application to the Space Shuttle main engine. II Turbulent flow

An implicit finite difference code cast in general curvilinear coordinates is further developed for three-dimensional incompressible turbulent flows. The code is based on the method of pseudocompressibility and utilizes the Beam and Warming implicit approximate factorization algorithm to achieve computational efficiency. A multiple-zone method is further extended to include composite-grids to overcome the excessive computer memory required for solving turbulent flows in complex three-dimensional geometries. A simple turbulence model is proposed for internal flows. The code is being used for the Space Shuttle Main Engine (SSME) internal flow analyses.

Chang, J. L. C.↗

Some Navier-Stokes calculations for the CAST-10 airfoil

Results of the simulation of the viscous flow past the CAST-10 airfoil were shown for different flow conditions. Since the experiments provide only surface pressures and force coefficients, the comparison to the numerical results relies on these. Good agreement of these results is found for the lower Mach number cases except for the shock position. As numerical experiments indicate, this seems to be due to the turbulent shock boundary layer interaction which is not correctly modelled by the algebraic turbulence model employed. For the lower Mach number case the influence of the transition location was investigated, too. Changing the transition location at the lower surface has much more influence on the pressure distribution than changing it on the upper side. For the higher Mach number case the double shock structure found in the experiment for the lower Reynolds numbers was not reproduced by the numerical solutions. The reason for this is unknown though it may be due to the turbulence modelling. For the higher Reynolds number a better resolution of the boundary layer is needed in the computation in order to recover the experimental pressure plateau; but then the shock position is still found downstream of the experimental one.

Schwamborn, D.↗

Feedback control for unsteady flow and its application to the stochastic Burgers equation

The study applies mathematical methods of control theory to the problem of control of fluid flow with the long-range objective of developing effective methods for the control of turbulent flows. Model problems are employed through the formalism and language of control theory to present the procedure of how to cast the problem of controlling turbulence into a problem in optimal control theory. Methods of calculus of variations through the adjoint state and gradient algorithms are used to present a suboptimal control and feedback procedure for stationary and time-dependent problems. Two types of controls are investigated: distributed and boundary controls. Several cases of both controls are numerically simulated to investigate the performances of the control algorithm. Most cases considered show significant reductions of the costs to be minimized. The dependence of the control algorithm on the time-descretization method is discussed.

Choi, Haecheon↗

Effect of mild atherosclerosis on flow resistance in a coronary artery casting of man

An in-vitro flow study was conducted in a mildly atherosclerotic main coronary artery casting of man using sugar-water solutions simulating blood viscosity. Steady flow results indicated substantial increases in pressure drop, and thus flow resistance at the same Reynolds number, above those for Poiseuille flow by 30 to 100 percent in the physiological Reynolds number range from about 100 to 400. Time-averaged pulsatile flow data showed additional 5 percent increases in flow resistance above the steady flow results. Both pulsatile and steady flow data from the casting were found to be nearly equal to those from a straight, axisymmetric model of the casting up to a Reynolds number of about 200, above which the flow resistance of the casting became gradually larger than the corresponding values from the axisymmetric model.

Back, L. H.↗

Method of Modeling and Simulation of Shaped External Occulters

The present invention relates to modeling an external occulter including: providing at least one processor executing program code to implement a simulation system, the program code including: providing an external occulter having a plurality of petals, the occulter being coupled to a telescope; and propagating light from the occulter to a telescope aperture of the telescope by scalar Fresnel propagation, by: obtaining an incident field strength at a predetermined wavelength at an occulter surface; obtaining a field propagation from the occulter to the telescope aperture using a Fresnel integral; modeling a celestial object at differing field angles by shifting a location of a shadow cast by the occulter on the telescope aperture; calculating an intensity of the occulter shadow on the telescope aperture; and applying a telescope aperture mask to a field of the occulter shadow, and propagating the light to a focal plane of the telescope via FFT techniques.

Lyon, Richard G.↗

A Study of the Role of Clouds in the Relationship Between Land Use/Land Cover and the Climate and Air Quality of the Atlanta Area

The goal of Project ATLANTA is to derive a better scientific understanding of how land cover changes associated with urbanization affect local and regional climate and air quality. Clouds play a significant role in this relationship. Using GOES images, we found that in a 63-day period (5 July-5 September 1996) there were zero days which were clear for the entire daylight period. Days which are cloud-free in the morning become partly cloudy with small cumulus clouds in the afternoon in response to solar heating. This result casts doubt on the applicability of California-style air quality models which run in perpetual clear skies. Days which are clear in the morning have higher ozone than those which are cloudy in the morning. Using the RAMS model, we found that urbanization increases the skin surface temperature by about 1.0-1.5 C on average under cloudy conditions, with an extreme of +3.5 C. Clouds cool the surface due to their shading effect by 1.5-2.0 C on average, with an extreme of 5.0 C. RAMS simulates well the building stage of the cumulus cloud field, but does poorly in the decaying phase. Next year's work: doing a detailed cloud climatology and developing improved RAMS cloud simulations.

Kidder, Stanley Q.↗

Transonic aeroelastic analysis of launch vehicle configurations

A numerical study of the aeroelastic stability of typical launch vehicle configurations in transonic flight is performed. Recent computational fluid dynamics techniques are used to simulate the transonic aerodynamic flow fields, as opposed to relying on experimental data for the unsteady aerodynamic pressures. The flow solver is coupled to an appropriate structural representation of the vehicle. The aerodynamic formulation is based on the thin layer approximation to the Reynolds-Averaged Navier-Stokes equations, where the account for turbulent mixing is done by the two-layer Baldwin and Lomax algebraic eddy viscosity model. The structural-dynamic equations are developed considering free-free flexural vibration of an elongated beam with variable properties and are cast in modal form. Aeroelastic analyses are performed by integrating simultaneously in the two sets of equations. By tracing the growth or decay of a perturbed oscillation, the aeroelastic stability of a given constant configuration can be ascertained. The method is described in detail, and results that indicate its application are presented. Applications include some validation cases for the algorithm developed, as well as the study of configurations known to have presented flutter programs in the past.

Filgueirasdeazevedo, Joao Luiz↗

Approximate Co-Location of Precipitation and Low-Level Westerlies in Tropical Monthly Means

In summer monsoon regions the monthly mean precipitation regions coincide approximately well with regions of westerlies at low-levels. An included chart shows a 15-year (1979-1993) mean August 850 hPa zonal wind from NCEP/NCAR reanalysis dataset and Xie-Arkin precipitation. It shows a region of westerlies covering most of Northern Indian Ocean and extending to northwestern Pacific. This region coincides well with the region with precipitation greater than 6 mm/day. Obviously the coincidence is not exact; the region of larges; zonal wind in the Arabian Sea is in a region of relatively low precipitation and is far from the region of maximum precipitation in Bay of Bengal. Also, in a zonally averaged sense between 40E and 140E, the latitude of maximum precipitation is slightly higher than that of the maximum zonal wind. Low-level westerlies are also found in regions west of Central America and in western Africa north of the equator. These regions are also closely associated with precipitation centers. Across equator from these westerlies regions there are regions of strong easterlies. Also, on their poleward side the westerly regions are flanked by weaker easterly regions. In February, similar observation can be found in the Australian monsoon area and in South America monsoon region; again the regions of westerlies coincide well with regions of maximum precipitation. As in the northern hemisphere, the maximum precipitation is found to the cast of the maximum zonal wind. The two maxima lie almost at the same latitude with that of the westerlies slightly closer to the equator. In the non-monsoon seasons the low- level westerlies can also be found in the tropical precipitation regions, the longitudinal range of the westerlies is undiminished and the speed of the westerlies is not much weaker than that found in February. The interpretation of these observational facts is the goal of this investigation. The approach taken is numerical simulation with the Goddard Earth Observation Systems atmospheric general circulation model, and its aqua-planet version, combined with theoretical arguments.

Choa, Winston C.↗

Multiple Damage Progression Paths in Model-Based Prognostics

Model-based prognostics approaches employ domain knowledge about a system, its components, and how they fail through the use of physics-based models. Component wear is driven by several different degradation phenomena, each resulting in their own damage progression path, overlapping to contribute to the overall degradation of the component. We develop a model-based prognostics methodology using particle filters, in which the problem of characterizing multiple damage progression paths is cast as a joint state-parameter estimation problem. The estimate is represented as a probability distribution, allowing the prediction of end of life and remaining useful life within a probabilistic framework that supports uncertainty management. We also develop a novel variance control mechanism that maintains an uncertainty bound around the hidden parameters to limit the amount of estimation uncertainty and, consequently, reduce prediction uncertainty. We construct a detailed physics-based model of a centrifugal pump, to which we apply our model-based prognostics algorithms. We illustrate the operation of the prognostic solution with a number of simulation-based experiments and demonstrate the performance of the chosen approach when multiple damage mechanisms are active

Daigle, Matthew↗