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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

A numerical study of transient, thermally-conductive solar wind

A numerical analysis of transient solar wind starting at the solar surface and arriving at 1 AU is performed by an implicit numerical method. The model hydrodynamic equations include thermal conduction terms for both steady and unsteady simulations. Simulation results show significant influence of thermal conduction on both steady and time-dependent solar wind. Higher thermal conduction results in higher solar wind speed, higher temperature, but lower plasma density at 1 AU. Higher base temperature at the solar surface gives lower plasma speed, lower temperature, but higher density at 1 AU. Higher base density, on the other hand, gives lower velocity, lower temperature, but higher density at 1 AU.

Han, S. M.↗
NASA NTRS1989

Helicopter flight control system design and evaluation for NOE operations using controller inversion techniques

A flight control system design technique is described that allows rapid evaluation of candidate control laws using computer simulated trajectories. The two-step technique first computes the optimum trajectory for the basic aircraft, and then reconstructs the pilot's control displacements necessary to fly this trajectory, given a particular control system architecture. Since the flight control system itself does not add to the acceleration potential of the helicopter, one may evaluate several proposed control laws in terms of the resulting stick inputs for the same maneuver. The method is illustrated through application to implicit and explicit model-following control laws designed for a battlefield helicopter with and without an auxiliary propulsor. It is shown that the technique is a valuable aid for selection of promising control laws for further simulator studies and flight investigations.

Mckillip, Robert M., Jr.↗
NASA NTRS1991

The effect of bleed configuration on shock/boundary layer interactions

The flow field characteristics are simulated numerically in an oblique shock wave/laminar boundary layer interaction for three different bleed configurations. The numerical solution for the flow field is obtained for the strong conservation-law form of the two-dimensional compressible Navier-Stokes equations using an implicit scheme. The computations model the flow in the interaction region and inside the bleed slot for an impinging oblique shock of sufficient strength to cause boundary layer separation on a flat plate boundary layer. The computed results are presented for a normal bleed slot at three different locations, upstream, across and downstream of the shock impingement point. The detailed flow characteristics in the interaction zone and inside the bleed slot are compared for the three cases. The resulting surface pressure and shear stress distributions are also compared for the three bleed cases and for the shock wave boundary layer interaction without bleed.

Hamed, A.↗
NASA NTRS1993

Parameterization of sub-grid scale convection

The following topics are discussed: an overview of the cumulus parameterization problem; interactions between explicit and implicit processes in mesoscale models; effects of model grid size on the cumulus parameterization problem; parameterizing convective effects on momentum fields in mesoscale models; differences between slantwise and vertical cumulus parameterization; experiments with different closure hypotheses; and coupling cumulus parameterizations to boundary layer, stable cloud, and radiation schemes.

Frank, William↗
NASA NTRS1996

Modeling, Monitoring and Fault Diagnosis of Spacecraft Air Contaminants

Progress and results in the development of an integrated air quality modeling, monitoring, fault detection, and isolation system are presented. The focus was on development of distributed models of the air contaminants transport, the study of air quality monitoring techniques based on the model of transport process and on-line contaminant concentration measurements, and sensor placement. Different approaches to the modeling of spacecraft air contamination are discussed, and a three-dimensional distributed parameter air contaminant dispersion model applicable to both laminar and turbulent transport is proposed. A two-dimensional approximation of a full scale transport model is also proposed based on the spatial averaging of the three dimensional model over the least important space coordinate. A computer implementation of the transport model is considered and a detailed development of two- and three-dimensional models illustrated by contaminant transport simulation results is presented. The use of a well established Kalman filtering approach is suggested as a method for generating on-line contaminant concentration estimates based on both real time measurements and the model of contaminant transport process. It is shown that high computational requirements of the traditional Kalman filter can render difficult its real-time implementation for high-dimensional transport model and a novel implicit Kalman filtering algorithm is proposed which is shown to lead to an order of magnitude faster computer implementation in the case of air quality monitoring.

Ramirez, W. Fred↗
NASA NTRS2008

Preliminary Design of a Ramjet for Integration with Ground-Based Launch Assist

This viewgraph presentation reviews the preliminary design of a ramjet for integration with a ground based launch assist. The reasons for the use of ground-based launch assist and the proposed mechanism for a system are reviewed. The use of a Optimal Trajectory by Implicit Simulation (OTIS), to model the flight and comparison with an actual rocket trajectory is given. The OTIS system is reviewed, The benefits of a launch assist system are analyzed concluding that a launch assist can provide supersonic speeds thus allowing ignition of ramjet without an onboard compressor. This means a further reduction in total launch weight. The Ramjet study is reviewed next. This included a review of the ONX simulations, the verification of the ONX results with the use of Holloman Sled experiment data as derived from the Feasibility of Ramjet Engine Test Capability on The Holloman AFB Sled Track. The conclusion was that the ONX system was not sufficient to meet the needs for the modeling required. The GECAT (Graphical Engine Cycle Analysis Tool) is examined. The results of the GECAT simulations was verified with data from Stataltex and D21 flights. The Next steps are: to create a GECAT Model of a launch assist ramjet, to adjust the geometry to produce the desired thrust, and to survey the ramjet's performance over a range of Mach numbers. The assumptions and requirements of a launch assist ramjet are given, and the acceptable flight regimes are reviewed.

Sayles, Emily L.↗
NASA NTRS2015

Damage Instability and Transition From Quasi-Static to Dynamic Fracture

In a typical mechanical test, the loading phase is intended to be a quasi-static process, while the failure and collapse is usually a dynamic event. The structural strength and modes of damage can seldom be predicted without accounting for these two aspects of the response. For a proper prediction, it is therefore essential to use tools and methodologies that are capable of addressing both aspects of responses. In some cases, implicit quasi-static models have been shown to be able to predict the entire response of a structure, including the unstable path that leads to fracture. However, is it acceptable to ignore the effect of inertial forces in the formation of damage? In this presentation we examine aspects of the damage processes that must be simulated for an accurate prediction of structural strength and modes of failure.

Davila, Carlos G.↗
DOE OSTI2021

Geodyn Material Library: Pseudocap models for dry porous tocks

This report describes the second edition of the Pseudocap Strength models for porous rocks implemented in GEODYN material library. The first model was developed in 2007 and calibrated for concrete. Then, the model parameters were calibrated based on triaxial tests reported for limestones and sandstones of various porosities. In these models some key parameters were chosen as functions of the reference porosity,Φ. Since then multiple modifications were implemented in the model, therefore, it has been recalibrated for some common porous materials such as limestones, sandstones, alluvium, tuffs and granite. Two types of models are described in this repot. The first type (called Pseudocap Model or PM) is for rocks from a specific location. Parameters were calibrated for several specific geologic materials. The second type (called Generic Pseudocap Model or GPM) is useful for the sites where only basic information (rock type, porosity) is available. Generic models include built-in correlations between porosities and other mechanical properties observed for certain rock types. The models of both types were validated by comparing not only to quasi-static triaxial tests for these materials but also to shock Hugoniot data and spherical explosion data for some materials. All models were derived in the frame of isotropic plasticity. They are designed to be used in explicit finite element/difference codes. Tangent stiffness tensor is not provided but can be calculated numerically for the model to be used in implicit finite element codes. For an isotropic material the stress can be decomposed into volumetric and deviatoric parts. The volumetric part is modeled using an Equation of state (EOS) which calculates the pressure and the bulk sound speed as functions of the internal specific energy and density. Here a simple, Mie-Gruneisen EOS is presented, but tabulated EOS (LEOS) provided by the library can be used as well. The stress is limited by the yield surface which depends on three invariants of the stress tensor and specific internal energy. In addition, to capture the strain-rate dependence a simple multiplier is used for the yield surface which depends on the equivalent plastic strain rate. The failure surface (the ultimate yield, Y f , defined later) is chosen in the Hoek-Brown form, commonly used in rock mechanics. It includes measurable parameters such as Unconfined Compressive Strength (UCS) as well as scale parameters characterizing the quality of the rock such as GSI (Geologic Strength Index). Thus, even though the model is calibrated for small samples it offers a way to extrapolate the strength to the field scale using geological characterization of the rock mass. The model captures effects of brittle-ductile transition in rocks by introducing a cap multiplier to the yield function. The rate of dilatancy (bulking) is proportional to the slope of the yield surface affected by the cap. Therefore, it takes place only at low confinements when the pressure is less than the brittle-ductile transition pressure, P BD . On the contrary, the porous compaction takes place at pressures higher than P BD . The cap moves as the porosity is compacted or new porosity is generated due to dilatancy. The porous compaction is modeled using an evolution equation which includes deviatoric stress so that the onset of compaction corresponds to the cap surface. The model captures effects of shear-enhanced compaction which is an important for porous rocks. Section 2 describes the modeling framework and Section 3 presents the model calibration procedure. Section 4 compares experimental data for various rocks versus model predictions. The model parameters used for this comparison are given in Appendix. The files with material constants are available with the latest GEODYN material library distribution.

58 GEOSCIENCES↗
DOE OSTI2021

Graphical Gaussian Process Regression Model for Aqueous Solvation Free Energy Prediction of Organic Molecules in Redox Flow Battery

The solvation free energy of organic molecules is a critical parameter in determining emergent properties such as solubility, liquid-phase equilibrium constants, and pKa and redox potentials in an organic redox flow battery. In this work, we present a machine learning (ML) model that can learn and predict the aqueous solvation free energy of an organic molecule using Gaussian process regression method based on a new molecular graph kernel. To investigate the performance of the ML model on electrostatic interaction, the nonpolar interaction contribution of solvent and the conformational entropy of solute in solvation free energy, three data sets with implicit or explicit water solvent models, and contribution of conformational entropy of solute are tested. We demonstrate that our ML model can predict the solvation free energy of molecules at chemical accuracy with a mean absolute error of less than 1 kcal/mol for subsets of the QM9 dataset and the Freesolv database. To solve the general data scarcity problem for a graph-based ML model, we propose a dimension reduction algorithm based on the distance between molecular graphs, which can be used to examine the diversity of the molecular data set. It provides a promising way to build a minimum training set to improve prediction for certain test sets where the space of molecular structures is predetermined.

25 ENERGY STORAGE↗
NASA NTRS1993

Understanding land surface evapotranspiration with satellite multispectral measurements

Quantitative use of remote multispectral measurements to study and map land surface evapotranspiration has been a challenging issue for the past 20 years. Past work is reviewed against process physics. A simple two-layer combination-type model is used which is applicable to both vegetation and bare soil. The theoretic analysis is done to show which land surface properties are implicitly defined by such evaporation models and to assess whether they are measurable as a matter of principle. Conceptual implications of the spatial correlation of land surface properties, as observed by means of remote multispectral measurements, are illustrated with results of work done in arid zones. A normalization of spatial variability of land surface evaporation is proposed by defining a location-dependent potential evaporation and surface temperature range. Examples of the application of remote based estimates of evaporation to hydrological modeling studies in Egypt and Argentina are presented.

Menenti, M.↗
NASA NTRS1983

Effects of expansion ratio on the calculation of parallel-walled backward-facing step flows - Comparison of four models of turbulence

This paper presents a numerical study of the effects of expansion ratio on two-dimensional separating and reattaching flows in plane backward-facing step geometries with parallel walls. Closure of the Reynolds equations was achieved by four different turbulence models: k-epsilon, 'modified' k-epsilon, algebraic stress, and 'modified' algebraic stress models. The k-epsilon model relates the Reynolds stresses to the mean rate of strain through the definition of an isotropic turbulent viscosity. The more advanced algebraic stress model calculates the stresses from implicit algebraic relationships containing the stresses themselves, the mean rate of strain, and the turbulent kinetic energy and its dissipation rate. 'Modified' versions of the models employ a new dissipation rate equation whose production term was made more sensitive to streamwise curvature effects. A new nonequilibrium wall function treatment proposed by Chieng and Launder (1980) was also incorporated into each model.

Sindir, M. M.↗
DOE OSTI2021

Numerical analyses of exponential time-differencing schemes for the solution of atmospheric models

In high-resolution numerical weather prediction models, fast-moving acoustic waves must be treated in a stable manner. Among the implicit–explicit (IMEX) schemes used for the solution of these models, horizontally explicit, vertically implicit (HEVI) methods show good stability and scalability on massively parallel machines. In this work, we present two classes of exponential time-differencing (ETD) methods for atmospheric models that use a HEVI splitting strategy, one being a three-stage method with the addition of artificial diffusion, the second based on a Strang splitting approach. Overall, the stability properties of the methods are analyzed and numerical examples are provided, which compare time-step restrictions and cost-to-accuracy ratios of the new methods with those for existing approaches.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1995

Use of midlatitude soil moisture and meteorological observations to validate soil moisture simulations with biosphere and bucket models

Soil moisture observations in sites with natural vegetation were made for several decades in the former Soviet Union at hundreds of stations. In this paper, the authors use data from six of these stations from different climatic regimes, along with ancillary meteorological and actinometric data, to demonstrate a method to validate soil moisture simulations with biosphere and bucket models. Some early and current general circulation models (GCMs) use bucket models for soil hydrology calculations. More recently, the Simple Biosphere Model (SiB) was developed to incorporate the effects of vegetation on fluxes of moisture, momentum, and energy at the earth's surface into soil hydrology models. Until now, the bucket and SiB have been verified by comparison with actual soil moisture data only on a limited basis. In this study, a Simplified SiB (SSiB) soil hydrology model and a 15-cm bucket model are forced by observed meteorological and actinometric data every 3 h for 6-yr simulations at the six stations. The model calculations of soil moisture are compared to observations of soil moisture, literally 'ground truth,' snow cover, surface albedo, and net radiation, and with each other. For three of the stations, the SSiB and 15-cm bucket models produce good simulations of seasonal cycles and interannual variations of soil moisture. For the other three stations, there are large errors in the simulations by both models. Inconsistencies in specification of field capacity may be partly responsible. There is no evidence that the SSiB simulations are superior in simulating soil moisture variations. In fact, the models are quite similar since SSiB implicitly has a bucket embedded in it. One of the main differences between the models is in the treatment of runoff due to melting snow in the spring -- SSiB incorrectly puts all the snowmelt into runoff. While producing similar soil moisture simulations, the models produce very different surface latent and sensible heat fluxes, which would have large effects on GCM simulations.

Robock, Alan↗
NASA NTRS1989

A locally implicit finite element analysis algorithm for 3-D compressible flows

The paper outlines the research results on the development of a finite-element analysis procedure for three-dimensional compressible viscous flows. Some new algorithmic features of the computational model, such as a locally implicit solution scheme which avoids the necessity of computation and assembly of global matrices. The new algorithm makes it applicable to large-scale three-dimensional flow problems of practical interest, for example, the analysis of flows in Space Shuttle Main Engine (SSME).

Reddy, K. C.↗
NASA NTRS1976

Comments on the role of diagonal dominance in implicit difference methods

Numerical tests were made for a model of the Navier-Stokes equations using a second-order accurate implicit scheme which guarantees diagonal dominance. The results suggest that the failure of implicit methods using large marching steps may not always be attributed to the lack of diagonal dominance in the coefficient matrix. In some cases the failure may be caused by a nonlinear instability associated with the solution method.

Rudy, D. H.↗
NASA NTRS1979

A structural model of the adaptive human pilot

A compensatory tracking model of the human pilot is offered which attempts to provide a more realistic representation of the human's signal processing structure than that which is exhibited by pilot models currently in use. Two features of the model distinguish it from other representations of the human pilot. First, proprioceptive information from the control stick or manipulator constitutes one of the major feedback paths in the model, providing feedback of vehicle output rate due to control activity. Implicit in this feedback loop is a model of the vehicle dynamics which is valid in and beyond the region of crossover. Second, error-rate information is continuously derived and independently but intermittently controlled. An output injected remnant model is offered and qualitatively justified on the basis of providing a measure of the effect of inaccuracies such as time variations in the pilot's internal model of the controlled element dynamics. The data from experimental tracking tasks involving five different controlled element dynamics and one nonideal viewing condition were matched with model generated describing functions and remnant power spectral densities.

Hess, R. A.↗
NASA NTRS2004

Improved Ionospheric Electrodynamic Models and Application to Calculating Joule Heating Rates

Improved techniques have been developed for empirical modeling of the high-latitude electric potentials and magnetic field aligned currents (FAC) as a function of the solar wind parameters. The FAC model is constructed using scalar magnetic Euler potentials, and functions as a twin to the electric potential model. The improved models have more accurate field values as well as more accurate boundary locations. Non-linear saturation effects in the solar wind-magnetosphere coupling are also better reproduced. The models are constructed using a hybrid technique, which has spherical harmonic functions only within a small area at the pole. At lower latitudes the potentials are constructed from multiple Fourier series functions of longitude, at discrete latitudinal steps. It is shown that the two models can be used together in order to calculate the total Poynting flux and Joule heating in the ionosphere. An additional model of the ionospheric conductivity is not required in order to obtain the ionospheric currents and Joule heating, as the conductivity variations as a function of the solar inclination are implicitly contained within the FAC model's data. The models outputs are shown for various input conditions, as well as compared with satellite measurements. The calculations of the total Joule heating are compared with results obtained by the inversion of ground-based magnetometer measurements. Like their predecessors, these empirical models should continue to be a useful research and forecast tools.

Weimer, D. R.↗
DOE OSTI2026

Wake-Resolving Acoustic Tomography: Advances through Numerical Covariance Methods

Acoustic tomography offers path-integrated measurements of atmospheric velocity and temperature fluctuations with high spatial resolution. Classical implementations of time-dependent stochastic inversion rely on homogeneous, isotropic covariance models that are poorly suited to the anisotropic structure of wind turbine wakes. By directly estimating heterogeneous covariances from large-eddy simulations (LESs) into the time-dependent stochastic inversion operator, we relax implicit assumptions in the analytical models used historically. Retrievals using these LES-informed models improve agreement with true fields in variance, turbulent kinetic energy, and spectral content compared to analytical and precursor-based covariance models. The results indicate that LES-informed covariance models can enhance the accuracy of acoustic tomography retrievals in complex, anisotropic flows such as wind turbine wakes in some cases and highlight instances where analytical models still offer competitive performance, despite their simplifying assumptions.

17 WIND ENERGY↗