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 847 records · Page 47

Experimental facilities for system identification

Future space vehicles will differ significantly from the space systems used in the past. The planned spacecraft configurations will include extremely large structures, up to 100-200 meters across. Because the allowable launch weights are limited, the large space structures must be constructed of lightweight, flexible elements, and active control of the shape and attitude of the spacecraft will be required. The behavior of large structures is characterized by many closely spaced natural modes, and some applications may also include large on-board disturbances. Consequently, the control and disturbance forces will invariably spill over to a large number of modes. Structural identification will be necessary for precision pointing and shape controls to be effective. An accurate mathematical model of a structure is essential for the success of precision control system design. However, the currently available analytical modeling codes, such as NASTRAN, are incapable of producing numerical models of the required precision. The favored approach for model determination is to refine the mathematical model based on experiments in the orbital environment. It is to this effect that we must develop techniques to reliably generate the accurate models required for future missions, and the arena in which these methods will be developed and proved is through ground-based system identification experiments and demonstrations. The Air Force Astronautics Laboratory (AFAL) has expanded the on-site ground-test facilities in recent years, and additional sites are planned for the immediate future. We shall review the plans for laboratory growth at AFAL with regard to the type of experiments proposed and the availability of the new facilities.

Das, Alok↗

On the relationship between tectonic plates and thermal mantle plume morphology

Models incorporating plate-like behavior, i.e., near uniform surface velocity and deformation concentrated at plate boundaries, into a convective system, heated by a mix of internal and basal heating and allowing for temperature dependent viscosity, were constructed and compared to similar models not possessing plate-like behavior. The simplified numerical models are used to explore how plate-like behavior in a convective system can effect the lower boundary layer from which thermal plumes form. A principal conclusion is that plate-like behavior can significantly increase the temperature drop across the lower thermal boundary layer. This temperature drop affects the morphology of plumes by determining the viscosity drop across the boundary layer. Model results suggest that plumes on planets possessing plate-like behavior, e.g., the Earth, may differ in morphologic type from plumes on planets not possessing plate-like behavior, e.g., Venus and Mars.

Lenardic, A.↗

Massively Parallel Assimilation of TOGA/TAO and Topex/Poseidon Measurements into a Quasi Isopycnal Ocean General Circulation Model Using an Ensemble Kalman Filter

A massively parallel ensemble Kalman filter (EnKF)is used to assimilate temperature data from the TOGA/TAO array and altimetry from TOPEX/POSEIDON into a Pacific basin version of the NASA Seasonal to Interannual Prediction Project (NSIPP)ls quasi-isopycnal ocean general circulation model. The EnKF is an approximate Kalman filter in which the error-covariance propagation step is modeled by the integration of multiple instances of a numerical model. An estimate of the true error covariances is then inferred from the distribution of the ensemble of model state vectors. This inplementation of the filter takes advantage of the inherent parallelism in the EnKF algorithm by running all the model instances concurrently. The Kalman filter update step also occurs in parallel by having each processor process the observations that occur in the region of physical space for which it is responsible. The massively parallel data assimilation system is validated by withholding some of the data and then quantifying the extent to which the withheld information can be inferred from the assimilation of the remaining data. The distributions of the forecast and analysis error covariances predicted by the ENKF are also examined.

Keppenne, Christian L.↗

Evaluation of Skin Temperatures Retrieved from GOES-8

Skin temperatures derived from geostationary satellites have the potential of providing the temporal and spatial resolution needed for model assimilation. To adequately assess the potential improvements in numerical model forecasts that can be made by assimilating satellite data, an estimate of the accuracy of the skin temperature product is necessary. A particular skin temperature algorithm, the Physical Split Window Technique, that uses the longwave infrared channels of the GOES Imager has shown promise in recent model assimilation studies to provide land surface temperatures with reasonable accuracy. A comparison of retrieved GOES-8 skin temperatures from this algorithm with in situ measurements is presented. Various retrieval algorithm issues are addressed including surface emissivity

Suggs, Ronnie, J.↗

The End of Comet Shoemaker-Levy 9

The pieces of the former Comet Shoemaker-Levy 9 struck Jupiter during the week of July 16-22 1994. We divide the general impact phenomena into four phases: (1) Atmospheric entry should have been accompanied by an entry flash, probably only detectable by the Galileo spacecraft. (2) Aerodynamic forces cause the impactor to explode. We show that the apparent altitude of the explosions is consistent with the analytical model proposed by Zahnle and Chyba et al, with numerical models calculated by Mac Low and Zahnle and Yabe et al., and with dynamically-based estimates of fragment masses Suggested by Scotti and Melosh and Asphaug and Benz. (3) The explosion produces a fireball that propagates most quickly up the wake, as predicted by Takata et al. and Crawford et al. The fireball's luminosity fades over a timescale of 30-60 seconds. (4) Each explosion produced a plume that reached thousands of kilometers above the cloud-tops. The plumes appear to have been dusty and SO illuminated when they reached sunlight. The plumes rose and fell over a period of ten or twenty minutes, governed by the simple laws of ballistics. When plume material (dust) fell back on the atmosphere it was reheated and glowed brightly in the thermal infrared for a period of order ten to twenty minutes; a large fraction of the impact energy was spent thus. Geometry introduces additional factors that complicate interpretation of the observed light curves.

Zahnle, K. J.↗

Combustion Of Metals In Reduced Gravity And Extraterrestrial Environments

The recent focus of this research project has been to model the combustion of isolated metal droplets and, in particular, to couple the existing theories and formulations of phenomena such as condensation, reaction kinetics, radiation, and surface reactions to formulate a more complete combustion model. A fully transient, one-dimensional (spherical symmetry) numerical model that uses detailed chemical kinetics, multi-component molecular transport mechanisms, condensation kinetics, and gas phase radiation heat transfer was developed. A coagulation model was used to simulate the particulate formation of MgO. The model was used to simulate the combustion of an Mg droplet in pure O2 and CO2. Methanol droplet combustion is considered as a test case for the solution method for both quasi-steady and fully transient simulations. Although some important processes unique to methanol combustion, such as water absorption at the surface, are not included in the model, the results are in sufficient agreement with the published data. Since the major part of the heat released in combustion of Mg, and in combustion of metals in general, is due to the condensation of the metal oxide, it is very important to capture the condensation processes correctly. Using the modified nucleation theory, an Arrhenius type rate expression is derived to calculate the condensation rate of MgO. This expression can be easily included in the CHEMKIN reaction mechanism format. Although very little property data is available for MgO, the condensation rate expression derived using the existing data is able to capture the condensation of MgO. An appropriate choice of the reference temperature to calculate the rate coefficients allows the model to correctly predict the subsequent heat release and hence the flame temperature.

Abbud-Madrid, A.↗

Laboratory for Atmospheres 2008 Technical Highlights

The 2008 Technical Highlights describes the efforts of all members of the Laboratory for Atmospheres. Their dedication to advancing Earth Science through conducting research, developing and running models, designing instruments, managing projects, running field campaigns, and numerous other activities, is highlighted in this report. The Laboratory for Atmospheres (Code 613) is part of the Earth Sciences Division (Code 610), formerly the Earth Sun Exploration Division, under the Sciences and Exploration Directorate (Code 600) based at NASA s Goddard Space Flight Center in Greenbelt, Maryland. In line with NASA s Exploration Initiative, the Laboratory executes a comprehensive research and technology development program dedicated to advancing knowledge and understanding of the atmospheres of Earth and other planets. The research program is aimed at understanding the influence of solar variability on the Earth s climate; predicting the weather and climate of Earth; understanding the structure, dynamics, and radiative properties of precipitation, clouds, and aerosols; understanding atmospheric chemistry, especially the role of natural and anthropogenic trace species on the ozone balance in the stratosphere and the troposphere; and advancing our understanding of physical properties of Earth s atmosphere. The research program identifies problems and requirements for atmospheric observations via satellite missions. Laboratory scientists conceive, design, develop, and implement ultraviolet, infrared, optical, radar, laser, and lidar technology for remote sensing of the atmosphere. Laboratory members conduct field measurements for satellite data calibration and validation, and carry out numerous modeling activities. These modeling activities include climate model simulations, modeling the chemistry and transport of trace species on regional-to-global scales, cloud-resolving models, and development of next-generation Earth system models. Interdisciplinary research is carried out in collaboration with other laboratories and research groups within the Earth Sciences Division, across the Sciences and Exploration Directorate, and with partners in universities and other Government agencies. The Laboratory for Atmospheres is a vital participant in NASA s research agenda. Our Laboratory often has relatively large programs, sizable satellite missions, and observational campaigns that require the cooperative and collaborative efforts of many scientists. We ensure an appropriate balance between our scientists responsibility for these large collaborative projects and their need for an active individual research agenda. This balance allows members of the Laboratory to continuously improve their scientific credentials. Members of the Laboratory interact with the general public to support a wide range of interests in the atmospheric sciences. Among other activities, the Laboratory raises the public s awareness of atmospheric science by presenting public lectures and demonstrations, by making scientific data available to wide audiences, by teaching, and by mentoring students and teachers. The Laboratory makes substantial efforts to attract new scientists to the various areas of atmospheric research. We strongly encourage the establishment of partnerships with Federal and state agencies that have operational responsibilities to promote the societal application of our science products. This report describes our role in NASA s mission, gives a broad description of our research, and summarizes our scientists major accomplishments during calendar year 2008. The report also contains useful information on human resources, scientific interactions, and outreach activities.

Cote, Charles E.↗

Predictive Modeling of Fast-Curing Thermosets in Nozzle-Based Extrusion

This work presents an approach to modeling the dynamic spreading and curing behavior of thermosets in nozzle-based extrusions. Thermosets cover a wide range of materials, some of which permit low-temperature processing with subsequent high-temperature and high-strength working properties. Extruding thermosets may overcome the limited working temperatures and strengths of conventional thermoplastic materials used in additive manufacturing. This project aims to produce technology for the fabrication of thermoset-based structures leveraging advances made in nozzle-based extrusion, such as fused deposition modeling (FDM), material jetting, and direct writing. Understanding the synergistic interactions between spreading and fast curing of extruded thermosetting materials will provide essential insights for applications that require accurate dimensional controls, such as additive manufacturing [1], [2] and centrifugal coating/forming [3]. Two types of thermally curing thermosets -- one being a soft silicone (Ecoflex 0050) and the other being a toughened epoxy (G/Flex) -- served as the test materials in this work to obtain models for cure kinetics and viscosity. The developed models align with extensive measurements made with differential scanning calorimetry (DSC) and rheology. DSC monitors the change in the heat of reaction, which reflects the rate and degree of cure at different crosslinking stages. Rheology measures the change in complex viscosity, shear moduli, yield stress, and other properties dictated by chemical composition. By combining DSC and rheological measurements, it is possible to establish a set of models profiling the cure kinetics and chemorheology without prior knowledge of chemical composition, which is usually necessary for sophisticated mechanistic modeling. In this work, we conducted both isothermal and dynamic measurements with both DSC and rheology. With the developed models, numerical simulations yielded predictions of diameter and height of droplets, along with width and height of extruded lines cured at varied temperatures. Experimental results carried out on a goniometric platform and a nozzle-based 3D printer showed agreement with the numerical simulations. Finally, this presentation will show how the models are adaptable to the planning of tool paths and designs in additive manufacturing.

Xie, Jingjin↗

Modeling the near-field ionospheric disturbances during earthquakes

We present the latest development of the Wave Perturbation-Global Ionosphere-Thermosphere Model (WP-GITM), a threedimensional physics-based numerical model for seismic/tsunami-ionosphere coupling via atmospheric acoustic-gravity waves. WP-GITM was previously applied to simulate the ionospheric perturbations resulted from the epicentral crustal movement by assuming spherical acoustic-gravity waves originated from a point source, which was specified by the seismic measurement at a single location nearby the epicenter. In this work, we extend WP-GITM to include the effects of Rayleigh surface waves and adapt WP-GITM to utilize seismic measurements from more than one location in the attempt to capture the radiation pattern of the seismic source. We apply the new WP-GITM to model the near-field co-seismic ionospheric perturbations during the 16 September 2015 Illapel earthquake. The comparison between the simulated ionospheric total electron content perturbations and the GPS observations shows promising results.

Ravanelli, Michela↗

A radiation model for calculating atmospheric corrections to remotely sensed infrared measurements, version 2

A numerical model is developed which calculates the atmospheric corrections to infrared radiometric measurements due to absorption and emission by water vapor, carbon dioxide, and ozone. The corrections due to aerosols are not accounted for. The transmissions functions for water vapor, carbon dioxide, and water are given. The model requires as input the vertical distribution of temperature and water vapor as determined by a standard radiosonde. The vertical distribution of carbon dioxide is assumed to be constant. The vertical distribution of ozone is an average of observed values. The model also requires as input the spectral response function of the radiometer and the nadir angle at which the measurements were made. A listing of the FORTRAN program is given with details for its use and examples of input and output listings. Calculations for four model atmospheres are presented.

Boudreau, R. D.↗

A simple model for the interaction between vertical eddy heat fluxes and static stability

A numerical model for studying the interaction of vertical eddy heat fluxes and vertical temperature structure in midlatitude regions is described. The temperature profile for the model was derived from calculations of the equilibrium among heating rates in simplified representations of large-scale vertical eddy heat flux, moist convection and radiation. An eddy flux profile is calculated based on the quasi-geostrophic, liner baroclinic instability of a single wave. Model equilibrium states for summer and winter conditions are compared with observations, and the results are discussed in detail.

Gutowski, W. J., Jr.↗

A Model of Thermal Conductivity for Planetary Soils: Theory for Cemented Soils - 2

A numerical model of heat conduction through particulate media made of spherical grains cemented by various bonding agents is presented. The pore-filling gas conductivity, volume fraction, and thermal conductivity of the cementing phase are tunable parameters. Cement fractions <0.001-0.01% in volume have small effects on the soil bulk thermal conductivity. A significant conductivity increase (factor 3-8) is observed for bond fractions of 0.01 to 1% in volume. In the 1 to 15% bond fraction domain, the conductivity increases continuously but less intensely (25-100% conductivity increase compared to a 1% bond system). Beyond 15% of cements, the conductivity increases vigorously and the bulk conductivity rapidly approaches that of bedrock. The composition of the cements (i.e. conductivity) has little influence on the bulk thermal inertia of the soil, especially if the volume of bond <10%. These results indicate that temperature measurements are sufficient to detect cemented soils and quantify the amount of cementing phase, but the mineralogical nature of the bonds and the typical grain size are unlikely to be determined from orbit. On Mars, a widespread surface unit characterized by a medium albedo (0.19-0.26) and medium/high thermal inertia (200-600 J s(0.5)/sq m/K) has long been hypothesized to be associated with a duricrust. The fraction of cement required to fit the thermal data is less than approx.1-5% by volume. This small amount of material is consistent with orbital observations, confirming that soil cementation is an important factor controlling the thermal inertia of the Martian surface

Piqueux, S.↗

Modeling and Observations of Outlet Canyons from Lake Overflow Floods on Early Mars

Numerous observations from both orbital remote sensing [1-3] and Mars Curiosity [4] suggest that lakes were once part of the martian landscape. From orbital data, one of the key lines of evidence for past paleolakes is the existence of several hundred valley network-fed basins – usually craters – that have outlet valleys that remain perched above their floors. The existence of outlets requires that water ponded to the point that it overflowed confining topography. Beyond recognizing these landforms, there has been only limited work reconstructing the morphometry, formative hydrology, and incision history for these outlets. Here, we describe our recently published observations of outlets and ongoing numerical modeling looking at these factors.

Fassett, Caleb I.↗

A numerical study of the NOL large-scale test

The large-scale gap test was modeled numerically using the reactive hydrodynamics code 2DE with forest fire burn rates. The model gives good agreement between calculated and experimental critical gap values. The calibration of peak pressure in the gap versus gap length is obtained from these calculations, and is in good agreement with the published experiment curve. The two dimensional nature of the gap test is evident from the curvature of the shock wave in the gap and test sample, and from the observed distance of run to detonation in the test samples.

Bowan, A. L.↗

Finite Element Modeling of Magnetically-Damped Convection during Solidification

A fully 3-D, transient finite element model is developed to represent the magnetic damping effects on complex fluid flow, heat transfer and electromagnetic field distributions in a Sn- 35.5%Pb melt undergoing unidirectional solidification. The model is developed based on our in- house finite element code for the fluid flow, heat transfer and electromagnetic field calculations. The numerical model is tested against numerical and experimental results for water as reported in literature. Various numerical simulations are carried out for the melt convection and temperature distribution with and without the presence of a transverse magnetic field. Numerical results show that magnetic damping can be effectively applied to stabilize melt flow, reduce turbulence and flow levels in the melt and over a certain threshold value a higher magnetic field resulted in a greater reduction in velocity. Also, for the study of melt flow instability, a long enough running time is needed to ensure the final fluid flow recirculation pattern. Moreover, numerical results suggest that there seems to exist a threshold value of applied magnetic field, above which magnetic damping becomes possible and below which the 0 convection in the melt is actually enhanced.

deGroh, H. C.↗

Seasonal-scale Observational Data Analysis and Atmospheric Phenomenology for the Cold Land Processes Experiment

The Cold Land Processes Experiment (CLPX) experiment emphasized the development of a strong synergism between process-oriented understanding, land surface models and microwave remote sensing. Our work sought to investigate which topographically- generated atmospheric phenomena are most relevant to the CLPX MSA's for the purpose of evaluating their climatic importance to net local moisture fluxes and snow transport through the use of high-resolution data assimilation/atmospheric numerical modeling techniques. Our task was to create three long-term, scientific quality atmospheric datasets for quantitative analysis (for all CLPX researchers) and provide a summary of the meteorologically-relevant phenomena of the three MSAs (see Figure) over northern Colorado. Our efforts required the ingest of a variety of CLPX datasets and the execution an atmospheric and land surface data assimilation system based on the Navier-Stokes equations (the Local Analysis and Prediction System, LAPS, and an atmospheric numerical weather prediction model, as required) at topographically- relevant grid spacing (approx. 500 m). The resulting dataset will be analyzed by the CLPX community as a part of their larger research goals to determine the relative influence of various atmospheric phenomena on processes relevant to CLPX scientific goals.

Poulos, Gregory S.↗

Subduction in an Eddy-Resolving State Estimate of the Northeast Atlantic Ocean

Are eddies an important contributor to subduction in the eastern subtropical gyre? Here, an adjoint model is used to combine a regional, eddy-resolving numerical model with observations to produce a state estimate of the ocean circulation. The estimate is a synthesis of a variety of in- situ observations from the Subduction Experiment, TOPEX/POSEIDON altimetry, and the MTI General Circulation Model. The adjoint method is successful because the Northeast Atlantic Ocean is only weakly nonlinear. The state estimate provides a physically-interpretable, eddy-resolving information source to diagnose subduction. Estimates of eddy subduction for the eastern subtropical gyre of the North Atlantic are larger than previously calculated from parameterizations in coarse-resolution models. Furthermore, eddy subduction rates have typical magnitudes of 15% of the total subduction rate. Eddies contribute as much as 1 Sverdrup to water-mass transformation, and hence subduction, in the North Equatorial Current and the Azores Current. The findings of this thesis imply that the inability to resolve or accurately parameterize eddy subduction in climate models would lead to an accumulation of error in the structure of the main thermocline, even in the relatively-quiescent eastern subtropical gyre.

WATER MASS TRANSFORMATIONS↗

Sensitivity of a climatologically-driven sea ice model to the ocean heat flux

Ocean heat flux sensitivity was studied on a numerical model of sea ice covering the Weddell Sea region of the southern ocean. The model is driven by mean monthly climatological atmospheric variables. For each model run, the ocean heat flux is uniform in both space and time. Ocean heat fluxes below 20 W m to the minus 2 power do not provide sufficient energy to allow the ice to melt to its summertime thicknesses and concentrations by the end of the 14 month simulation, whereas ocean heat fluxes of 30 W m to the minus 2 power and above result in too much ice melt, producing the almost total disappearance of ice in the Weddell Sea by the end of the 14 months. These results are dependent on the atmospheric forcing fields.

Parkinson, C. L.↗