Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model initialization”

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 1,351 records · Page 75
All sourcesNASA NTRSDOE OSTI
NASA NTRS1980

Whole planet cooling and the radiogenic heat source contents of the earth and moon

Thermal evolution models based on subsolidus whole mantle convection which indicate that the surface heat flows of the earth and the moon do not necessarily provide good measures of the total amounts of radioactives in these bodies have been constructed. These models assume an initially hot state, but with a wide variety of choices for the parameters characterizing the rheology and convective vigor. All models are constrained to be consistent with present-day surface heat fluxes, and many of the terrestrial models are consistent with the mantle viscosities indicated by postglacial rebound. In the lunar models, heat generation is typically only 70-80% of the surface heat flow, even with allowance for the strong near-surface enhancement of radioactives. Despite the simplicity of these models, the persistence of a significant difference between heat generation and heat output indicates that this difference is real and should be incorporated in geochemical modeling of planets.

Schubert, G.↗
NASA NTRS1984

Stability of finite difference models containing two boundaries or interfaces

The stability of finite difference models of hyperbolic initial boundary value problems is connected with the propagation and reflection of parasitic waves. Wave propagation ideas are applied to models containing two boundaires or interfaces, where repeated reflection of trapped wave packets is a potential new source of instability. Various known instability phenomena are accounted for in a unified way. Results show: (1) dissipativity does not ensure stability when three or more formulas are concatenated at a boundary or internal interface; (2) algebraic GKS instabilities can be converted by a second boundary to exponential instabilities only when an infinite numerical reflection coefficient is present; and (3) GKS-stability and P-stability can be established in certain problems by showing that all numerical reflection coefficients have modulus less than 1.

Trefethen, L. N.↗
NASA NTRS1985

Stability of finite-difference models containing two boundaries or interfaces

The stability of finite difference models of hyperbolic initial boundary value problems is connected with the propagation and reflection of parasitic waves. Wave propagation ideas are applied to models containing two boundaries or interfaces, where repeated reflection of trapped wave packets is a potential new source of instability. Various known instability phenomena are accounted for in a unified way. Results show: (1) dissipativity does not ensure stability when three or more formulas are concatenated at a boundary or internal interface; (2) algebraic GKS instabilities can be converted by a second boundary to exponential instabilities only when an infinite numerical reflection coefficient is present; and (3) GKS-stability and P-stability can be established in certain problems by showing that all numerical reflection coefficients have modulus less than 1.

Trefethen, L. N.↗
NASA NTRS1996

The Aggregate Representation of Terrestrial Land Covers Within Global Climate Models (GCM)

This project had four initial objectives: (1) to create a realistic coupled surface-atmosphere model to investigate the aggregate description of heterogeneous surfaces; (2) to develop a simple heuristic model of surface-atmosphere interactions; (3) using the above models, to test aggregation rules for a variety of realistic cover and meteorological conditions; and (4) to reconcile biosphere-atmosphere transfer scheme (BATS) land covers with those that can be recognized from space; Our progress in meeting these objectives can be summarized as follows. Objective 1: The first objective was achieved in the first year of the project by coupling the Biosphere-Atmosphere Transfer Scheme (BATS) with a proven two-dimensional model of the atmospheric boundary layer. The resulting model, BATS-ABL, is described in detail in a Masters thesis and reported in a paper in the Journal of Hydrology Objective 2: The potential value of the heuristic model was re-evaluated early in the project and a decision was made to focus subsequent research around modeling studies with the BATS-ABL model. The value of using such coupled surface-atmosphere models in this research area was further confirmed by the success of the Tucson Aggregation Workshop. Objective 3: There was excellent progress in using the BATS-ABL model to test aggregation rules for a variety of realistic covers. The foci of attention have been the site of the First International Satellite Land Surface Climatology Project Field Experiment (FIFE) in Kansas and one of the study sites of the Anglo-Brazilian Amazonian Climate Observational Study (ABRACOS) near the city of Manaus, Amazonas, Brazil. These two sites were selected because of the ready availability of relevant field data to validate and initiate the BATS-ABL model. The results of these tests are given in a Masters thesis, and reported in two papers. Objective 4: Progress far exceeded original expectations not only in reconciling BATS land covers with those that can be recognized from space, but also in then applying remotely-sensed land cover data to map aggregate values of BATS parameters for heterogeneous covers and interpreting these parameters in terms of surface-atmosphere exchanges.

Shuttleworth, W. James↗
NASA NTRS1985

Vibration response of geometrically nonlinear elastic beams subjected to pulse and impact loading

An understanding of the response of composite laminates to foreign object impact is of fundamental importance to the design and the accurate evaluation of aircraft and aerospace structures. A characterization of the basic deformation response of the system is needed for the accurate and consistent modeling of the damage caused by the impact. The response to low velocity impact of foreign objects (less than 10 m/sec at 10 J of impact energy or less) is of particular interest because it corresponds to tool drop problems. An investigation of the vibration response of a simple square beam was conducted, taking into account both sharp initial velocity pulse conditions and simulated impact conditions. Attention is given to nonlinear beam theory, boundary conditions and materials modeling, initial conditions and impact modeling, and a solution of the obtained system of partial differential equations by the finite difference approach.

Moyer, E. T., Jr.↗
DOE OSTI2020

Anomalous Transport in Three-Dimensional Discrete Fracture Networks: Interplay Between Aperture Heterogeneity and Injection Modes

In this work, we study how the interplay between fracture aperture heterogeneity and tracer injection mode controls fluid flow and tracer transport in three-dimensional (3D) discrete fracture networks (DFNs). The direct 3-D DFN simulations show that tracer injection mode has substantial effects on tracer spreading across all levels of aperture heterogeneity. The key controlling factor for effective transport is the initial Lagrangian velocity distribution, which is determined by the interplay between injection mode and aperture heterogeneity. The fundamental difference between initial Lagrangian velocity distribution and domain-scale Eulerian velocity distribution plays a vital role in determining anomalous transport. We effectively capture the observed anomalous transport using an upscaled transport model that incorporates initial velocity distribution, stationary velocity distribution, velocity correlation length, and average advective tortuosity. With the upscaled transport model, we accurately capture the evolution of Lagrangian velocity distribution and predict longitudinal spreading in 3-D DFN.

58 GEOSCIENCES↗
NASA NTRS1989

Molecular modeling of calmodulin: a comparison with crystallographic data

Two methods of side-chain placement on a modeled protein have been examined. Two molecular models of calmodulin were constructed that differ in the treatment of side chains prior to optimization of the molecule. A virtual bond analysis program developed by Purisima and Scheraga was used to determine the backbone conformation based on 2.2 angstroms resolution C alpha coordinates for the molecules. In the first model, side chains were initially constructed in an extended conformation. In the second model, a conformational grid search technique was employed. Calcium ions were treated explicitly during energy optimization using CHARMM. The models are compared to a recently published refined crystal structure of calmodulin. The results indicate that the initial choices for side-chains, but also significant effects on the main-chain conformation and supersecondary structure. The conformational differences are discussed. Analysis of these and other methods makes possible the formulation of a methodology for more appropriate side-chain placement in modeled proteins.

NASA Discipline Exobiology↗
DOE OSTI2020

Smart-PGSim: Using Neural Network to Accelerate AC-OPF Power Grid Simulation

In this work we address the problem of accelerating complex power-grid simulation through machine learning ( ML). Specifically, we develop a framework, Smart-PGSim,which generates multitask-learning (MTL) neural network (NN)models to predict the initial values of variables critical to the problem convergence. MTL models allow information sharing when predicting multiple dependent variables while including customized layers to predict individual variables. We show that,to achieve the required accuracy, it is paramount to embed domain-specific constraints derived from the specific power-grid components in the MTL model. Smart-PGSim then employs the predicted initial values as a high-quality initial condition for the power-grid numerical solver (warm start), resulting in both higher performance compared to state-of-the-art solutions while maintaining the required accuracy. Smart-PGSim brings 2.60×speedup on average (up to 3.28×) computed over 10,000 problems, without losing solution optimality.

machine learning, neural networks↗
DOE OSTI2020

A new soil mixing layer model for simulating conservative solute loss from initially saturated soil to surface runoff

Surface water pollution due to solute loss from soil to surface runoff is a serious threat to the environment, and it is of great significance to predict the solute loss. While soil mixing layer theory has been widely used to predict the solute loss, it has been observed that the conventional soil mixing layer models with constant mixing layer depth underestimate solute concentration in surface runoff, and thus underestimate the solute loss after surface runoff starts. This study introduced a new soil mixing layer model with a new concept that the mixing layer depth increases downward with time after surface runoff occurs. This new concept was implemented by adding a mixing front between the soil mixing layer and its underlying soil layer, which is conceptually similar to the wetting front of soil infiltration given by the Green-Ampt equation. This is equivalent to treating the soil mixing layer as a deforming control volume. A mechanistic explanation for the new model was presented. This study also developed a closed form solution of the new model. The new model was evaluated by simulating five laboratory experiments conducted by two groups of researchers. The first two experiments have ponding water, and consider different drainage and runoff conditions at both early and late runoff times. The other three experiments do not have ponding water, and consider different infiltration conditions at the bottom of soil boxes. Simulation results of the five experiments show that the new model resolves the underestimation problem of conventional soil mixing layer models. The new model is subject to a total of six limits such as that the new model only considers conservative solute. Three limits are related to the assumptions used to derive the closed form solution, and they may be alleviated in a future research as briefly discussed at the end of the paper.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A Note on the Probability of Initiation Problem

The purpose of this note is to draw attention to Pál for his derivation of the Pál-Bell equation which is used by LLNL and LANL to model Probability of Initiation (POI) problems. Both laboratories credit the equation to Bell and Lee. Although the equation formulated by Bell and Lee and the equation derived by Pál appear to be the same, the equation by Bell and Lee is flawed. Scattering, which was absent in Bell’s original formulation of a POI problem, was later introduced by Bell and Lee by replacing a fission term with a scattering term. Such a replacement, however, leads to an incorrect average number $\overline{v}$ of neutrons that is produced in a fission event. By approaching a POI problem from probabilistic point of view, Pál formulates the problem as a branching process that i) describes the growth of a population of neutrons in a fissile medium, and ii) predicts $\overline{v}$ correctly. There are five reasons for writing this report. The first is to clarify the failing in Bell and Lee’s derivation. The second is to establish the origin of the equation used by the laboratories. The third is to elucidate the fixed point technique formulated by Mokhtar-Kharroubi to solve the Pál-Bell equation. The fourth is to derive the α eigenvalue that is associated with a super-critical situation. The fifth is to describe in the Appendix a second order positivity preserving technique that accelerates the convergence of the first order fixed point method developed by Mokhtar-Kharroubi.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1987

Comments on the challenge of using mesoscale data in mesoscale numerical models

The dangers of addressing the initialization issues for limited-area mesoscale models by extending the lessons learned during the development of global analysis and prediction systems are discussed. Lack of impact with data inserts at one time suggests that the lateral boundary conditions imposed on the limited-area models might force the model simulation toward a preferred solution, work against the new data being inserted into the model and, therefore, limit the potential impact that this data can have on the model system. The second potential pitfall involves the imposition of balance constraints on the data that are being inserted into the model to compute winds from temperature data and/or temperature from wind data.

Uccellini, Louis W.↗
NASA NTRS1989

Incorporation of a progressive failure analysis method in the CSM testbed software system

Analysis of the postbuckling behavior of composite shell structures pose many difficult and challenging problems in the field of structural mechanics. Current analysis methods perform well for most cases in predicting the postbuckling response of undamaged components. To predict component behavior accurately at higher load levels, the analysis must include the effects of local material failures. The CSM testbed software system is a highly modular structural analysis system currently under development at Langley Research Center. One of the primary goals of the CSM testbed is to provide a software environment for the development of advanced structural analysis methods and modern numerical methods which will exploit advanced computer architecture such as parallel-vector processors. Development of a progressive failure analysis method consists of the design and implementation of a processor which will perform the ply-level progressive failure analysis and the development of a geometrically nonlinear analysis procedure which incorporates the progressive failure processor. Regarding the development of the progressive failure processor, two components are required: failure criteria and a degradation model. For the initial implementation, the failure criteria of Hashin will be used. For a matrix failure which typically indicates the development of transverse matrix cracks, the ply properties will be degraded. Work to date includes the design of the progressive failure analysis processor and initial plans for the controlling geometrically nonlinear analysis procedure. The implementation of the progressive failure analysis has begun. Access to the model database and the Hashin failure criteria are completed. Work is in progress on the input/output operations for the processor related data and the finite element model updating procedures. In total the progressive failure processor is approximately one-third complete.

Arenburg, Robert T.↗
NASA NTRS2004

Advances in Scientific Balloon Thermal Modeling

The National Aeronautics and Space Administration's Balloon Program office has long acknowledged that the accurate modeling of balloon performance and flight prediction is dependant on how well the balloon is thermally modeled. This ongoing effort is focused on developing accurate balloon thermal models that can be used to quickly predict balloon temperatures and balloon performance. The ability to model parametric changes is also a driver for this effort. This paper will present the most recent advances made in this area. This research effort continues to utilize the "Thrmal Desktop" addition to AUTO CAD for the modeling. Recent advances have been made by using this analytical tool. A number of analyses have been completed to test the applicability of this tool to the problem with very positive results. Progressively detailed models have been developed to explore the capabilities of the tool as well as to provide guidance in model formulation. A number of parametric studies have been completed. These studies have varied the shape of the structure, material properties, environmental inputs, and model geometry. These studies have concentrated on spherical "proxy models" for the initial development stages and then to transition to the natural shaped zero pressure and super pressure balloons. An assessment of required model resolution has also been determined. Model solutions have been cross checked with known solutions via hand calculations. The comparison of these cases will also be presented. One goal is to develop analysis guidelines and an approach for modeling balloons for both simple first order estimates and detailed full models. This papa presents the step by step advances made as part of this effort, capabilities, limitations, and the lessons learned. Also presented are the plans for further thermal modeling work.

Bohaboj, T.↗
DOE OSTI2024

Complementing the CCS Class VI Well Permit Process with DOE-NETL's SMART Initiative Tools and Workflows

This is a presentation on model explorer developed under SMART initiative Task 2. Our team will present the current progress of the model explorer in using machine learning models to accelerate CCS project at GWPC meeting. Model explorer bring new capabilities, (fast, Realtime, and accurate) that can help CCS stakeholders including regulatory agencies, public and site operators make faster decisions and process information and data.

Hosseini, Seyyed↗
DOE OSTI2021

Effects of finite element discretization and model simplification on calculations of ductile failure initiation

The finite element method is a scheme to discretize the infinite number of degrees of freedom in continuum-level problems down to a finite number of degrees of freedom. This discretization is done in conjunction with methods that also reduce the field differential equations to sets of algebraic ones that can be solved by arithmetical operations. Therefore, solutions attained by finite element models are approximations to the exact solutions of the field equations.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1991

Constraints on cosmological models from the IRAS QDOT survey results

The results of the recent IRAS survey are used to explore constraints on the nature of the Galaxy density field. Realistic simulations of the survey made from N-body simulations of the n = -1 power-law power spectrum and the Omega = 1 CDM models show that the power-law model is basically consistent with all QDOT variance results, while the CDM model has a smaller variance than the observation at the 20/h Mpc Gaussian smoothing scale. It is concluded that the QDOT results are best explained if the power spectrum of the Galaxy density field has a slope of about -1 near a wavelength of about 100/h Mpc, while the CDM model has a slope of about 0 at that scale. The possibility of marginal non-Gaussian behavior at the 20/h Mpc scale reported by Saunders et al. (1990) is examined for gravitational instability models with Gaussian initial conditions. It is argued that the observed (delta exp 3) versus (delta exp 2) in redshift space is not what these models would predict.

Park, Changbom↗
NASA NTRSDate not supplied

Heliospheric energetic particle transport: Analysis of near-field-aligned particle propagation for SEP events observed by Wind and MAVEN

Magnetic field alignments of spacecraft over large distances in the heliosphere are rare and are usually very limited in duration. Cruise phases of planetary transfers, however, are an exception to this rule, given the Hohmann-Parker effect. The transfer of the MAVEN s/c in 2014 is one such example. Multiple (~10) solar particle events occurred and were detected at both MAVEN and the Wind s/c, originating from solar activity near the foot points of both s/c. We show initial analysis results of the data collected by the solar wind and energetic particle instruments on both s/c while they were more than 0.2 AU apart, but practically Parker field aligned. Using a 1D model, we present initial simulations in qualitative agreement with energetic electron measurements. Next step is to implement the 2D model with approximate particle release from the Sun, and transport durations between Sun, Earth, and MAVEN. We will discuss implications for this data-model comparison, including the possibility to constrain particle scattering inside 1 AU, as well as its radial dependence between Earth and MAVEN.

Solar Energetic Particles↗