Search NASA⌕ Search

SEARCH · Search NASA

Results for “fully linear model”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS1978

Stochastic partial differential equations in turbulence related problems

The theory of stochastic partial differential equations (PDEs) and problems relating to turbulence are discussed by employing the theories of Brownian motion and diffusion in infinite dimensions, functional differential equations, and functional integration. Relevant results in probablistic analysis, especially Gaussian measures in function spaces and the theory of stochastic PDEs of Ito type, are taken into account. Linear stochastic PDEs are analyzed through linearized Navier-Stokes equations with a random forcing. Stochastic equations for waves in random media as well as model equations in turbulent transport theory are considered. Markovian models in fully developed turbulence are discussed from a stochastic equation viewpoint.

Chow, P.-L.↗
NASA NTRS1991

Computed stress fields in GaAs during vertical Bridgman growth

The thermoelastic stress field of GaAs crystals grown by the vertical Bridgman method has been approximated using a linear elastic, axisymmetric stress model. The model included the effects due to the thermal field and the gravitational field, as well as the interaction with the ampoule, but did not fully account for the elastic anisotropy of the crystal. Realistic thermal fields for different furnace temperatures and insulation zone thicknesses were computed with a similar model. The difference between the Von Mises stress and the critical resolved shear stress (CRSS) was used as a measure of the number of dislocations present. Different sets of growth parameters were examined to determine which produced the lowest amount of stress in the crystal. Results show that the crystal-ampoule interaction is the most important parameter in dislocation generation. In addition, for the same sticking boundary conditions, concave interfaces have lower stress than convex interfaces.

Rosch, William↗
NASA NTRS2005

Development of a Stirling System Dynamic Model With Enhanced Thermodynamics

The Stirling Convertor System Dynamic Model developed at NASA Glenn Research Center is a software model developed from first principles that includes the mechanical and mounting dynamics, the thermodynamics, the linear alternator, and the controller of a free-piston Stirling power convertor, along with the end user load. As such it represents the first detailed modeling tool for fully integrated Stirling convertor-based power systems. The thermodynamics of the model were originally a form of the isothermal Stirling cycle. In some situations it may be desirable to improve the accuracy of the Stirling cycle portion of the model. An option under consideration is to enhance the SDM thermodynamics by coupling the model with Gedeon Associates Sage simulation code. The result will be a model that gives a more accurate prediction of the performance and dynamics of the free-piston Stirling convertor. A method of integrating the Sage simulation code with the System Dynamic Model is described. Results of SDM and Sage simulation are compared to test data. Model parameter estimation and model validation are discussed.

Regan, Timothy F.↗
DOE OSTI2021

Cross-Code verification and sensitivity analysis to effectively model the electrothermal instability

This manuscript presents verification cases that are developed to study the electrothermal instability (ETI). Specific verification cases are included to ensure that the unit physics components necessary to model the ETI are accurate, providing a path for fluid-based codes to effectively simulate ETI in the linear and nonlinear growth regimes. Two software frameworks with different algorithmic approaches are compared for accuracy in their ability to simulate diffusion of a magnetic field, linear growth of the ETI, and a fully nonlinear ETI evolution. The nonlinear ETI simulations show early time agreement, with some differences emerging, as noted in the wavenumber spectrum, late into the nonlinear development of ETI. Here, a sensitivity study explores the role of equation-of-state (EOS), vacuum density, and vacuum resistivity. EOS and vacuum resistivity are found to be the most critical factors in the modeling of nonlinear ETI development.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2025

Robust Trajectory Optimization for NRHO Rendezvous Using SPICE Kernel Relative Motion

In this paper, robust optimization is performed on trajectory correction maneuvers during the lunar lander return phase of an Artemis mission, treating the trajectory from one hour after low lunar orbit departure to arrival in the vicinity of the lunar Gateway as a relative motion problem. To enable rapid stochastic optimization techniques requiring many candidate trajectories, SPICE kernel relative motion as implemented by the Quadratic Interpolated State Transition (QIST) system is used as the underlying dynamics propagation. The optimization is performed with a genetic optimizer using linear covariance (LinCov) software in a simplified operational context, taking into account the availability of navigation sensors with varying measurement models, ranges, and accuracies. No numerical integration is used, since the relative motion around Gateway is fully characterized with the a priori computation of the QIST coefficients. Maneuver placements are computed to optimize the minimum 3σ delta-v of the trajectory, the position dispersion at a target point, and a convex combination of these two metrics. An order of magnitude runtime improvement is provided over legacy methods with less than 10% error introduced. All QIST results are shown to be in-family with legacy methods. The tradespace for optimal delta-v design is found to range from 77.0 to 93.9 m/s, while the range of optimal dispersion is between 1.4 and 11.7 km.

Relative Motion↗
DOE OSTI2014

Accelerating dissipative particle dynamics simulations on GPUs: Algorithms, numerics and applications

We present a scalable dissipative particle dynamics simulation code, fully implemented on the Graphics Processing Units (GPUs) using a hybrid CUDA/MPI programming model, which achieves 10-30 times speedup on a single GPU over 16 CPU cores and almost linear weak scaling across a thousand nodes. A unified framework is developed within which the efficient generation of the neighbor list and maintaining particle data locality are addressed. Our algorithm generates strictly ordered neighbor lists in parallel, while the construction is deterministic and makes no use of atomic operations or sorting. Such neighbor list leads to optimal data loading efficiency when combined with a two-level particle reordering scheme. A faster in situ generation scheme for Gaussian random numbers is proposed using precomputed binary signatures. We designed custom transcendental functions that are fast and accurate for evaluating the pairwise interaction. The correctness and accuracy of the code is verified through a set of test cases simulating Poiseuille flow and spontaneous vesicle formation. Computer benchmarks demonstrate the speedup of our implementation over the CPU implementation as well as strong and weak scalability. A large-scale simulation of spontaneous vesicle formation consisting of 128 million particles was conducted to further illustrate the practicality of our code in real-world applications.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2002

Covariance Propagation and Partial Eigendecomposition Filtering on the Continuum: The Cases of Advective Dynamics

As a motivation for this lecture, we begin by stating a paradox that challenges our fundamental understanding of covariance evolution (at least it challenged my own). Attempting to resolve this 'divergence paradox' leads us to introduce the continuum fundamental solution operator for the dynamics under consideration, which will be advection dynamics in this lecture. This operator is the object that is approximated by the discrete 'tangent linear model. We then show how the fundamental solution operator can be used to describe the solution of the continuum covariance evolution equation. This description is complete enough to resolve fully the divergence paradox.

Cohn, S.↗
DOE OSTI2024

Phase-space entropy cascade and irreversibility of stochastic heating in nearly collisionless plasma turbulence

We consider a nearly collisionless plasma consisting of a species of “test particles” in one spatial and one velocity dimension, stirred by an externally imposed stochastic electric field—a kinetic analog of the Kraichnan model of passive advection. The mean effect on the particle distribution function is turbulent diffusion in velocity space—known as stochastic heating. Accompanying this heating is the generation of fine-scale structure in the distribution function, which we characterize with the collisionless (Casimir) invariant C 2 ∝ ∫ ∫ d x d v 〈 f 2 〉 —a quantity that here plays the role of (negative) entropy of the distribution function. We find that C 2 is transferred from large scales to small scales in both position and velocity space via a phase-space cascade enabled by both particle streaming and nonlinear interactions between particles and the stochastic electric field. We compute the steady-state fluxes and spectrum of C 2 in Fourier space, with k and s denoting spatial and velocity wave numbers, respectively. In our model, the nonlinearity in the evolution equation for the spectrum turns into a fractional Laplacian operator in k space, leading to anomalous diffusion. Whereas even the linear phase mixing alone would lead to a constant flux of C 2 to high s (towards the collisional dissipation range) at every k , the nonlinearity accelerates this cascade by intertwining velocity and position space so that the flux of C 2 is to both high k and high s simultaneously. Integrating over velocity (spatial) wave numbers, the k -space ( s -space) flux of C 2 is constant down to a dissipation length (velocity) scale that tends to zero as the collision frequency does, even though the rate of collisional dissipation remains finite. The resulting spectrum in the inertial range is a self-similar function in the ( k , s ) plane, with power-law asymptotics at large k and s . Our model is fully analytically solvable, but the asymptotic scalings of the spectrum can also be found via a simple phenomenological theory whose key assumption is that the cascade is governed by a “critical balance” in phase space between the linear and nonlinear timescales. We argue that stochastic heating is made irreversible by this entropy cascade and that, while collisional dissipation accessed via phase mixing occurs only at small spatial scales rather than at every scale as it would in a linear system, the cascade makes phase mixing even more effective overall in the nonlinear regime than in the linear one. Published by the American Physical Society 2024

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

FSISPH: An SPH formulation for impacts between dissimilar materials

Here, we present an SPH formulation with several new features designed to better model the fully-compressible interaction of dissimilar materials. We developed the new method to simulate the atmospheric entry and break-up of small celestial bodies in planetary atmospheres. The formulation uses a unity-based, density-energy discretization of the hydrodynamic conservation laws with linear-corrected kernel gradients. To account for variations in compressibility, we use an HLLC approximate Riemann solver to adjust the velocity gradient at material interfaces. To handle large transverse velocity discontinuities, we introduce a simple slip interface model that limits the artificial viscosity at material interfaces. Diffusion is optionally applied through the velocity gradient and this allows the density and specific thermal energy to evolve in a manner more consistent with the first law of thermodynamics in comparison to other more direct diffusion schemes. We also introduce a material-local second-order artificial conduction scheme used to smooth the specific thermal energy field. Material damage fits neatly under this framework by treating the damage front as a material interface. The method has been implemented as a solver, FSISPH, within the code, Spheral++, and is publicly available on github. We test our new solver on a number of classic shock, mixing, and multi-material problem. The components we outline can significantly improve accuracy of SPH for problems with sharp contact discontinuities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

A multi-dimensional Vlasov-Fokker-Planck code for arbitrarily anisotropic high-energy-density plasmas

To study the kinetic physics of High-Energy-Density Laboratory Plasmas, we have developed the parallel relativistic 2D3P Vlasov-Fokker-Planck code Oshun. The numerical scheme uses a Cartesian mesh in configuration-space and incorporates a spherical harmonic expansion of the electron distribution function in momentum-space. The expansion is truncated such that the necessary angular resolution of the distribution function is retained for a given problem. Finite collisionality causes rapid decay of the high-order harmonics, thereby providing a natural truncation mechanism for the expansion. The code has both fully explicit and implicit field-solvers and employs a linearized Fokker-Planck collision operator. Oshun has been benchmarked against well-known problems, in the highly kinetic limit to model collisionless relativistic instabilities, and in the hydrodynamic limit to recover transport coefficients. The performance of the code, its applicability, and its limitations are discussed in the context of simple problems with relevance to inertial fusion energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2003

Role of genetic background in induced instability

Genomic instability is effectively induced by ionizing radiation. Recently, evidence has accumulated supporting a relationship between genetic background and the radiation-induced genomic instability phenotype. This is possibly due to alterations in proteins responsible for maintenance of genomic integrity or altered oxidative metabolism. Studies in human cell lines, human primary cells, and mouse models have been performed predominantly using high linear energy transfer (LET) radiation, or high doses of low LET radiation. The interplay between genetics, radiation response, and genomic instability has not been fully determined at low doses of low LET radiation. However, recent studies using low doses of low LET radiation suggest that the relationship between genetic background and radiation-induced genomic instability may be more complicated than these same relationships at high LET or high doses of low LET radiation. The complexity of this relationship at low doses of low LET radiation suggests that more of the population may be at risk than previously recognized and may have implications for radiation risk assessment.

NASA Discipline Radiation Health↗
DOE OSTI2025

Predicting Flow in Fracture Networks With Quantum Algorithms

Uncertainty quantification plays a crucial role in the modeling of subsurface flow. For instance, uncertainties in the properties of geologic fracture networks significantly impact flow, requiring numerous simulations to accurately estimate quantities of interest. However, each simulation is computationally expensive because it requires solving a large linear system to capture features that involve both small and large fractures. An example is in percolation, where the interaction of many small fractures (which cumulatively can have a large surface area) with the rock matrix must be modeled precisely. Quantum computing is an emerging tool with the potential to address this issue. Quantum algorithms offer a significant speedup in solving linear systems, achieving efficiencies that are challenging to match with classical approaches. These classical approaches include direct solvers, such as LU decomposition, and iterative methods, notably preconditioned conjugate gradient, commonly used in subsurface modeling to solve large sparse systems. However, applying quantum algorithms to geologic fracture flow requires careful attention to algorithmic and problem-specific constraints to fully realize this quantum advantage. In this work we describe a quantum algorithm for generalized Monte Carlo applications with a quadratic speedup over the classical approaches which can be combined with the quantum speedup, currently under investigation, for solving quantum linear systems for subsurface flow. We show that for quantum algorithms the computational cost of estimating a quantity of interest for a statistical ensemble of networks is roughly the same as that of a single realization, essentially implying that one can get uncertainty quantification for free.

58 GEOSCIENCES↗
DOE OSTI2015

Large-scale 3D geoelectromagnetic modeling using parallel adaptive high-order finite element method

Here, we have investigated the use of the adaptive high-order finite-element method (FEM) for geoelectromagnetic modeling. Because high-order FEM is challenging from the numerical and computational points of view, most published finite-element studies in geoelectromagnetics use the lowest order formulation. Solution of the resulting large system of linear equations poses the main practical challenge. We have developed a fully parallel and distributed robust and scalable linear solver based on the optimal block-diagonal and auxiliary space preconditioners. The solver was found to be efficient for high finite element orders, unstructured and nonconforming locally refined meshes, a wide range of frequencies, large conductivity contrasts, and number of degrees of freedom (DoFs). Furthermore, the presented linear solver is in essence algebraic; i.e., it acts on the matrix-vector level and thus requires no information about the discretization, boundary conditions, or physical source used, making it readily efficient for a wide range of electromagnetic modeling problems. To get accurate solutions at reduced computational cost, we have also implemented goal-oriented adaptive mesh refinement. The numerical tests indicated that if highly accurate modeling results were required, the high-order FEM in combination with the goal-oriented local mesh refinement required less computational time and DoFs than the lowest order adaptive FEM.

58 GEOSCIENCES↗
NASA NTRS2004

Autonomous Guidance of Agile Small-scale Rotorcraft

This report describes a guidance system for agile vehicles based on a hybrid closed-loop model of the vehicle dynamics. The hybrid model represents the vehicle dynamics through a combination of linear-time-invariant control modes and pre-programmed, finite-duration maneuvers. This particular hybrid structure can be realized through a control system that combines trim controllers and a maneuvering control logic. The former enable precise trajectory tracking, and the latter enables trajectories at the edge of the vehicle capabilities. The closed-loop model is much simpler than the full vehicle equations of motion, yet it can capture a broad range of dynamic behaviors. It also supports a consistent link between the physical layer and the decision-making layer. The trajectory generation was formulated as an optimization problem using mixed-integer-linear-programming. The optimization is solved in a receding horizon fashion. Several techniques to improve the computational tractability were investigate. Simulation experiments using NASA Ames 'R-50 model show that this approach fully exploits the vehicle's agility.

Mettler, Bernard↗
DOE OSTI2016

Broadly Applicable Dose Model for Radiotherapy Beams from 6-25 MV

The goal of modern external beam radiotherapy is to deliver a highly targeted radiation dose to a specified treatment volume while sparing the healthy surrounding tissue as much as possible. However, some healthy tissue radiation exposure is unavoidable. Healthy tissue in the margin of the treatment field will receive absorbed doses on the order of the prescribed dose, and tissue far out of field will receive stray dose from scattered and leakage radiation. These doses can lead to late effects such as secondary cancers both near the treatment field and far out of field in tissues only receiving stray radiation dose. Many researchers have investigated algorithms to model absorbed dose from external beam photon radiation therapy; however, none of these works have fully addressed the stray dose far from the treatment field. The AAPM Task Group Report No. 36 described three main sources of stray radiation from external beam radiation therapy delivered with electron linear accelerators. Radiation scattered from the treatment head, known as head scatter, is primarily important within about 10 cm from the field edge. Patient scatter is an important source up to around 30 cm from the field edge. Finally, leakage radiation emanates from the treatment enclosure and predominates the stray radiation dose beyond about 30 cm. Monte Carlo simulations have been a useful research tool for modeling stray dose, but these methods have not found use in clinical settings due to their complexity and long computational times. The feasibility of analytical models to predict stray dose from radiation therapy has been supported by several works, but few attempts have been made at developing models accurate for both in-field and out-of-field doses. Jagetic and Newhauser reported on one such model that predicts absorbed dose from therapeutic, scatter, and leakage radiation with promising results. This model was validated at 6 MV for Conformal Radiation Therapy (CRT) beams from an Elekta linac. Intensity Modulated Radiation Therapy (IMRT) was not considered, and the model has not been tested at beam energies greater than 6 MV. Additionally, this model relies on proprietary, machine specific parameters, and so adapting the model for use on treatment units from other manufacturers is not straightforward. In this study we aim to test whether a physics-based analytical modeling approach is broadly applicable and if sufficient accuracy can be achieved without the need for proprietary, machine specific parameters. In order to accomplish this, a new analytical model was developed that does not rely on machine specific parameters, so that it can easily be trained on measurements from whatever treatment machine is available at a given clinic. This model was tested using measured data from a variety of treatment machines and techniques in the 6 MV to 25 MV energy interval. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2014

Deformation Behavior of Solid-Solution-Strengthened Mg-9wt%Al Alloy: In-Situ Neutron Diffraction and Elastic-Viscoplastic Self-Consistent Modeling

In situ neutron diffraction and elastic viscoplastic self-consistent (EVPSC) modeling have been employed to understand the deformation mechanisms of the loading unloading process under uniaxial tension in a solid-solution-strengthened extruded Mg 9 wt.% Al alloy. The initial texture measured by neutron diffraction shows that the {00.2} basal planes in most grains are tilted around 20 30 from the extrusion axis, indicating that basal slip should be easily activated in a majority of grains under tension. Non-linear stress strain responses are observed during unloading and reloading after the material is fully plastically deformed under tension. In situ neutron diffraction measurements have also demonstrated the non-linear behavior of lattice strains during unloading and reloading, revealing that load redistribution continuously occurs between soft and hard grain orientations. The predicted macroscopic stress strain curve and the lattice strain evolution by the EVPSC model are in good agreement with the experimental data. The EVPSC model provides the relative activities of the available slip and twinning modes, as well as the elastic and plastic strains of the various grain families. It is suggested that the non-linear phenomena in the macroscopic stress strain responses and microscopic lattice strains during unloading and reloading are due to plastic deformation by the operation of basal a slip in the soft grain orientations (e.g. {10.1}, {11.2} and {10.2} grain families).

Lee, Sooyeol↗
NASA NTRS2017

Crack Opening Displacement Behavior in Ceramic Matrix Composites

Ceramic Matrix Composites (CMC) modeling and life prediction strongly depend on oxidation, and therefore require a thorough understanding of when matrix cracks occur, the extent of cracking for given conditions (time-temperature-environment-stress), and the interactions of matrix cracks with fibers and interfaces. In this work, the evolution of matrix cracks in a melt-infiltrated Silicon Carbide/Silicon Carbide (SiC/SiC) CMC under uniaxial tension was examined using scanning electron microscopy (SEM) combined with digital image correlation (DIC) and manual crack opening displacement (COD) measurements. Strain relaxation due to matrix cracking, the relationship between COD's and applied stress, and damage evolution at stresses below the proportional limit were assessed. Direct experimental observation of strain relaxation adjacent to regions of matrix cracking is presented and discussed. Additionally, crack openings were found to increase linearly with increasing applied stress, and no crack was found to pass fully through the gage cross-section. This observation is discussed in the context of the assumption of through-cracks for all loading conditions and fiber architectures in oxidation modeling. Finally, the combination of SEM with DIC is demonstrated throughout to be a powerful means for damage identification and quantification in CMC's at stresses well below the proportional limit.

acoustic emission↗
DOE OSTI2018

A Generalized Backward Euler algorithm for the numerical integration of a viscous breakage model

This paper discusses the formulation and the numerical performance of a fully implicit algorithm used to integrate a rate-dependent model defined within a breakage mechanics framework. For this purpose, a Generalized Backward Euler (GBE) algorithm has been implemented according to two different linearization strategies: The former is derived by a direct linearization of the constitutive equations, while the latter introduces rate effects through a consistency parameter. Additionally, the accuracy and efficiency of the GBE algorithm have been investigated by (1) performing material point analyses and (2) solving initial boundary value problems. In both cases, the overall performance of the underlying algorithm is inspected for a range of loading rates, thus simulating comminution from slow to fast dynamic problems. As the viscous response of the breakage model can be recast through a viscous nucleus function, the presented algorithm can be considered as a general framework to integrate constitutive equations relying on the overstress approach typical of Perzyna-like viscoplastic models.

36 MATERIALS SCIENCE↗