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 883 records · Page 49

Experience with a vectorized general circulation weather model on Star-100

A version of an atmospheric general circulation model was vectorized to run on a CDC STAR 100. The numerical model was coded and run in two different vector languages, CDC and LRLTRAN. A factor of 10 speed improvement over an IBM 360/95 was realized. Efficient use of the STAR machine required some redesigning of algorithms and logic. This precludes the application of vectorizing compilers on the original scalar code to achieve the same results. Vector languages permit a more natural and efficient formulation for such numerical codes.

Soll, D. B.↗

Can Machine Learning Reveal Precursors of Reversals of the Geomagnetic Axial Dipole Field?

It is well known that the axial dipole part of Earth’s magnetic field reverses polarity, so that the magnetic North Pole becomes the South Pole and vice versa. The timing of reversals is well documented for the past 160 Myr, but the conditions that lead to a reversal are still not well understood. It is not known if there are reliable ‘precursors’ of reversals (events that indicate that a reversal is upcoming) or what they might be. We investigate if machine learning (ML) techniques can reliably identify precursors of reversals based on time-series of the axial magnetic dipole field. The basic idea is to train a classifier using segments of time-series of the axial magnetic dipole. This training step requires modification of standard ML techniques to account for the fact that we are interested in rare events—a reversal is unusual, while a non-reversing field is the norm. Without our tweak, the ML classifiers lead to useless predictions. Perhaps even more importantly, the usable observational record is limited to 0–2 Ma and contains only five reversals, necessitating that we determine if the data are even sufficient to reliably train and validate an ML algorithm. To answer these questions we use several ML classifiers (linear/non-linear support vector machines and long short-term memory networks), invoke a hierarchy of numerical models (from simplified models to 3-D geodynamo simulations), and two palaeomagnetic reconstructions (PADM2M and Sint-2000). The performance of the ML classifiers varies across the models and the observational record and we provide evidence that this is not an artefact of the numerics, but rather reflects how ‘predictable’ a model or observational record is. Studying models of Earth’s magnetic field via ML classifiers thus can help with identifying shortcomings or advantages of the various models. For Earth’s magnetic field, we conclude that the ability of ML to identify precursors of reversals is limited, largely due to the small amount and low frequency resolution of data, which makes training and subsequent validation nearly impossible. Put simply: the ML techniques we tried are not currently capable of reliably identifying an axial dipole moment (ADM) precursor for geomagnetic reversals. This does not necessarily imply that such a precursor does not exist, and improvements in temporal resolution and length of ADM records may well offer better prospects in the future.

Kyle Gwirtz↗

A Numerical Climate Observing Network Design Study

This project was concerned with three related questions of an optimal design of a climate observing system: 1. The spatial sampling characteristics required from an ARGO system. 2. The degree to which surface observations from ARGO can be used to calibrate and test satellite remote sensing observations of sea surface salinity (SSS) as it is anticipated now. 3. The more general design of an climate observing system as it is required in the near future for CLIVAR in the Atlantic. An important question in implementing an observing system is that of the sampling density required to observe climate-related variations in the ocean. For that purpose this project was concerned with the sampling requirements for the ARGO float system, but investigated also other elements of a climate observing system. As part of this project we studied the horizontal and vertical sampling characteristics of a global ARGO system which is required to make it fully complementary to altimeter data with the goal to capture climate related variations on large spatial scales (less thanAttachment: 1000 km). We addressed this question in the framework of a numerical model study in the North Atlantic with an 1/6 horizontal resolution. The advantage of a numerical design study is the knowledge of the full model state. Sampled by a synthetic float array, model results will therefore allow to test and improve existing deployment strategies with the goal to make the system as optimal and cost-efficient as possible. Attachment: "Optimal observations for variational data assimilation".

Stammer, Detlef↗

Analysis of the early spectra and light curve of SN 1987A

Numerical modeling of supernova spectra, light curves, and hydrodynamics requires physical inputs, numerical techniques, approximations, and assumptions which must be thoroughly understood in order to study the details of supernova explosions. Here, we discuss some of these in the context of the early evolution of supernova 1987A. Gray radiation-hydrodynamics is used to calculate the bolometric light curve and the hydrodynamic evolution of the supernova. Synthetic spectra are then obtained for the resulting density and velocity structure. The spectrum calculations are performed using a special-relativistic treatment of the radiative transfer equation in the comoving frame, line blanketing by about 10(exp 5) spectral lines, and departures from local thermodynamic equilibrium (LTE) for H I, He I, Mg II, and Ca II. We find that we are able to simultaneously fit the early light curve and spectra reasonably well, using a progenitor model from Arnett (1991a), without fine-tuning the free parameters. Temperature structures and radiative equilibrium, non-LTE effects, homologous expansion, and mean opacities are discussed.

Hauschildt, Peter H.↗

Confinement of hot, hard X-ray producing electrons in solar flares

Possible thermal models for solar, hard X-ray emission, consisting of small volumes in which the electrons are rapidly heated to 400 million K, are examined to determine under what conditions such models can be more efficient than nonthermal models. The analysis is restricted to one dimension by supposing that the hot electrons are confined in a flux tube and that conduction and convection act only along the magnetic field. A conduction front, thin in relation to other length scales involved, separates the hot electrons from the relatively cold plasma. The range of conduction-front velocity is studied by systematizing the present knowledge of heat conduction and applying it first to a simple model in which convection is neglected, and then to a detailed numerical model in which convection is included. Convection in which mass motion goes to zero at the conduction front does not lead to a change in the conduction-front velocity compared to conduction alone, but does constitute another nonradiative energy loss which reduces the source efficiency.

Smith, D. F.↗

A general computer model for predicting the performance of gas sorption refrigerators

Projected performance requirements for cryogenic spacecraft sensor cooling systems which demand higher reliability and longer lifetimes are outlined. The gas/solid sorption refrigerator is viewed as a potential solution to cryogenic cooling needs. A software model of an entire gas sorption refrigerator system was developed. The numerical model, evaluates almost any combination and order of refrigerator components and any sorbent-sorbate pair or which the sorption isotherm data are available. Parametric curves for predicting system performance were generated for two types of refrigerators, a LaNi5-H2 absorption cooler and a Charcoal-N2 adsorption cooler. It is found that precooling temperature and heat exchanger effectiveness affect the refrigerator performance. It is indicated that gas sorption refrigerators are feasible for a number of space applications.

Sigurdson, K. B.↗

A General-applications Direct Global Matrix Algorithm for Rapid Seismo-acoustic Wavefield Computations

A new matrix method for rapid wave propagation modeling in generalized stratified media, which has recently been applied to numerical simulations in diverse areas of underwater acoustics, solid earth seismology, and nondestructive ultrasonic scattering is explained and illustrated. A portion of recent efforts jointly undertaken at NATOSACLANT and NORDA Numerical Modeling groups in developing, implementing, and testing a new fast general-applications wave propagation algorithm, SAFARI, formulated at SACLANT is summarized. The present general-applications SAFARI program uses a Direct Global Matrix Approach to multilayer Green's function calculation. A rapid and unconditionally stable solution is readily obtained via simple Gaussian ellimination on the resulting sparsely banded block system, precisely analogous to that arising in the Finite Element Method. The resulting gains in accuracy and computational speed allow consideration of much larger multilayered air/ocean/Earth/engineering material media models, for many more source-receiver configurations than previously possible. The validity and versatility of the SAFARI-DGM method is demonstrated by reviewing three practical examples of engineering interest, drawn from ocean acoustics, engineering seismology and ultrasonic scattering.

Schmidt, H.↗

On the validation of a code and a turbulence model appropriate to circulation control airfoils

A computer code for calculating flow about a circulation control airfoil within a wind tunnel test section has been developed. This code is being validated for eventual use as an aid to design such airfoils. The concept of code validation being used is explained. The initial stages of the process have been accomplished. The present code has been applied to a low-subsonic, 2-D flow about a circulation control airfoil for which extensive data exist. Two basic turbulence models and variants thereof have been successfully introduced into the algorithm, the Baldwin-Lomax algebraic and the Jones-Launder two-equation models of turbulence. The variants include adding a history of the jet development for the algebraic model and adding streamwise curvature effects for both models. Numerical difficulties and difficulties in the validation process are discussed. Turbulence model and code improvements to proceed with the validation process are also discussed.

Viegas, J. R.↗

On the validation of a code and a turbulence model appropriate to circular control airfoils

A computer code for calculating flow about a circulation control airfoil within a wind tunnel test section was developed. This code is being validated for eventual use as an aid to design such airfoils. The concept of code validation being used is explained. The initial stages of the process were accomplished. The present code was applied to a low subsonic, 2-D flow about a circulation control airfoil for which extensive data exist. Two basic turbulence models and variants thereof were successfully introduced into the algorithm, the Baldwin-Lomax algebraic and the Jones-Launder two equation models of turbulence. The variants include adding a history of the jet development for the algebraic model and adding streamwise curvature effects for both models. Numerical difficulties and difficulties in the validation process are discussed. Turbulence model and code improvements to proceed with the validation process are also discussed.

Viegas, J. R.↗

The role of radiation-dynamics interaction in regional numerical weather prediction

The role of radiation-dynamics interaction in regional numerical weather prediction of severe storm environment and mesoscale convective systems over the United States is researched. Based upon the earlier numerical model simulation experiments, it is believed that such interaction can have a profound impact on the dynamics and thermodynamics of regional weather systems. The research will be carried out using real-data model forecast experiments performed on the Cray-X/MP computer. The forecasting system to be used is a comprehensive mesoscale prediction system which includes analysis and initialization, the dynamic model, and the post-forecast diagnosis codes. The model physics are currently undergoing many improvements in parameterizing radiation processes in the model atmosphere. The forecast experiments in conjunction with in-depth model verification and diagnosis are aimed at a quantitative understanding of the interaction between atmospheric radiation and regional dynamical processes in mesoscale models as well as in nature. Thus, significant advances in regional numerical weather prediction can be made. Results shall also provide valuable information for observational designs in the area of remote sensing techniques to study the characteristics of air-land thermal interaction and moist processes under various atmospheric conditions.

Chang, Chia-Bo↗

Global ocean tide mapping using TOPEX/Poseidon altimetry

The investigation's main goals are to produce accurate tidal maps of the main diurnal, semidiurnal, and long-period tidal components in the world's deep oceans. This will be done by the application of statistical estimation techniques to long time series of altimeter data provided by the TOPEX/POSEIDON mission, with additional information provided by satellite tracking data. In the prelaunch phase, we will use in our simulations and preliminary work data supplied by previous oceanographic missions, such as Seasat and Geosat. These results will be of scientific interest in themselves. The investigation will also be concerned with the estimation of new values, and their uncertainties, for tidal currents and for the physical parameters appearing in the Laplace tidal equations, such as bottom friction coefficients and eddy viscosity coefficients. This will be done by incorporating the altimetry-derived charts of vertical tides as boundary conditions in the integration of those equations. The methodology of the tidal representation will include the use of appropriate series expansions such as ocean-basin normal modes and spherical harmonics. The results of the investigation will be space-determined tidal models of coverage and accuracy superior to that of the present numerical models of the ocean tides, with the concomitant benefits to oceanography and associated disciplinary fields.

Sanchez, Braulio V.↗

Radiation Transport Modeling and Assessment to Better Predict Radiation Exposure, Dose, and Toxicological Effects to Human Organs on Long Duration Space Flights

NASA's long-range plans include possible human exploratory missions to the moon and Mars within the next quarter century. Such missions beyond low Earth orbit will expose crews to transient radiation from solar particle events which include high-energy galactic cosmic rays and high-energy protons. Because the radiation levels in space are high and the missions long, adequate shielding is needed to minimize the deleterious health effects of exposure to radiation. The focus of this study is radiation exposure to the blood-forming organs of the NASA astronauts. NASA/JSC developed the Phantom Torso Experiment for Organ Dose Measurements which housed active and passive dosimeters that would monitor and record absorbed radiation levels at vital organ locations. This experiment was conducted during the STS-9 I mission in May '98 and provided the necessary space radiation data for correlation to results obtained from the current analytical models used to predict exposure to the blood-forming organs. Numerous models (i.e., BRYNTRN and HZETRN) have been developed and used to predict radiation exposure. However, new models are continually being developed and evaluated. The Space Environment Information Systems (SPENVIS) modeling program, developed by the Belgian Institute for Space Aeronomy, is to be used and evaluated as a part of the research activity. It is the intent of this research effort to compare the modeled data to the findings from the STS-9 I mission; assess the accuracy and efficiency of this model; and to determine its usefulness for predicting radiation exposure and developing better guidelines for shielding requirements for long duration manned missions.

Denkins, Pamela↗

Leveraging Pre-storm Soil Moisture Estimates for Enhanced Land Surface Model Calibration in Ungauged Hydrologic Basins

Despite long-standing efforts, hydrologists still lack robust tools for calibrating land surface model (LSM) streamflow estimates within ungauged basins. Using surface soil moisture estimates from the Soil Moisture Active Passive Level 4 Soil Moisture (L4_SM) product, precipitation observations, and streamflow gauge measurements for 617 medium-scale (200-10,000 km2) basins in the contiguous United States, we measure the temporal (Spearman) rank correlation between antecedent (i.e., pre-storm) surface soil moisture (ASM) and the storm-scale runoff coefficient (RC; the fraction of storm-scale precipitation accumulation converted into streamflow). In humid and semi-humid basins, this rank correlation is shown to be sufficiently strong to allow for the substitution of storm-scale RC observations (available only in basins that are both lightly regulated and gauged) with high-quality ASM values (available quasi-globally from L4_SM) in streamflow calibration procedures. Using this principle, we define a new, basin-wise LSM streamflow calibration approach based on L4_SM alone and successfully apply it to identify LSM configurations that produce a high rank correlation with observed RC. However, since the approach cannot detect RC bias, it is less successful in identifying LSM configurations with low mean-absolute error. Plain Text Summary Accurately forecasting the fraction of rainfall that runs off into streams, as opposed to infiltrates into the soil, is critical for flash-flood prediction, water-resource monitoring, and tracking the transport of nutrients from agricultural fields into local waterways. Such forecasting is typically performed by hydrologic models that attempt to represent the physical processes responsible for surface runoff generation. However, to provide accurate streamflow forecasts, these models typically need to be calibrated against actual streamflow observations. This is problematic given the relatively poor, and declining, global availability of stream gauges. This paper presents a novel model calibration strategy that uses soil moisture from remote sensing and numerical modeling in place of streamflow observations during calibration. This transition has significant practical advantages because, unlike streamflow observations, the soil moisture data are continuously available across space. Our results demonstrate that this new approach can significantly improve hydrologic models within humid and semi-humid basins lacking sufficient ground-based instrumentation for traditional streamflow calibration.

hydrologic models↗

Martian B Storm Evolution: Modeling Dust Activity over the Receding South Polar CO 2 Ice Cap at Southern Hemisphere Summer Solstice

Observations of annually recurring, elevated (25 km) atmospheric warming and associated dust concentrations over the south polar cap near Southern Hemisphere (SH) summer solstice are the defining features of the regional dust storm known as the B storm [1]. Our observational analysis of MCS and TES temperature and dust retrievals include all Mars Years (MY) lacking a global dust storm (GDS) between MY 26 and MY 33 (Figure 1). The data indicate the B storm is initiated around perihelion (L s =252°) in the southern midlatitudes and reaches maximum height and intensity over the south pole near SH summer solstice at L s ~267°. At peak intensity, mid-level atmospheric temperatures are elevated by as much as 30 K and mid-level atmospheric dust mixing ratios are more than triple pre-B storm levels (1.5-2 ppm to 7-8 ppm). Cap edge lifting may contribute significantly to initial dust mobilization in the B storm since surface wind stress along the receding cap edge is likely strengthened by the combined effects of the sea breeze circulation and katabatic winds at this time of year [2, 3, 4]. We are investigating dust lifting and lofting near the south pole of Mars during SH summer using the NASA Ames Mars Global Climate Model (MGCM). The NASA Ames MGCM is a finite-volume numerical model with a cubed sphere grid ideal for modeling processes at high latitudes. Initial simulations at low resolution (4x4 degrees) capture some elevated dust and associated warming at the south pole around L s =265°, but the model fails to capture the intensity of the B storm. We will examine and present the model-predicted surface stress and vertical velocities near the receding seasonal CO 2 cap edge and at the south pole during solstice. Finally, we will present higher resolution simulations to better capture small-scale dynamics and features at the pole.

Courtney M. L. Batterson↗

Use of an adjoint model for finding triggers for Alpine lee cyclogenesis

The authors propose a new procedure, designated the adjoint-based genesis diagnostic (AGD) procedure, for studying triggering mechanisms and the subsequent genesis of the synoptic phenomena of interest. This procedure makes use of a numerical model sensitivity to initial conditions and the nonlinear evolution of the initial perturbations that are designed using this sensitivity. The model sensitivity is evaluated using the associated adjoint model. This study uses the dry version of the National Center for Atmospheric Research Mesoscale Adjoint Modeling System (MAMS) for the numerical experiments. The authors apply the AGD procedure to two cases of Alpine lee cyclogenesis that were observed during the Alpine Experiment special observations period. The results show that the sensitivity fields that are produced by the adjoint model and the associated initial perturbations are readily related to the probable triggering mechanisms for these cyclones. Additionally, the nonlinear evolution of these initial perturbations points toward the physical processes involved in the lee cyclone formation. The AGD experiments for a weak cyclone case indicate that the MAMS forecast model has an underrepresented topographic forcing due to the sigma vertical coordinate and that this model error can be compensated by adjustments in the initial conditions that are related to the triggering mechanisms, which is not associated with the topographic blocking mechanism.

Vukicevic, Tomislava↗

A coupled ice-ocean model of ice breakup and banding in the marginal ice zone

A coupled ice-ocean numerical model for the marginal ice zone is considered. The model consists of a nonlinear sea ice model and a two-layer (reduced gravity) ocean model. The dependence of the upwelling response on wind stress direction is discussed. The results confirm earlier analytical work. It is shown that there exist directions for which there is no upwelling, while other directions give maximum upwelling in terms of the volume of uplifted water. The ice and ocean is coupled directly through the stress at the ice-ocean interface. An interesting consequence of the coupling is found in cases when the ice edge is almost stationary. In these cases the ice tends to break up a few tenths of kilometers inside of the ice edge.

Smedstad, O. M.↗

Assimilation of Altimeter Data into a Quasigeostrophic Model of the Gulf Stream System: Dynamical Considerations - Part 1

The dynamical consequences of constraining a numerical model with sea surface height data have been investigated. The model used for this study is a quasigeostrophic model of the Gulf Stream region. The data that have been assimilated are maps of sea surface height obtained as the superposition of sea surface height variability deduced from the Geosat altimeter measurements and a mean field constructed from historical hydrographic data. The method used for assimilating the data is the nudging technique. Nudging has been implemented in such a way as to achieve a high degree of convergence of the surface model fields toward the observations. The assimilation of the surface data is thus equivalent to the prescription of a surface pressure boundary condition. The authors analyzed the mechanisms of the model adjustment and the characteristics of the resultant equilibrium state when the surface data are assimilated. Since the surface data are the superposition of a mean component and an eddy component, in order to understand the relative role of these two components in determining the characteristics of the final equilibrium state, two different experiments have been considered: in the first experiment only the climatological mean field is assimilated, while in the second experiment the total surface streamfunction field (mean plus eddies) has been used. It is shown that the model behavior in the presence of the surface data constraint can be conveniently described in terms of baroclinic Fofonoff modes. The prescribed mean component of the surface data acts as a 'surface topography' in this problem. Its presence determines a distortion of the geostrophic contours in the subsurface layers, thus constraining the mean circulation in those layers. The intensity of the mean flow is determined by the inflow/outflow conditions at the open boundaries, as well as by eddy forcing and dissipation.

Capotondi, Antonietta↗

Modelling the Phanerozoic carbon cycle and climate: constraints from the 87Sr/86Sr isotopic ratio of seawater

A numerical model describing the coupled evolution of the biogeochemical cycles of carbon, sulfur, calcium, magnesium, phosphorus, and strontium has been developed to describe the long-term changes of atmospheric carbon dioxide and climate during the Phanerozoic. The emphasis is on the effects of coupling the cycles of carbon and strontium. Various interpretations of the observed Phanerozoic history of the seawater 87Sr/86Sr ratio are investigated with the model. More specifically, the abilities of continental weathering, volcanism, and surface lithology in generating that signal are tested and compared. It is suggested that the observed fluctuations are mostly due to a changing weatherability over time. It is shown that such a conclusion is very important for the modelling of the carbon cycle. Indeed, it implies that the conventional belief that the evolution of atmospheric carbon dioxide and climate on a long time scale is governed by the balance between the volcanic input of CO2 and the rate of silicate weathering is not true. Rather carbon exchanges between the mantle and the exogenic system are likely to have played a key role too. Further, the increase of the global weathering rates with increasing surface temperature and/or atmospheric CO2 pressure usually postulated in long-term carbon cycle and climate modelling is also inconsistent with the new model. Other factors appear to have modulated the weatherability of the continents through time, such as mountain building and the existence of glaciers and ice sheets. Based on these observations, a history of atmospheric carbon dioxide and climate during Phanerozoic time, consistent with the strontium isotopic data, is reconstructed with the model and is shown to be compatible with paleoclimatic indicators, such as the timing of glaciation and the estimates of Cretaceous paleotemperatures.

NASA Discipline Number 52-30↗