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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
NASA NTRS1980

Acid rain: Mesoscale model

A mesoscale numerical model of the Florida peninsula was formulated and applied to a dry, neutral atmosphere. The prospective use of the STAR-100 computer for the submesoscale model is discussed. The numerical model presented is tested under synoptically undisturbed conditions. Two cases, differing only in the direction of the prevailing geostrophic wind, are examined: a prevailing southwest wind and a prevailing southeast wind, both 6 m/sec at all levels initially.

Hsu, H. M.↗
NASA NTRS1981

Advantages of multigrid methods for certifying the accuracy of PDE modeling

Numerical techniques for assessing and certifying the accuracy of the modeling of partial differential equations (PDE) to the user's specifications are analyzed. Examples of the certification process with conventional techniques are summarized for the three dimensional steady state full potential and the two dimensional steady Navier-Stokes equations using fixed grid methods (FG). The advantages of the Full Approximation Storage (FAS) scheme of the multigrid technique of A. Brandt compared with the conventional certification process of modeling PDE are illustrated in one dimension with the transformed potential equation. Inferences are drawn for how MG will improve the certification process of the numerical modeling of two and three dimensional PDE systems. Elements of the error assessment process that are common to FG and MG are analyzed.

Forester, C. K.↗
DOE OSTI2024

Code-to-code comparison between FLASH and HYDRA in gas-puff Z-pinch modeling

The numerical modeling of gas-puff Z pinches involves the nonlinear coupling of a broad range of complex, multi-physics phenomena that makes such simulations challenging. The challenge is further compounded by nonlinear instabilities that can impact the dynamics of imploding gas-puff Z pinches, such as the magneto Rayleigh–Taylor instability (MRTI). If the growth rate and amplitude of the latter is comparable to the relevant timescales and properties of the imploding plasma, the MRTI can dramatically alter implosion dynamics, dictate pinch stability, and govern the plasma properties achievable in pulsed-power-driven laboratory experiments. National Laboratories and academic teams have developed numerical tools that can accurately model Z-pinch configurations and provide reliable design capabilities that can guide experimental choices and assist in interpreting experimental results. Most such tools, however, are not broadly available. Here, we present newly developed Z-pinch simulation capabilities of the publicly available FLASH code, applied in the study of MRTI growth and dynamical effects in gas-puff implosions. To verify the new implementations, we perform a comparison of FLASH gas-puff implosion simulations with previously published calculations with the HYDRA code from Lawrence Livermore National Laboratory, which have been validated with experimental data from the CESZAR pulsed-power driver at the University of California, San Diego. The experiments involved double- and triple-nozzle configurations, in an experimental attempt to stabilize the pinch to the MRTI. The code-to-code comparison shows similar results between the FLASH and HYDRA simulations, supporting the use of FLASH in the modeling of future gas-puff Z-pinch experiments at CESZAR.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2014

Assessing the Risk of Crew Injury Due to Dynamic Loads During Spaceflight

Spaceflight requires tremendous amounts of energy to achieve Earth orbit and to attain escape velocity for interplanetary missions. Although the majority of the energy is managed in such a way as to limit the accelerations on the crew, several mission phases may result in crew exposure to dynamic loads. In the automotive industry, risk of serious injury can be tolerated because the probability of a crash is remote each time a person enters a vehicle, resulting in a low total risk of injury. For spaceflight, the level of acceptable injury risk must be lower to achieve a low total risk of injury because the dynamic loads are expected on each flight. To mitigate the risk of injury due to dynamic loads, the NASA Human Research Program has developed a research plan to inform the knowledge gaps and develop relevant tools for assessing injury risk. The risk of injury due to dynamic loads can be further subdivided into extrinsic and intrinsic risk factors. Extrinsic risk factors include the vehicle dynamic profile, seat and restraint design, and spacesuit design. Human tolerance to loads varies considerably depending on the direction, amplitude, and rise-time of acceleration therefore the orientation of the body to the dynamic vector is critical to determining crew risk of injury. Although a particular vehicle dynamic profile may be safe for a particular design, the seat, restraint, and suit designs can affect the risk of injury due to localized effects. In addition, characteristics intrinsic to the crewmember may also contribute to the risk of injury, such as crewmember sex, age, anthropometry, and deconditioning due to spaceflight, and each astronaut may have a different risk profile because of these factors. The purpose of the research plan is to address any knowledge gaps in the risk factors to mitigate injury risk. Methods for assessing injury risk have been well documented in other analogous industries and include human volunteer testing, human exposure to dynamic environments, post-mortem human subject (PMHS) testing, animal testing, anthropomorphic test devices (ATD), dynamic models of the human, numerical models of ATDs, and numerical models of the human. Each has inherent strengths and limitations. For example, human volunteer testing is advantageous because a population can be selected that is similar to the astronaut corps; however, because of the inherent ethical limitations, only sub-injurious conditions can be tested. PMHSs can be tested in a variety of conditions including injurious levels, but the responses are not completely analogous to living human subjects. In addition, it is exceedingly difficult to select a PMHS population that is similar to the astronaut corps. ATDs are currently widely used in the automotive industry and military because they are highly repeatable and durable. Unfortunately, because they are mechanical models of the human body, the biofidelity of the responses are limited to dynamic conditions used to validate the ATD. Numerical models of the ATD, in addition to the strengths and limitations for ATDs, are easy to use for a variety of designs before a design is fabricated, but also have additional limitations for ATDs, are easy to use for a variety of designs before a design is fabricated, but also have additional uncertainty. Dynamic models are simple and easy to use, but do not account for localized effects of the seat and suit. Finally, numerical models of the human have the potential to have the most advantages; however, the current models are not validated for the conditions expected during spaceflight. To properly assess spaceflight conditions with numerical human models, human data would be needed to optimize the model responses for those conditions. Using the appropriate assessment method with the knowledge gained for each risk factor, an appropriate approach for mitigating the risk of injury due to dynamic loads can be developed ensuring crew safety in future NASA vehicles.

Somers, J. T.↗
NASA NTRS1982

On the stability of Schwarzschild's triaxial galaxy model

The numerical model for a triaxial galaxy constructed by Schwartzschild (1979) has been tested by means of fully self-consistent three-dimensional n-body numerical experiments that use a representation of this model to start the integration. The integration, study models, initial load, and correcting potential of the experiment are detailed. Preliminary checks included density values, potential values, and integration in potential of an ideal model. The results showed some particle loss and no significant shape change. The model is rugged and robust, and no growing disturbances with growth rate in excess of 0.5 per crossing time could be detected. Runs were done at 50 percent excess halo and with tube orbits reversed.

Smith, B. F.↗
NASA NTRS1991

Measurement of the performance of a spiral wound polyimide regenerator in a pulse tube refrigerator

A regenerator for use in a pulse tube refrigerator has been constructed from a polyimide (polypyromellitimide or PPMI) whose small ratio of thermal conductivity to heat capacity make it a good candidate for a regenerator material in cryocoolers. The regenerator was fabricated using 25 micron thick photoresist strips bonded to a 50 micron thick sheet of PPMI. This composite sheet was wound in jelly-roll fashion around a mandrel and inserted into the regenerator housing. The photoresist strips, formed using a photolithographic technique, provided a 25 micron spacing for the axial flow of gas between each layer of PPMI. Ineffectiveness results are presented for this material under actual operating conditions in a pulse tube refrigerator and compared with a numerical model. The numerical model indicated that a polyimide regenerator would perform much better than one constructed of stainless steel screen, but the experimental results showed the opposite behavior. Measured values for the ineffectiveness were 0.003 for the stainless steel screen and 0.017 for the polyimide.

Rawlins, Wayne↗
DOE OSTI2024

System and method for wave prediction

A method and system for prediction of wave properties include collecting time-series data streams from one or more wave measurement devices and processing the data to identify data parameters to establish boundary conditions of a numerical model. The numerical model may be used to compute a predicted wave field of time-series data for a variety of wave properties at a target location.

Previsic, Mirko↗
NASA NTRS1988

A study of the dynamics of the Intertropical Convergence Zone (ITCZ) in a symmetric atmosphere-ocean model

A numerical model of the circulation of a coupled axisymmetric atmosphere-ocean system was constructed to investigate the physical factors governing the location and intensity of the Intertropical Convergence Zone (ITCZ) over oceans and over land. The results of several numerical integrations are presented to illustrate the interaction of the individual atmospheric and oceanic circulations. It is shown that the ITCA cannot be located at the equator because the atmosphere-ocean system is unstable for lateral displacements of the ITCA from an equilibrium position at the equator.

Charney, J. G.↗
NASA NTRS1993

The generation and use of numerical shape models for irregular Solar System objects

We describe a procedure that allows the efficient generation of numerical shape models for irregular Solar System objects, where a numerical model is simply a table of evenly spaced body-centered latitudes and longitudes and their associated radii. This modeling technique uses a combination of data from limbs, terminators, and control points, and produces shape models that have some important advantages over analytical shape models. Accurate numerical shape models make it feasible to study irregular objects with a wide range of standard scientific analysis techniques. These applications include the determination of moments of inertia and surface gravity, the mapping of surface locations and structural orientations, photometric measurement and analysis, the reprojection and mosaicking of digital images, and the generation of albedo maps. The capabilities of our modeling procedure are illustrated through the development of an accurate numerical shape model for Phobos and the production of a global, high-resolution, high-pass-filtered digital image mosaic of this Martian moon. Other irregular objects that have been modeled, or are being modeled, include the asteroid Gaspra and the satellites Deimos, Amalthea, Epimetheus, Janus, Hyperion, and Proteus.

Simonelli, Damon P.↗
NASA NTRS1995

Recent improvements to and validation of the one dimensional NASA wave rotor model

A numerical model has been developed at the NASA Lewis Research Center which can predict both the unsteady flow quantities within a wave rotor passage and the steady averaged flows in the ports. The model is based on the assumptions of one-dimensional, unsteady, perfect gas flow. The model assesses not only the dominant wave behavior, but the loss effects of finite passage opening time, leakage from the passage ends, viscosity, and heat transfer to and from the passages. The model operates in the rotor reference frame; however, until recently no account was made for the often significant effect of the rotor circumferential velocity component. The present model accounts for this by modifying the passage boundary conditions, allowing the internal computational scheme to remain the rotor reference frame, while quantities such as inlet duct stagnation properties may be specified in the fixed or absolute reference frame. Accurate modeling of this effect is critical to successful wave rotor analysis and design, particularly in off-design predictions where the flows in the inlet ducts are mismatched with the rotor passages and significant turning may take place (i.e., work is done on the gas). The relative simplicity of the model makes it useful for design and optimization, as well as analysis, of wave rotor cycles for many applications. This report, building on several earlier papers, describes the most recent modifications to the model. These include accounting for the relative/absolute transition at the passage boundaries and refinements to the viscous source term correlation which resulted from this accounting. Comparison of model predictions with measured data is then presented and discussed.

Paxson, Daniel E.↗
NASA NTRS1979

A nested numerical tidal model of the southern New England bight

Efforts were focused on the development and application of a three-dimensional numerical model for predicting pollutant and sediment transport in estuarine and coastal environments. To successfully apply the pollutant and sediment transport model to Rhode Island coastal waters, it was determined that the flow field in this region had to be better described through the use of existing numerical circulation models. A nested, barotropic numerical tidal model was applied to the southern New England Bight (Long Island, Block Island, Rhode Island Sounds, Buzzards Bay, and the shelf south of Block Island). Forward time and centered spatial differences were employed with the bottom friction term evaluated at both time levels. Using existing tide records on the New England shelf, adequate information was available to specify the tide height boundary condition further out on the shelf. Preliminary results are within the accuracy of the National Ocean Survey tide table data.

Gordon, R. B.↗
NASA NTRS1989

Numerical simulations - Some results for the 2- and 3-D Hubbard models and a 2-D electron phonon model

Numerical simulations on the half-filled three-dimensional Hubbard model clearly show the onset of Neel order. Simulations of the two-dimensional electron-phonon Holstein model show the competition between the formation of a Peierls-CDW state and a superconducting state. However, the behavior of the partly filled two-dimensional Hubbard model is more difficult to determine. At half-filling, the antiferromagnetic correlations grow as T is reduced. Doping away from half-filling suppresses these correlations, and it is found that there is a weak attractive pairing interaction in the d-wave channel. However, the strength of the pair field susceptibility is weak at the temperatures and lattice sizes that have been simulated, and the nature of the low-temperature state of the nearly half-filled Hubbard model remains open.

Scalapino, D. J.↗
NASA NTRS1972

A numerical study of a three dimensional spherical thermospheric density and wind model

Numerical calculations of the generation and propagation of the two important fundamental symmetric tidal wave modes - the diurnal mode (1, 1, 1,) and the semidiurnal mode (2, 2, 2) - were performed applying a realistic model thermosphere and taking into account heat conduction and the temporally and spatially varying ion-neutral collision number. Both wave modes are predominantly generated by the solar EUV heat input. It is shown that the latitude structure of the (1, 1, 1)-mode which is identical with the Hough function(1, -1) within the lower non-dissipative atmosphere degenerates into the spherical function P sub 1, 1 at thermospheric heights. The pressure field of this mode constitutes the observed pressure bulge of the thermosphere, the diurnal component of which peaks at 15 h L. T. The electric polarization field of the geomagnetic Sq current generates a significant fraction of this wave mode at F layer heights. This wave component shifts the total horizontal wind system to earlier times by about 1 hour in agreement with ionospheric observations. The latitude structure of the (2, 2, 2) mode is identical with the Hough function (2, 2) within the lower non-dissipative atmosphere. It degenerates to the spherical function P sub 2, 2 at thermospheric heights.

Volland, H.↗
NASA NTRSDate not supplied

Using XR for Improving Scientific Discovery With Numerical Weather Models

Earth science (ES) digital twins will help us understand the complex interactions and interrelationships that make up our Earth system and the impacts of earth science phenomena on it. Our work addresses two underdeveloped areas in current ES digital twin work: improving the understanding and interaction with ES model outputs by using Virtual and Mixed Reality (XR) tools and improving the non-intuitive mapping of continuous ES natural phenomena to gridded reference frames in current numerical models. Traditionally, scientists working on ES view and analyze the results of calculated or measured observables with static 1-dimensional (1D), 2D or 3D plots displayed on flat computer screens or paper. Using such limited mediums, it can be very difficult to identify, track and understand the evolution of key features due to poor viewing angles and the nature of flat computer screens. In addition, numerical models, such as the NASA Goddard Earth Observing System (GEOS) ES model, are almost exclusively formulated, visualized and analyzed in an Eulerian reference frame with fixed grid points in space and time. However, ES phenomena such as convective clouds, hurricanes and wildfire smoke plumes are visualized and analyzed in a Lagrangian reference frame: therefore it is often difficult and unnatural to understand these phenomena in relation to each other, visualized either in an Eulerian or Lagrangian context. In 3D visualizations, data generally takes one of three forms: gridded (e.g., voxelized) data, where space is divided into regions; point clouds, where data is represented as a set of points; and meshes, where objects are rendered as surfaces composed of small polygons (usually triangles). A gridded, Eulerian reference frame has been the default representation for the 2D visual analysis of atmospheric data in part because the numerical methods used to generate atmospheric model data in the first place use a gridded approach, with equations defining the relationships between the physical variables in each of a grid's cells across successive timesteps. In our work, we are particularly interested in data from GEOS. Another reason why gridded representations tend to be used for visualizing data from such models is because trajectories are difficult to interpret from representations on 2D surfaces, due to line-of-sight ambiguity. Instead of a fixed grid from GEOS, we embed a trajectory model to simulate particles' movement throughout a GEOS run. We then ingest these particle trajectories as animated point clouds with a NASA open source XR toolkit, the Mixed Reality Exploration Toolkit (MRET), and merge GEOS data with ES phenomena data onto one combined visualization that the user can intuitively interact with. Efficient rendering of arbitrarily large point clouds is an ongoing challenge being addressed by the computer science community, with the GPU-based optimizations and efficient GPU memory utilization a common theme of recent advances, especially for XR, where sustained high frame rate is mandatory to save the user from suffering due to simulation sickness. In this work, we describe and evaluate our progress in choosing and implementing appropriate methods for rendering arbitrarily large point clouds within MRET for XR. While tracking the XR headset enables the immersion of a user within a 3D scene of a data visualization, tracking of XR handheld controllers or user’s hands enables us to implement intuitive user interactions with the visualized datasets. Conventional tools require a user working with an ES visualization to conduct many interactions to commit their intended selections or manipulations with a visualized dataset; for example to specify a set of points in 3D space. Doing so in a 2D flat screen interface has traditionally required specifying a set of points in three distinct 2D coordinate systems (XY, XZ, and YZ), which is cumbersome. In other scientific domains, it has been shown that specifying or selecting a location or volume in XR using handheld controllers or tracked hands allows for greater speed and accuracy. We anticipate the same will hold true for atmospheric data, and we will share initial results of measuring the utility of such an interface. Notably, as the data being visualized is generated by GEOS as a prediction based on initial conditions, an intended application of our tool is to serve as part of an iterative feedback loop. Through XR, a scientist will review and manipulate a GEOS model run, modifying the conditions as needed to do subsequent runs of GEOS. Thereby, XR-based improvements to speed and accuracy of 3D tagging of points minimizes the effort required by both the scientist and the computer cluster conducting the necessary calculations.

Thomas Grubb↗
NASA NTRS1985

Cirrus clouds. I - A cirrus cloud model. II - Numerical experiments on the formation and maintenance of cirrus

A simplified cirrus cloud model is presented which may be used to investigate the role of various physical processes in the life cycle of a cirrus cloud. The model is a two-dimensional, time-dependent, Eulerian numerical model where the focus is on cloud-scale processes. Parametrizations are developed to account for phase changes of water, radiative processes, and the effects of microphysical structure on the vertical flux of ice water. The results of a simulation of a thin cirrostratus cloud are given. The results of numerical experiments performed with the model are described in order to demonstrate the important role of cloud-scale processes in determining the cloud properties maintained in response to larger scale forcing. The effects of microphysical composition and radiative processes are considered, as well as their interaction with thermodynamic and dynamic processes within the cloud. It is shown that cirrus clouds operate in an entirely different manner than liquid phase stratiform clouds.

Starr, D. OC.↗
NASA NTRS1978

Treatment of cloudiness in regional scale - Numerical weather prediction models

A numerical weather prediction model is considered in which cloud water and precipitation are explicitly predicted as functions of three space dimensions and time. Approximations in the thermodynamic and continuity equations are discussed in detail, as well as the cumulus parameterization and the cloud microphysics parameterization. An example of a real-data forecast indicates that very accurate quantitative precipitation and cloud cover forecasts are possible using initial observations from only the relatively sparse rawinsonde network in the U.S.

Kreitzberg, C. W.↗
NASA NTRS2006

Constraining Numerical Geodynamo Modeling with Surface Observations

Numerical dynamo solutions have traditionally been generated entirely by a set of self-consistent differential equations that govern the spatial-temporal variation of the magnetic field, velocity field and other fields related to dynamo processes. In particular, those solutions are obtained with parameters very different from those appropriate for the Earth s core. Geophysical application of the numerical results therefore depends on correct understanding of the differences (errors) between the model outputs and the true states (truth) in the outer core. Part of the truth can be observed at the surface in the form of poloidal magnetic field. To understand these differences, or errors, we generate new initial model state (analysis) by assimilating sequentially the model outputs with the surface geomagnetic observations using an optimal interpolation scheme. The time evolution of the core state is then controlled by our MoSST core dynamics model. The final outputs (forecasts) are then compared with the surface observations as a means to test the success of the assimilation. We use the surface geomagnetic data back to year 1900 for our studies, with 5-year forecast and 20-year analysis periods. We intend to use the result; to understand time variation of the errors with the assimilation sequences, and the impact of the assimilation on other unobservable quantities, such as the toroidal field and the fluid velocity in the core.

Kuang, Weijia↗