Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Direct numerical simulation”

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,351 records · Page 75
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Reliable Devices Yield Stable Quantum Computations

Stable quantum computation requires noisy results to remain bounded even in the presence of noise fluctuations. Yet non-stationary noise processes lead to drift in the varying characteristics of a quantum device that can greatly influence the circuit outcomes. Here we address how temporal and spatial variations in noise relate device reliability to quantum computing stability. First, our approach quantifies the differences in statistical distributions of characterization metrics collected at different times and locations using Hellinger distance. We then validate an analytical bound that relates this distance directly to the stability of a computed expectation value. Our demonstration uses numerical simulations with models informed by the washington superconducting transmon device. We find that the stability metric is consistently bounded from above by the corresponding Hellinger distance, which can be cast as a specified tolerance level. These results underscore the significance of reliable quantum computing devices and the impact for stable quantum computation.

Dasgupta, Samudra↗
NASA NTRS1986

Lift augmentation via spanwise tip blowing - A numerical study

Numerical simulations of a low aspect ratio wing with and without a spanwise directed jet issuing from the wing tip have been performed. The results show that the tip vortex is displaced outward and upward by the blowing. This gives rise to a local lift augmentation mechanism, vortex lift caused by the vortex core being above the wing, and a global mechanism, the reduction of induced velocities due to greater apparent spin.

Childs, R. E.↗
NASA NTRS1989

A multi-domain spectral computation of three-dimensional laminar horseshoe vortex flow using incompressible Navier-Stokes equations

Multidomain spectral methods are presently used to numerically simulate a strut-wall intersection's laminar horseshoe vortex flow through direct solution of the three-dimensional, incompressible, time-dependent Navier-Stokes equations. Direct expansion in Chebyshev polynomials and spectral element method spatial discretization of flow dependence are used to achieve high-order accuracy, and minimal dispersion errors. Low and moderate Reynolds number results are presented to illustrate the method application.

Tan, C. S.↗
NASA NTRS1990

Mathematical modeling of the flow field and particle motion in a rotating bioreactor at unit gravity and microgravity

The biotechnology group at NASA Johnson Space Center is developing systems for culturing mammalian cells that stimulate some aspect of microgravity and provide a low shear environment for microgravity-based studies on suspension and anchorage dependent cells. The design of these vessels for culturing cells is based on the need to suspend cells and aggregates of cells and microcarrier beads continually in the culturing medium. The design must also provide sufficient circulation for adequate mass transfer of nutrients to the cells and minimize the total force on the cells. Forces, resulting from sources such as hydrodynamic fluid shear and collisions of cells and walls of the vessels, may damage delicate cells and degrade the formation of three dimensional structures. This study examines one particular design in both unit gravity and microgravity based on two concentric cylinders rotating in the same direction at different speeds to create a Couette flow between them. A numerical simulation for the flow field and the trajectories of particles in the vessel. The flow field for the circulation of the culturing medium is modeled by the Navier-Stokes equations. The forces on a particle are assumed to be drag from the fluid's circulation, buoyancy from the gravitational force and centrifugal force from the rotation of the vessel. The problem requires first solving the system of partial differential equations for the fluid flow by a finite difference method and then solving the system of ordinary differential equations for the trajectories by Gear's stiff method. Results of the study indicate that the trajectories in unit gravity and microgravity are very similar except for small spatial deviations on the fast time scale in unit gravity. The total force per unit cross sectional area on a particle in microgravity, however, is significantly smaller than the corresponding value in unit gravity, which is also smaller than anticipated. Hence, this study indicates that this design for a bioreactor with optimal rates of rotation can provide a good environment for culturing cells in microgravity with adequate circulation and minimal force on the cells.

Boyd, Ernest J.↗
NASA NTRS1991

Testing continuum descriptions of low-Mach-number shock structures

Numerical experiments have been performed on normal shock waves with Monte Carlo Direct Simulations (MCDS's) to investigate the validity of continuum theories at very low Mach numbers. Results from the Navier-Stokes and the Burnett equations are compared to MCDS's for both hard-sphere and Maxwell gases. It is found that the maximum-slope shock thicknesses are described equally well (within the MCDS computational scatter) by either of the continuum formulations for Mach numbers smaller than about 1.2. For Mach numbers greater that 1.2, the Burnett predictions are more accurate than the Navier-Stokes results. Temperature-density profile separations are best described by the Burnett equations for Mach numbers greater than about 1.3. At lower Mach numbers the MCDS scatter is too great to differentiate between the two continuum theories. For all Mach numbers above one, the shock shapes are more accurately described by the Burnett equations.

Pham-Van-diep, Gerald C.↗
NASA NTRS1992

Preliminary results of a Monte Carlo simulation of burst directions by BATSE

The effect of varying BATSE trigger parameters on the modeled BATSE detection rate is investigated numerically. Quantitative estimates are obtained of the systematic effects generated by the complex interplay of trigger times, variations in the signal-to-noise ratio for burst detection as well as the time scale over which backgrounds are determined, and the number of required detections above threshold. It is found, in particular that employing a third detection, at a reduced signal-to-noise ratio, to resolve a burst does not increase the detection rate. It is also shown that the BATSE detection rate is insensitive to variations in the time scale over which the background is determined.

Higdon, J. C.↗
NASA NTRS1995

Theoretical properties of Omega-loops in the convective zone of the Sun. 3: Extended updrafts

It was pointed out in an earlier paper that the continuing emergence of Omega-loops at localized sites on the surface of the Sun indicates a continuing updraft at those sites. The updraft evidently extends all the way from the base of the convective zone to somewhat near (approximately 10(exp 9)cm) the surface. We pointed out that such updrafts enhance the convective heat transport to the surface, accounting for a major part of the increased solar brightness or irradiance during times of solar activity. The problem is to work out, as much as possible, the dynamical nature of the extended updrafts, initiated as the wakes of successive rising Omega-loops and driven therafter by the convective forces. The question is, does the updraft take on a long-lived columnar form of its own, or does it never devlop beyond a sequence of rising wakes, resembling beads on a string? The dynamics of a columnar updraft is complicated by both the large Reynolds number and the strong stratication of the atmosphere, and by a total lack of direct observational information. Extended slender updrafts are not a spontaneous occurrence in numerical simulations of thermal convection in a stratified atmosphere, although slender concentrated downdrafts commonly occur. This paper examines several aspects of a columnar updrft in a convective atmosphere under various idealized circumstances to investigate to what extent that state can be maintained against the diminishing vorticity and expansion in the updraft. It appears that the successive passage of Omega-loops from the bottom to the top of the convective zone is an essential feature of the continuing existence of the updraft.

Parker, E. N.↗
NASA NTRS2016

Dynamic Beam Solutions for Real-Time Simulation and Control Development of Flexible Rockets

In this study, flexible rockets are structurally represented by linear beams. Both direct and indirect solutions of beam dynamic equations are sought to facilitate real-time simulation and control development for flexible rockets. The direct solution is completed by numerically integrate the beam structural dynamic equation using an explicit Newmark-based scheme, which allows for stable and fast transient solutions to the dynamics of flexile rockets. Furthermore, in the real-time operation, the bending strain of the beam is measured by fiber optical sensors (FOS) at intermittent locations along the span, while both angular velocity and translational acceleration are measured at a single point by the inertial measurement unit (IMU). Another study in this paper is to find the analytical and numerical solutions of the beam dynamics based on the limited measurement data to facilitate the real-time control development. Numerical studies demonstrate the accuracy of these real-time solutions to the beam dynamics. Such analytical and numerical solutions, when integrated with data processing and control algorithms and mechanisms, have the potential to increase launch availability by processing flight data into the flexible launch vehicle's control system.

flexible body dyanmics fiber optic strain sensor p↗
DOE OSTI2021

On the Use of Smart Meter Data to Estimate the Voltage Magnitude on the Primary Side of Distribution Service Transformers: Preprint

This paper develops a novel method to estimate the voltage magnitude on the primary side of distribution service transformers. The proposed method relies exclusively on smart meters, and therefore it is fully data-driven. This is an important feature because electric utilities have detailed models of only the primary network--that is, the network between the distribution substation and the primary side of service transformers that are installed closer to end-customer sites. The network that connects the secondary side of service transformers to end-customer sites, referred to as the secondary network, is simply represented by a lumped load. For each secondary network, the proposed method uses data acquired from only 2 smart meters: the closest and the farthest--in the sense of electrical distance--from the service transformer. As a reference to this feature, the proposed method is named SM2Vp. To our knowledge, this is the first time a method is shown to provide actionable information for real-time operation and control of power distribution grids using only two smart meters per secondary network. This is important because utilities have experienced barriers in managing and using large data sets for real-time operation and control. SM2Vp is primarily intended to provide pseudo-measurements for distribution system state estimation, but it can also be used directly for voltage control schemes. The performance of SM2Vp is demonstrated by numerical simulations carried out on three secondary network synthetic models and by using field data provided by a utility partner serving customers in southwestern California. A maximum relative error of approximately 3.9% or less is observed for the primary voltage magnitude estimates in all numerical experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2020

Multibody Simulation of REMORA CubeSat Docking to and Pushing a Spent Rocket Booster

This paper details the multibody simulation of three phases of critical importance to the feasibility of the REMORA CubeSat space debris mitigation concept: the final approach of the CubeSat to a spent rocket booster; the grappling of the spent rocket booster using a robotic arm; and the pushing of the spent rocket booster to divert its course from another on-orbit asset. The extension of a robotic mobility and manipulation modeling toolkit (M3TK) from multibody dynamics simulation of manipulators and ground vehicles to simulation of orbital robotics is outlined. This includes the identification of the appropriate parameters required to concisely and generically describe thruster loads, thruster mixing, spacecraft control, and spacecraft navigation for the purpose of on-orbit robotics simulation. A high-level spacecraft navigator commands maneuvers to target spacecraft states. A PID spacecraft controller takes the target states and calculates desired forces and torques. A thruster mixer solves a quadratic program to determine the optimal thruster firing times for the propulsion system. Pulse width modulated actuation of eight canted cold gas thrusters is used in the simulated approach to a rocket nozzle from a distance of 200 m. A five degree of freedom robotic arm is controlled to position a pair of pincers to grasp the rocket nozzle. Contact dynamics are used to accurately simulate the grasping of the rocket nozzle by the pincers. A similar, but separate simulation is performed to assess the ability of the REMORA CubeSat to push the large spent rocket booster. This diversion maneuver makes use of an additional, larger thruster, and pushes the rocket booster in excess of 400 m. Appropriate motor control gains on the robotic arm are found to be higher during the pushing phase than those which are appropriate during free motion; this increase promotes rigidity of the arm and allows it to properly direct the pushing force. Challenges encountered in time step selection for numerical stability of the simulation are also discussed.

Mukherjee, Rudranarayan↗
DOE OSTI2020

Will LISA Detect Harmonic Gravitational Waves from Galactic Cosmic String Loops?

Rapid advancement in the observation of cosmic strings has been made in recent years placing increasingly stringent constraints on their properties, with $G\mu\lesssim 10^{-11}$ from Pulsar Timing Array (PTA). Cosmic string loops with low string tension clump in the Galaxy due to slow loop decay and low gravitational recoil, resulting in great enhancement to loop abundance in the Galaxy. With an average separation of down to just a fraction of a kpc, and the total power of gravitational wave (GW) emission dominated by harmonic modes spanning a wide angular scale, resolved loops located in proximity are powerful, persistent, and highly monochromatic sources of GW with a harmonic signature not replicated by any other sources, making them prime targets for direct detection by the upcoming Laser Interferometer Space Antenna (LISA), whose frequency range is well-matched. Unlike detection of bursts where the detection rate scales with loop abundance, the detection rate for harmonic signal is the result of a complex interplay between the strength of GW emission, loop abundance, and other sources of noise, and is most suitably studied through numerical simulations. We develop a robust and flexible framework for simulating loops in the Galaxy for predicting direct detection of harmonic signal from resolved loops by LISA. Our simulation reveals that the most accessible region in the parameter space consists of large loops $\alpha=0.1$ with low tension $10^{-21}\lesssim G\mu\lesssim 10^{-19}$. Direct detection of field theory cosmic strings is unlikely, with the detection probability $p_{\mathrm{det}}\lesssim 2\%$ for a 1-year mission. An extension suggests that direct detection of cosmic superstrings with a low intercommutation probability is very promising. Searching for harmonic GW signal from resolved loops through LISA observations will potentially place physical constraints on string theory.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Modeling and Experimental Demonstration of Flux Spreading in Light Trapping Planar-Cavity Solar Enclosed Particle Receivers

This study experimentally validates and numerically models the flux-spreading effect in a light-trapping planar-cavity solar receiver) for particle-based concentrating solar power systems. The receiver's shallow cavity with vertical planar walls redistributes concentrated solar flux, reducing peak intensity and achieving uniform heat flux. On-sun tests at National Renewable Energy Laboratory's High-Flux Solar Furnace under flux up to 1500 kW/m2 measured cavity wall temperatures, which were compared with Monte Carlo ray-tracing (SolTrace) and computational numerical simulations. Three angular absorptance models were evaluated: constant absorptance, a Pyromark-based directional model, and a Fresnel-based Cr2O3 model. The Fresnel-derived model showed the best agreement with experiments, achieving high correlation (PC > 0.85), structural similarity (SSIM > 0.98), and signal-to-noise ratios (PSNR > 40 dB), with temperature prediction errors of 1-11%. Results confirm that flux spreading mitigates local overheating and validate the integrated modeling approach, supporting the solar receiver scalability for high-efficiency, high-temperature concentrating solar power applications.

14 SOLAR ENERGY↗
NASA NTRS1990

Comparison of thermal and nonthermal hard X-ray emission in electron-heated solar flares

Direct nonthermal bremsstrahlung and electronbeam-energized hot plasma thermal bremsstrahlung calculations are presently conducted on the basis of numerical simulation results for solar atmospheric response to heating during solar flares by nonthermal electron beams. It is established that while emission is dominated by the thermal component below a certain 'cross-over' energy, the direct brehmsstrahlung component becomes more important at higher energies. The dependence of the cross-over energy on the parameters of the electron energy input is also investigated, and the minimum photon energy which should be observable with future instruments in order that the nonthermal component not be swamped by the self-consistent thermal counterpart created by beam heating.

Li, Peng↗
DOE OSTI2019

Deep Direct-Use Feasibility Study Numerical Modeling and Uncertainty Analysis using iTOUGH2 for West Virginia University

To reduce the geothermal exploration risk, a feasibility study is performed for a deep direct-use system proposed at the West Virginia University (WVU) Morgantown campus. This study applies numerical simulations to investigate reservoir impedance and thermal production. Because of the great depth of the geothermal reservoir, few data are available to characterize reservoir features and properties. As a result, the study focuses on the following three aspects: 1. model choice for predicting reservoir impedance and thermal breakthrough: after investigating three potential models (one single permeability model and two dual permeability models) for flow through fractured rock, it is decided to use single permeability model for further analysis; 2. well placement (horizontal vs. vertical) options: horizontal well placement seems to be more robust to heterogeneity and the impedance is more acceptable; 3. Prediction uncertainty: the most influential parameters are identified using a First-Order-Second-Moment uncertainty propagation analysis, and the uncertain range of the model predictions is estimated by performing a Monte Carlo simulation. Heterogeneity has a large impact on the perdition, therefore, is considered in the predictive model and uncertainty analysis. The numerical model results and uncertainty analysis are used for economic analysis. The dataset submitted here support the described study. Manuscript is submitted to Geothermics, will be linked once paper is accepted.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Performance estimation of resonance-enhanced dual-buoy wave energy converter using coupled time-domain simulation

This article presents the modeling methodology and performance evaluation of the resonance-enhanced dual-buoy WEC (Wave Energy Converter) by HEM (hydrodynamic & electro-magnetic) fully-coupled-dynamics time-domain-simulation program. The numerical results are systematically compared with the authors’ 1/6-scale experiment. With a direct-drive linear generator, the WEC consists of dual floating cylinders and a moon-pool between the cylinders, which can utilize three resonance phenomena from moon-pool dynamics as well as heave motions of inner and outer buoys. The contact and friction between the two buoys observed in the experiment are also properly modeled in the time-domain simulation by the Coulomb-friction model. Moon-pool resonance peaks significantly exaggerated in linear potential theory are empirically adjusted through comparisons with measured values. A systematic comparative study between the simulations and experiments with and without PTO (power-take-off) is conducted, and the relative heave displacements/velocities and power outputs are well matched. Then, parametric studies are carried out with the simulation program to determine optimum generator parameters. The performance with various wave conditions is also assessed.

16 TIDAL AND WAVE POWER↗
NASA NTRSDate not supplied

Computational Study of Noise Characteristics for 30P30N Leading-edge Slat and Krueger Flap

In this study, numerical simulations are performed to investigate the flowfield, the noise sources, and their detailed characteristics around high-lift configurations featuring either a conventional leading-edge slat or a Krueger flap. While the flowfields and potential noise sources are similar for the two types of configurations, the noise directivity patterns exhibit some differences. The influence of the cavity where the Krueger flap is stowed on the lower surface is also examined via numerical simulations with and without the cavity. Results show that the cavity does not have a significant influence on the flowfield around the Krueger flap and the gap region. However, when the cavity is present, additional noise sources are observed around its rear part, although they are not as significant as the noise sources on the lower side of the Krueger flap. These additional noise sources contribute to the directivity of high noise levels, especially in the upstream direction.

Aeroacoustics↗
DOE OSTI2021

Large Eddy Simulation of Lean Mixed-Mode Combustion Assisted by Partial Fuel Stratification in a Spark-Ignition Engine

Partial fuel stratification (PFS) is a promising fuel injection strategy to improve the stability of lean combustion by applying a small amount of pilot injection right before spark timing. Mixed-mode combustion, which makes use of end-gas autoignition following conventional deflagration-based combustion, can be further utilized to speed up the overall combustion. In this study, PFS-assisted mixed-mode combustion in a lean-burn direct injection spark-ignition (DISI) engine is numerically investigated using multi-cycle large eddy simulation (LES). A previously developed hybrid G-equation/well-stirred reactor combustion model for the well-mixed operation is extended to the PFS-assisted operation. The experimental spray morphology is employed to derive spray model parameters for the pilot injection. The LES-based model is validated against experimental data and is further compared with the Reynolds-averaged Navier–Stokes (RANS)-based model. Overall, both RANS and LES predict the mean pressure and heat release rate traces well, while LES outperforms RANS in capturing the cycle-to-cycle variation (CCV) and the combustion phasing in the mass burned space. Liquid and vapor penetrations obtained from the simulations agree reasonably well with the experiment. Detailed flame structures predicted from the simulations reveal the transition from a sooting diffusion flame to a lean premixed flame, which is consistent with experimental findings. LES captures more wrinkled and stretched flames than RANS. Finally, the LES model is employed to investigate the impacts of fuel properties, including heat of vaporization (HoV) and laminar burning speed (SL). Combustion phasing is found more sensitive to SL than to HoV, with a larger fuel property sensitivity of the heat release rate from autoignition than that from deflagration. Moreover, the combustion phasing in the PFS-assisted operation is shown to be less sensitive to SL compared with the well-mixed operation.

fuel property sensitivity↗
DOE OSTI2024

Steady-state properties of multi-orbital systems using quantum Monte Carlo

A precise dynamical characterization of quantum impurity models with multiple interacting orbitals is challenging. In quantum Monte Carlo methods, this is embodied by sign problems. A dynamical sign problem makes it exponentially difficult to simulate long times. A multi-orbital sign problem generally results in a prohibitive computational cost for systems with multiple impurity degrees of freedom even in static equilibrium calculations. Here, we present a numerically exact inchworm method that simultaneously alleviates both sign problems, enabling simulation of multi-orbital systems directly in the equilibrium or nonequilibrium steady-state. The method combines ideas from the recently developed steady-state inchworm Monte Carlo framework [Erpenbeck et al., Phys. Rev. Lett. 130, 186301 (2023)] with other ideas from the equilibrium multi-orbital inchworm algorithm [Eidelstein et al., Phys. Rev. Lett. 124, 206405 (2020)]. We verify our method by comparison with analytical limits and numerical results from previous methods.

Chemistry↗