Search NASA⌕ Search

SEARCH · Search NASA

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 829 records · Page 46

Vertical transport of water in the Martian boundary layer

We are continuing our examination of the transport of H2O through the martian boundary layer, and we have written a one-dimensional numerical model of the exchange of H2O between the atmosphere and subsurface of Mars through the planetary boundary layer (PBL). Our goal is to explore the mechanisms of H2O exchange, and to elucidate the role played by the regolith in the local H2O budget. The atmospheric model includes effects of Coriolis, pressure gradient, and frictional forces for momentum, as well as radiation, sensible heat flux, and advection for heat. The model differs from Flasar and Goody by use of appropriate Viking-based physical constants and inclusion of the radiative effects of atmospheric dust. We specify the pressure gradient force or compute it from a simple slope model. The subsurface model accounts for conduction of heat and diffusion of H2O through a porous adsorbing medium in response to diurnal forcing. The model is initialized with depth-independent H2O concentrations (2 kg M(exp -3)) in the regolith, and a dry atmosphere. The model terminates when the atmospheric H2O column abundance stabilizes at 0.1 percent per sol.

Zent, Aaron P.↗

What Went Wrong: Explaining Counterexamples

Model checking, initially successful in the field of hardware design, has recently been applied to software. One of the chief advantages of model checking is the production of counterexamples demonstrating that a system does not satisfy a speci cation. However, it may require a great deal of human effort to extract the essence of an error from even a detailed source-level trace of a failing run. We use an automated method for nding multiple versions of an error (and similar executions that do not produce an error), and analyze these executions to produce a more succinct description of the key elements of the error. The description produced includes identi cation of portions of the source code crucial to distinguishing failing and succeeding runs, di erences in invariants between failing and non-failing runs, and information on the necessary changes in scheduling and environmental actions needed to cause suc- cessful runs to fail. In addition, this analysis allows a classi cation of errors by features such as whether they are purely concurrent (i.e. can be induced by changing only thread scheduling).

Groce, Alex↗

Finite element analyses of single particle crushing tests incorporating computed tomography imaging and damage mechanics

Fracture of individual grains within a granular assembly affect the strength and deformation behavior of granular materials under large stress states. Interparticle forces, leading to contact stresses and ultimately fracture initiation, are influenced by particle morphology. A numerical method that incorporates both particle fracture and morphology can provide a more accurate model of a granular material’s macro-scale response. This research investigates an approach to modeling fracture initiation and propagation within individual grains, utilizing high resolution X-ray computed tomography (CT) imaging to consider grain morphology and an explicit finite element code optimized for high performance computing which incorporates damage mechanics. Using precisely measured force-displacement curves from single particle crushing tests of manufactured quartz spheres, the method is validated, and the effects of fracture properties on grain fragmentation are determined. This approach is then used to simulate single particle crushing of Ottawa sand, with grain morphology incorporated through discretized CT images. Here the effect of grain orientation on the force and applied displacement at catastrophic splitting is also explored. The proposed discrete particle based finite element framework can be applied to modeling granular assemblies, considering the effect of individual particle shape and fracture on the assembly’s deformation response through well calibrated numerical simulations.

42 ENGINEERING↗

A comparison of limited-area energetic processes between observations and primitive equation model predictions

Energetic analyses of the NMC initial conditions and NMC six-layer primitive equation operational prediction model 12-hr forecast for a developing cyclone are presented. Consideration is given to the total kinetic energy, the energetics of the divergent and nondivergent flows and the baroclinic (vertical shear flow) and barotropic (vertical mean flow) components of the kinetic energy. It is found that the model initial conditions lose 10-15% of the kinetic energy at various levels compared to a limited-area multivariate statistical analysis of the observational data, leading to a decrease in the horizontal kinetic energy flux, a misrepresentation of the synoptic scale wave system in the 12-hr forecast. Similar results are obtained for the nondivergent flow, while the divergent flow energetics are not reproduced accurately by the model. The horizontal flux terms of the vertical mean and vertical shear energetics are also not found to be reproduced in the upper levels, although horizontal flux contributions to the baroclinic component are improved at middle and lower levels. Finally, vertical shear kinetic energy generation is found to be well represented in the model prediction, however kinetic energy conversion between vertical shear and mean flow is not reproduced in the lower layer.

Alpert, J. C.↗

Durability and Damage Tolerance of High Temperature Polymeric Composites

Modern durability and damage tolerance predictions for composite material systems rely on accurate estimates of the local stress and material states for each of the constituents, as well as the manner in which the constituents interact. In this work, an number of approaches to estimating the stress states and interactions are developed. First, an elasticity solution is presented for the problem of a penny-shaped crack in an N-phase composite material system opened by a prescribed normal pressure. The stress state around such a crack is then used to estimate the stress concentrations due to adjacent fiber fractures in composite materials. The resulting stress concentrations are then used to estimate the tensile strength of the composite. The predicted results are compared with experimental values. In addition, a cumulative damage model for fatigue is presented. Modifications to the model are made to include the effects of variable amplitude loading. These modifications are based upon the use of remaining strength as a damage metric and the definition of an equivalent generalized time. The model is initially validated using results from the literature. Also, experimental data from APC-2 laminates and IM7/K3B laminates are used in the model. The use of such data for notched laminates requires the use of an effective hole size, which is calculated based upon strain distribution measurements. Measured remaining strengths after fatigue loading are compared with the predicted values for specimens fatigued at room temperature and 350 F (177 C).

Case, Scott W.↗

Understanding Solar Eruptions with SDO/HMI Measuring Photospheric Flows, Testing Models, and Steps Towards Forecasting Solar Eruptions

The Solar Dynamics Observatory (SDO) is carrying the first full-disk imaging vector magnetograph, the Helioseismic and Magnetic Imager (HMI), into an inclined geosynchronous orbit. This magnetograph will provide nearly continuous measurements of photospheric vector magnetic fields at cadences of 90 seconds to 12 minutes with 1" resolution, precise pointing, and unfettered by atmospheric seeing. The enormous data stream of 1.5 Terabytes per day from SAO will provide an unprecedented opportunity to understand the mysteries of solar eruptions. These ground-breaking observations will permit the application of a new technique, the differential affine velocity estimator for vector magnetograms (DAVE4VM), to measure photospheric plasma flows in active regions. These measurements will permit, for the first time, accurate assessments of the coronal free energy available for driving CMEs and flares. The details of photospheric plasma flows, particularly along magnetic neutral-lines, are critical to testing models for initiating coronal mass ejections (CMEs) and flares. Assimilating flows and fields into state-of-the art 3D MHD simulations that model the highly stratified solar atmosphere from the convection zone to the corona represents the next step towards achieving NASA's Living with a Star forecasting goals of predicting "when a solar eruption leading to a CME will occur." Our presentation will describe these major science and predictive advances that will be delivered by SDO/HMI.

Schuck, Peter W.↗

Understanding Solar Eruptions with SDO/HMI Measuring Photospheric Flows, Testing Models, and Steps Towards Forecasting Solar Eruptions

The imminent launch of Solar Dynamics Observatory (SDO) will carry the first full-disk imaging vector magnetograph, the Helioseismic and Magnetic Imager (HMI), into an inclined geosynchronous orbit. This magnetograph will provide nearly continuous measurements of photospheric vector magnetic fields at cadences of 90 seconds to 12 minutes with I" resolution, precise pointing, and unfettered by atmospheric seeing. The enormous data stream of ~1.5 Terabytes per day from SDO will provide an unprecedented opportunity to understand the mysteries of solar eruptions. These ground-breaking observations will permit the application of a new technique, the differential affine velocity estimator for vector magnetograms (DAVE4VM), to measure photospheric plasma flows in active regions. These measurements will permit, for the first time, accurate assessments of the coronal free energy available for driving CMEs and flares. The details of photospheric plasma flows, particularly along magnetic neutral-lines, are critical to testing models for initiating coronal mass ejections (CMEs) and flares. Assimilating flows and fields into state-of-the art 3D MHD simulations that model the highly stratified solar atmosphere from the convection zone to the corona represents the next step towards achieving NASA's Living with a Star forecasting goals of predicting "when a solar eruption leading to a CME will occur." This talk will describe these major science and predictive advances that will be delivered by SDO /HMI.

Schuck, Peter W.↗

Operational and Technical Updates to the Object Reentry Survival Analysis Tool

The Object Reentry Survival Analysis Tool (ORSAT) has been used in the NASA Orbital Debris Program Office for over 25 years to estimate risk due to uncontrolled reentry of spacecraft and rocket bodies. Development over the last 3 years has included: a major change to the treatment of carbon fiber- and glass fiber-reinforced plastics (CFRP and GFRP, respectively); an updated atmospheric model; a new model for computing casualty area around an impacting debris object; and a newly-implemented scheme to determine the breakup altitude of a reentry object. Software also was written to automatically perform parameter sweeps in ORSAT to allow for uncertainty quantification and sensitivity analysis for components with borderline demisability. These updates have improved the speed and fidelity of the reentry analysis performed using ORSAT, and have allowed for improved engineering understanding by estimating the uncertainty for each component’s survivability. A statistical model for initial conditions captures the latitude bias in population density, a large improvement over the previous inclination-based latitude-averaged models. A sample spacecraft has been analyzed with standard techniques using ORSAT 6.2.1 and again using all the updated models; we will demonstrate the variation in the total debris casualty area and overall expectation of casualty.

Ostrom, C.↗

Primordial inhomogeneities in the expanding universe. I - Density and velocity distributions of galaxies in the vicinities of rich clusters

The density profiles and Hubble flow deviations in the vicinities of rich galaxy clusters are derived for a variety of models of initial density and velocity perturbations at the recombination epoch. The galaxy correlation function, measured with respect to the Abell clusters, is used to normalize the theoretical models. The angular scales of the required primordial inhomogeneities are calculated. It is found that the resulting density profiles around rich clusters are surprisingly insensitive to the shape of the initial perturbations and also to the cosmological density parameter, Omega. However, it is shown that the distribution of galaxy radial velocities can provide a possible means of deriving Omega.

Silk, J.↗

Downstream Separation of Formic Acid with Anion-Exchange Resin from Electrocatalytic Carbon Dioxide (CO 2 ) Conversion: Adsorption, Kinetics, and Equilibrium Modeling

The electrocatalytic carbon dioxide reduction reaction (CO 2 RR) that produces liquid formic acid (FA) is one potential route to minimize the CO 2 concentration in the atmosphere by utilizing CO 2 as a feedstock. However, one of the main disadvantages of this route is the high energy demand required for purification through conventional distillation. This method of purification has been shown to result in a large carbon footprint in life-cycle analysis (LCA). Therefore, finding an alternative purification technique that can integrate with an electrochemical CO 2 RR process with the potential to lower the carbon footprint would be of great interest. Under the experimental conditions of the CO 2 RR, FA remains in its basic form as the formate anion in the catholyte. Anion exchange resins with their basic properties can separate the formate anion from the catholyte mixture as an alternative to energy intensive distillation. In this work, three anion exchange resins with increasing basicity order: Amberlite IRA-96 < Amberlite IRA-910 < Ambersep 900 are tested to separate FA in the presence of aqueous solution of potassium bicarbonate catholyte with varying levels of resin and FA. Kinetics and equilibrium studies data for the FA adsorption are interpreted using several kinetics and isotherm models. The kinetics data fit well with a pseudo-first-order model at high initial FA concentrations and a pseudo-second-order model at low initial FA concentrations. The experimental data can be best explained with the Freundlich isotherm model. The fitted results show that strongly basic Ambersep 900 is more effective at separating FA, with a maximum adsorption capacity of 336.7 mg/g compared to Amberlite IRA-96 (275.2 mg/g) and Amberlite IRA-910 (209.2 mg/g) in neutral water. However, weakly basic Amberlite IRA-96 (110.8 mg/g) can separate FA more efficiently from potassium bicarbonate catholyte than Amberlite IRA-910 (99.9 mg/g) and Ambersep 900 (40.0 mg/g). As a result, competitive adsorption occurs in the presence of potassium bicarbonate catholyte, but the formic acid adsorption capacity is still sufficient for adequate separation.

20 FOSSIL-FUELED POWER PLANTS↗

Turbulence modeling for non-equilibrium flow

The work performed during this year has involved further assessment and extension of the k-epsilon-v(exp 2) model, and initiation of work on scalar transport. The latter is introduced by the contribution of Y. Shabany to this volume. Flexible, computationally tractable models are needed for engineering CFD. As computational technology has progressed, the ability and need to use elaborate turbulence closure models has increased. The objective of our work is to explore and develop new analytical frameworks that might extend the applicability of the modeling techniques. In past years the development of a method for near-wall modeling was described. The method has been implemented into a CFD code and its viability has been demonstrated by various test cases. Further tests are reported herein. Non-equilibrium near-wall models are needed for some heat transfer applications. Scalar transport seems generally to be more sensitive to non-equilibrium effects than is momentum transport. For some applications turbulence anisotropy plays a role and an estimate of the full Reynolds stress tensor is needed. We have begun work on scalar transport per se, but in this brief I will only report on an extension of the k-epsilon-v(exp 2) model to predict the Reynolds stress tensor.

Durbin, P. A.↗

Use of A-Train Aerosol Observations to Constrain Direct Aerosol Radiative Effects (DARE) Comparisons with Aerocom Models and Uncertainty Assessments

We describe a technique for combining multiple A-Train aerosol data sets, namely MODIS spectral AOD (aerosol optical depth), OMI AAOD (absorption aerosol optical depth) and CALIOP aerosol backscatter retrievals (hereafter referred to as MOC retrievals) to estimate full spectral sets of aerosol radiative properties, and ultimately to calculate the 3-D distribution of direct aerosol radiative effects (DARE). We present MOC results using almost two years of data collected in 2007 and 2008, and show comparisons of the aerosol radiative property estimates to collocated AERONET retrievals. Use of the MODIS Collection 6 AOD data derived with the dark target and deep blue algorithms has extended the coverage of the MOC retrievals towards higher latitudes. The MOC aerosol retrievals agree better with AERONET in terms of the single scattering albedo (ssa) at 441 nm than ssa calculated from OMI and MODIS data alone, indicating that CALIOP aerosol backscatter data contains information on aerosol absorption. We compare the spatio-temporal distribution of the MOC retrievals and MOC-based calculations of seasonal clear-sky DARE to values derived from four models that participated in the Phase II AeroCom model intercomparison initiative. Overall, the MOC-based calculations of clear-sky DARE at TOA over land are smaller (less negative) than previous model or observational estimates due to the inclusion of more absorbing aerosol retrievals over brighter surfaces, not previously available for observationally-based estimates of DARE. MOC-based DARE estimates at the surface over land and total (land and ocean) DARE estimates at TOA are in between previous model and observational results. Comparisons of seasonal aerosol property to AeroCom Phase II results show generally good agreement best agreement with forcing results at TOA is found with GMI-MerraV3. We discuss sampling issues that affect the comparisons and the major challenges in extending our clear-sky DARE results to all-sky conditions. We present estimates of clear-sky and all-sky DARE and show uncertainties that stem from the assumptions in the spatial extrapolation and accuracy of aerosol and cloud properties, in the diurnal evolution of these properties, and in the radiative transfer calculations.

aerosol↗

Comparison Between Surf and Multi-Shock Forest Fire High Explosive Burn Models

PAGOSA1 has several different burn models used to model high explosive detonation. Two of these, Multi-Shock Forest Fire and Surf, are capable of modeling shock initiation. Accurately calculating shock initiation of a high explosive is important because it is a mechanism for detonation in many accident scenarios (i.e. fragment impact). Comparing the models to pop-plot data give confidence that the models are accurately calculating detonation or lack thereof. To compare the performance of these models, pop-plots2 were created from simulations where one two cm block of PBX 9502 collides with another block of PBX 9502.

54 ENVIRONMENTAL SCIENCES↗

Development of Microreactor Automated Control System (MACS): Surrogate Plant-level Modeling and Control Algorithms Integration

This report discusses progress on the modeling and integration task as part of the development of Microreactor Automated Control System (MACS). What follows is a discussion of the software tools used to develop a surrogate microreactor plant-level model, MACS module development, and a summary of the findings from integration with a hardware-in-the-loop (HIL) setup at Idaho National Laboratory (INL). Oak Ridge National Laboratory worked with INL to understand available hardware at INL (such as existing control platforms). With this information and using a preliminary framework for MACS with software interfaces defined, surrogate models and initial automated control algorithms were created for integration into the hardware to provide a HIL demonstration platform for MACS. Available reduced order reactor models that leverage existing microreactor neutronics and thermal hydraulics models have been integrated in the surrogate plant-level model. The surrogate model has been evaluated for sensitivity of all parameters and the MACS software modules (including the surrogate model) have been demonstrated to interact with the hardware. The integration, however, pointed to the need for further improvements to the framework—and specifically improvements to the interfaces to ensure that the HIL simulator is capable of longer-term stable operation with MACS in the loop. These improvements will be the focus of future research.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

The formation of protostellar disks. I - 1 M(solar)

Hydrodynamical calculations of the collapse of an axisymmetric, rotating one solar mass protostellar cloud, including the effects of radiative transfer and radiative acceleration but without magnetic fields, are presented. The results include calculations of infrared protostellar spectra as a function of time and viewing angle. A numerical algorithm involving explicit nested grids is used to resolve the region of initial disk formation and at the same time to include the outer regions in the calculation. The central part of the protostar is modeled approximately. Initial conditions are systematically varied to investigate their influence on the evolution and final configuration of central star plus circumstellar disk. The initial state for the standard case is a centrally condensed molecular cloud core of one solar mass with a mean density of 8 x 10 exp -18 g/cu cm and a specific angular momentum at the outer edge of 7 x 10 exp 20 sq cm/s. The collapse is followed for 8 x 10 exp 4 yr, at which point 0.45 solar mass is contained in a rapidly rotating central object and most of the remainder in a surrounding equilibrium disk. The stability of this final structure is qualitatively analyzed.

Yorke, Harold W.↗

Real-time Adaptive Control Using Neural Generalized Predictive Control

The objective of this paper is to demonstrate the feasibility of a Nonlinear Generalized Predictive Control algorithm by showing real-time adaptive control on a plant with relatively fast time-constants. Generalized Predictive Control has classically been used in process control where linear control laws were formulated for plants with relatively slow time-constants. The plant of interest for this paper is a magnetic levitation device that is nonlinear and open-loop unstable. In this application, the reference model of the plant is a neural network that has an embedded nominal linear model in the network weights. The control based on the linear model provides initial stability at the beginning of network training. In using a neural network the control laws are nonlinear and online adaptation of the model is possible to capture unmodeled or time-varying dynamics. Newton-Raphson is the minimization algorithm. Newton-Raphson requires the calculation of the Hessian, but even with this computational expense the low iteration rate make this a viable algorithm for real-time control.

Haley, Pam↗

Mid-Air Range-Visual-Inertial Estimator Initialization for Micro Air Vehicles

Monocular Visual-Inertial Odometry (VIO) has become ubiquitous for navigation of autonomous Micro Air Vehicles (MAVs). Yet, state-of-the-art VIO is still very failure- prone, which can have dramatic consequences. To prevent this, VIO must be able to re-initialize in mid-air, either during a free fall or on a constant horizontal velocity trajectory after attitude control has been re-established. However, for both of these trajectories, the visual scale cannot be observed with VIO batch initializers because of the absence of acceleration change. We propose to use a small and lightweight laser range finder and a scene facet model to initialize vision-based navigation at the right scale under any motion condition and over any scene structure. This new range constraint is integrated into a visual-inertial bundle adjustment initializer. We demonstrate in simulation and real data how this approach can address mid-air state estimation failure in real-time and analyze the sensitivity of estimators to initial state errors.

Brockers, Roland↗