Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time integrator”

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

Numerical investigations of low-density nozzle flow by solving the Boltzmann equation

A two-dimensional finite-difference code to solve the BGK-Boltzmann equation has been developed. The solution procedure consists of three steps: (1) transforming the BGK-Boltzmann equation into two simultaneous partial differential equations by taking moments of the distribution function with respect to the molecular velocity u(sub z), with weighting factors 1 and u(sub z)(sup 2); (2) solving the transformed equations in the physical space based on the time-marching technique and the four-stage Runge-Kutta time integration, for a given discrete-ordinate. The Roe's second-order upwind difference scheme is used to discretize the convective terms and the collision terms are treated as source terms; and (3) using the newly calculated distribution functions at each point in the physical space to calculate the macroscopic flow parameters by the modified Gaussian quadrature formula. Repeating steps 2 and 3, the time-marching procedure stops when the convergent criteria is reached. A low-density nozzle flow field has been calculated by this newly developed code. The BGK Boltzmann solution and experimental data show excellent agreement. It demonstrated that numerical solutions of the BGK-Boltzmann equation are ready to be experimentally validated.

Deng, Zheng-Tao↗
NASA NTRS2000

Progress in the Development of a Class of Efficient Low Dissipative High Order Shock-capturing Methods

In a series of papers, Olsson (1994, 1995), Olsson & Oliger (1994), Strand (1994), Gerritsen Olsson (1996), Yee et al. (1999a,b, 2000) and Sandham & Yee (2000), the issue of nonlinear stability of the compressible Euler and Navier-Stokes Equations, including physical boundaries, and the corresponding development of the discrete analogue of nonlinear stable high order schemes, including boundary schemes, were developed, extended and evaluated for various fluid flows. High order here refers to spatial schemes that are essentially fourth-order or higher away from shock and shear regions. The objective of this paper is to give an overview of the progress of the low dissipative high order shock-capturing schemes proposed by Yee et al. (1999a,b, 2000). This class of schemes consists of simple non-dissipative high order compact or non-compact central spatial differencings and adaptive nonlinear numerical dissipation operators to minimize the use of numerical dissipation. The amount of numerical dissipation is further minimized by applying the scheme to the entropy splitting form of the inviscid flux derivatives, and by rewriting the viscous terms to minimize odd-even decoupling before the application of the central scheme (Sandham & Yee). The efficiency and accuracy of these scheme are compared with spectral, TVD and fifth- order WENO schemes. A new approach of Sjogreen & Yee (2000) utilizing non-orthogonal multi-resolution wavelet basis functions as sensors to dynamically determine the appropriate amount of numerical dissipation to be added to the non-dissipative high order spatial scheme at each grid point will be discussed. Numerical experiments of long time integration of smooth flows, shock-turbulence interactions, direct numerical simulations of a 3-D compressible turbulent plane channel flow, and various mixing layer problems indicate that these schemes are especially suitable for practical complex problems in nonlinear aeroacoustics, rotorcraft dynamics, direct numerical simulation or large eddy simulation of compressible turbulent flows at various speeds including high-speed shock-turbulence interactions, and general long time wave propagation problems. These schemes, including entropy splitting, have also been extended to freestream preserving schemes on curvilinear moving grids for a thermally perfect gas (Vinokur & Yee 2000).

Yee, H. C.↗
DOE OSTI2021

Search for Multi-flare Neutrino Emissions in 10 yr of IceCube Data from a Catalog of Sources

A recent time-integrated analysis of a catalog of 110 candidate neutrino sources revealed a cumulative neutrino excess in the data collected by IceCube between 2008 April 6 and 2018 July 10. This excess, inconsistent with the background hypothesis in the Northern Hemisphere at the 3.3σ level, is associated with four sources: NGC 1068, TXS 0506+056, PKS 1424+240, and GB6 J1542+6129. Here, this Letter presents two time-dependent neutrino emission searches on the same data sample and catalog: a point-source search that looks for the most significant time-dependent source of the catalog by combining space, energy, and time information of the events, and a population test based on binomial statistics that looks for a cumulative time-dependent neutrino excess from a subset of sources. Compared to previous time-dependent searches, these analyses enable a feature to possibly find multiple flares from a single direction with an unbinned maximum-likelihood method. M87 is found to be the most significant time-dependent source of this catalog at the level of 1.7σ post-trial, and TXS 0506+056 is the only source for which two flares are reconstructed. The binomial test reports a cumulative time-dependent neutrino excess in the Northern Hemisphere at the level of 3.0σ associated with four sources: M87, TXS 0506+056, GB6 J1542+6129, and NGC 1068.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1982

An analysis of the multibeam altimeter

The design concept and computer simulations of a multibeam altimeter (MA) to provide satellite-based oceanographic-feature data for map construction within relatively short times are presented. The limitations of the pulse-limited radar altimeter used on GEOS-3 and SEASAT-A, in particular the narrow swath of each measurement pass, are characterized. The MA uses interferometry with two 1.2-m dish antennas separated by an 11-m boom to cover a total swath width of 100 km at satellite altitude 800 km, nominal repetition rate 6 kp/sec, and nominal integration time 2.5 sec. Simulation analysis of antenna element pattern, echo shape, tracking algorithms, tracking noise, and measurement accuracy was performed, and the general feasibility of ocean-current mapping with a one-week time delay, using two MA-equipped satellites, was demonstrated.

Bush, G. B.↗
NASA NTRS2007

The EBIT Calorimeter Spectrometer: A New, Permanent User Facility at the LLNL EBIT

The EBIT Calorimeter Spectrometer (ECS) has recently been completed and is currently being installed at the EBIT facility at the Lawrence Livermore National Laboratory. The ECS will replace the smaller XRS/EBIT spectrometer that has been in almost continuous operation since 2000. The XRS/EBIT was based on a spare laboratory cryostat and an engineering model detector system from the Suzaku/XRS observatory. The new ECS spectrometer was built from the ground up to be a low maintenance, high performance microcalorimeter spectrometer with 4 eV resolution at 6 keV, 32 detector channels, 10 us event timing, and capable of uninterrupted acquisition sessions of over 70 hours at 50 mK. The XRSIEBIT program has been extremely successful, producing over two-dozen refereed publications on topics such as laboratory astrophysics, atomic physics, nuclear physics, and calibration of the spectrometers for the National Ignition Facility, with many more publications in preparation. The ECS spectrometer will continue this work into the future with improved spectral resolution, integration times, and ease-of-use. We designed the ECS instrument with TES detectors in mind by using the same highly successful magnetic shielding as our laboratory TES cryostats. This design will lead to a future TES instrument at the LLNL EBIT. This proposed future instrument would include a hybrid detector system with 0.8 eV resolution in the band from 0.1-1.0 keV, 2 eV from 0.1-10 keV, and 30 eV from 0.5-100 keV, with high quantum efficiency in each band. Here we discuss the legacy of the XRS/EBIT program, the performance of the new ECS spectrometer, and plans for a future TES spectrometer.

Porter, S.↗
DOE OSTI2021

Discrete Green’s functions and spectral graph theory for computationally efficient thermal modeling

Here, this work concerns solutions of the heat equation with the spectral graph method, for which the temperature is defined at discrete points in the domain and the spatial relationship among the points is described by a graph. The heat equation on the graph is solved using matrix techniques involving the eigenvectors and eigenvalues of the Laplacian matrix. The spectral graph approach precludes the computationally intensive meshing and numerous time-integration steps of the finite element method. In the present work, the spectral graph method is extended to include heat loss at the boundaries with a generalized boundary condition, and physics-based edge weights are introduced which simplify the calibration process. From this approach a discrete Green’s function is defined which allows for solutions under a variety of heating conditions including: space-varying initial conditions; time-and-space varying internal heating; and, time-and-space-varying heating at boundaries of type 1 (Dirichlet), type 2 (Neumann) and type 3 (Robin). Results are provided for benchmark heat transfer problems in one spatial dimension and in three spatial dimensions, and verification is provided by comparison with exact analytical solutions and finite difference solutions. The spectral graph method converges within 0.4% error of the analytical solution. The practical utility of the approach is demonstrated by thermal simulation of a multilayer additive manufacturing process. The spectral graph results are compared to experimentally-obtained temperature data for two metal parts, with error less than 5% of the experimental measurements, with computation time less than one minute on a desktop computer.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Simulation of 24,000 Electron Dynamics: Real-Time Time-Dependent Density Functional Theory (TDDFT) with the Real-Space Multigrids (RMG)

Here, we present the theory, implementation, and benchmarking of a real-time time-dependent density functional theory (RT-TDDFT) module within the RMG code, designed to simulate the electronic response of molecular systems to external perturbations. Our method offers insights into nonequilibrium dynamics and excited states across a diverse range of systems, from small organic molecules to large metallic nanoparticles. Benchmarking results demonstrate excellent agreement with established TDDFT implementations and showcase the superior stability of our time integration algorithm, enabling long-term simulations with minimal energy drift. The scalability and efficiency of RMG on massively parallel architectures allow for simulations of complex systems, such as plasmonic nanoparticles with thousands of atoms. Future extensions, including nuclear and spin dynamics, will broaden the applicability of this RT-TDDFT implementation, providing a powerful toolset for studies of photoactive materials, nanoscale devices, and other systems where real-time electronic dynamics is essential.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

A higher-order finite-element implementation of the nonlinear Fokker–Planck collision operator for charged particle collisions in a low density plasma

Collisions between particles in a low density plasma are described by the Fokker–Planck collision operator. In applications, this nonlinear integro-differential operator is often approximated by linearised or ad-hoc model operators due to computational cost and complexity. In this work, we present an implementation of the nonlinear Fokker–Planck collision operator written in terms of Rosenbluth potentials in the Rosenbluth–MacDonald–Judd (RMJ) form. The Rosenbluth potentials may be obtained either by direct integration or by solving partial differential equations (PDEs) similar to Poisson's equation: we optimise for performance and scalability by using sparse matrices to solve the relevant PDEs. We represent the distribution function using a tensor-product continuous-Galerkin finite-element representation and we derive and describe the implementation of the weak form of the collision operator. We present tests demonstrating a successful implementation using an explicit time integrator and we comment on the speed and accuracy of the operator. Finally, we speculate on the potential for applications in the current and next generation of kinetic plasma models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

An immersed interface method for the 2D vorticity-velocity Navier-Stokes equations with multiple bodies

We present an immersed interface method for the vorticity-velocity form of the 2D Navier Stokes equations that directly addresses challenges posed by nonconvex immersed bodies, multiply connected domains, and the calculation of force distributions on immersed surfaces. The immersed interface method is re-interpreted as a polynomial extrapolation of flow quantities and boundary conditions into the immersed solid bodies, reducing computational cost and enabling simulations with nonconvex bodies that could not be discretized with previous immersed interface methods. In the flow, the vorticity transport equation is discretized using a conservative finite difference scheme and explicit Runge-Kutta time integration. The velocity reconstruction problem is transformed to a scalar Poisson equation that is discretized with conservative finite differences, and solved using an FFT-accelerated iterative algorithm. The use of conservative differencing throughout leads to exact enforcement of a discrete Kelvin's theorem, allowing for simulations with multiply connected domains and outflow boundaries that have challenged other immersed interface vortex methods. We also explore novel methods for recovering time-dependent pressure distributions on immersed bodies within a vorticity-based method and present a novel control volume formulation for recovering aerodynamic moments from only the vorticity and velocity fields. The method achieves second order spatial accuracy and third order temporal accuracy, and is validated on a variety of 2D flows in internal and free-space domains.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2010

Confidence-Based Feature Acquisition

Confidence-based Feature Acquisition (CFA) is a novel, supervised learning method for acquiring missing feature values when there is missing data at both training (learning) and test (deployment) time. To train a machine learning classifier, data is encoded with a series of input features describing each item. In some applications, the training data may have missing values for some of the features, which can be acquired at a given cost. A relevant JPL example is that of the Mars rover exploration in which the features are obtained from a variety of different instruments, with different power consumption and integration time costs. The challenge is to decide which features will lead to increased classification performance and are therefore worth acquiring (paying the cost). To solve this problem, CFA, which is made up of two algorithms (CFA-train and CFA-predict), has been designed to greedily minimize total acquisition cost (during training and testing) while aiming for a specific accuracy level (specified as a confidence threshold). With this method, it is assumed that there is a nonempty subset of features that are free; that is, every instance in the data set includes these features initially for zero cost. It is also assumed that the feature acquisition (FA) cost associated with each feature is known in advance, and that the FA cost for a given feature is the same for all instances. Finally, CFA requires that the base-level classifiers produce not only a classification, but also a confidence (or posterior probability).

Wagstaff, Kiri L.↗
NASA NTRS1988

A NASTRAN DMAP alter for linear buckling analysis under dynamic loading

A unique modification to the NASTRAN solution sequence for transient analysis with direct time integration (COSMIC NASTRAN rigid format 9) was developed and incorporated into a DMAP alter. This DMAP alter calculates the buckling stability of a dynamically loaded structure, and is used to predict the onset of structural buckling under stress wave loading conditions. The modified solution sequence incorporates the linear buckling analysis capability (rigid format 5) of NASTRAN into the existing Transient solution rigid format in such a way as to provide a time dependent eigensolution which is used to assess the buckling stability of the structure as it responds to the impulsive load. As a demonstration of the validity of this modified solution procedure, the dynamic buckling of a prismatic bar subjected to an impulsive longitudinal compression is analyzed and compared to the known theoretical solution. In addition, a dynamic buckling analysis is performed for the analytically less tractable problem of the localized dynamic buckling of an initially flawed composite laminate under transverse impact loading. The addition of this DMAP alter to the transient solution sequence in NASTRAN facilitates the prediction of both time and mode of buckling.

Aiello, Robert A.↗
DOE OSTI2025

Using quantum noise correlation analysis to measure ion temperature

Quantum noise correlation analysis was tested at the proof-of-concept level as a technique to measure ion temperature in a plasma. If eventually successful, this technique could enable a compact, inexpensive, and robust ion temperature diagnostic suitable for a burning plasma environment. Ion temperature is a key parameter determining the fusion performance of a burning plasma, as the fusion cross-section has a strong dependence on ion temperature. This ion temperature diagnostic would require only a small optical view of the plasma through a port to passively record impurity line emission. The instrumentation would be remote from the reactor behind the neutron and bio-shielding. The technique relies solely on quantum correlations of the photons emitted by a plasma impurity to measure ion temperature; there is no grating dispersion of an emission line-width or pulse-height analysis of photon energy. This measurement innovation was tested with instrumentation consisting of two single-photon detectors with high timing resolution, a time-tagging unit, and simple light collection optics. This instrumentation measures the second-order correlation between the light intensity falling on the two detectors. The next steps beyond the proof-of-concept level will be development of diagnostic designs for application of this technique to high-temperature and burning plasmas. Arrays of single-photon avalanche detectors to multiplex measurements of photon correlation will be required to reduce signal integration time to an acceptable duration.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1989

A NASTRAN DMAP alter for linear buckling analysis under dynamic loading

A modification to the NASTRAN solution sequence for transient analysis with direct time integration (COSMIC NASTRAN rigid format 9) was developed and incorporated into a DMAP alter. This DMAP alter calculates the buckling stability of a dynamically loaded structure, and is used to predict the onset of structural buckling under stress-wave loading conditions. The modified solution sequence incorporates the linear buckling analysis capability (rigid format 5) of NASTRAN into the existing Transient solution rigid format in such a way as to provide a time dependent eigensolution which is used to assess the buckling stability of the structure as it responds to the impulsive load. As a demonstration of the validity of this modified solution procedure, the dynamic buckling of a prismatic bar subjected to an impulsive longitudinal compression is analyzed and compared to the known theoretical solution. In addition, a dynamic buckling analysis is performed for the analytically less tractable problem of the localized dynamic buckling of an initially flawed composite laminate under transverse impact loading. The addition of this DMAP alter to the transient solution sequence in NASTRAN facilitates the computational prediction of both the time at which the onset of dynamic buckling occurs in an impulsively loaded structure, and the dynamic buckling mode shapes of that structure.

Aiello, Robert A.↗
NASA NTRS2021

PatchworkWave: A Multipatch Infrastructure for Multiphysics/Multiscale/Multiframe/Multimethod Simulations at Arbitrary Order

We present an extension of the PatchworkMHD code [1], itself an MHD-capable extension of thePatch-workcode [2], for which several algorithms presented here were co-developed. Its purpose is to create a“multipatch” scheme compatible with numerical simulations of arbitrary equations of motion at any dis-cretization order in space and time. In thePatchworkframework, the global simulation is comprised of anarbitrary number of moving, local meshes, or “patches”, which are free to employ their own resolution, co-ordinate system/topology, physics equations, reference frame, and in our new approach, numerical method.Each local patch exchanges boundary data with a single global patch on which all other patches residethrough a client-router-server parallelization model. In generalizingPatchworkto be compatible witharbitrary order time integration,PatchworkMHDandPatchworkWavehave significantly improved theinterpatch interpolation accuracy by removing an interpolation of interpolated data feedback present in theoriginalPatchworkcode. Furthermore, we extendPatchworkto bemultimethodby allowing multiplestate vectors to be updated simultaneously, with each state vector providing its own interpatch interpolationand transformation procedures. As such, our scheme is compatible with nearly any set of hyperbolic partialdifferential equations. We demonstrate our changes through the implementation of a scalar wave toy-modelthat is evolved on arbitrary, time dependent patch configurations at 4th order accuracy.

Dennis B Bowen↗
NASA NTRS2005

Computations of Flow over a Hump Model Using Higher Order Method with Turbulence Modeling

Turbulent separated flow over a two-dimensional hump is computed by solving the RANS equations with k - omega (SST) turbulence model for the baseline, steady suction and oscillatory blowing/suction flow control cases. The flow equations and the turbulent model equations are solved using a fifth-order accurate weighted essentially. nonoscillatory (WENO) scheme for space discretization and a third order, total variation diminishing (TVD) Runge-Kutta scheme for time integration. Qualitatively the computed pressure distributions exhibit the same behavior as those observed in the experiments. The computed separation regions are much longer than those observed experimentally. However, the percentage reduction in the separation region in the steady suction case is closer to what was measured in the experiment. The computations did not predict the expected reduction in the separation length in the oscillatory case. The predicted turbulent quantities are two to three times smaller than the measured values pointing towards the deficiencies in the existing turbulent models when they are applied to strong steady/unsteady separated flows.

Balakumar, P.↗
NASA NTRS2019

Deep Mapping of Small Solar System Bodies with Galactic Cosmic Ray Secondary Particle Showers

We will investigate the use of galactic cosmic ray (GCR) secondary particles to probe the deep interiors of small solar system bodies (SSBs), including comets, asteroids, and geologic structures on the surfaces of airless bodies. Applications include solar system science, planetary defense, and resource utilization. Our Phase I study demonstrated that muons, the long-range charged component of GCR showers, can penetrate SSBs up to a km in diameter, providing information on their interior structure. Muons produced in Earth’s atmosphere have been applied to image the interior of large objects for science and engineering. In Phase I, we found that the production of muons in the solid surfaces of airless bodies is much smaller than in Earth’s atmosphere. Nevertheless, the flux of transmitted muons is sufficient to detect inclusions within an asteroid or comet in a reasonable amount of time, ranging from hours to weeks, depending on the size of the SSB and the density contrast, position and size of the inclusion. For asteroids and comets, large density variations (e.g., porous soil or ice versus solid rock) are relatively easy to detect. The intrinsic spatial resolution of muon radiography (“muography”) is on the scale of a few meters. The spatial resolution that can be achieved in practice depends on signal intensity and integration time (counting statistics), the angular resolution of the muon tracker (hodoscope) and details of data reduction and analysis methodology. Our Phase II project will assess remaining unknowns for the application of muography to determining the interior structure of SSBs, assess risks for implementation, and provide a roadmap for development of SSB muography beyond the NIAC program. To achieve our objectives, we will focus on four interrelated tasks: Task1) Signal and background characterization: Characterize the production and transmission of muons and secondary particle backgrounds made by cosmic ray showers in SSBs; and near-surface features from radiographic and tomographic data; Task2) Imaging studies: Develop methods to determine the density structure of SSB interiors and near-surface features from radiographic and tomographic data; Task3) Instrument design: Using simulations and bench-top laboratory experiments, investigate specific concepts for the design of compact hodoscopes and components; Task4) Synthesis: Combine the results of the first three tasks to determine the range of applicability of the method, identify the steps needed for maturation of the concept, and explore concepts for a pilot muography mission.

Prettyman, Thomas H.↗
NASA NTRS2004

Computations of Flow Over a Hump Model Using Higher Order Method With Turbulence Modeling

Turbulent separated flow over a two-dimensional hump is computed by solving the RANS equations with k-omega (SST) turbulence model for the baseline, steady suction and oscillatory blowing/suction flow control cases. The flow equations and the turbulent model equations are solved using a fifth-order accurate weighted essentially nonoscillatory (WENO) scheme for space discretization and a third order, total variation diminishing (TVD) Runge-Kutta scheme for time integration. Qualitatively the computed pressure distributions exhibit the same behavior as they are observed in the experiments. The computed separation regions are much longer than that are observed. However, the percentage reduction in the separation region in the steady suction case is closer to that was measured in the experiment. The computations did not predict the expected reduction in the separation length in the oscillatory case. The predicted turbulent quantities are two to three times smaller than that are measured and it points towards the deficiencies in the existing turbulent models when they are applied to strong steady/unsteady separated flows.

Balakumar, Ponnampalam↗
NASA NTRS2017

Single Object & Time Series Spectroscopy with JWST NIRCam

JWST will enable high signal-to-noise spectroscopic observations of the atmospheres of transiting planets with high sensitivity at wavelengths that are inaccessible with HST or other existing facilities. We plan to exploit this by measuring abundances, chemical compositions, cloud properties, and temperature-pressure parameters of a set of mostly warm (T 600 - 1200 K) and low mass (14 -200 Earth mass) planets in our guaranteed time program. These planets are expected to have significant molecular absorptions of H2O, CH4, CO2, CO, and other molecules that are key for determining these parameters and illuminating how and where the planets formed. We describe how we will use the NIRCam grisms to observe slitless transmission and emission spectra of these planets over 2.4 - 5.0 microns wavelength and how well these observations can measure our desired parameters. This will include how we set integration times, exposure parameters, and obtain simultaneous shorter wavelength images to track telescope pointing and stellar variability. We will illustrate this with specific examples showing model spectra, simulated observations, expected information retrieval results, completed Astronomer's Proposal Tools observing templates, target visibility, and other considerations.

NIRCam↗