Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical modelling”

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 127 records · Page 7

Using Numerical Modeling to Simulate Space Capsule Ground Landings

Experimental work is being conducted at the National Aeronautics and Space Administration s (NASA) Langley Research Center (LaRC) to investigate ground landing capabilities of the Orion crew exploration vehicle (CEV). The Orion capsule is NASA s replacement for the Space Shuttle. The Orion capsule will service the International Space Station and be used for future space missions to the Moon and to Mars. To evaluate the feasibility of Orion ground landings, a series of capsule impact tests are being performed at the NASA Langley Landing and Impact Research Facility (LandIR). The experimental results derived at LandIR provide means to validate and calibrate nonlinear dynamic finite element models, which are also being developed during this study. Because of the high cost and time involvement intrinsic to full-scale testing, numerical simulations are favored over experimental work. Subsequent to a numerical model validated by actual test responses, impact simulations will be conducted to study multiple impact scenarios not practical to test. Twenty-one swing tests using the LandIR gantry were conducted during the June 07 through October 07 time period to evaluate the Orion s impact response. Results for two capsule initial pitch angles, 0deg and -15deg , along with their computer simulations using LS-DYNA are presented in this article. A soil-vehicle friction coefficient of 0.45 was determined by comparing the test stopping distance with computer simulations. In addition, soil modeling accuracy is presented by comparing vertical penetrometer impact tests with computer simulations for the soil model used during the swing tests.

Heymsfield, Ernie↗

Analytical and Numerical Modeling of Sensor Port Acoustics

The purpose of this Technical Publication (TP) is to provide background, relevant theory, and examples for acoustic response analysis of a sensor port. A new approach is devised theoretically and computationally that captures the true acoustic response of a sensor port. This TP summarizes the acoustics back-ground, the port response theoretical development, and provides comparisons of a port acoustic response using an analytical model and computational acoustics. The effects of nonlinear acoustics and acoustic propagation in liquids is also examined. Additionally, this TP describes the design of a specialized filter using the predicted sensor port response that can be applied to data for correction.

sensor port↗

Numerical Modeling of the Transient Chilldown Process of a Cryogenic Propellant Transfer Line

Before cryogenic fuel depots can be fully realized, efficient methods with which to chill down the spacecraft transfer line and receiver tank are required. This paper presents numerical modeling of the chilldown of a liquid hydrogen tank-to-tank propellant transfer line using the Generalized Fluid System Simulation Program (GFSSP). To compare with data from recently concluded turbulent LH2 chill down experiments, seven different cases were run across a range of inlet liquid temperatures and mass flow rates. Both trickle and pulse chill down methods were simulated. The GFSSP model qualitatively matches external skin mounted temperature readings, but large differences are shown between measured and predicted internal stream temperatures. Discrepancies are attributed to the simplified model correlation used to compute two-phase flow boiling heat transfer. Flow visualization from testing shows that the initial bottoming out of skin mounted sensors corresponds to annular flow, but that considerable time is required for the stream sensor to achieve steady state as the system moves through annular, churn, and bubbly flow. The GFSSP model does adequately well in tracking trends in the data but further work is needed to refine the two-phase flow modeling to better match observed test data.

Chilldown↗

Numerical modeling of enclosure convection

A numerical study on the steady and unsteady natural convection in two-dimensional rectangular enclosures has been performed by a time-accurate ADI finite difference scheme. The study covered a range of Rayleigh numbers between 1000 and 10 to the 7th, aspect ratios between 0.2 and 10.0, and tilt angles between -90 (heating from bottom) and +90 deg (heating from top). Various Prandtl numbers have been studied, but only the results of water (Pr = 7.0) are reported here due to space limitations. The physics revealed, however, includes the convection phenomena and the Rayleigh-Benard stability, as well as the combined mechanism of these two. The onset of secondary cells is determined by using a velocity map, which is simpler and cleaner, instead of a streamline plot. The critical Ra number for the occurrence of these secondary cells is shown to be lower than can be detected by experimental studies. On the Rayleigh-Benard stability part, a second transition from stable single-cell convection to periodic multicellular convection is disclosed.

Duh, J. C.↗

Numerical model of long-lived Jovian vortices

The extension of the measured zonal velocity profile into the adiabatic interior of Jupiter, while eddies and large oval structures are confined to a shallow stably-stratified upper layer, are assumed in a nonlinear numerical model of long-lived Jovian vortices. In agreement of the observed flows of Jupiter, each vortex is stationary with respect to the shear flow at a critical latitude that is close to the latitude of the vortex center. The solutions obtained are strongly nonlinear, in contrast to the solitary wave solutions that are the weakly nonlinear extensions of ultralong linear waves. The merging of two stable vortices upon collision, rather than the non-interaction predicted by solitary wave theory, is in keeping with Jovian vortex observations. It is suggested that long-lived vortices maintain themselves against dissipation by absorbing smaller vortices produced by convection.

Ingersoll, A. P.↗

Final results of an experiment in operational forecasting of sea breeze thunderstorms using a mesoscale numerical model

Sea breeze thunderstorms during quiescent synoptic conductions account for 40 percent of Florida rainfall, and are the dominant feature of April-October weather at the Kennedy Space Center (KSC). An effort is presently made to assess the feasibility of a mesoscale numerical model in improving the point-specific thunderstorm forecasting accuracy at the KSC, in the 2-12 hour time frame. Attention is given to the Applied Regional Atmospheric Modeling System.

Lyons, Walter A.↗

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.↗

Numerical modeling of troposphere-induced gravity wave propagation

Sources of internal gravity waves (IGW) in the upper atmosphere are assumed to be meteorological processes in the troposphere. These sources are vertically and horizontally inhomogeneous and time dependent. In order to describe the IGW propagation from such sources, a numerical solution of a system of hydrodynamical equations is required. In addition, it is necessary to take into account the influence of the altitude latitude inhomogeneity of the temperature and wind fields on the IGW propagation as well as the processes of dissipation. An algorithm is proposed for numerical modelling of the IGW propagation over a limited area from tropospheric local sources to the upper atmosphere. The algorithm takes into account all the above features. A spectral grid method is used with the expansion of wave fields into the Fourier series over longitude. The upper limit conditions were obtained from the requirement of a limited energy dissipation rate in an atmospheric column. The no slip (zero velocity) condition was used at the Earth's surface.

Gavrilov, N. M.↗

Numerical Modeling of Thermofluid Transients During Chilldown of Cryogenic Transfer Lines

The chilldown of fluid transfer lines is an important part of using cryogenic systems such as those found in both ground and space based applications. The chilldown process is a complex combination of both thermal and fluid transient phenomena. A cryogenic liquid flows through a transfer line that is initially at a much higher temperature than the cryogen. Transient heat transfer processes between the liquid and transfer line cause vaporization of the liquid, and this phase change can cause transient pressure and flow surges in the liquid. As the transfer line is cooled, these effects diminish until the liquid reaches a steady flow condition in the chilled transfer line. If these transient phenomena are not properly accounted for in the design process of a cryogenic system, it can lead to damage or failure of system components during operation. For such cases, analytical modeling is desirable for ensuring that a cryogenic system transfer line design is adequate for handling the effects of a chilldown process. The purpose of this paper is to present the results of a numerical model developed using Generalized Fluid System Simulation Program (GFSSP)'s new fluid transient capability in combination with its previously developed thermal transient capability to predict pressure and flow surge in cryogenic transfer lines during a chilldown process. An experiment performed by the National Bureau of Standards (NBS) in 1966 has been chosen as the baseline comparison case for this work. NBS s experimental set-up consisted of a 10.59 cubic foot supply dewar, an inlet valve, and a 200 foot long, in Outside Diameter (OD) vacuum jacketed copper transfer line that exhausted to atmosphere. Three different inlet valves, an in-port ball valve, a 1-in-port globe valve and a 1-in-port gate valve, were used in NBS's experiments. Experiments were performed using both liquid hydrogen and liquid nitrogen as the fluids. The proposed paper will include detailed comparisons of GFSSP's predictions with NBS's experimental results.

Majumdar, Alok↗

An Open-Source Numerical Model for Mitigating Refractory Alloy Hot Cracking Susceptibility

Refractory alloys are susceptible to solidification cracking during welding and 3D printing. Composition control is an effective method of controlling solidification cracking. This work evaluates the effect of compositional variation in refractory metal systems on a computed solidification cracking susceptibility. A numerical model has been developed using Python code and open-source CALPHAD software to calculate Kou’s crack susceptibility index. The model is validated against past weldability studies performed on several refractory alloy systems. The approach is extended towards the development of new alloys with improved 3D printability and weldability and is shown to have utility in defining compositional limits for existing alloys and feedstocks. Furthermore, the model will aid in determining process controls for powder reuse and recycling.

pycalphad↗

Effects of the pressure perturbation field in numerical models of unidirectionally sheared thunderstorm convection - Two versus three dimensions

The physical roles of 'buoyant' and 'dynamic' pressure components, and the distinction between buoyant and hydrostatic pressure perturbations, are aspects of the pressure perturbation field in strongly sheared convective storms studied by means of two- and three-dimensional anelastic numerical modeling experiments with common environmental profiles. The pressure analysis clarifies the differences between two- and three-dimensional storms. In the main updraft, strong midlevel thermal buoyancy is partly opposed by a downward-perturbed vertical pressure gradient force. This, however, occurs to a much greater extent in two dimensions than in three, contributing to smaller net upward accelerations. While the buoyant and hydrostatic fields are intimately related to the total buoyancy distribution, the buoyant pressure perturbation is smoother and of lower amplitude than its hydrostatic counterpart. For the model experiments, this distinction is far greater in three dimensions than in two, in association with the smaller scale of the active convection in three dimensions.

Schlesinger, R. E.↗

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 Modeling of Extensional Necking Instabilities: Application to Ganymede's Grooved Terrain

Ganymede s pervasive 5-10 km-wavelength grooves have been suggested to result from a necking instability during an epoch of lithospheric extension, but to date few quantitative studies of groove formation have been performed. We present two-dimensional numerical models of necking instabilities under conditions that are appropriate to Ganymede at the time of groove formation. Preliminary simulations indicate that extensional necking instabilities can occur under a range of conditions, many of which may be relevant to Ganymede. The form of the surface topography produced by these instabilities varies as a function of the strain rate, amount of extension, initial topographic perturbation, and rheological parameters.

Bland, M. T.↗

Numerical model of solar dynamic radiator for parametric analysis

Growth power requirements for Space Station Freedom will be met through addition of 25 kW solar dynamic (SD) power modules. Extensive thermal and power cycle modeling capabilities have been developed which are powerful tools in Station design and analysis, but which prove cumbersome and costly for simple component preliminary design studies. In order to aid in refining the SD radiator to the mature design stage, a simple and flexible numerical model was developed. The model simulates heat transfer and fluid flow performance of the radiator and calculates area mass and impact survivability for many combinations of flow tube and panel configurations, fluid and material properties, and environmental and cycle variations.

Rhatigan, Jennifer L.↗

Optimum employment of satellite indirect soundings as numerical model input

The characteristics of satellite-derived temperature soundings that would significantly affect their use as input for numerical weather prediction models were examined. Independent evaluations of satellite soundings were emphasized to better define error characteristics. Results of a Nimbus-6 sounding study reveal an underestimation of the strength of synoptic scale troughs and ridges, and associated gradients in isobaric height and temperature fields. The most significant errors occurred near the Earth's surface and the tropopause. Soundings from the TIROS-N and NOAA-6 satellites were also evaluated. Results again showed an underestimation of upper level trough amplitudes leading to weaker thermal gradient depictions in satellite-only fields. These errors show a definite correlation to the synoptic flow patterns. In a satellite-only analysis used to initialize a numerical model forecast, it was found that these synoptically correlated errors were retained in the forecast sequence.

Horn, L. H.↗

Transient Numerical Modeling of Catalytic Channels

This paper presents a transient model of catalytic combustion suitable for isolated channels and monolith reactors. The model is a lumped two-phase (gas and solid) model where the gas phase is quasi-steady relative to the transient solid. Axial diffusion is neglected in the gas phase; lateral diffusion, however, is accounted for using transfer coefficients. The solid phase includes axial heat conduction and external heat loss due to convection and radiation. The combustion process utilizes detailed gas and surface reaction models. The gas-phase model becomes a system of stiff ordinary differential equations while the solid phase reduces, after discretization, into a system of stiff ordinary differential-algebraic equations. The time evolution of the system came from alternating integrations of the quasi-steady gas and transient solid. This work outlines the numerical model and presents some sensitivity studies on important parameters including internal transfer coefficients, catalytic surface site density, and external heat-loss (if applicable). The model is compared to two experiments using CO fuel: (1) steady-state conversion through an isothermal platinum (Pt) tube and (2) transient propagation of a catalytic reaction inside a small Pt tube. The model requires internal mass-transfer resistance to match the experiments at lower residence times. Under mass-transport limited conditions, the model reasonably predicted exit conversion using global mass-transfer coefficients. Near light-off, the model results did not match the experiment precisely even after adjustment of mass-transfer coefficients. Agreement improved for the first case after adjusting the surface kinetics such that the net rate of CO adsorption increased compared to O2. The CO / O2 surface mechanism came from a sub-set of reactions in a popular CH4 / O2 mechanism. For the second case, predictions improved for lean conditions with increased external heat loss or adjustment of the kinetics as in the first case. Finally, the results show that different initial surface-species distribution leads to different steady-states under certain conditions. These results demonstrate the utility of a lumped two-phase model of a transient catalytic combustor with detailed chemistry.

Struk, Peter M.↗

Development and Implementation of Large-Scale Numerical Models for Shape Memory Mars Spring Tires

The current work focuses on implementation of the user-defined shape memory alloy (SMA) model in the finite element analysis program ABAQUS for large-scale simulations of the Spring Tires made of SMA developed at the NASA Glenn Research Center. The main objective of this study was to improve and optimize the structural design of the SMA tires through in-depth numerical analysis and sensitivity studies. Various design variables (wire diam., coil diam., pitch, spring length, and bead angle) were varied to study their influence on the global load-displacement response of the tire construct. A detailed investigation of the three-dimensional stress states was also carried out to enhance understanding of the local changes as the tire goes through global deformation. It was concluded that a robust numerical model with a good predictive capability, together with a thoughtfully crafted sensitivity study, can result in reduced design iterations required to reach a desired tire performance. This will, in turn, lead to significantly reduced manufacturing time, required labor, and testing expense.

shape memory alloys↗

Regenerator performance in a Vuilleumier refrigerator compared with a third-order numerical model

A three-stage Vuilleumier refrigerator was used to measure the performance of various third stage regenerators. The refrigerator operates between 2.5 and 5.0 Hz and, depending on the material used in the third stage regenerator, achieves temperatures of 8 to 20 K at the cold end of the third stage. This paper presents a comparison of regenerator performance for four regenerator materials: 229 micron diameter spheres of Pb(+)5 pct Sb, 229 micron diameter spheres of brass, 216 micron irregularly-shaped GdRh powder, and a mixture of 229 microns and 762 microns diameter spheres of Pb(+5) pct Sb. The experimental results are compared with a first-order model that neglects the void volume within the regenerator and with a third-order model that considers the effect of pressure oscillations in the regenerator void volume. Experimental results indicate that regenerator losses are dominated by the pressure oscillation in the void volume rather than the mass flow through the temperature gradient in the regenerator. These results are consistent with the third-order numerical model. This model shows that the heat capacity of the gas in the void space as well as the heat capacity of the matrix influences the regenerator performance.

Bradley, P. E.↗