Search NASA⌕ Search

SEARCH · Search NASA

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,765 records · Page 98

On the physical size of the Milky Way globular cluster NGC 7089 (M2)

ABSTRACT We study the outer regions of the Milky Way globular cluster NGC 7089 based on new Dark Energy Camera observations. The resulting background-cleaned stellar density profile reveals the existence of an extended envelope. We confirm previous results that cluster stars are found out up to ∼1° from the cluster’s centre, which is nearly three times the value of the most robust tidal radii estimations. We also used results from direct N-body simulations in order to compare with the observations. We found a fairly good agreement between the observed and numerically generated stellar density profiles. Because of the existence of gaps and substructures along globular cluster tidal tails, we closely examined the structure of the outer cluster region beyond the Jacobi radius. We extended the analysis to a sample of 35 globular clusters, 20 of them with observed tidal tails. We found that if the stellar density profile follows a power law ∝ r−α, the α slope correlates with the globular cluster present mass, in the sense that, the more massive the globular cluster, the smaller the α value. This trend is not found in globular clusters without observed tidal tails. The origin of such a phenomenon could be related, among other reasons, to the proposed so-called potential escapers or to the formation of globular clusters within dark matter minihaloes.

79 ASTRONOMY AND ASTROPHYSICS↗

Physics-constrained, low-dimensional models for magnetohydrodynamics: First-principles and data-driven approaches

Plasmas are highly nonlinear and multiscale, motivating a hierarchy of models to understand and describe their behavior. However, there is a scarcity of plasma models of lower fidelity than magnetohydrodynamics (MHD), although these reduced models hold promise for understanding key physical mechanisms, efficient computation, and real-time optimization and control. Galerkin models, obtained by projection of the MHD equations onto a truncated modal basis, and data-driven models, obtained by modern machine learning and system identification, can furnish this gap in the lower levels of the model hierarchy. This work develops a reduced-order modeling framework for compressible plasmas, leveraging decades of progress in projection-based and data-driven modeling of fluids. We begin by formalizing projection-based model reduction for nonlinear MHD systems. To avoid separate modal decompositions for the magnetic, velocity, and pressure fields, we introduce an energy inner product to synthesize all of the fields into a dimensionally consistent, reduced-order basis. Next, we obtain an analytic model by Galerkin projection of the Hall-MHD equations onto these modes. We illustrate how global conservation laws constrain the model parameters, revealing symmetries that can be enforced in data-driven models, directly connecting these models to the underlying physics. We demonstrate the effectiveness of this approach on data from high-fidelity numerical simulations of a three-dimensional spheromak experiment. Finally, this manuscript builds a bridge to the extensive Galerkin literature in fluid mechanics and facilitates future principled development of projection-based and data-driven models for plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Enabling Predictive Scale-Bridging Simulations through Active Learning (Institutional Computing Annual Report (Project w20_alscalebridging)) [Slides]

The goal of this project was to develop, demonstrate, and provide a new capability to achieve greater physical fidelity in large-scale simulations, rather than the usual brute-force increases in the number of mesh elements or particles. This was done by using machine learning (ML) techniques to develop emulators for subscale physics that can be used in coarse-scale continuum simulations, trained on fine-scale molecular dynamics (MD) simulations that are launched on-the-fly via active learning. raditional inertial confinement fusion (ICF) simulations rely on numerical diffusion to simulate molecular effects such as non-local transport and mixing without truly accounting for molecular interactions; our approach directly accounts for this physics.

74 ATOMIC AND MOLECULAR PHYSICS↗

Numerical simulation of the flow field and fuel sprays in an IC engine

A two-dimensional model for axisymmetric piston-cylinder configurations is developed to study the flow field in two-stroke direct-injection Diesel engines under motored conditions. The model accounts for turbulence by a two-equation model for the turbulence kinetic energy and its rate of dissipation. A discrete droplet model is used to simulate the fuel spray, and the effects of the gas phase turbulence on the droplets is considered. It is shown that a fluctuating velocity can be added to the mean droplet velocity every time step if the step is small enough. Good agreement with experimental data is found for a range of ambient pressures in Diesel engine-type microenvironments. The effects of the intake swirl angle in the spray penetration, vaporization, and mixing in a uniflow-scavenged two-stroke Diesel engine are analyzed. It is found that the swirl increases the gas phase turbulence levels and the rates of vaporization.

Nguyen, H. L.↗

Transonic Symposium: Theory, Application, and Experiment, volume 1, part 2

In order to assess the state of the art in transonic flow disciplines and to glimpse at future directions, NASA-Langley held a Transonic Symposium. Emphasis was placed on steady, three dimensional external, transonic flow and its simulation, both numerically and experimentally. The symposium included technical sessions on wind tunnel and flight experiments; computational fluid dynamic applications; inviscid methods and grid generation; viscous methods and boundary layer stability; and wind tunnel techniques and wall interference. This, being volume 1, is unclassified.

Foughner, Jerome T., Jr.↗

Lateral charge transport from heavy-ion tracks in integrated circuit chips

A 256K DRAM has been used to study the lateral transport of charge (electron-hole pairs) induced by direct ionization from heavy-ion tracks in an IC. The qualitative charge transport has been simulated using a two-dimensional numerical code in cylindrical coordinates. The experimental bit-map data clearly show the manifestation of lateral charge transport in the creation of adjacent multiple-bit errors from a single heavy-ion track. The heavy-ion data further demonstrate the occurrence of multiple-bit errors from single ion tracks with sufficient stopping power. The qualitative numerical simulation results suggest that electric-field-funnel-aided (drift) collection accounts for single error generated by an ion passing through a charge-collecting junction, while multiple errors from a single ion track are due to lateral diffusion of ion-generated charge.

Zoutendyk, J. A.↗

Numerical experiments in supersonic boundary-layer stability

The three-dimensional (3-D) time-dependent compressible Navier-Stokes equations are numerically solved by a Fourier-Chebyshev collocation method to study the stability of supersonic flows over a flat plate. Several direct simulations carried out in this study suggest the existence of a secondary instability that might provide a route to transition. The interaction of the modes involved in the secondary instability is possibly amenable to a Floquet-type analysis. Pertinent differences between this instability and the analogous incompressible K-type secondary instability are pointed out. Some preliminary results of a 2-D direct simulation of the nonlinear evolution of a second mode perturbation wave are also discussed.

Erlebacher, Gordon↗

Comet Shoemaker-Levy 9, Jupiter, and Impact Shock Chemistry

Four years ago this month, a hitherto unknown comet in loose orbit around Jupiter passed so near the giant planet that it was torn apart into 20 fragments by tides. One orbit later, two years ago this month, the fragments of doomed comet P/Shoemaker Levy (SL9) fell into Jupiter. The enormous energies of these impacts (the largest fragments were nearly 1 km across and, hitting at 60 km/s, released some 2-4 x 10(exp 27) ergs) produced enormous explosions. Several of the ejecta plumes were imaged towering 3000 km above Jupiter's limb. The heat released when the plumes fell was considerable and easily observed on Earth. The impacts produced strong shocks, both promptly at the impact site and again, later, and over thousands of kilometers, when the ejecta plume reentered the atmosphere. The focus of this talk will be to discuss what the SL9 impacts taught us about impact shock chemistry - the processes, the ingredients, the results - and what inferences we may draw for impacts on early Earth. Shock chemistry generates a suite of molecules not usually seen on Jupiter. The most surprising report was of a huge amount of diatomic sulfur S2 at the site of the G impact. Other unusual products include CS, CS2, OCS, H2S, SO2, HCN, CO, and H2O; although H2S and H2O are doubtless abundant below the visible clouds. Hot or enhanced CH4 and NH3 were also detected. A general rule of shock chemistry is that CO forms until either C or O is exhausted. If O greater than C, the other products are oxidized, and excess O goes to H2O. If C greater than O, the other products are reduced, and excess C goes to HCN, C2H2, and a wide variety of more complicated organics. Ultimately, given time, the carbon would react all the way to graphite, but in practice the reactions are incomplete. The dark ejecta debris were probably composed in part of carbonaceous particles generated by the shocks. In a sense, the SL9 impacts performed the famous Miller-Urey experiment on a grand scale, with one result being the production of a lot of complex brown organic solids (called "tholins"). We use, a straightforward chemical kinetics model for the H, N, C, O, S system to follow the nonequilibrium chemistry behind the shocks. The model traces the evolving chemical composition of a parcel of gas by directly integrating the web of chemical reactions. Pressure and temperature histories of the parcels are patterned after those calculated by numerical hydrodynamic simulations of the ejecta plume. A given plume parcel is generally shocked twice; t.e a parcel shocked near the impact site is ejected at high velocity and is shocked again when it reenters the atmosphere. The final state of the gas depends mostly on the second shock, provided that the latter is hot enough. The chemical evidence is ambiguous, but most indications are that C greater than O in the shocked, reacting gas. Telltale signatures of abundant oxygen - SO2, SO, CO2, O2 - were not seen, while signatures of abundant carbon - CS, CS2, and HCN - were. On the other hand, abundant H2O would appear to require O greater than C, and two other observed sulfur species, S2 and OCS, appear to form more easily in a somewhat oxidized gas, presumable vaporized from the comet itself. Since on general principles one expects the -comet to have had a more-or-less cosmic composition, i.e. O greater than C, the production of CS, CS2, and HCN probably requires C greater than O in the shocked jovian air. This in turn implies that even the largest fragments released the bulk of their energy above the jovian water table, in all likelihood above 5 bars . There is no evidence in favor of the proposition that a significant amount of wet jovian air was shocked strongly enough to coax water to react; i.e. wet jovian air saw only temperatures significantly below 2000 K.

Zahnle, Kevin↗

Hypersonic Shock Interactions About a 25 deg/65 deg Sharp Double Cone

This paper presents the results of a numerical study of shock interactions resulting from Mach 10 air flow about a sharp double cone. Computations are made with the direct simulation Monte Carlo (DSMC) method by using two different codes: the G2 code of Bird and the DAC (DSMC Analysis Code) code of LeBeau. The flow conditions are the pretest nominal free-stream conditions specified for the ONERA R5Ch low-density wind tunnel. The focus is on the sensitivity of the interactions to grid resolution while providing information concerning the flow structure and surface results for the extent of separation, heating, pressure, and skin friction.

Moss, James N.↗

Adaptive Control of Small Outboard-Powered Boats for Survey Applications

Four autopilot controllers have been developed in this work that can both hold a desired heading and follow a straight line. These PID, adaptive PID, neuro-adaptive, and adaptive augmenting control algorithms have all been implemented into a numerical simulation of a 33-foot center console vessel with wind, waves, and current disturbances acting in the perpendicular (across-track) direction of the boat s desired trajectory. Each controller is tested for its ability to follow a desired heading in the presence of these disturbances and then to follow a straight line at two different throttle settings for the same disturbances. These controllers were tuned for an input thrust of 2000 N and all four controllers showed good performance with none of the controllers significantly outperforming the others when holding a constant heading and following a straight line at this engine thrust. Each controller was then tested for a reduced engine thrust of 1200 N per engine where each of the three adaptive controllers reduced heading error and across-track error by approximately 50% after a 300 second tuning period when compared to the fixed gain PID, showing that significant robustness to changes in throttle setting was gained by using an adaptive algorithm.

VanZwieten, T.S.↗

Modal Analysis and Correlation of International Space Station Early Configurations

This paper will summarize the modal analysis and model refinement results for the International Space Station (ISS) Flight-2A mated configuration. To maintain structural integrity of the ISS, structural loading distributions have been rigorously analyzed through numerical simulation and taken into account during the design of the structure and mission operations. The accuracy of the analysis results is directly affected by the precision of mathematical models and estimated input forces. A mathematical model of an ISS configuration is composed of individual component math models. Each component model is required to be correlated with ground test data.

Kaouk, Mohamed↗

Preferrential Concentration of Particles in Protoplanetary Nebula Turbulence

Preferential concentration in turbulence is a process that causes inertial particles to cluster in regions of high strain (in-between high vorticity regions), with specifics depending on their stopping time or Stokes number. This process is thought to be of importance in various problems including cloud droplet formation and aerosol transport in the atmosphere, sprays, and also in the formation of asteroids and comets in protoplanetary nebulae. In protoplanetary nebulae, the initial accretion of primitive bodies from freely-floating particles remains a problematic subject. Traditional growth-by-sticking models encounter a formidable "meter-size barrier" [1] in turbulent nebulae. One scenario that can lead directly from independent nebula particulates to large objects, avoiding the problematic m-km size range, involves formation of dense clumps of aerodynamically selected, typically mm-size particles in protoplanetary turbulence. There is evidence that at least the ordinary chondrite parent bodies were initially composed entirely of a homogeneous mix of such particles generally known as "chondrules" [2]. Thus, while it is arcane, turbulent preferential concentration acting directly on chondrule size particles are worthy of deeper study. Here, we present the statistical determination of particle multiplier distributions from numerical simulations of particle-laden isotopic turbulence, and a cascade model for modeling turbulent concentration at lengthscales and Reynolds numbers not accessible by numerical simulations. We find that the multiplier distributions are scale dependent at the very largest scales but have scale-invariant properties under a particular variable normalization at smaller scales.

protoplanetary nebulae↗

Using Co-Optimized Machine Learned Manifolds for Modeling Chemically Reacting Flows

Chemically reacting flows play a key role in a wide range of engineered systems, from chemical and polymer processing to combustion-based energy conversion technologies. Simulations of these flows involve solving a coupled set of partial differential equations for mass, momentum, energy, and all relevant chemical species in the system. Chemical reaction pathways may be extremely complex and involve hundreds or more intermediate species, with reactions that occur over timescales varying by several orders of magnitude - presenting a significant numerical stiffness challenge. The combination of these factors makes simulation of chemically reacting flows vastly more expensive than nonreactive simulations, and often makes direct solution of the governing equations intractable. It is necessary to apply lower-fidelity models in place of the detailed governing equations in order to reduce computational cost to enable reacting flow simulation tools to be used in the engineering design process. Many of the models employed for this purpose are based on reducing the dimension of the thermochemical state, motivated by the observation that the observed thermochemical states in a system lie on a low-dimensional manifold in thermochemical state space. This behavior occurs due to the fast equilibration of certain reactive and transport processes, and physics-based manifold models rely on idealized assumptions about the balance of timescales and the way in which chemistry and transport are coupled. In this work, we apply a novel method for data-driven manifold-based modeling that can leverage data from high-fidelity reacting flow simulations to improve model accuracy in cases where the physics-based modeling assumptions break down. The approach is designed to be broadly applicable across chemically reacting flow systems but is applied here to turbulent combustion modeling.

machine learning↗

Impact of the Earthquake Rupture on Ground-Motion Variability of the 24 August 2016 M w 6.2 Amatrice, Italy, Earthquake

Here, the devastating 24 August 2016 M w 6.2 earthquake that struck Amatrice, Italy, marked the beginning of a prolonged seismic sequence dominated by three subsequent M w ≥6.0 events in the central Apennines region. The earthquake destroyed Amatrice’s historic center, claiming the lives of 299 individuals and causing widespread damage in the neighboring villages. The severity of the ground shaking, with a recorded maximum acceleration of 850 cm/s 2 on the east–west component at the Amatrice station, was far greater than the predicted acceleration based on the Italian ground-motion model (GMM). As pointed out by several investigations, the observed ground-motion amplitude and its spatial variability during the earthquake can be linked to specific rupture characteristics, including slip distribution and rupture directivity effects revealed by the observed data (Tinti et al., 2016; Pischiutta et al., 2021). In this study, we conducted physics-based 3D numerical simulations of ground motion for the Amatrice earthquake for frequencies up to 3 Hz. We employed a series of kinematic rupture models and a well-constrained local 3D velocity model incorporating surface topography. The kinematic rupture realizations were generated using multiscale hybrid and fully stochastic models, following the technique proposed by Graves and Pitarka (2016). We focused on assessing the sensitivity of near-fault ground-motion amplitudes to earthquake rupture characteristics, in particular, the spatial slip pattern. To evaluate the quality of our simulations, we employed goodness-of-fit measurements performed in comparisons of simulated and recorded ground motions. The simulated ground motions compare well with the recorded data and predictions from GMMs for Italy, ITA18 (Lanzano et al., 2019). However, we found that the simulated interevent ground-motion variability (randomness in the source process) of peak ground velocity, σ (PGV) is higher than the constant σ (PGV) predicted by conventional GMMs. Our simulations using several rupture scenarios demonstrate that the near-fault ground-motion amplification pattern is directly related to the slip distribution pattern.

3D Ground Motions Simulations↗

Light hydrocarbon and noble gas migration as an analogue for potential CO 2 leakage: numerical simulations and field data from three hydrocarbon systems

Abstract Sites where anthropogenic CO 2 captured from industrial sources is stored in deep geological formations for climate change mitigation are required to show secure retention of the injected CO 2 . Monitoring, reporting, and verification (MRV) plans are needed to indicate that no CO 2 release has occurred. We explored the degree to which direct comparison between a surface anomaly and reservoir geochemistry using various geochemical parameters can be used for attribution. We used data collected on light hydrocarbons and noble gases throughout the sedimentary column at three CO 2 enhanced oil recovery (EOR) sites to understand the processes that may cause fluid evolution. Light hydrocarbon and noble gases were sampled from reservoirs, gas‐bearing intervals above reservoirs, and groundwater. Vertical profiles indicated that lighter components are relatively enriched during migration (i.e. chromatographic separation). Static and numerical models were designed to simulate episodic gas migration and geochemical alteration of these geochemical parameters from solubility and sorption. The effects of hydrocarbon solubility were minimal (Bernard ratio changes within 5.2%) although field data were within the range of expected alterations from sorption. Forward models of CO 2 migration and noble gas interactions showed that CO 2 stripping causes an enrichment of crustal noble gases. In areas where natural fluxes of CO 2 from depth are non‐existent, the occurrence of crustal noble gas signature may distinguish fugitive CO 2 from the reservoir from natural near‐surface sources, and could be considered to explain apparent fluid anomalies. © 2019 Society of Chemical Industry and John Wiley & Sons, Ltd.

Anderson, Jacob↗

Numerical study of core formation of asymmetrically driven cone-guided targets

Compression of a directly driven fast ignition cone-sphere target with a finite number of laser beams is numerically studied using a three-dimensional hydrodynamics code IMPACT-3D. The formation of a dense plasma core is simulated for 12-, 9-, 6-, and 4-beam configurations of the GEKKO XII laser. The complex 3D shapes of the cores are analyzed by elucidating synthetic 2D x-ray radiographic images in two orthogonal directions. Finally, the simulated x-ray images show significant differences in the core shape between the two viewing directions and rotation of the stagnating core axis in the top view for the axisymmetric 9- and 6-beam configurations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Supersonic Jet Noise Reductions Predicted with Increased Jet Spreading Rate

In this paper, predictions are made of noise radiation from single, supersonic, axisymmetric jets. We examine the effects of changes in operating conditions and the effects of simulated enhanced mixing that would increase the spreading rate of the jet shear layer on radiated noise levels. The radiated noise in the downstream direction is dominated by mixing noise and it is well described by the instability wave noise radiation analysis. A numerical prediction scheme is used for the mean flow providing an efficient method to obtain the mean flow development for various operating conditions and to simulate the enhanced mixing. Using far field radiated noise measurements as a reference, the calculations predict that enhanced jet spreading results in a reduction of radiated noise.

Dahl, Milo D.↗