Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Validated Numerical Models for the Convective Extinction of Fuel Droplets (CEFD)

The NASA Nebraska Space Grant (NSGC) & EPSCoR programs have continued their effort to support outstanding research endeavors by funding the Numerical Simulation of the Combustion of Fuel Droplets study at the University of Nebraska at Lincoln (UNL). This team of researchers has developed a transient numerical model to study the combustion of suspended and moving droplets. The engines that propel missiles, jets, and many other devices are dependent upon combustion. Therefore, data concerning the combustion of fuel droplets is of immediate relevance to aviation and aeronautical personnel, especially those involved in flight operations. The experiments being conducted by Dr. Gogos and Dr. Nayagam s research teams, allow investigators to gather data for comparison with theoretical predictions of burning rates, flame structures, and extinction conditions. The consequent improved fundamental understanding of droplet combustion may contribute to the clean and safe utilization of fossil fuels (Williams, Dryer, Haggard & Nayagam, 1997, f 2). The present state of knowledge on convective extinction of fuel droplets derives from experiments conducted under normal gravity conditions. However, any data obtained with suspended droplets under normal gravity are grossly affected by gravity. The need to obtain experimental data under microgravity conditions is therefore well justified and addresses one of the goals of NASA's Human Exploration and Development of Space (HEDS) microgravity combustion experiment.

Gogos, George↗
NASA NTRS1992

A general numerical model for wave rotor analysis

Wave rotors represent one of the promising technologies for achieving very high core temperatures and pressures in future gas turbine engines. Their operation depends upon unsteady gas dynamics and as such, their analysis is quite difficult. This report describes a numerical model which has been developed to perform such an analysis. Following a brief introduction, a summary of the wave rotor concept is given. The governing equations are then presented, along with a summary of the assumptions used to obtain them. Next, the numerical integration technique is described. This is an explicit finite volume technique based on the method of Roe. The discussion then focuses on the implementation of appropriate boundary conditions. Following this, some results are presented which first compare the numerical approximation to the governing differential equations and then compare the overall model to an actual wave rotor experiment. Finally, some concluding remarks are presented concerning the limitations of the simplifying assumptions and areas where the model may be improved.

Paxson, Daniel W.↗
NASA NTRS1982

A preliminary study of the application of HCMM satellite data to define initial and boundary conditions for numerical models: A case study in St. Louis, Missouri

Infrared and visible HCMM data were used to examine the potential application of these data to define initial and boundary conditions for mesoscale numerical models. Various boundary layer models were used to calculate the distribution of the surface heat flux, specific humidity depression (the difference between the specific humidity in the air at approxmately the 10 m level and the specific humidity at the ground), and the eddy vicosity in a 72 km by 72 km area centered about St. Louis, Missouri. Various aspects of the implications of the results on the meteorology of St. Louis are discussed. Overall, the results indicated that a reasonable estimate of the surface heat flux, urban albedo, ground temperature, and specific humidity depression can be obtained using HCMM satellite data. Values of the ground-specific humidity can be obtained if the distribution of the air-specific humidity is available. More research is required in estimating the absolute magnitude of the specific humidity depression because calculations may be sensitive to model parameters.

Vukovich, F. M.↗
NASA NTRS1992

Preconditioning for stratospheric sudden warmings - Sensitivity studies with a numerical model

Results are presented of numerical experiments performed to determine what constitutes preconditioning for stratospheric sudden warmings. A number of 20-d runs were performed with a 3D primitive equation model. Runs varied in that they had different (zonally symmetric) initial fields and different lower boundary planetary wave fields, both based on NMC data for particular days. Initial and boundary data were used from the two minor and the sole major sudden warmings of January-February 1979 and some minor and the sole major sudden warmings of January-February 1989. To isolate the properties of the initial conditions from preconditioned cases and those with no warming, hybrid initial fields were created by combining data from two different days, one preconditioned and the other not. Runs with the hybrid initial conditions indicate that it is the lower stratospheric winds that are important in characterizing a preconditioned flow. It is suggested that it would be possible for a preconditioned state to exist for at least a short time without a warming occurring.

Smith, Anne K.↗
NASA NTRS1994

Transient Numerical Modeling of the Combustion of Bi-Component Liquid Droplets: Methanol/Water Mixture

This study shows that liquid mixtures of methanol and water are attractive candidates for microgravity droplet combustion experiments and associated numerical modeling. The gas phase chemistry for these droplet mixtures is conceptually simple, well understood and substantially validated. In addition, the thermodynamic and transport properties of the liquid mixture have also been well characterized. Furthermore, the results obtained in this study predict that the extinction of these droplets may be observable in ground-based drop to tower experiments. Such experiments will be conducted shortly followed by space-based experiments utilizing the NASA FSDC and DCE experiments.

Marchese, A. J.↗
NASA NTRS2009

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↗
NASA NTRS2021

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↗
DOE OSTI2020

A Digital Twin Approach to Study Additive Manufacturing Processing Using Embedded Optical Fiber Sensors and Numerical Modeling

One of the major challenges for metal-powder-based additive manufacturing is measuring and mitigating residual strain induced during the manufacturing processes. This article reports distributed fiber optic sensors embedded in Inconel alloy components as experimental means to validate numerical models of additive manufacturing process. Electroplating was used to deposit a metal protective jacket onto standard telecom single-mode fibers for strain measurements, Fiber sensors were embedded in an Inconel alloy substrate using the laser engineered net shaping (LENS) process. Here by using a Rayleigh-scattering optical frequency domain reflectometer (OFDR), temperature changes, and residual strain in the metal substrate were monitored with 5 mm spatial resolution during the LENS process. Using finite element analysis, temperature and strain profiles induced by the LENS deposition processes were also numerically studied. Discrepancies between the simulated temperature and strain profiles and those measured directly were less than 10%. Results presented in this article demonstrates a digital twin approach to fuse modeling results with distributed fiber sensor measurement data to study additive manufacturing process toward design and fabrication process optimization.

36 MATERIALS SCIENCE↗
NASA NTRS2015

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↗
NASA NTRS1989

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

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

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

OC6 Project Phase IV: Validation of Numerical Models for Novel Floating Offshore Wind Support Structures

The fourth phase of the Offshore Code Comparison, Collaboration, Continued, with Correlation and unCertainty (OC6) aims to assess the ability of the tools to model a novel floating wind design. The design that was selected for this project is the TetraSpar designed by Stiesdal Offshore Technologies. The platform is a unique design that features a spar like stability with a very low center of gravity, but does not require deep-water ports. This combination is achieved with a system made of two separate steel structures, a hull at the water’s surface that supports the wind turbine tower, and a keel that is suspended below the main body. This lower keel has a very high density and drives the total platform’s low center of gravity. Scaled model tests of the TetraSpar design were performed in 2018-2019 by the University of Maine. The tests were performed with 1:43 Froude scaling. The data collected in the model tests provide a detailed description of the response of the platform in various controlled and relatively well understood conditions. Comparison with these data provides a good opportunity for validation of numerical models. The OC6 Phase IV work will focus on the specific parameters used in these scaled model tests. The load cases (LC) considered in this project are as follows: LC 1.X - Equilibrium and Static Offset Tests LC 2.X - Platform Free Decays LC 3.X - Turbine Thrust (Wind Only) LC 4.X - Regular and Irregular Waves (Wave only) LC 5.X - Combined Wind and Wave Loading

17 WIND ENERGY↗
NASA NTRS1986

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

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

Numerical Modeling of Air-Blast Suppression as a Function of Explosive-Charge Burial Depth

As a chemical explosion is buried, the mechanism for acoustic wave generation transitions from fully gas-generated at the surface to completely spall-induced at full containment depth. The fully gas-generated and completely spall-induced signals in the acoustic waveform are well described; however, the transition between these two end-members eludes numerical modeling because of the complex phenomena that are involved. The phenomena of crater formation and explosive cloud evolution are simulated using an Eulerian hydrocode that incorporates geomaterials with strength and porosity. Having accurately modeled these phenomena, we can confidently predict the propagation and relative strength of the gas-generated and spall-induced pulses in the recorded acoustic waveform. The numerical predictions agree with observations from the historical Stagecoach experiment as well as modern recordings from the Source Physics Experiment. In particular, the peak pressure p generated by an explosion is initially due to the gas-generated mechanism and decays with scaled depth of burial d s (depth d scaled by the cube-root of explosive yield w 1/3 ) as exp(-d s ) but then transitions near a scaled depth of 6 m/ton 1/3 to the spall-generated mechanism in which the decay is d$_{s}^{-7/4}$. This decay form is related to the strong ground-motion attenuation relationship that affects spall strength. So these results can improve seismoacoustic inverse models for the explosive source that need to account for the gas-generated and spall-induced signals and their effect on peak pressures and other acoustic signal features.

58 GEOSCIENCES↗
NASA NTRS2003

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↗
DOE OSTI2025

Resolution requirements for numerical modeling of neutrino quantum kinetics

Neutrino quantum kinetics is a rapidly evolving field in computational astrophysics, with a primary focus on collective neutrino oscillations in core-collapse supernovae and postmerger phases of binary neutron star mergers. In recent years, there has been considerable debate concerning resolution dependence in numerical simulations. In this paper, we conduct a comprehensive resolution study in both angular- and spatial directions by using two independent schemes of quantum kinetic neutrino transport: finite volume and pseudospectral methods. We complement our discussion by linear stability analysis including inhomogeneous modes. Our result suggests that decreasing spatial resolutions underestimates the growth of flavor instability, and then leads to wrong asymptotic states of flavor conversions, which potentially has a critical impact on astrophysical consequences. We further delve into numerical results of low resolution simulations, that reveals the underlying mechanism responsible for numerical artifacts caused by insufficient resolutions. Finally, this study settles the debate on requirements of resolutions and serves as a guideline for numerical modeling of quantum kinetic neutrino transport.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗