Search NASA⌕ Search

SEARCH · Search NASA

Results for “Initial conditions problem”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Development of an integrated BEM approach for hot fluid structure interaction

The development of a boundary element formulation for the study of hot fluid-structure interaction in earth-to-orbit engine hot section components is described. The initial primary thrust of the program to date was directed quite naturally toward the examination of fluid flow, since boundary element methods for fluids are at a much less developed state. This required the development of integral formulations for both the solid and fluid, and some preliminary infrastructural enhancements to a boundary element code to permit coupling of the fluid-structure problem. Boundary element formulations are implemented in two dimensions for both the solid and the fluid. The solid is modeled as an uncoupled thermoelastic medium under plane strain conditions, while several formulations are investigated for the fluid. For example, both vorticity and primitive variable approaches are implemented for viscous, incompressible flow, and a compressible version is developed. All of the above boundary element implementations are incorporated in a general purpose two-dimensional code. Thus, problems involving intricate geometry, multiple generic modeling regions, and arbitrary boundary conditions are all supported.

Dargush, Gary F.↗
NASA NTRS1995

Existence of frozen-in coordinate systems

The 'frozen-in' coordinate systems were first introduced in the works on 'reconnection' and 'magnetic barrier' theories (see review by M.l.Pudovkin and V.S.Semenov, Space Sci. Rev. 41,1 1985). The idea was to utilize the mathematical apparatus developed for 'general relativity' theory to simplify obtaining solutions to the ideal MHD equations set. Magnetic field (B), plasma velocity (v), and their vector product were used as coordinate vectors. But there exist no stationary solutions of ideal MHD set that satisfies the required boundary conditions at infinity (A.D.Chertkov, Solar Wind Seven Conf.,Pergamon Press,1992,165) having non-zero vector product of v and B where v and B originate from the same sphere. The existence of a solution is the hidden mine of the mentioned theories. The solution is constructed in the coordinate system, which is unknown and indeterminate before obtaining this solution. A substitution of the final solution must be done directly into the initial MHD set in order to check the method. One can demonstrate that 'solutions' of Petschek's problem, obtained by 'frozen-in' coordinate systems, does not satisfy just the 'frozen-in' equation, i.e. induction equation. It stems from the fact that Petschek's 're-connection' model, treated as a boundary problem, is over determined. This problem was incorrectly formulated.

Chertkov, A. D.↗
NASA NTRS2003

Real-Time Attitude Independent Three Axis Magnetometer Calibration

In this paper new real-time approaches for three-axis magnetometer sensor calibration are derived. These approaches rely on a conversion of the magnetometer-body and geomagnetic-reference vectors into an attitude independent observation by using scalar checking. The goal of the full calibration problem involves the determination of the magnetometer bias vector, scale factors and non-orthogonality corrections. Although the actual solution to this full calibration problem involves the minimization of a quartic loss function, the problem can be converted into a quadratic loss function by a centering approximation. This leads to a simple batch linear least squares solution. In this paper we develop alternative real-time algorithms based on both the extended Kalman filter and Unscented filter. With these real-time algorithms, a full magnetometer calibration can now be performed on-orbit during typical spacecraft mission-mode operations. Simulation results indicate that both algorithms provide accurate integer resolution in real time, but the Unscented filter is more robust to large initial condition errors than the extended Kalman filter. The algorithms are also tested using actual data from the Transition Region and Coronal Explorer (TRACE).

Crassidis, John L.↗
DOE OSTI2022

Fast Reactor Physics Model Verification Studies using ARC and PyARC Workflows

PyARC was recently developed at Argonne National Laboratory to automate many of the tasks required in the ARC (Argonne Reactor Computation) fast reactor simulation workflow, from input file generation, code execution, data transfer between ARC codes, and output postprocessing. PyARC will likely be the path forward to train new users of the ARC codes with the goal of wide adoption by the national laboratories, academia, and industry. In particular, for the ANL-JAEA collaboration under the Civil Nuclear Working Group (CNWG) project agreement NE-01, PyARC will be used to model the Joyo and EBR-II reactors for comparisons with measured data and calculated results from JAEA (Task 3: Fast Reactor Fuel and Core). As an additional avenue for verification and validation, this report investigates the use of PyARC towards a variety of existing ARC-based reactor models, in order to understand its efficacy in replicating the behavior of base ARC codes and better understand any limitations within modeling realistic fast reactor problems. To this end, PyARC was used to model the Joyo MKI, RBEC Benchmark-M, PRISM Mod-B, and EBR-II Run 138B cores, and its results were compared to those from existing ARC-based models. It was found that for hexagonal-based geometries PyARC was able to replicate the behavior of ARC codes to within 10 pcm for small reactor cores, and ~150pcm difference in eigenvalue for larger cores. These discrepancies are attributed primarily to differences in local mesh refinement options between ARC and PyARC, which currently cannot be resolved with PyARC’s latest version (1.6.0). In some of these cases, PyARC was used to model steady-state problems with initial core compositions originating from a prior REBUS depletion calculation. While PyARC was not designed to support such steady-state calculations, workarounds were applied to replicate the behavior of ARC-based calculations as closely as possible. Thus, these results demonstrate the wide extent to which they can be applied to fast reactor problems while still providing immense benefit to the user in terms of automating and standardizing common routines within the fast reactor analysis workflow. This study concluded that PyARC will be suitable for modeling the steady-state conditions of the EBR-II and Joyo fast reactors as part of the CNWG project agreement.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

New Virtual Test Bed Capabilities: Virtual DOME Model and New Updates to Repository

The Department of Energy (DOE) Office of Nuclear Energy National Reactor Innovation Center accelerates the deployment of novel reactor concepts by establishing both physical and virtual spaces for building and testing various components, systems, and complete pilot plants. The Virtual Test Bed represents the virtual arm of the National Reactor Innovation Center and is a joint effort with the DOE Nuclear Energy Advanced Modeling and Simulation Program. The Virtual Test Bed mission is to accelerate the deployment of advanced reactors by facilitating the adoption of cutting-edge DOE advanced modeling and simulation tools to design, evaluate, and license reactors. This is primarily achieved by storing example challenge problems in an externally available repository and by developing models to fill the M&S gaps needed for potential demonstrators. Activities conducted this fiscal year focused on developing of a Demonstration of Microreactor Experiments shield model to help accelerate the confirmatory analysis required for the reactor demonstration. This model and workflow will allow developers to leverage advanced modeling and simulation tools to ensure their reactor demonstration concept will meet dose requirements and that the surrounding shield will stay within concrete temperature limits during steady-state and transient operation conditions. An initial model has been developed to evaluate the temperature distribution in the concrete shield during steady-state operation, including neutron and gamma heating effects. Various modeling strategies have been examined to understand their applicability and limitations with different reactor designs to make the workflow as reactor-agnostic as possible and computationally effective to maximize its usability. In addition to describing the Demonstration of Microreactor Experiments shield model and associated results, this report summarizes other accomplishments regarding repository maintenance and improvement and new external models hosted on the repository.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS2006

On the design of a science orbit about Europa

A science mission about Europa requires high-inclination low-altitude orbits. However, perturbations of Jupiter on the orbiter result in instability. Previous approaches to maximize the lifetime of the orbiter use the doubly averaged problem. We work with the unaveraged equations and find unstable periodic orbits with long lifetimes. These low-altitude repeat ground track solutions exist at all inclinations, making them suitable for mapping missions. The governing dynamics include Hill's model and a Europa gravity field based on synchronous moon theory. Inclusion of additional gravity terms is trivial to the solution method, and for the case of J3, we find a marginal impact on orbit lifetime. The science orbits are found to last on the order of 1 year when the initial conditions are achieved to 11 significant digits and 4 months when only 3 significant digits are achieved. Finally, we demonstrate that the solutions are robust in a realistic ephemeris model, finding average lifetimes of 3 to 4 months for wide range of initial conditions with peak lifetimes of up to 6 months.

unstable↗
NASA NTRS2008

Predicting Chandra CCD Degradation with the Chandra Radiation Model

Not long after launch of the Chandra X-Ray Observatory, it was discovered that the Advanced CCD Imaging Spectrometer (ACIS) detector was rapidly degrading due to radiation. Analysis by Chandra personnel showed that this degradation was due to 10w energy protons (100 - 200 keV) that scattered down the optical path onto the focal plane. In response to this unexpected problem, the Chandra Team developed a radiation-protection program that has been used to manage the radiation damage to the CCDs. This program consists of multiple approaches - scheduled sating of the ACIS detector from the radiation environment during passage through radiation belts, real-time monitoring of space weather conditions, on-board monitoring of radiation environment levels, and the creation of a radiation environment model for use in computing proton flux and fluence at energies that damage the ACIS detector. This radiation mitigation program has been very successful. The initial precipitous increase in the CCDs' charge transfer inefficiency (CTI) resulting from proton damage has been slowed dramatically, with the front-illuminated CCDS having an increase in CTI of only 2.3% per year, allowing the ASIS detector's expected lifetime to exceed requirements. This paper concentrates on one aspect of the Chandra radiation mitigation program, the creation of the Chandra Radiation Model (CRM). Because of Chandra's highly elliptical orbit, the spacecraft spends most of its time outside of the trapped radiation belts that present the severest risks to the ACIS detector. However, there is still a proton flux environment that must be accounted for in all parts of Chandra's orbit. At the time of Chandra's launch there was no engineering model of the radiation environment that could be used in the outer regions of the spacecraft's orbit, so the CRM was developed to provide the flux environment of 100 - 200 keV protons in the outer magnetosphere, magnetosheath, and solar wind regions of geospace. This presentation describes CRM, its role in Chandra operations, and its prediction of the ACIS CTI increase.

Minow, Joseph I.↗
DOE OSTI2021

Real time trajectory optimization for hybrid energy management utilizing connected information technologies

The disclosure is directed to solving a full trajectory optimization problem in real-time for a hybrid electric vehicle (HEV) such that future driving conditions and energy usage may be fully considered in determining optimal engine energy usage and battery energy usage in real-time during a trip. An electronic control unit of the HEV may be configured to: receive route information for a route to be driven by the HEV; and after receiving the route information, iterating the operations of: measuring a current state of charge (SOC) of the battery; using at least the measured SOC and an initial co-state value stored in a memory, performing a process to iteratively update the co-state value to obtain an updated co-state value; using at least the updated co-state value, computing an updated control value; and applying the updated control value to control a usage of the battery and the internal combustion engine.

Huang, Mike X.↗
NASA NTRS1985

Mass Transport Phenomena Between Bubbles and Dissolved Gases in Liquids Under Reduced Gravity Conditions

The long term objective of the experiment is to observe the dissolution of isolated, immobile gas bubbles of specified size and composition in a solvent liquid of known concentration in the reduced gravity environment of earth orbit. Preliminary bubble dissolution experiment conducted both in the NASA Lewis 2.2 sec drop tower and in normal gravity using SO2 - Toluene system were not completely successful in their objective. The method of gas injection and lack of bubble interface stabiliy experienced due to the extreme solubility of SO in Toluene has the effects of changing the problem from that of bubble dissolution to one of bubble formation stability and subsequent dissolution in a liquid of unknown initial solute concentration. Current work involves further experimentation in order to refine the bubble injection system and to investigate the concept of having a bubble with a critical radius in a state of unstable equilibrium.

Dewitt, K. J.↗
DOE OSTI2022

Segmenting water and shadow regions within WorldView imagery using local binary patterns

Many shadow detection algorithms pertaining to remotely sensed imagery exist. Several of these algorithms exploit the spectral characteristics of shadows within imagery to identify shaded regions. However, these algorithms can have problems when water is also present within the imagery because water shares similar spectral characteristics with shadows. Many of these algorithms are applied to small image subsets instead of the image as a whole and frequently are applied to urban environments that require additional use of the normalized difference water index or other features to affect the removal of water from the shadow mask. This diversity of the image scene content coupled with the complexity and wide variety of environmental conditions that satellite imagery can acquire makes reliably separating shadow and water within larger images a complex problem. Thresholding various spectral indices to produce segmentation maps can be challenging when large, complex, and often imbalanced scenes are captured, which may require manual adjustment of algorithmic parameters for different image areas. Here, we present a methodology that makes use of the near-infrared channel using a watershed segmentation algorithm, local binary pattern measurements, and a support vector machine to classify shadow and water within a full WorldView scene. Results are promising with initial accuracies above 97% for both shadow and water.

47 OTHER INSTRUMENTATION↗
NASA NTRS1985

The aerodynamic challenges of SRB recovery

Recovery and reuse of the Space Shuttle solid rocket boosters was baselined to support the primary goal to develop a low cost space transportation system. The recovery system required for the 170,000-lb boosters was for the largest and heaviest object yet to be retrieved from exoatmospheric conditions. State-of-the-art design procedures were ground-ruled and development testing minimized to produce both a reliable and cost effective system. The ability to utilize the inherent drag of the boosters during the initial phase of reentry was a key factor in minimizing the parachute loads, size and weight. A wind tunnel test program was devised to enable the accurate prediction of booster aerodynamic characteristics. Concurrently, wind tunnel, rocket sled and air drop tests were performed to develop and verify the performance of the parachute decelerator subsystem. Aerodynamic problems encountered during the overall recovery system development and the respective solutions are emphasized.

Bacchus, D. L.↗
NASA NTRS1985

Foundations for the numerical solution of the Euler equations

The Navier-Stokes equations represent an extremely good model of the physical phenomena encountered in most aeronautical problems. However, the computational resource needed to solve the Navier-Stokes equations are so large that even with today's supercomputers, it is necessary to make use of simpler models. A large number of external aerodynamic problems can be accurately described by a simpler model. This model consists of an outer inviscid flow plus a boundary-layer thickness correction for the vehicle shape. The outer inviscid model may be represented by the potential equation or by the Euler equation. The present paper provides the foundations for the numerical solution of the Euler equations. The governing equations are considered, taking into account conservation laws, the medium, the differential form of the conservation laws, generalized solutions, shock-fitting, and characteristics. Attention is also given to initial and boundary conditions, existence and uniqueness, and rotational phenomena.

Salas, M. D.↗
NASA NTRS1997

Regularization for Atmospheric Temperature Retrieval Problems

Passive remote sensing of the atmosphere is used to determine the atmospheric state. A radiometer measures microwave emissions from earth's atmosphere and surface. The radiance measured by the radiometer is proportional to the brightness temperature. This brightness temperature can be used to estimate atmospheric parameters such as temperature and water vapor content. These quantities are of primary importance for different applications in meteorology, oceanography, and geophysical sciences. Depending on the range in the electromagnetic spectrum being measured by the radiometer and the atmospheric quantities to be estimated, the retrieval or inverse problem of determining atmospheric parameters from brightness temperature might be linear or nonlinear. In most applications, the retrieval problem requires the inversion of a Fredholm integral equation of the first kind making this an ill-posed problem. The numerical solution of the retrieval problem requires the transformation of the continuous problem into a discrete problem. The ill-posedness of the continuous problem translates into ill-conditioning or ill-posedness of the discrete problem. Regularization methods are used to convert the ill-posed problem into a well-posed one. In this paper, we present some results of our work in applying different regularization techniques to atmospheric temperature retrievals using brightness temperatures measured with the SSM/T-1 sensor. Simulation results are presented which show the potential of these techniques to improve temperature retrievals. In particular, no statistical assumptions are needed and the algorithms were capable of correctly estimating the temperature profile corner at the tropopause independent of the initial guess.

Velez-Reyes, Miguel↗
DOE OSTI2026

Structure-aware Initialization via Numerical Continuation and Informed Priors

Scientific machine learning (SciML) often operates in ill-conditioned, weakly identifiable regimes due to limited data or indirect observations. In such settings, optimization and inference are highly sensitive to the starting point, making initialization--often under-reported--a consequential degree of freedom. Random initialization is not a neutral default as it induces an implicit prior over candidate solutions and can systematically bias the result, producing large run-to-run variability. Here, we formalize this view by treating initialization as a hidden confounder in SciML and develop a unifying theory for structure-aware initialization via numerical continuation, constructing warm starts from related problem instances. Across representative tasks, including physics-informed neural networks, maximum likelihood estimation, and variational inference, warm starts have been shown to consistently reduce optimization effort and improve reliability.

Data integrity↗
NASA NTRS1971

Film stability in a vertical rotating tube with a core-gas flow.

The linear hydrodynamic stability of a thin-liquid layer flowing along the inside wall of a vertical tube rotating about its axis in the presence of a core-gas flow is examined. The stability problem is formulated under the conditions that the liquid film is thin, the density and viscosity ratios of gas to liquid are small and the relative (axial) pressure gradient in the gas is of the same order as gravity. The resulting eigenvalue problem is first solved by a perturbation method appropriate to axisymmetric long-wave disturbances. The damped nature (to within the thin-film and other approximations made) of the nonaxisymmetric and short-wave disturbances is noted. In view of the limitations on a truncated perturbation solution when the disturbance wavenumber is not small, an initial value method using digital computer is presented. Stability characteristics of neutral, growing, and damped modes are presented showing the influences of rotation, surface tension, and the core-gas flow. Energy balance in a neutral mode is also illustrated.

Sarma, G. S. R.↗
NASA NTRS2017

Incorporation of Failure Into an Orthotropic Three-Dimensional Model with Tabulated Input Suitable for Use in Composite Impact Problems

The need for accurate material models to simulate the deformation, damage and failure of polymer matrix composites under impact conditions is becoming critical as these materials are gaining increased use in the aerospace and automotive communities. The aerospace community has identified several key capabilities which are currently lacking in the available material models in commercial transient dynamic finite element codes. To attempt to improve the predictive capability of composite impact simulations, a next generation material model is being developed for incorporation within the commercial transient dynamic finite element code LS-DYNA. The material model, which incorporates plasticity, damage and failure, utilizes experimentally based tabulated input to define the evolution of plasticity and damage and the initiation of failure as opposed to specifying discrete input parameters such as modulus and strength. The plasticity portion of the orthotropic, three-dimensional, macroscopic composite constitutive model is based on an extension of the Tsai-Wu composite failure model into a generalized yield function with a non-associative flow rule. For the damage model, a strain equivalent formulation is used to allow for the uncoupling of the deformation and damage analyses. In the damage model, a semi-coupled approach is employed where the overall damage in a particular coordinate direction is assumed to be a multiplicative combination of the damage in that direction resulting from the applied loads in various coordinate directions. For the failure model, a tabulated approach is utilized in which a stress or strain based invariant is defined as a function of the location of the current stress state in stress space to define the initiation of failure. Failure surfaces can be defined with any arbitrary shape, unlike traditional failure models where the mathematical functions used to define the failure surface impose a specific shape on the failure surface. In the current paper, the complete development of the failure model is described and the generation of a tabulated failure surface for a representative composite material is discussed.

Impact↗
NASA NTRS2017

Implementation of a Tabulated Failure Model Into a Generalized Composite Material Model Suitable for Use in Impact Problems

The need for accurate material models to simulate the deformation, damage and failure of polymer matrix composites under impact conditions is becoming critical as these materials are gaining increased use in the aerospace and automotive communities. The aerospace community has identified several key capabilities which are currently lacking in the available material models in commercial transient dynamic finite element codes. To attempt to improve the predictive capability of composite impact simulations, a next generation material model is being developed for incorporation within the commercial transient dynamic finite element code LS-DYNA. The material model, which incorporates plasticity, damage and failure, utilizes experimentally based tabulated input to define the evolution of plasticity and damage and the initiation of failure as opposed to specifying discrete input parameters such as modulus and strength. The plasticity portion of the orthotropic, three-dimensional, macroscopic composite constitutive model is based on an extension of the Tsai-Wu composite failure model into a generalized yield function with a non-associative flow rule. For the damage model, a strain equivalent formulation is used to allow for the uncoupling of the deformation and damage analyses. For the failure model, a tabulated approach is utilized in which a stress or strain based invariant is defined as a function of the location of the current stress state in stress space to define the initiation of failure. Failure surfaces can be defined with any arbitrary shape, unlike traditional failure models where the mathematical functions used to define the failure surface impose a specific shape on the failure surface. In the current paper, the complete development of the failure model is described and the generation of a tabulated failure surface for a representative composite material is discussed.

polymer matrix composites↗
DOE OSTI2026

Nonlinear causality and strong hyperbolicity of baryon-rich Israel-Stewart hydrodynamics

We present the first set of fully-nonlinear, necessary and sufficient conditions guaranteeing causal evolution of the initial data for the Israel-Stewart hydrodynamic equations with shear and bulk viscosity coupled to a nonzero baryon current. These constraints not only provide nonlinear causality: they also (a) guarantee the existence of a locally well-posed evolution of the initial data (they enforce strong hyperbolicity) when excluding the endpoints of the bounds, (b) arise from purely algebraic constraints that make no underlying symmetry assumptions on the degrees of freedom and (c) propagate the relevant symmetries of the degrees of freedom over the entire evolution of the problem. Our work enforces a mathematically rigorous foundation for future studies of viscous relativistic hydrodynamics with baryon-rich matter including neutron star mergers and heavy-ion collisions.

Cordeiro, Ian [University of Illinois Urbana-Champ↗