Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “formulation systems”

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,765 records · Page 98
All sourcesNASA NTRSDOE OSTI
NASA NTRS1973

Space shuttle post-entry and landing analysis. Volume 2: Appendices

Four candidate navigation systems for the space shuttle orbiter approach and landing phase are evaluated in detail. These include three conventional navaid systems and a single-station one-way Doppler system. In each case, a Kalman filter is assumed to be mechanized in the onboard computer, blending the navaid data with IMU and altimeter data. Filter state dimensions ranging from 6 to 24 are involved in the candidate systems. Comprehensive truth models with state dimensions ranging from 63 to 82 are formulated and used to generate detailed error budgets and sensitivity curves illustrating the effect of variations in the size of individual error sources on touchdown accuracy. The projected overall performance of each system is shown in the form of time histories of position and velocity error components.

Crawford, B. S.↗
DOE OSTI2024

Scalable Multiphysics Block Preconditioning for Low Mach Number Compressible Resistive MHD with Application to Magnetic Confinement Fusion

This study investigates multiphysics block preconditioners that are critical in devising scalable Newton–Krylov iterative solvers for longer time-scale fully implicit fluid plasma models. The specific model of interest is the visco-resistive, low Mach number, compressible magnetohydrodynamics (MHD) model. This model describes the dynamics of conducting fluids in the presence of electromagnetic fields and can be used to study aspects of astrophysical phenomena, important science and technology applications, and basic plasma physics. The specific application of interest that motivates this study is the macroscopic simulation of longer time-scale stability and disruptions of magnetic confinement fusion devices, specifically the ITER Tokamak. The computational solution of the governing balance equations for mass, momentum, heat transfer, and magnetic induction for resistive MHD systems can be extremely challenging. These difficulties arise from both the strong nonlinear, nonsymmetric coupling of fluid and electromagnetic phenomena as well as the significant range of time and length scales that the interactions of these physical mechanisms produce. To handle the range of time and spatial scales of interest, a fully implicit unstructured variational multiscale finite element formulation is employed. For the scalable solution of the Newton linearized systems, fully coupled block preconditioners are designed to leverage algebraic multigrid subsolves. In conclusion, results are presented for the strong and weak scaling of the method as well as the robustness of these techniques for a large range of Lundquist numbers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Physics and the choice of regulators in functional renormalisation group flows

The Renormalisation Group is a versatile tool for the study of many systems where scale-dependent behaviour is important. Its functional formulation can be cast into the form of an exact flow equation for the scale-dependent effective action in the presence of an infrared regularisation. The functional RG flow for the scale-dependent effective action depends explicitly on the choice of regulator, while the physics does not. In this work, we systematically investigate three key aspects of how the regulator choice affects RG flows: (i) We study flow trajectories along closed loops in the space of action functionals varying both, the regulator scale and shape function. Such a flow does not vanish in the presence of truncations. Based on a definition of the length of an RG trajectory, we suggest a constructive procedure for devising optimised regularisation schemes within a truncation. (ii) In systems with various field variables, a choice of relative cutoff scales is required. At the example of relativistic bosonic two-field models, we study the impact of this choice as well as its truncation dependence. We show that a crossover between different universality classes can be induced and conclude that the relative cutoff scale has to be chosen carefully for a reliable description of a physical system. (iii) Non-relativistic continuum models of coupled fermionic and bosonic fields exhibit also dependencies on relative cutoff scales and regulator shapes. At the example of the Fermi polaron problem in three spatial dimensions, we illustrate such dependencies and show how they can be interpreted in physical terms. - Highlights: • We show the regulator dependence of functional renormalisation (FRG) group equations. • Solving FRG equations along closed loops provides a distance between two regulators. • We suggest a practical procedure for devising optimal regulators. • A universality class crossover can be induced by changing relative cutoff scales. • We examine the regulator dependence explicitly in the Fermi polaron problem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1989

A finite element approach for large motion dynamic analysis of multibody structures in space

A three-dimensional finite element formulation for modeling the transient dynamics of constrained multibody space sructures with truss-like configurations is presented. Convected coordinate systems are used to define rigid-body motion of individual elements in the system. These systems are located at one end of each element and are oriented such that one axis passes through the other end of the element. Deformation of each element, relative to its convected coordinate system, is defined by cubic flexural shape functions as used in finite element methods of structural analysis. The formulation is oriented toward joint dominated structures and places the generalized coordinates at the joint. A transformation matrix is derived to integrate joint degree-of-freedom into the equations of motion of the element. Based on the derivation, a general-purpose code LATDYN (Large Angle Transient DYNamics) was developed. Two examples are presented to illustrate the application of the code. For the spin-up of a flexible beam, results are compared with existing solutions available in the literature. For the deployment of one bay of a deployable space truss (the Minimast), results are verified by the geometric knowledge of the system and converged solution of a successively refined model.

Chang, Che-Wei↗
DOE OSTI2014

Estimating Residual Solids Volume In Underground Storage Tanks

The Savannah River Site liquid waste system consists of multiple facilities to safely receive and store legacy radioactive waste, treat, and permanently dispose waste. The large underground storage tanks and associated equipment, known as the 'tank farms', include a complex interconnected transfer system which includes underground transfer pipelines and ancillary equipment to direct the flow of waste. The waste in the tanks is present in three forms: supernatant, sludge, and salt. The supernatant is a multi-component aqueous mixture, while sludge is a gel-like substance which consists of insoluble solids and entrapped supernatant. The waste from these tanks is retrieved and treated as sludge or salt. The high level (radioactive) fraction of the waste is vitrified into a glass waste form, while the low-level waste is immobilized in a cementitious grout waste form called saltstone. Once the waste is retrieved and processed, the tanks are closed via removing the bulk of the waste, chemical cleaning, heel removal, stabilizing remaining residuals with tailored grout formulations and severing/sealing external penetrations. The comprehensive liquid waste disposition system, currently managed by Savannah River Remediation, consists of 1) safe storage and retrieval of the waste as it is prepared for permanent disposition; (2) definition of the waste processing techniques utilized to separate the high-level waste fraction/low-level waste fraction; (3) disposition of LLW in saltstone; (4) disposition of the HLW in glass; and (5) closure state of the facilities, including tanks. This paper focuses on determining the effectiveness of waste removal campaigns through monitoring the volume of residual solids in the waste tanks. Volume estimates of the residual solids are performed by creating a map of the residual solids on the waste tank bottom using video and still digital images. The map is then used to calculate the volume of solids remaining in the waste tank. The ability to accurately determine a volume is a function of the quantity and quality of the waste tank images. Currently, mapping is performed remotely with closed circuit video cameras and still photograph cameras due to the hazardous environment. There are two methods that can be used to create a solids volume map. These methods are: liquid transfer mapping / post transfer mapping and final residual solids mapping. The task is performed during a transfer because the liquid level (which is a known value determined by a level measurement device) is used as a landmark to indicate solids accumulation heights. The post transfer method is primarily utilized after the majority of waste has been removed. This method relies on video and still digital images of the waste tank after the liquid transfer is complete to obtain the relative height of solids across a waste tank in relation to known and usable landmarks within the waste tank (cooling coils, column base plates, etc.). In order to accurately monitor solids over time across various cleaning campaigns, and provide a technical basis to support final waste tank closure, a consistent methodology for volume determination has been developed and implemented at SRS.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Electromagnetic full-$f$ gyrokinetics in the tokamak edge with discontinuous Galerkin methods

Here, we present an energy-conserving discontinuous Galerkin scheme for the full-f electromagnetic gyrokinetic system in the long-wavelength limit. We use the symplectic formulation and solve directly for ∂A II /∂t, the inductive component of the parallel electric field, using a generalized Ohm’s law derived directly from the gyrokinetic equation. Linear benchmarks are performed to verify the implementation and show that the scheme avoids the Amp`ere cancellation problem. We perform a nonlinear electromagnetic simulation in a helical open-field-line system as a rough model of the tokamak scrape-off layer using parameters from the National Spherical Torus Experiment (NSTX). This is the first published nonlinear electromagnetic gyrokinetic simulation on open field lines. Comparisons are made to a corresponding electrostatic simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Accelerating and stabilizing the convergence of vibrational self-consistent field calculations via the direct inversion of the iterative subspace (vDIIS) algorithm

The vibrational self-consistent field (VSCF) method yields anharmonic states and spectra for molecular vibrations, and it serves as the starting point for more sophisticated correlated-vibration methods. Convergence of the iterative, non-linear optimization in VSCF calculations can be erratic or altogether unsuccessful, particularly for chemical systems involving low-frequency motions. In this work, a vibrational formulation of the Direct Inversion of the Iterative Subspace method of Pulay is presented and investigated. This formulation accounts for distinct attributes of the vibrational and electronic cases, including the expansion of each single-mode vibrational wavefunction in its own basis set. The resulting Direct Inversion of the Iterative Subspace method is shown to substantially accelerate VSCF convergence in all convergent cases as well as rectify many cases where Roothaan-based methods fail. Performance across systems ranging from small, rigid molecules to weakly bound molecular clusters is investigated in this analysis.

Chemistry↗
DOE OSTI2018

Coupled elasticity, plastic slip, and twinning in single crystal titanium loaded by split-Hopkinson pressure bar

Here coupled elasticity, plastic slip, and twinning in high-purity single crystal titanium loaded by split-Hopkinson pressure bar (SHPB) are investigated, in the framework of large deformation and by utilizing the finite element method (FEM). A thermodynamically consistent system of equations for combined plastic slip and twinning is formulated. Novel kinematics is proposed to develop the driving forces for slip and twinning processes and to consider plastic slip in twins. Dislocation based crystal plasticity is proposed, with an emphasis on the interactions among all slip modes in multi-variant/multiphase heterogeneous materials. A mechanism for dislocation density evolution during twinning is proposed. The impact loading along the [0001] and [1$0\bar{1}$1] directions of single crystal high purity titanium is investigated. The evolution of stress-strain field, dislocation density, and volume fraction for each variant in the sample during the entire loading process are obtained and discussed in detail. Results in simulations show that for the [0001] specimen due to the existence of six symmetry primary twin variants with respect to the loading direction and fifteen different slip systems in each variant, the fields of stress and dislocation density are very homogenous in the sample. The simulation results are compared with experimental data for the [0001] and [1$0\bar{1}$1] specimens. The experimental results are explained and interpreted for both specimens, which include the pole figures, the stress-strain curves, and volume fraction of variants. The simulation results provide an important insight into mechanical responses of high-purity single crystal titanium under high rate loading, and this paper advances the mesoscale modeling and simulations in coupled plastic slip and twinning.

36 MATERIALS SCIENCE↗
NASA NTRS1988

Automated seed manipulation and planting

Activities for the Fall Semester, 1987 focused on investigating the mechanical/electrical properties of wheat seeds and forming various Seed Planting System (SPS) concepts based on those properties. The Electrical Division of the design group was formed to devise an SPS using electrostatic charge fields for seeding operations. Experiments concerning seed separation using electrical induction (rearranging of the charges within the seed) were conducted with promising results. The seeds, when exposed to the high voltage and low current field produced by a Van de Graff generator, were observed to move back and forth between two electrodes. An SPS concept has been developed based on this phenomena, and will be developed throughout the Spring Semester, 1988. The Mechanical Division centered on SPS concepts involving valves, pumps, and fluids to separate and deliver seeds. An SPS idea utilizing the pressure difference caused by air as it rushes out of holes drilled in the wall of a closed container has been formulated and will be considered for future development. Also, a system of seed separation and delivery employing a combination of centrifugal force, friction, and air flow was considered.

Garcia, Ray↗
NASA NTRS2014

Development Of FIAT-Based Parametric Thermal Protection System Mass Estimating Relationships For NASA's Multi-Mission Earth Entry Concept

An integrated tool called the Multi Mission System Analysis for Planetary Entry Descent and Landing (M-SAPE) is being developed as part of NASAs In-Space Propulsion Technology (ISPT) program. Part of M-SAPEs development requires the formulation of mass estimating relationships (MERs) to determine the vehicle's Thermal Protection System (TPS) material and required thickness for safe Earth entry. The objective of this study was to develop MERs using simple correlations that were non-ITAR and matched as accurately as possible NASAs high-fidelity ablation modeling tool, the Fully Implicit Ablation and Thermal Analysis Program (FIAT ). These MERs would be a first-estimate for feasibility studies; it is understood that higher-fidelity modeling like FIAT would be necessary once a proposed trajectory was down-selected. The trajectory space for these MERS consisted of 840 different trajectories, and a materials heating limit was the main constraint for an allowable trajectory. MERs for the vehicle fore body included the ablating materials Phenolic Impregnated Carbon Ablator (PICA ) and Carbon Phenolic atop Advanced Carbon-Carbon. For the backshell the materials were Silicone Impregnated Reusable Ceramic Ablator (SIRCA ), Acusil II, SLA-561V, and LI-900. The MERFIAT ratio indicates MERs are accurate to within 14 percent (at one standard deviation) of FIAT prediction, and the most any MER can under-predict TPS thickness is 18.7 percent of FIAT prediction. This poster focuses on the development of these MERs, the resulting equations, model limitations, and model accuracy.

thermal protection system↗
DOE OSTI2018

Metal Hydride Material Development for High Efficiency and Low Cost Hydrogen Compressors

To make the hydrogen economy realistic at a large scale, one of the main technical issues to be overcome is relative to the high-pressure hydrogen delivery. Traditional hydrogen compressors suffer from high operating and maintenance costs and high capital investment. Metal hydride based compressors represent a realistic alternative that can potentially overcome all of the technical hurdles of the traditional mechanical systems. However, to date, the metal hydrides, available for this type of application (e.g. rare earth based materials), are generally based on expensive compounds and elements. A novel materials-genome based high throughput technique was established to define and develop a new class of materials for hydrogen compression. Specific material properties, required to meet the techno-economic targets for hydrogen compression systems, were identified. After validation of the high throughput methodology by comparing the results with the available material data, the technique was applied to identify preliminary new material formulations that meet the targets. A novel techno-economic analysis constrained machine learning system was developed and applied to identify the new set of metal hydride materials.

08 HYDROGEN↗
NASA NTRS1974

A Computing Method for Sound Propagation Through a Nonuniform Jet Stream

Understanding the principles of jet noise propagation is an essential ingredient of systematic noise reduction research. High speed computer methods offer a unique potential for dealing with complex real life physical systems whereas analytical solutions are restricted to sophisticated idealized models. The classical formulation of sound propagation through a jet flow was found to be inadequate for computer solutions and a more suitable approach was needed. Previous investigations selected the phase and amplitude of the acoustic pressure as dependent variables requiring the solution of a system of nonlinear algebraic equations. The nonlinearities complicated both the analysis and the computation. A reformulation of the convective wave equation in terms of a new set of dependent variables is developed with a special emphasis on its suitability for numerical solutions on fast computers. The technique is very attractive because the resulting equations are linear in nonwaving variables. The computer solution to such a linear system of algebraic equations may be obtained by well-defined and direct means which are conservative of computer time and storage space. Typical examples are illustrated and computational results are compared with available numerical and experimental data.

Padula, S. L.↗
NASA NTRS2012

Engineering a Multimission Approach to Navigation Ground Data System Operations

The Mission Design and Navigation (MDNAV) Section at the Jet Propulsion Laboratory (JPL) supports many deep space and earth orbiting missions from formulation to end of mission operations. The requirements of these missions are met with a multimission approach to MDNAV ground data system (GDS) infrastructure capable of being shared and allocated in a seamless and consistent manner across missions. The MDNAV computing infrastructure consists of compute clusters, network attached storage, mission support area facilities, and desktop hardware. The multimission architecture allows these assets, and even personnel, to be leveraged effectively across the project lifecycle and across multiple missions simultaneously. It provides a more robust and capable infrastructure to each mission than might be possible if each constructed its own. It also enables a consistent interface and environment within which teams can conduct all mission analysis and navigation functions including: trajectory design; ephemeris generation; orbit determination; maneuver design; and entry, descent, and landing analysis. The savings of these efficiencies more than offset the costs of increased complexity and other challenges that had to be addressed: configuration management, scheduling conflicts, and competition for resources. This paper examines the benefits of the multimission MDNAV ground data system infrastructure, focusing on the hardware and software architecture. The result is an efficient, robust, scalable MDNAV ground data system capable of supporting more than a dozen active missions at once.

Mission Design and Navigation (MDNAV)↗
DOE OSTI2024

A Particle-in-Cell Method for Plasmas with a Generalized Momentum Formulation, Part II: Enforcing the Lorenz Gauge Condition

In a previous paper Christlieb et al. (A particle-in-cell method for plasmas with a generalized momentum formulation, part I: Model formulation, 2024), we developed a new particle-in-cell (PIC) method for the relativistic Vlasov–Maxwell system in which the electromagnetic fields and the equations of motion for the particles were cast in terms of scalar and vector potentials through a Hamiltonian formulation. This new method evolved the potentials under the Lorenz gauge using integral equation methods. New methods to construct spatial derivatives of the potentials that converge at the same rates as the fields were also presented. The new particle method was compared against standard explicit discretizations, including the well-known FDTD-PIC method, for a range of applications involving sheaths and particle beams. Here, this paper extends this new class of methods by focusing on the enforcement the Lorenz gauge condition in both exact and approximate forms using co-located meshes. A time-consistency property of the proposed field solver for the vector potential form of Maxwell’s equations is established, which is shown to preserve the equivalence between the semi-discrete Lorenz gauge condition and the analogous semi-discrete continuity equation. Using this property, we present three methods to enforce a semi-discrete gauge condition. The first method introduces an update for the continuity equation that is consistent with the discretization of the Lorenz gauge condition. Both the finite difference and spectral implementations satisfy this discrete gauge condition to machine precision. The second approach we propose enforces a semi-discrete continuity equation using the boundary integral solution to the field equations. The potential benefit of this approach is that it eliminates spatial derivatives that appear on the particle data, namely the current density, which is often calculated by linear combinations of low-order spline basis functions. This method is ideally suited to boundary integral equation methods that invert multi-dimensional operators without dimensional splitting techniques and will be the subject of future work. The third approach introduces a gauge correcting method that makes direct use of the gauge condition to modify the scalar potential and uses local maps for both the charge and current densities. This results in a gauge error, as the maps do not enforce the continuity equation. The vector potential coming from the current density is taken to be exact, and using the Lorenz gauge, we compute a correction to the scalar potential that makes the two potentials satisfy the gauge condition. This method also enforces the gauge condition to machine precision. We demonstrate two of the proposed methods in the context of periodic domains. Problems defined on bounded domains, including those with complex geometric features remain an ongoing effort. However, this work shows that it is possible to design computationally efficient methods that can effectively enforce the Lorenz gauge condition in a non-staggered PIC formulation.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1993

Forecasting the impact of virtual environment technology on maintenance training

To assist NASA and the Air Force in determining how and when to invest in virtual environment (VE) technology for maintenance training, we identified possible roles for VE technology in such training, assessed its cost-effectiveness relative to existing technologies, and formulated recommendations for a research agenda that would address instructional and system development issues involved in fielding a VE training system. In the first phase of the study, we surveyed VE developers to forecast capabilities, maturity, and estimated costs for VE component technologies. We then identified maintenance tasks and their training costs through interviews with maintenance technicians, instructors, and training developers. Ten candidate tasks were selected from two classes of maintenance tasks (seven aircraft maintenance and three space maintenance) using five criteria developed to identify types of tasks most likely to benefit from VE training. Three tasks were used as specific cases for cost-benefit analysis. In formulating research recommendations, we considered three aspects of feasibility: technological considerations, cost-effectiveness, and anticipated R&D efforts. In this paper, we describe the major findings in each of these areas and suggest research efforts that we believe will help achieve the goal of a cost-effective VE maintenance training system by the next decade.

Schlager, Mark S.↗
DOE OSTI2026

Efficient Computation of Doppler-Broadened Elastic Scattering Kernel Moments Using Ladder-Operator Formulation

Anefficient routine for computing Legendre moments of the Doppler-broadened elastic scattering kernel, including resonance scattering effects, has been implemented in the ISOXML module of Griffin. Isotropic scattering in the center-of-mass system and the ideal gas model for target motion are assumed. A ladder-operator formulation is introduced to compute all Legendre moments from order 0 to N simultaneously, enabling near-linear scaling of computational cost with respect to the maximum Legendre order. A physics-based strategy for constructing outgoing energy grids has also been developed, in which a tailored base grid is combined with adaptive refinement to maintain accuracy while limiting the number of outgoing energy points. For energies between resonances, a constant cross-section model is employed to further reduce computational cost. In addition, a quantitative criterion is derived to determine isotope-wise cut-off incident energies based on a prescribed up-scattering probability coverage. For 238U, up to incident energies of approximately 75, 230, and 661 eV at 294, 900, and 2500 K (corresponding to a 2% up-scattering probability threshold), computation of P0 kernels requires 1–8 s and computation of P0–P5 kernels requires 0.4–4 min using a single thread, while maintaining 1–3% relative error in up-scattering probability. These results demonstrate that the proposed formulation enables accurate and computationally practical Doppler-broadened kernel generation for online multigroup cross-section production in Griffin.

Doppler-broadening↗
DOE OSTI2016

Coherent states, quantum gravity, and the Born-Oppenheimer approximation. I. General considerations

This article, as the first of three, aims at establishing the (time-dependent) Born-Oppenheimer approximation, in the sense of space adiabatic perturbation theory, for quantum systems constructed by techniques of the loop quantum gravity framework, especially the canonical formulation of the latter. The analysis presented here fits into a rather general framework and offers a solution to the problem of applying the usual Born-Oppenheimer ansatz for molecular (or structurally analogous) systems to more general quantum systems (e.g., spin-orbit models) by means of space adiabatic perturbation theory. The proposed solution is applied to a simple, finite dimensional model of interacting spin systems, which serves as a non-trivial, minimal model of the aforesaid problem. Furthermore, it is explained how the content of this article and its companion affect the possible extraction of quantum field theory on curved spacetime from loop quantum gravity (including matter fields).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1994

General analytic results for nonlinear waves and solitons in molecular clouds

We study nonlinear wave phenomena in self-gravitating fluid systems, with a particular emphasis on applications to molecular clouds. This paper presents analytical results for one spatial dimension. We show that a large class of physical systems can be described by theories with a 'charge density' q(rho); this quantity replaces the density on the right-hand side of the Poisson equation for the gravitational potential. We use this formulation to prove general results about nonlinear wave motions in self-gravitating systems. We show that in order for stationary waves to exist, the total charge (the integral of the charge density over the wave profile) must vanish. This 'no-charge' property for solitary waves is related to the capability of a system to be stable to gravitational perturbations for arbitrarily long wavelengths. We find necessary and sufficient conditions on the charge density for the existence of solitary waves and stationary waves. We study nonlinear wave motions for Jeans-type theories (where q(rho) = rho-rho(sub 0)) and find that nonlinear waves of large amplitude are confined to a rather narrow range of wavelengths. We also study wave motions for molecular clouds threaded by magnetic fields and show how the allowed range of wavelengths is affected by the field strength. Since the gravitational force in one spatial dimension does not fall off with distance, we consider two classes of models with more realistic gravity: Yukawa potentials and a pseudo two-dimensional treatment. We study the allowed types of wave behavior for these models. Finally, we discuss the implications of this work for molecular cloud structure. We argue that molecular clouds can support a wide variety of wave motions and suggest that stationary waves (such as those considered in this paper) may have already been observed.

Adams, Fred C.↗