Search NASA⌕ Search

SEARCH · Search NASA

Results for “Implicit modeling”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

MARIKA - A model revision system using qualitative analysis of simulations

This paper describes portions of a novel system called MARIKA (Model Analysis and Revision of Implicit Key Assumptions) to automatically revise a model of the normal human orientation system. The revision is based on analysis of discrepancies between experimental results and computer simulations. The discrepancies are calculated from qualitative analysis of quantitative simulations. The experimental and simulated time series are first discretized in time segments. Each segment is then approximated by linear combinations of simple shapes. The domain theory and knowledge are represented as a constraint network. Incompatibilities detected during constraint propagation within the network yield both parameter and structural model alterations. Interestingly, MARIKA diagnosed a data set from the Massachusetts Eye and Ear Infirmary Vestibular Laboratory as abnormal though the data was tagged as normal. Published results from other laboratories confirmed the finding. These encouraging results could lead to a useful clinical vestibular tool and to a scientific discovery system for space vestibular adaptation.

Groleau, Nicolas↗
NASA NTRS1977

A solution to the surface intersection problem

An application-independent geometric model within a data base framework should support the use of Boolean operators which allow the user to construct a complex model by appropriately combining a series of simple models. The use of these operators leads to the concept of implicitly and explicitly defined surfaces. With an explicitly defined model, the surface area may be computed by simply summing the surface areas of the bounding surfaces. For an implicitly defined model, the surface area computation must deal with active and inactive regions. Because the surface intersection problem involves four unknowns and its solution is a space curve, the parametric coordinates of each surface must be determined as a function of the arc length. Various subproblems involved in the general intersection problem are discussed, and the mathematical basis for their solution is presented along with a program written in FORTRAN IV for implementation on the IBM 370 TSO system.

Timer, H. G.↗
NASA NTRS1995

Diabatic forcing and intialization with assimilation of cloud water and rainwater in a forecast model

In this study, diabatic forcing, and liquid water assimilation techniques are tested in a semi-implicit hydrostatic regional forecast model containing explicit representations of grid-scale cloud water and rainwater. Diabatic forcing, in conjunction with diabatic contributions in the initialization, is found to help the forecast retain the diabatic signal found in the liquid water or heating rate data, consequently reducing the spinup time associated with grid-scale precipitation processes. Both observational Special Sensor Microwave/Imager (SSM/I) and model-generated data are used. A physical retrieval method incorporating SSM/I radiance data is utilized to estimate the 3D distribution of precipitating storms. In the retrieval method the relationship between precipitation distributions and upwelling microwave radiances is parameterized, based upon cloud ensemble-radiative model simulations. Regression formulae relating vertically integrated liquid and ice-phase precipitation amounts to latent heating rates are also derived from the cloud ensemble simulations. Thus, retrieved SSM/I precipitation structures can be used in conjunction with the regression-formulas to infer the 3D distribution of latent heating rates. These heating rates are used directly in the forecast model to help initiate Tropical Storm Emily (21 September 1987). The 14-h forecast of Emily's development yields atmospheric precipitation water contents that compare favorably with coincident SSM/I estimates.

Raymond, William H.↗
NASA NTRS1988

Development of the general interpolants method for the CYBER 200 series of supercomputers

The General Interpolants Method (GIM) is a 3-D, time-dependent, hybrid procedure for generating numerical analogs of the conservation laws. This study is directed toward the development and application of the GIM computer code for fluid dynamic research applications as implemented for the Cyber 200 series of supercomputers. An elliptic and quasi-parabolic version of the GIM code are discussed. Turbulence models, algebraic and differential equations, were added to the basic viscous code. An equilibrium reacting chemistry model and an implicit finite difference scheme are also included.

Stalnaker, J. F.↗
NASA NTRSDate not supplied

Validation and Sensitivity Analyses of Arc-Jet Performance of Woven Thermal Protection Entry Systems

Woven thermal protection materials for the Adaptable, Deployable, Entry and Placement Technology (ADEPT) entry system are simulated using a dual layer formalism in a finite volume implicit framework. The model developed and demonstrated here includes multi-dimensional and time-varying thermal response of the material to arc-jet test conditions. This work further demonstrates the sensitivities of the macro-scale thermal response to the through-thickness and in-plane conductivities of the carbon fabric, and to other experimental measurements such as cloth emissivity and thickness. The modeling can be extended to the ADEPT system that would encompass a single-piece woven heat shield configuration

Pratibha Raghunandan↗
NASA NTRS1984

Forecasts of the 500 mb height using a dynamically oriented statistical model

The forecast skill of a simple dynamically inspired statistical model of the Northern Hemisphere 500 mb height field is evaluated in spectral and physical space for a variety of forecast lead times (1-32 days) and predictand averaging times (1-32 days). The model includes viscous damping, wave propagation, climatology and implicit stochastic forcing. The largest model skill was found for forecasts of the zonal flow and the largest waves. In general, the largest forecast skills were also associated with the largest forecast error, there being a slight geographic phase shift of the skill with respect to the error. Model skills for climate (time-averaged) forecasts are greater when using instantaneous rather than time averages to forecast time averages. Analysis of model errors suggests areas for improvement in representing forcing terms and model physics. However, the model error fields are largely 'white noise' which suggests that global forecast skills substantially larger than those obtained here are unlikely to be achieved by more sophisticated models.

Roads, J. O.↗
NASA NTRS1987

Molecular bridge effects on distant charge tunneling

The question arises as to whether different hydrocarbon bridges should give different electron-trransfer rates. This question is answered on the basis of energetic and topological (interference) effects that can be gleaned from one-electron models. A discussion of model compound electron-transfer rates based on this interpretation is given. The approximations implicit in the periodic models used here (and in previous work) are carefully justified.

Onuchic, Jose Nelson↗
NASA NTRS1993

Flexible language constructs for large parallel programs

The goal of the research described is to develop flexible language constructs for writing large data parallel numerical programs for distributed memory (MIMD) multiprocessors. Previously, several models have been developed to support synchronization and communication. Models for global synchronization include SIMD (Single Instruction Multiple Data), SPMD (Single Program Multiple Data), and sequential programs annotated with data distribution statements. The two primary models for communication include implicit communication based on shared memory and explicit communication based on messages. None of these models by themselves seem sufficient to permit the natural and efficient expression of the variety of algorithms that occur in large scientific computations. An overview of a new language that combines many of these programming models in a clean manner is given. This is done in a modular fashion such that different models can be combined to support large programs. Within a module, the selection of a model depends on the algorithm and its efficiency requirements. An overview of the language and discussion of some of the critical implementation details is given.

Rosing, Matthew↗
NASA NTRS1984

Implicit Numerical Methods in Meteorology

The development of a fully implicit finite-difference model, whose time step is chosen solely to resolve accurately the physical flow of interest is discussed. The method is based on an operator factorization which reduces the dimensionality of the implicit approach: at each time step only (spatially) one-dimensional block-tridiagonal linear systems must be solved. The scheme uses two time levels and is second-order accurate in time. Compact implicit spatial differences are used, yielding fourth-order accuracy both vertically and horizontally. In addition, the development of a fully interactive computer code is discussed. With this code the user will have a choice of models, with various levels of accuracy and sophistication, which are imbedded, as subsets of the fully implicit 3D code.

Augenbaum, J.↗
NASA NTRS1998

Sensitivity of Cirrus Simulations in Idealized Situations: The WG2 Test Cases

GCSS Cirrus Cloud Systems Working Group (WG2) is presently conducting a comparison of cirrus cloud models for idealized initial conditions. The experiments involve binary (off/on) tests of model sensitivity to infrared radiative processes, and thermal stratification, and vertical wind shear for situations of weakly forced (3 cm/s uplift) cold (-60 to -70 C) and warm (-35 to -50 C) cirrus clouds. A range of model types are involved including parcel, SCM, 2-D CRM, 3-D CRM and LES models. The test cases will be described and results from 2-dimensional cirrus cloud models with bulk microphysics (implicit second moment scheme) and explicit bin microphysics will be compared. Vertical ice mass flux (particle fall speed) is a critical model component leading to significant intermodel differences. Efforts are ongoing to better quantify this aspect. Future plans of WG2 will also be briefly described and include model comparisons for a well-observed case of cold (ARM IOP) cirrus and of warm (EUCREX) cirrus, as well as, a joint activity with WG4 to consider the treatment of anvil cirrus in a variety of models.

Starr, David OC.↗
NASA NTRS1991

A practical implementation of turbulence models for the computation of three-dimensional separated flows

An upwind MUSCL-type implicit scheme for the three-dimensional Navier-Stokes equations is presented and details on the implementation for three-dimensional flows of a 'diagonal' upwind implicit operator are developed. Turbulence models for separated flows are also described with an emphasis on the numerical specificities of the Johnson-King nonequilibrium model. Good predictions of separated two- and three-dimensional flows are demonstrated.

Marx, Yves P.↗
NASA NTRS1991

Fatigue life prediction of an intermetallic matrix composite at elevated temperatures

A strain-based fatigue life prediction method is proposed for an intermetallic matrix composite (IMC) under tensile cyclic loadings at elevated temperatures. Styled after the Universal Slopes method, the model utilizes the composite's tensile properties to estimate fatigue life. Factors such as fiber volume ratio, number of plys and temperature dependence are implicitly incorporated into the model through these properties. The model constants are determined by using unidirectional fatigue data at temperatures of 425 and 815 C. Fatigue lives from two independent sources are used to verify the model at temperatures of 650 and 760 C. Cross-ply lives at 760 C are also predicted. It is demonstrated that the correlation between experimental and predicted lives is within a factor of two.

Bartolotta, Paul A.↗
NASA NTRS1991

Fatigue life prediction of an intermetallic matrix composite at elevated temperatures

A strain-based fatigue life prediction method is proposed for an intermetallic matrix composite (IMC) under tensile cyclic loadings at elevated temperatures. Styled after the 'Universal Slopes' method, the model utilizes the composite's tensile properties to estimate fatigue life. Factors such as fiber volume ratio (Vf), number of plys and temperature dependence are implicitly incorporated into the model through these properties. The model constants are determined by using unidirectional fatigue data at temperatures of 425 and 815 C. Fatigue lives from two independent sources are used to verify the model at temperatures of 650 and 760 C. Cross-ply lives at 760 C are also predicted. It is demonstrated that the correlation between experimental and predicted lives is within a factor of two.

Bartolotta, P. A.↗
NASA NTRS2014

Development Of FIAT-Based Parametric Thermal Protection System Mass Estimating Relationships For NASA's Multi-Mission Earth Entry Concept

An integrated tool called the Multi Mission System Analysis for Planetary Entry Descent and Landing (M-SAPE) is being developed as part of NASAs In-Space Propulsion Technology (ISPT) program. Part of M-SAPEs development requires the formulation of mass estimating relationships (MERs) to determine the vehicle's Thermal Protection System (TPS) material and required thickness for safe Earth entry. The objective of this study was to develop MERs using simple correlations that were non-ITAR and matched as accurately as possible NASAs high-fidelity ablation modeling tool, the Fully Implicit Ablation and Thermal Analysis Program (FIAT ). These MERs would be a first-estimate for feasibility studies; it is understood that higher-fidelity modeling like FIAT would be necessary once a proposed trajectory was down-selected. The trajectory space for these MERS consisted of 840 different trajectories, and a materials heating limit was the main constraint for an allowable trajectory. MERs for the vehicle fore body included the ablating materials Phenolic Impregnated Carbon Ablator (PICA ) and Carbon Phenolic atop Advanced Carbon-Carbon. For the backshell the materials were Silicone Impregnated Reusable Ceramic Ablator (SIRCA ), Acusil II, SLA-561V, and LI-900. The MERFIAT ratio indicates MERs are accurate to within 14 percent (at one standard deviation) of FIAT prediction, and the most any MER can under-predict TPS thickness is 18.7 percent of FIAT prediction. This poster focuses on the development of these MERs, the resulting equations, model limitations, and model accuracy.

thermal protection system↗
NASA NTRS1998

The Stratigraphy and Evolution of the Lunar Crust

Reconstruction of stratigraphic relationships in the ancient lunar crust has proved to be a formidable task. The intense bombardment during the first 700 m.y. of lunar history has severely perturbed the original stratigraphy and destroyed the primary textures of all but a few nonmare rocks. However, a knowledge of the crustal stratigraphy as it existed prior to the cataclysmic bombardment about 3.9 Ga is essential to test the major models proposed for crustal origin, i.e., crystal fractionation in a global magmasphere or serial magmatism in a large number of smaller bodies. Despite the large difference in scale implicit in these two models, both require an efficient separation of plagioclase and mafic minerals to form the anorthositic crust and the mafic mantle. Despite the havoc wreaked by the large body impactors, these same impact processes have brought to the lunar surface crystalline samples derived from at least the upper half of the lunar crust, thereby providing an opportunity to reconstruct the stratigraphy in areas sampled by the Apollo missions. As noted, ejecta from the large multiring basins are dominantly, or even exclusively, of crustal origin. Given the most recent determinations of crustal thicknesses, this implies an upper limit to the depth of excavation of about 60 km. Of all the lunar samples studied, a small set has been recognized as "pristine", and within this pristine group, a small fraction have retained some vestiges of primary features formed during the earliest stages of crystallization or recrystallization prior to 4.0 Ga. We have examined a number of these samples that have retained some record of primary crystallization to deduce thermal histories from an analysis of structural, textural, and compositional features in minerals from these samples. Specifically, by quantitative modeling of (1) the growth rate and development of compositional profiles of exsolution lamellae in pyroxenes and (2) the rate of Fe-Mg ordering in orthopyroxenes, we can constrain the cooling rates of appropriate lunar samples. These cooling rates are used to compute depths of burial at the time of crystallization, which enable us to reconstruct parts of the crustal stratigraphy as it existed during the earliest stages of lunar history.

McCallum, I. Stewart↗
NASA NTRS2022

The Simulation and Subseasonal Forecasting of Hydrological Variables: Insights from a Simple Water Balance Model

Past work has shown that a land surface model’s (LSM’s) implicit (not explicitly coded) relationships between soil moisture and both evapotranspiration (ET) and runoff largely determine the LSM’s hydrological behavior. Here we estimate the relationships that appear to be operating in the real world and compare them to those of the LSM component of a state-of-the-art Earth system model (ESM). The two sets of relationships are determined by calibrating them within a simple water balance model (WBM): once using stream gauge observations from small, unregulated rivers over the eastern half of the U.S., and once using the runoffs generated by the LSM as part of a state-of-the-art atmospheric reanalysis. Hydrological simulations and subseasonal hydrological forecasts performed with the two calibrated versions of the WBM provide two key results. First, the version calibrated to the LSM-generated runoffs does successfully reproduce, to first order, the hydrological behavior of the full LSM within its ESM environment. Second, of the two WBM versions, the one calibrated to the observations reproduces more accurately a broad collection of fully independent streamflow observations as well as a similarly broad collection of in-situ soil moisture measurements. Taken together, the two results suggest that the observations-calibrated ET and runoff efficiency functions do successfully represent, at least to some degree, soil moisture controls over hydrological variability in Nature and can serve as potentially useful targets for further LSM development.

Water Balance Model↗
NASA NTRS1979

In-flight simulation with pilot-center of gravity offset and velocity mismatch

Similarity transformations which preserve modal characteristics and pilot's acceleration cues in in-flight simulation are presented. The model transformation for lateral acceleration matching is developed. A velocity-mismatch example, based on a VRA simulation of the Space Shuttle, illustrates that acceleration matching is achieved at the expense of mismatching in cues which are secondary to the simulated piloting task, while primarily cues are preserved. The approach is applicable for both implicit and explicit model-following, and it can easily be extended to the longitudinal case.

Stengel, R. F.↗
NASA NTRS1995

Numerical Simulations of High Enthalpy Pulse Facilities

Axisymmetric flows within shock tubes and expansion tubes are simulated including the effects of finite rate chemistry and both laminar and turbulent boundary layers. The simulations demonstrate the usefulness of computational fluid dynamics for characterizing the flows in high enthalpy pulse facilities. The modeling and numerical requirements necessary to simulate these flows accurately are also discussed. Although there is a large body of analysis which explains and quantifies the boundary layer growth between the shock and the interface in a shock tube, there is a need for more detailed solutions. Phenomena such as thermochemical nonequilibrium. or turbulent transition behind the shock are excluded in the assumptions of Mirels' analysis. Additionally there is inadequate capability to predict the influence of the boundary layer on the expanded gas behind the interface. Quantifying the gas in this region is particularly important in expansion tubes because it is the location of the test gas. Unsteady simulations of the viscous flow in shock tubes are computationally expensive because they must follow features such as a shock wave over the length of the facility and simultaneously resolve the small length scales within the boundary layer. As a result, efficient numerical algorithms are required. The numerical approach of the present work is to solve the axisymmetric gas dynamic equations using an finite-volume formulation where the inviscid fluxes are computed with a upwind TVD scheme. Multiple species equations are included in the formulation so that finite-rate chemistry can be modeled. The simulations cluster grid points at the shock and interface and translate this clustered grid with these features to minimize numerical errors. The solutions are advanced at a CFL number of less than one based on the inviscid gas dynamics. To avoid limitations on the time step due to the viscous terms, these terms are treated implicitly. This requires a block tri-diagonal matrix inversion along each line of cells normal to the wall. The cost of this inversion is more than offset by the larger allowable time step. The source terms representing the finite-rate chemical kinetics are also treated implicitly. An algebraic turbulence model for compressible flow is used. The flow in a low pressure shock tube is computed and the results are compared with Mirels'analysis. The driven gas is nitrogen at 70 Pa, and the incident shock speed is approximately 2.9 km/sec so that there is little dissociation. The simulations include a laminar boundary layer and are run until the limiting flow regime is achieved. At this limit, the shock and interface travel at the same velocity because the amount of driven gas between these two features remains the same: the mass flow across the shock is equal to the mass of gas being entrained at the interface by the boundary layer. Simulations with several grids are presented to establish the grid independence of the solution, Good agreement is achieved between Mirels' correlations and the computations. This is expected since the flow conditions are chosen to be consistent with the assumptions used in Mirels' analysis. This comparison adds credibility to the numerical approach and highlights some of the differences between the theory and the detailed simulations. In addition, simulations of the HYPULSE expansion tube are presented for two operating conditions and the computations are compared to experimental data. The operating gas for both cases is nitrogen. One test condition is at a total enthalpy of 15.2 MJ/Kg and a relatively low pressure of 2 kPa. This case is characterized by a laminar boundary layer and significant chemical nonequilibrium. in the acceleration gas. The second test condition is at a total enthalpy of 10.2 MJ/Kg and a pressure of 38 kPa and is characterized by a turbulent boundary layer. The simulations compare well with experiment and reveal that the nonuniformity in pressure observed during the test time is related to variations in the boundary layer displacement thickness.

Wilson, Gregory J.↗