Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “consistent discretization”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Space Shuttle Main Engine Propellant Path Leak Detection Using Sequential Image Processing

Initial research in this study using theoretical radiation transport models established that the occurrence of a leak is accompanies by a sudden but sustained change in intensity in a given region of an image. In this phase, temporal processing of video images on a frame-by-frame basis was used to detect leaks within a given field of view. The leak detection algorithm developed in this study consists of a digital highpass filter cascaded with a moving average filter. The absolute value of the resulting discrete sequence is then taken and compared to a threshold value to produce the binary leak/no leak decision at each point in the image. Alternatively, averaging over the full frame of the output image produces a single time-varying mean value estimate that is indicative of the intensity and extent of a leak. Laboratory experiments were conducted in which artificially created leaks on a simulated SSME background were produced and recorded from a visible wavelength video camera. This data was processed frame-by-frame over the time interval of interest using an image processor implementation of the leak detection algorithm. In addition, a 20 second video sequence of an actual SSME failure was analyzed using this technique. The resulting output image sequences and plots of the full frame mean value versus time verify the effectiveness of the system.

Smith, L. Montgomery↗
NASA NTRS1999

Mass Movement and Landform Degradation on the Icy Galilean Satellites: Results of the Galileo Nominal Mission

The Galileo mission has revealed remarkable evidence of mass movement and landform degradation on the icy Galilean satellites of Jupiter. Weakening of surface materials coupled with mass movement reduces the topographic relief of landforms by moving surface materials down-slope. Throughout the Galileo orbiter nominal mission we have studied all known forms of mass movement and landform degradation of the icy galilean satellites, of which Callisto, by far, displays the most degraded surface. Callisto exhibits discrete mass movements that are larger and apparently more common than seen elsewhere. Most degradation on Ganymede appears consistent with sliding or slumping, impact erosion, and regolith evolution. Sliding or slumping is also observed at very small (100 m) scale on Europa. Sputter ablation, while probably playing some role in the evolution of Ganymede's and Callisto's debris layers, appears to be less important than other processes. Sputter ablation might play a significant role on Europa only if that satellite's surface is significantly older than 10(exp 8) years, far older than crater statistics indicate. Impact erosion and regolith formation on Europa are probably minimal, as implied by the low density of small craters there. Impact erosion and regolith formation may be important on the dark terrains of Ganymede, though some surfaces on this satellite may be modified by sublimation-degradation. While impact erosion and regolith formation are expected to operate with the same vigor on Callisto as on Ganymede, most of the areas examined at high resolution on Callisto have an appearance that implies that some additional process is at work, most likely sublimation-driven landform modification and mass wasting. The extent of surface degradation ascribed to sublimation on the outer two Galilean satellites implies that an ice more volatile than H2O is probably involved.

Moore, Jeffrey M.↗
NASA NTRS1993

An Upwind Multigrid Algorithm for Calculating Flows on Unstructured Grids

An algorithm is described that calculates inviscid, laminar, and turbulent flows on triangular meshes with an upwind discretization. A brief description of the base solver and the multigrid implementation is given, followed by results that consist mainly of convergence rates for inviscid and viscous flows over a NACA four-digit airfoil section. The results show that multigrid does accelerate convergence when the same relaxation parameters that yield good single-grid performance are used; however, larger gains in performance can be realized by doing less work in the relaxation scheme.

Bonhaus, Daryl L.↗
NASA NTRS1978

Design and fabrication of a skin stringer discrete tube actively cooled structural panel

The design optimization and practical implementation of actively cooled structural panel concepts was investigated. The desired actively cooled structural panel consisted of the cooled skin and a substructure. The primary load carrying components were fabricated from 2024-T3 aliminum alloy. The 3003-H14 coolant passage tubing was chosen because of its excellent corrosion resistance, workability needed to obtain the desired cross sectional shape, and strength. The Epon 951 adhesive was selected for its excellent structural properties and is the thinnest of available films, 0.064 mm. The Eccobond 58C silver filled epoxy was chosen because of its high thermal conductivity, and the alumina filled Epon 828 was chosen for structural and expansion characteristics.

Anthony, F. M.↗
DOE OSTI2026

Qubit lattice algorithm simulations of the scattering of a bounded two dimensional electromagnetic pulse from an infinite planar dielectric interface

Qubit lattice algorithm (QLA) simulations are performed for a two-dimensional spatially bounded pulse propagating onto a plane interface between two dielectric slabs. QLA is an initial value scheme that consists of a sequence of unitary collision and streaming operators, with appropriate potential operators, that recover Maxwell equations in inhomogeneous dielectric media to the second order in the lattice discreteness. For the case of total internal reflection, there is transient energy transfer into the second medium due to the evanescent fields as the Poynting unit vector of the pulse is rotated from its incident to reflected direction. Because of the finite spatial extent of the pulse, a self-consistent Goos–Hanchen-type displacement along the interface is found without imposing any explicit interface boundary conditions on the fields. For normal incidence, the standard Fresnel coefficients are recovered for appropriately averaged QLA fields. Energy is conserved at all times to seven significant figures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Efficient general method for numerically modeling laser pulse propagation, overlap, and lifetime effects in amplifiers

An efficient numerical time-dependent general method is developed to address incoherent pulse overlap and lifetime effects in laser amplifiers. The alternating propagation-population laser energetics method (APPLE) has been validated against a semi-discrete coupled rate equation numerical method (SDRE) and analytic formalisms in bounding cases. APPLE is based on decoupled rates applied to a time-dependent framework where both space-time-dependent populations and pulse energetics are consistently updated in each time step. A significant advantage of APPLE lies in its conceptual simplicity, ease of implementation, and relatively small computational cost. SDRE tracks the populations through coupled rates and uses the method of lines to discretize the hyperbolic partial differential transport equations allowing for use of ordinary differential equation solvers. With reasonably sized mesh, we report both energetic and power pulse shape relative differences on the order of one percent between the models over a large range of initial conditions.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Prediction and Analysis of Utah FORGE Injection Activities using a Coupled Thermo-hydro-mechanical and Earthquake (THM+E) Modeling Workflow

A coupled thermo-hydro-mechanical (THM) numerical workflow that is capable of modeling seismic slip is critical for the successful development of enhanced geothermal systems (EGS). By integrating key physical processes, this workflow enables accurate simulation of temperature and pressure diffusions, stress changes, and induced seismicity. As a result, it serves as a vital tool for predicting induced seismicity and optimizing reservoir stimulation strategies. The Utah FORGE (Frontier Observatory for Research in Geothermal Energy) project, located near Milford, Utah, is a U.S. Department of Energy initiative aimed at advancing EGS technology. In April 2024, eight new stimulation stages (Stages 3R-10) were conducted in well 16A (injection well) subsequent to the first series of stimulation (Stages 1-3) performed in April, 2022. To monitor the induced seismicity, geophones were deployed in wells 58-32, 56-32, and 78B-32, while fiber optic cables were also installed in wells 16B, 78-32, and 78B-32 to collect microseismic data and detect frac hits Preliminary analyses of microseismic catalogs and fiber optic data suggest that the stimulated fractures in Stages 3R–6 closely align with that generated during Stage 3, indicating that the new stimulations were likely reactivating the previously stimulated fracture. To better understand the underlying process, a comprehensive modeling approach that can accurately capture thermal, hydrological, mechanical, and seismic responses is essential. In this work, we propose and utilize a coupled thermo-hydro-mechanical and earthquake (THM+E) simulation workflow to numerically investigate the stimulation activities on well 16A. The specific objective is to confirm whether the new stimulation stages (Stages 3R–6) reactivated fractures previously stimulated during Stage 3. For this purpose, we perform THM+E simulations individually for Stages 3, 3R, 4, and 5, incorporating the discrete fracture networks (DFNs) created by the plane-fitting technique based on the microseismic catalogs. The simulation workflow consists of two separate models: a THM model and an earthquake model, coupled in a one-way manner. Detailed descriptions of the workflow are provided in Section 3. Simulation results are presented in terms of injection pressure, permeability evolution, and predicted seismic catalogs, which are then compared with field data for further analyses. This report is structured as follows. In Section 2, we present detailed analyses of the field data and propose the hypothesis that the new stimulation stages (Stages 3R–6) were probably reactivating the previously stimulated fractures in Stage 3. In Section 3, we introduce the coupled THM+E workflow and the problem setup to validate our hypothesis, followed by the simulation results for each stage in Section 4. Meanwhile, discussions are included to analyze the model predictions and their comparison with field data. Lastly, we conclude the report and outline future plans in Section 5.

15 GEOTHERMAL ENERGY↗
NASA NTRS1971

The effects of propagation and source distribution on cosmic ray composition and anisotropy.

We consider the propagation and source distribution of cosmic rays. The principal requirement for the various models we consider is that they should be capable of holding particles in dense regions of the galactic disk for periods of time sufficient to produce the observed fragmentation products of cosmic rays. This can be achieved by both simple and compound diffusion provided that suitable mean free paths and boundary conditions are chosen. The bulk of the anisotropy is caused by the discrete nature of the cosmic ray sources. However, models which reproduce the fragmentation products will in general yield anisotropies consistent with available upper limits.

Ramaty, R.↗
NASA NTRS1985

Towards developing robust algorithms for solving partial differential equations on MIMD machines

Methods for efficient computation of numerical algorithms on a wide variety of MIMD machines are proposed. These techniques reorganize the data dependency patterns to improve the processor utilization. The model problem finds the time-accurate solution to a parabolic partial differential equation discretized in space and implicitly marched forward in time. The algorithms are extensions of Jacobi and SOR. The extensions consist of iterating over a window of several timesteps, allowing efficient overlap of computation with communication. The methods increase the degree to which work can be performed while data are communicated between processors. The effect of the window size and of domain partitioning on the system performance is examined both by implementing the algorithm on a simulated multiprocessor system.

Saltz, J. H.↗
NASA NTRS1988

Towards developing robust algorithms for solving partial differential equations on MIMD machines

Methods for efficient computation of numerical algorithms on a wide variety of MIMD machines are proposed. These techniques reorganize the data dependency patterns to improve the processor utilization. The model problem finds the time-accurate solution to a parabolic partial differential equation discretized in space and implicitly marched forward in time. The algorithms are extensions of Jacobi and SOR. The extensions consist of iterating over a window of several timesteps, allowing efficient overlap of computation with communication. The methods increase the degree to which work can be performed while data are communicated between processors. The effect of the window size and of domain partitioning on the system performance is examined both by implementing the algorithm on a simulated multiprocessor system.

Saltz, Joel H.↗
NASA NTRS1994

Auroral activity associated with Kelvin-Helmholtz instability at the inner edge of the low-latitude boundary layer

Auroral activity occurred in the late afternoon sector (approx. 16 MLT) in the northern hemisphere during the passage at Earth of an interplanetary magnetic cloud on January 14, 1988. The auroral activity consisted of a very dynamic display which was preceded and followed by quiet auroral displays. During the quiet displays, discrete rayed arcs aligned along the geomagnetic L shells were observed. In the active stage, rapidly evolving spiral forms centered on magnetic zenith were evident. The activity persisted for many minutes and was characterized by the absence of directed motion. They were strongly suggestive of intense filaments of upward field-aligned currents embedded in the large-scale region 1 current system. Distortions of the flux ropes as they connect from the equatorial magnetosphere to the ionosphere were witnessed. We assess as possible generating mechanisms three nonlocal sources known to be associated with field-aligned currents. Of these, partial compressions of the magnetosphere due to variations of solar wind dynamic pressure seem an unlikely source. The possibility that the auroral forms are due to reconnection is investigated but is excluded because the active aurora were observed on the closed field line region just equatorward of the convection reversal boundary. To support this conclusion further, we apply recent results on the mapping of ionospheric regions to the equatorial plane based on the Tsyganenko 1989 model (Kaufmann et al., 1993). We find that for comparable magnetic activity the aurora map to the equatorial plane at X(sub GSM) = approx. 3 R(sub E) and approx. 2 R(sub E) inward of the magnetopause, that is, the inner edge of the boundary layer close to dusk. Since the auroral forms are manifestly associated with magnetic field shear, a vortical motion at the equatorial end of the flux rope is indicated, making the Kelvin-Helmholtz instability acting at the inner edge of the low-latitude boundary layer the most probable generating source.

Farrugia, C. J.↗
NASA NTRS1995

Physical properties of the OMC-2 and OMC-3 cores from CS and C(18)O observations

We have investigated the properties of the OMC-2 and OMC-3 cores in the Orion giant molecular cloud using high spatial spectral resolution observations of several transitions of the (13)CO, C(18)O, C(S-32) and C(S-34) molecules taken with the SEST telescope. The OMC-2 core consists of one clump (22 solar mass) with a radius of 0.11 pc surrounded by a cluster of 11 discrete infrared sources. The H2 column density and volume density in the center of this clump are 2 x 10(exp 22)/sq cm and 9 x 10(exp 5)/cu cm respectively. From a comparison between physical parameters derived from C(18)O and C(S-32) observations we conclude that the molecular envelope around the core has been completely removed by these sources and that only the very dense gas is left. OMC-3 shows a more complex elongated structure in C(18)O and CS than OMC-2. The C(S-32) and C(S-34) maps show that the denser region can be separated into at least sub-cores of roughly equal sizes (radius approximately equals 0.13 pc), with n(H2) = 6 x 10(exp 5)/cu cm, and a mass of 10 solar mass (from C(S-32)). The very different masses obtained for the central core from C(18)O and C(S-32) (55 and 12 solar mass respectively) indicate that a massive envelope is still present around the very dense sub-cores. We report the first detection of several molecular outflows in OMC-3. The presence of an IRAS source and the first detection of these outflows confirm that star formation is going on in OMC-3. Based on the different physical properties of these regions compared with OMC-1, OMC-2 appears to be in an intermediate evolutionary stage between OMC-1 and OMC-3.

Castets, A.↗
NASA NTRS2001

Designing Adaptive Low-Dissipative High Order Schemes for Long-Time Integrations

A general framework for the design of adaptive low-dissipative high order schemes is presented. It encompasses a rather complete treatment of the numerical approach based on four integrated design criteria: (1) For stability considerations, condition the governing equations before the application of the appropriate numerical scheme whenever it is possible; (2) For consistency, compatible schemes that possess stability properties, including physical and numerical boundary condition treatments, similar to those of the discrete analogue of the continuum are preferred; (3) For the minimization of numerical dissipation contamination, efficient and adaptive numerical dissipation control to further improve nonlinear stability and accuracy should be used; and (4) For practical considerations, the numerical approach should be efficient and applicable to general geometries, and an efficient and reliable dynamic grid adaptation should be used if necessary. These design criteria are, in general, very useful to a wide spectrum of flow simulations. However, the demand on the overall numerical approach for nonlinear stability and accuracy is much more stringent for long-time integration of complex multiscale viscous shock/shear/turbulence/acoustics interactions and numerical combustion. Robust classical numerical methods for less complex flow physics are not suitable or practical for such applications. The present approach is designed expressly to address such flow problems, especially unsteady flows. The minimization of employing very fine grids to overcome the production of spurious numerical solutions and/or instability due to under-resolved grids is also sought. The incremental studies to illustrate the performance of the approach are summarized. Extensive testing and full implementation of the approach is forthcoming. The results shown so far are very encouraging.

Yee, Helen C.↗
NASA NTRS2005

Designing Adaptive Low-Dissipative High Order Schemes for Long-Time Integrations

A general framework for the design of adaptive low-dissipative high order schemes is presented. It encompasses a rather complete treatment of the numerical approach based on four integrated design criteria: (1) For stability considerations, condition the governing equations before the application of the appropriate numerical scheme whenever it is possible. (2) For consistency, compatible schemes that possess stability properties, including physical and numerical boundary condition treatments, similar to those of the discrete analogue of the continuum are preferred. (3) For the minimization of numerical dissipation contamination, efficient and adaptive numerical dissipation control to further improve nonlinear stability and accuracy should be used. (4) For practical considerations, the numerical approach should be efficient and applicable to general geometries, and an efficient and reliable dynamic grid adaptation should be used if necessary. These design criteria are, in general, very useful to a wide spectrum of flow simulations. However, the demand on the overall numerical approach for non-linear stability and accuracy is much more stringent for long-time integration of complex multiscale viscous shock/shear/turbulence/acoustics interactions and numerical combustion. Robust classical numerical methods for less complex flow physics are not suitable or practical for such applications. The present approach is designed expressly to address such flow problems, especially unsteady flows.

Yee, H. C.↗
NASA NTRS2014

Application of Model-based Prognostics to a Pneumatic Valves Testbed

Pneumatic-actuated valves play an important role in many applications, including cryogenic propellant loading for space operations. Model-based prognostics emphasizes the importance of a model that describes the nominal and faulty behavior of a system, and how faulty behavior progresses in time, causing the end of useful life of the system. We describe the construction of a testbed consisting of a pneumatic valve that allows the injection of faulty behavior and controllable fault progression. The valve opens discretely, and is controlled through a solenoid valve. Controllable leaks of pneumatic gas in the testbed are introduced through proportional valves, allowing the testing and validation of prognostics algorithms for pneumatic valves. A new valve prognostics approach is developed that estimates fault progression and predicts remaining life based only on valve timing measurements. Simulation experiments demonstrate and validate the approach.

prediction↗
NASA NTRS2001

Mechanical behavior in living cells consistent with the tensegrity model

Alternative models of cell mechanics depict the living cell as a simple mechanical continuum, porous filament gel, tensed cortical membrane, or tensegrity network that maintains a stabilizing prestress through incorporation of discrete structural elements that bear compression. Real-time microscopic analysis of cells containing GFP-labeled microtubules and associated mitochondria revealed that living cells behave like discrete structures composed of an interconnected network of actin microfilaments and microtubules when mechanical stresses are applied to cell surface integrin receptors. Quantitation of cell tractional forces and cellular prestress by using traction force microscopy confirmed that microtubules bear compression and are responsible for a significant portion of the cytoskeletal prestress that determines cell shape stability under conditions in which myosin light chain phosphorylation and intracellular calcium remained unchanged. Quantitative measurements of both static and dynamic mechanical behaviors in cells also were consistent with specific a priori predictions of the tensegrity model. These findings suggest that tensegrity represents a unified model of cell mechanics that may help to explain how mechanical behaviors emerge through collective interactions among different cytoskeletal filaments and extracellular adhesions in living cells.

NASA Discipline Cell Biology↗
NASA NTRS1982

Ordered mixed-layer structures in the Mighei carbonaceous chondrite matrix

High resolution transmission electron microscopy of the Mighei carbonaceous chondrite matrix has revealed the presence of a new mixed layer structure material. This mixed-layer material consists of an ordered arrangement of serpentine-type (S) and brucite-type (B) layers in the sequence SBBSBB. Electron diffraction and imaging techniques show that the basal periodicity is approximately 17 A. Discrete crystals of SBB-type material are typically curved, of small size (less than 1 micron) and show structural variations similar to the serpentine group minerals. Mixed-layer material also occurs in association with planar serpentine. Characteristics of SBB-type material are not consistent with known terrestrial mixed-layer clay minerals. Evidence for formation by a condensation event or by subsequent alteration of pre-existing material is not yet apparent.

Mackinnon, I. D. R.↗
DOE OSTI2024

Towards dynamical low-rank approximation for neutrino kinetic equations. Part I: Analysis of an idealized relaxation model

Dynamical low-rank approximation (DLRA) is an emerging tool for reducing computational costs and provides memory savings when solving high-dimensional problems. Here, in this work, we propose and analyze a semi-implicit dynamical low-rank discontinuous Galerkin (DLR-DG) method for the space homogeneous kinetic equation with a relaxation operator, modeling the emission and absorption of particles by a background medium. Both DLRA and the discontinuous Galerkin (DG) scheme can be formulated as Galerkin equations. To ensure their consistency, a weighted DLRA is introduced so that the resulting DLR-DG solution is a solution to the fully discrete DG scheme in a subspace of the standard DG solution space. Similar to the standard DG method, we show that the proposed DLR-DG method is well-posed. We also identify conditions such that the DLR-DG solution converges to the equilibrium. Numerical results are presented to demonstrate the theoretical findings.

97 MATHEMATICS AND COMPUTING↗