Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial 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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

An adaptive mesh-moving and refinement procedure for one-dimensional conservation laws

We examine the performance of an adaptive mesh-moving and /or local mesh refinement procedure for the finite difference solution of one-dimensional hyperbolic systems of conservation laws. Adaptive motion of a base mesh is designed to isolate spatially distinct phenomena, and recursive local refinement of the time step and cells of the stationary or moving base mesh is performed in regions where a refinement indicator exceeds a prescribed tolerance. These adaptive procedures are incorporated into a computer code that includes a MacCormack finite difference scheme wih Davis' artificial viscosity model and a discretization error estimate based on Richardson's extrapolation. Experiments are conducted on three problems in order to qualify the advantages of adaptive techniques relative to uniform mesh computations and the relative benefits of mesh moving and refinement. Key results indicate that local mesh refinement, with and without mesh moving, can provide reliable solutions at much lower computational cost than possible on uniform meshes; that mesh motion can be used to improve the results of uniform mesh solutions for a modest computational effort; that the cost of managing the tree data structure associated with refinement is small; and that a combination of mesh motion and refinement reliably produces solutions for the least cost per unit accuracy.

Biswas, Rupak↗
NASA NTRS2010

Diffuse-Illumination Systems for Growing Plants

Agriculture in both terrestrial and space-controlled environments relies heavily on artificial illumination for efficient photosynthesis. Plant-growth illumination systems require high photon flux in the spectral range corresponding with plant photosynthetic active radiation (PAR) (400 700 nm), high spatial uniformity to promote uniform growth, and high energy efficiency to minimize electricity usage. The proposed plant-growth system takes advantage of the highly diffuse reflective surfaces on the interior of a sphere, hemisphere, or other nearly enclosed structure that is coated with highly reflective materials. This type of surface and structure uniformly mixes discrete light sources to produce highly uniform illumination. Multiple reflections from within the domelike structures are exploited to obtain diffuse illumination, which promotes the efficient reuse of photons that have not yet been absorbed by plants. The highly reflective surfaces encourage only the plant tissue (placed inside the sphere or enclosure) to absorb the light. Discrete light sources, such as light emitting diodes (LEDs), are typically used because of their high efficiency, wavelength selection, and electronically dimmable properties. The light sources are arranged to minimize shadowing and to improve uniformity. Different wavelengths of LEDs (typically blue, green, and red) are used for photosynthesis. Wavelengths outside the PAR range can be added for plant diagnostics or for growth regulation

May, George↗
NASA NTRS2014

Entropy Stable Wall Boundary Conditions for the Compressible Navier-Stokes Equations

Non-linear entropy stability and a summation-by-parts framework are used to derive entropy stable wall boundary conditions for the compressible Navier-Stokes equations. A semi-discrete entropy estimate for the entire domain is achieved when the new boundary conditions are coupled with an entropy stable discrete interior operator. The data at the boundary are weakly imposed using a penalty flux approach and a simultaneous-approximation-term penalty technique. Although discontinuous spectral collocation operators are used herein for the purpose of demonstrating their robustness and efficacy, the new boundary conditions are compatible with any diagonal norm summation-by-parts spatial operator, including finite element, finite volume, finite difference, discontinuous Galerkin, and flux reconstruction schemes. The proposed boundary treatment is tested for three-dimensional subsonic and supersonic flows. The numerical computations corroborate the non-linear stability (entropy stability) and accuracy of the boundary conditions.

Parsani, Matteo↗
DOE OSTI2021

Characterization of mechanical discontinuities based on data-driven classification of compressional-wave travel times

Wave propagation and diffusive transport phenomena are influenced by the mechanical discontinuities in material. Here this study shows that certain bulk properties of the network of low-velocity mechanical discontinuities (e.g. air-filled cracks) in a material can be characterized by processing compressional-wave travel times using traditional data-driven classification techniques. To that end, we perform three tasks in chronological order: (1) use the discrete fracture network (DFN) method to create two-dimensional (2D) numerical models of crack-bearing material embedded with various types of low-velocity mechanical discontinuities, (2) use the fast marching method (FMM) to simulate the propagation of the wave/diffusion front from a single source through the 2D crack-bearing material to multiple receivers placed on the boundary of the material, and (3) train 9 data-driven classifiers to characterize the crack-bearing materials (i.e. bulk properties of the network of mechanical discontinuities in the crack-bearing material) by learning from the simulations of travel times detected by multiple receivers placed around the crack-bearing material. The classifiers identified the orientation, spatial distribution, and dispersion of the low-velocity mechanical discontinuities. Voting classifier performs the best among the 9 classifiers. For the characterization of bulk dispersion and distribution of discontinuities, the sensors located on the adjacent boundaries are more important; whereas for the characterization of bulk orientation of discontinuities, the sensors located on the opposite side are more important.

42 ENGINEERING↗
DOE OSTI2026

Magnetic brightening and nanoscale imaging of spin-polarized helical edge modes in ZrTe 5

Efficient charge transport remains a fundamental challenge for nanoelectronic devices, as their performance is constrained by high dissipation and the impedance mismatch between high-frequency signal sources and nanoscale circuit interfaces. Although topologically protected helical edge modes offer a dissipationless and backscattering-resilient alternative, achieving nanoscale control over these modes requires direct visualization of their spatial electrodynamics, especially under high-frequency operation. Here, by utilizing cryogenic magneto-infrared scattering-type scanning near-field optical microscopy (cm-IR-sSNOM), we image spin-polarized helical edge channels in ZrTe 5 at the nanoscale, revealing magnetic-field-induced brightening as a high-frequency electrodynamic signature of robust topological edge modes. Operating at 1.8 K and under a magnetic field of up to 5 T, we observe the emergence of edge-state polarizability at infrared frequencies that is notably resilient to the magnetic gaps that typically quench d.c. and microwave edge transport. Our results reveal a topological ‘two-lane’ spatial reorganization in which an external magnetic field induces a spin-population imbalance between counterpropagating edge modes. Favoured helical branches are confined against physical boundaries to activate a net infrared near-field contrast. This electrodynamic response scales linearly with the number of atomic layers, which confirms that individual layers in ZrTe 5 preserve their discrete quantum spin Hall identities. These findings may be useful for developing topological spintronic devices, as the magnetic infrared tunability of helical edge channels provides a pathway for low-loss nanoscale interconnects and high-speed information processing.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Sublattice modulated superconductivity in the kagome Hubbard model

Here we identify a superconducting order featuring spatial pair modulations on the kagome lattice subject to onsite Hubbard U and nearest neighbor V interactions. Within our functional renormalization group analysis, this state appears with a concomitant d-wave superconducting (SC) instability at zero lattice momentum, where it distinguishes itself through intra-unit cell modulations of the pairing function thus breaking the discrete space group symmetry. The relative weight of the sublattice modulated superconductor (SMS) and d-wave SC is influenced by the absolute interaction strength and coupling ratio V/U. Parametrically adjacent to this domain at weak coupling, we find an intra-unit cell modulated vestigial charge density wave and an s-wave SC instability. Our study provides a microscopic setting and thorough description of this novel SMS arising within a translation symmetry broken background.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Relaxation of a single defect spin by the low-frequency gyrotropic mode of a magnetic vortex

We excite the gyrotropic mode of a magnetic vortex and observe the resulting effect on the spin state of a nearby nitrogen-vacancy (NV) defect in diamond. Thin permalloy disks fabricated on a diamond sample are magnetized in a vortex state in which the magnetization curls around a central core. The magnetization dynamics of this configuration are described by a discrete spectrum of confined magnon modes as well as a low-frequency gyrotropic mode in which the vortex core precesses about its equilibrium position. Despite the spin transition frequencies being far-detuned from the modes of the ferromagnet, we observe enhanced relaxation of the NV spin when driving the gyrotropic mode. Moreover, we map the spatial dependence of the interaction between the vortex and the spin by translating the vortex core within the disk with an applied magnetic field, resulting in steplike motion as the vortex is pinned and de-pinned. Strong spin relaxation is observed when the vortex core is within approximately 250 nm of the NV center defect. We attribute this effect to the higher frequencies in the spectrum of the magnetic fringe field arising from the soliton-like nature of the gyrotropic mode when driven with sufficiently large amplitude.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1978

Numerical method for solution of systems of non-stationary spatially one-dimensional nonlinear differential equations

A computational scheme and a standard program is proposed for solving systems of nonstationary spatially one-dimensional nonlinear differential equations using Newton's method. The proposed scheme is universal in its applicability and its reduces to a minimum the work of programming. The program is written in the FORTRAN language and can be used without change on electronic computers of type YeS and BESM-6. The standard program described permits the identification of nonstationary (or stationary) solutions to systems of spatially one-dimensional nonlinear (or linear) partial differential equations. The proposed method may be used to solve a series of geophysical problems which take chemical reactions, diffusion, and heat conductivity into account, to evaluate nonstationary thermal fields in two-dimensional structures when in one of the geometrical directions it can take a small number of discrete levels, and to solve problems in nonstationary gas dynamics.

Morozov, S. K.↗
NASA NTRS2015

Entropy Stable Wall Boundary Conditions for the Three-Dimensional Compressible Navier-Stokes Equations

Non-linear entropy stability and a summation-by-parts framework are used to derive entropy stable wall boundary conditions for the three-dimensional compressible Navier-Stokes equations. A semi-discrete entropy estimate for the entire domain is achieved when the new boundary conditions are coupled with an entropy stable discrete interior operator. The data at the boundary are weakly imposed using a penalty flux approach and a simultaneous-approximation-term penalty technique. Although discontinuous spectral collocation operators on unstructured grids are used herein for the purpose of demonstrating their robustness and efficacy, the new boundary conditions are compatible with any diagonal norm summation-by-parts spatial operator, including finite element, finite difference, finite volume, discontinuous Galerkin, and flux reconstruction/correction procedure via reconstruction schemes. The proposed boundary treatment is tested for three-dimensional subsonic and supersonic flows. The numerical computations corroborate the non-linear stability (entropy stability) and accuracy of the boundary conditions.

Parsani, Matteo↗
DOE OSTI2022

Continuous Emulation and Multiscale Visualization of Traffic Flow Using Stationary Roadside Sensor Data

With the advent of the next-generation traffic monitoring systems, there has been a significant increase in the spatial-temporal resolution of vehicle mobility data in many cities. Effective analysis and visualization of such data can provide transportation planners with data-driven insights, which can facilitate the understanding of multiscale traffic dynamics. In this paper, we present a web-based traffic emulator for emulating and visualizing near-real-time and historical traffic flows on highways using data from road-side sensors. To construct a continuous traffic flow, the emulator adopts an analytical pipeline that can (a) integrate traffic data collected from discrete road-side radar detection sensors, (b) interpolate traffic conditions (vehicle speed and volume) on unmeasured road segments based on traffic flow theory, and (c) generate lane-specific vehicle trajectories and movements using a mathematically optimized representation of the road network. Our app also provides an integrated visual workflow that allows users to explore the interconnected traffic dynamics using an appropriate traffic flow visualization selected based on the level of detail. We devise two innovative geo-visualization techniques that utilize an animated strips-network representation and a lane usage matrix to visualize lane performances. To ensure a smooth emulation of large-scale traffic flow in an easy-to-access web environment, we implement the emulator using client-side GPU-accelerated techniques. Lastly, we close with a case study that visualizes traffic dynamics of two scenarios - an afternoon peak hour and a traffic accident - in Chattanooga, Tennessee. Our app visualizes the responses of traffic dynamics during different traffic conditions, and to the presence of the traffic accident at different spatial scales.

42 ENGINEERING↗
NASA NTRS1996

Aerodynamic Characteristics of an Aerospace Vehicle During a Subsonic Pitch-Over Maneuver

Time-dependent CFD has been used to predict aerospace vehicle aerodynamics during a subsonic rotation maneuver. The inviscid 3D3U code is employed to solve the 3-D unsteady flow field using an unstructured grid of tetrahedra. As this application represents a challenge to time-dependent CFD, observations concerning spatial and temporal resolution are included. It is shown that even for a benign rotation rate, unsteady aerodynamic effects are significant during the maneuver. Possibly more significant, however, the rotation maneuver creates ow asymmetries leading to yawing moment, rolling moment, and side force which are not present in the quasi-steady case. A series of steady solutions at discrete points in the maneuver are also computed for comparison with wind tunnel measurements and as a means of quantifying unsteady effects.

Kleb, William L.↗
NASA NTRS2010

Sources and Losses of Ring Current Ions

During geomagnetic quiet times, in-situ measurements of ring current energetic ions (few to few tens of keVs) from THEMIS spacecraft often exhibit multiple ion populations at discrete energies that extend from the inner magnetosphere to the magnetopause at dayside or plasma sheet at nightside. During geomagnetic storm times, the levels of fluxes as well as the mean energies of these ions elevated dramatically and the more smooth distributions in energies and distances during quiet times are disrupted into clusters of ion populations with more confined spatial extends. This reveals local plasma heating processes that might have come into play. Several processes have been proposed. Magnetotail dipolarization, sudden enhancement of field-aligned current, local current disruptions, and plasma waves are possible mechanisms to heat the ions locally as well as strong convections of energetic ions directly from the magnetotail due to reconnections. We will examine two geomagnetic storms on October 11, 2008 and July 22, 2009 to reveal possible heating mechanisms. We will analyze in-situ plasma and magnetic field measurements from THEMIS, GOES, and DMSP for the events to study the ion pitch angle distributions and magnetic field perturbations in the auroral ionosphere and inner magnetosphere where the plasma heating processes occur.

Chen, Sheng-Hsien↗
DOE OSTI2020

Effect of the numerical discretization scheme in Shock-Driven turbulent mixing simulations

In this work, we evaluate the effects of distinct numerical strategies in simulations of shock-driven turbulent mixing. An air-SF 6 -air gas curtain subjected to Mach 1.2 shock-waves is computed with Implicit Large-Eddy Simulation based on three numerical schemes: directional-split with Harten-Lax-van Leer (HLL) solver; directional-unsplit using HLL-Contact (HLLC) solver; and directional-unsplit with HLLC solver and a Low Mach number Correction (LMC). The results illustrate the importance of the numerical strategy to the accuracy of the predictions. Whereas both split and unsplit schemes result in a similar spatial development of the initial shock-driven instability, only the unsplit schemes can predict the turbulent mixing transition after reshock observed by the laboratory experiments. Such feature increases the mixing rate of the two fluids, this being particularly pronounced when the LMC is active due to i) the reduced flow Mach number, and ii) the larger effective Reynolds number. Since the selected mixing problem is driven by the deposition of vorticity at the fluids’ interface, the resultant flow physics is analyzed by investigating the contribution of distinct inviscid mechanisms to the production of vorticity: baroclinicity, stretching, and dilatation. As expected, it is observed that the production of vorticity is initially dominated by the baroclinicity mechanism. Yet, the relevance of the remainder mechanisms is enhanced after the reshock and may even surpass the baroclinicity term.

42 ENGINEERING↗
DOE OSTI2021

StressNet - Deep learning to predict stress with fracture propagation in brittle materials

Abstract Catastrophic failure in brittle materials is often due to the rapid growth and coalescence of cracks aided by high internal stresses. Hence, accurate prediction of maximum internal stress is critical to predicting time to failure and improving the fracture resistance and reliability of materials. Existing high-fidelity methods, such as the Finite-Discrete Element Model (FDEM), are limited by their high computational cost. Therefore, to reduce computational cost while preserving accuracy, a deep learning model, StressNet, is proposed to predict the entire sequence of maximum internal stress based on fracture propagation and the initial stress data. More specifically, the Temporal Independent Convolutional Neural Network (TI-CNN) is designed to capture the spatial features of fractures like fracture path and spall regions, and the Bidirectional Long Short-term Memory (Bi-LSTM) Network is adapted to capture the temporal features. By fusing these features, the evolution in time of the maximum internal stress can be accurately predicted. Moreover, an adaptive loss function is designed by dynamically integrating the Mean Squared Error (MSE) and the Mean Absolute Percentage Error (MAPE), to reflect the fluctuations in maximum internal stress. After training, the proposed model is able to compute accurate multi-step predictions of maximum internal stress in approximately 20 seconds, as compared to the FDEM run time of 4 h, with an average MAPE of 2% relative to test data.

36 MATERIALS SCIENCE↗
NASA NTRS2010

Development and Verification of the Charring Ablating Thermal Protection Implicit System Solver

The development and verification of the Charring Ablating Thermal Protection Implicit System Solver is presented. This work concentrates on the derivation and verification of the stationary grid terms in the equations that govern three-dimensional heat and mass transfer for charring thermal protection systems including pyrolysis gas flow through the porous char layer. The governing equations are discretized according to the Galerkin finite element method with first and second order implicit time integrators. The governing equations are fully coupled and are solved in parallel via Newton's method, while the fully implicit linear system is solved with the Generalized Minimal Residual method. Verification results from exact solutions and the Method of Manufactured Solutions are presented to show spatial and temporal orders of accuracy as well as nonlinear convergence rates.

Amar, Adam J.↗
DOE OSTI2022

Massively distributed fiber strain sensing using Brillouin lasing

Brillouin based distributed fiber sensors present a unique set of characteristics amongst fiber sensing architectures. They are able to measure absolute strain and temperature over long distances, with high spatial resolution, and very large dynamic range in off-the-shelf fiber. However, Brillouin sensors traditionally provide only modest sensitivity due to the weak dependence of the Brillouin frequency on strain and the high signal to noise ratio required to identify the resonance’s peak frequency to within a small fraction of its linewidth. Recently, we introduced a technique which substantially improves the precision of Brillouin fiber sensors by exciting a series of lasing modes in a fiber loop cavity that experience Brillouin amplification at discrete locations in the fiber. The narrow-linewidth and high intensity of the lasing modes enabled ultra-low noise Brillouin sensors with large dynamic range. However, our initial demonstration was only modestly distributed: measuring strain at 40, non-contiguous positions along a 400 m fiber. In this work, we greatly extend this methodology to enable fully distributed sensing at 1000 contiguous locations along 3.5 km of fiber—an order of magnitude increase in sensor count and range. This highly-multiplexed Brillouin fiber laser sensor provides a strain noise as low as 34 nɛ/√Hz and we analyze the limiting factors in this approach.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Multiscale design of nonlinear materials using a Eulerian shape optimization scheme

Motivated by recent advances in manufacturing, the design of materials is the focal point of interest in the material research community. One of the critical challenges in this field is finding optimal material microstructure for a desired macroscopic response. This work presents a computational method for the mesoscale-level design of particulate composites for an optimal macroscale-level response. The method relies on a custom shape optimization scheme to find the extrema of a nonlinear cost function subject to a set of constraints. Three key “modules” constitute the method: multiscale modeling, sensitivity analysis, and optimization. Multiscale modeling relies on a classical homogenization method and a nonlinear NURBS-based generalized finite element scheme to efficiently and accurately compute the structural response of particulate composites using a nonconformal discretization. A three-parameter isotropic damage law is used to model microstructure-level failure. An analytical sensitivity method is developed to compute the derivatives of the cost/constraint functions with respect to the design variables that control the microstructure's geometry. The derivation uncovers subtle but essential new terms contributing to the sensitivity of finite element shape functions and their spatial derivatives. Several structural problems are solved to demonstrate the applicability, performance, and accuracy of the method for the design of particulate composites with a desired macroscopic nonlinear stress-strain response.

42 ENGINEERING↗
DOE OSTI2020

Plasticity of irradiated materials at the nano and micro-scales

Here, we review here our recent work on plastic deformation in irradiated materials at the nano- and micro-scales, as revealed by Discrete Dislocation Dynamics (DDD) simulations. Two methods of including irradiation effects in the DDD framework are presented. The first directly captures the atomistic interaction mechanisms, while the second can effectively study high-dose irradiation. Computer simulations lead to new understanding of the dynamics of collective dislocation-irradiation defect interactions, as well as the quantitative analysis of the temporal and spatial characteristics associated with plastic instabilities. Based on these insights, theoretical models are developed to predict the critical conditions for dislocation channel formation. A simple probability model is proposed and demonstrated to predict the width of dislocation channels in bulk irradiated materials with good agreement with experimental data. The fundamental understanding of the origins of plastic flow localization in irradiated materials sheds light on the design of future generations of radiation-resistant materials.

36 MATERIALS SCIENCE↗