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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,693 records · Page 94

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Multi-dimensional modeling of mixture preparation in a direct injection engine fueled with gaseous hydrogen

With the recent advances of direct injection (DI) technology, introducing hydrogen into the combustion chamber through DI is being considered as a viable approach to circumvent backfire and pre-ignition encountered in early generations of hydrogen engines. As part of a broader vision to develop a robust numerical model to study hydrogen spark ignition (SI) combustion in internal combustion (IC) engines, the present numerical investigation focuses on mixture preparation in a hydrogen DI SI engine. This study is carried out with a single hole injector with gaseous hydrogen injected at 100 bar injection pressure. Simulations are carried out for high and low tumble configurations and validated against optical data acquired from planar laser induced fluorescence (PLIF) measurements. Varying mesh configurations are investigated for the impact on in-cylinder mixture distribution. A particular emphasis is placed on the effect of nozzle geometry and mesh orientation near the wall. Overall, the computational model is found to predict the mixture distribution in the combustion cylinder reasonably well. The results showed that the alignment of mesh with the flow direction is important to achieve good agreement between numerical analysis and optical measurement data.

08 HYDROGEN↗

NASA NTRS2011

Direct Simulation of Multiple Scattering by Discrete Random Media Illuminated by Gaussian Beams

The conventional orientation-averaging procedure developed in the framework of the superposition T-matrix approach is generalized to include the case of illumination by a Gaussian beam (GB). The resulting computer code is parallelized and used to perform extensive numerically exact calculations of electromagnetic scattering by volumes of discrete random medium consisting of monodisperse spherical particles. The size parameters of the scattering volumes are 40, 50, and 60, while their packing density is fixed at 5%. We demonstrate that all scattering patterns observed in the far-field zone of a random multisphere target and their evolution with decreasing width of the incident GB can be interpreted in terms of idealized theoretical concepts such as forward-scattering interference, coherent backscattering (CB), and diffuse multiple scattering. It is shown that the increasing violation of electromagnetic reciprocity with decreasing GB width suppresses and eventually eradicates all observable manifestations of CB. This result supplements the previous demonstration of the effects of broken reciprocity in the case of magneto-optically active particles subjected to an external magnetic field.

Mackowski, Daniel W.↗

DOE OSTI2023

Creating large Fock states and massively squeezed states in optics using systems with nonlinear bound states in the continuum

The quantization of the electromagnetic field leads directly to the existence of quantum mechanical states, called Fock states, with an exact integer number of photons. Despite these fundamental states being long-understood, and despite their many potential applications, generating them is largely an open problem. For example, at optical frequencies, it is challenging to deterministically generate Fock states of order two and beyond. Here, we predict the existence of an effect in nonlinear optics, which enables the deterministic generation of large Fock states at arbitrary frequencies. The effect, which we call an n-photon bound state in the continuum, is one in which a photonic resonance (such as a cavity mode) becomes lossless when a precise number of photons n is inside the resonance. Based on analytical theory and numerical simulations, we show that these bound states enable a remarkable phenomenon in which a coherent state of light, when injected into a system supporting this bound state, can spontaneously evolve into a Fock state of a controllable photon number. This effect is also directly applicable for creating (highly) squeezed states of light, whose photon number fluctuations are (far) below the value expected from classical physics (i.e., shot noise). In conclusion, we suggest several examples of systems to experimentally realize the effects predicted here in nonlinear nanophotonic systems, showing examples of generating both optical Fock states with large n (n > 10), as well as more macroscopic photonic states with very large squeezing, with over 90% less noise (10 dB) than the classical value associated with shot noise.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1991

Quantitative knowledge acquisition for expert systems

A common problem in the design of expert systems is the definition of rules from data obtained in system operation or simulation. While it is relatively easy to collect data and to log the comments of human operators engaged in experiments, generalizing such information to a set of rules has not previously been a direct task. A statistical method is presented for generating rule bases from numerical data, motivated by an example based on aircraft navigation with multiple sensors. The specific objective is to design an expert system that selects a satisfactory suite of measurements from a dissimilar, redundant set, given an arbitrary navigation geometry and possible sensor failures. The systematic development is described of a Navigation Sensor Management (NSM) Expert System from Kalman Filter convariance data. The method invokes two statistical techniques: Analysis of Variance (ANOVA) and the ID3 Algorithm. The ANOVA technique indicates whether variations of problem parameters give statistically different covariance results, and the ID3 algorithms identifies the relationships between the problem parameters using probabilistic knowledge extracted from a simulation example set. Both are detailed.

Belkin, Brenda L.↗

DOE OSTI2023

Defect‐Free Sound Insulator Using Single Metal‐Based Friction Stir Process Array

Metals are excellent conductors for phonon transportation such as vibration, sound, and heat. Generally, metal sound insulators require multimaterial structure or defects and unimetal sound insulators are challenging. Therefore, a design of a defect-free sound insulator made by single alloys with multiple friction stir processes (FSPs) is proposed. Periodic friction stir processing can induce superlattice-like local mechanical properties’ modifications. By experimental acoustic characterization, it is observed that FSP can introduce clear acoustic–elastic property contrast on an aluminum plate by the presence of stir zone and heat-affected zones. In numerical simulations, the signature FSP-induced property profile is periodically and parallelly arranged on a long aluminum plate. The transmission gap frequencies are present on the frequency spectrum with the sound propagation direction perpendicular to the FSP paths. Disorder offsets on FSP periodicity are further introduced. Anderson localization is found on a resonance frequency, which provides −11 dB sound reduction by an exponential decay. Due to the finite design length, the slight disorder can also enhance sound insulation in the periodic transmission gap frequency. In conclusion, with analysis and comparison with different configurations, the best performance in the models can achieve −30 dB sound insulation in the 350 mm-long aluminum alloy plate with 14 parallel FSPs.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Multiscale numerical investigation of ratchet growth damage effects in PBX 9502

This paper presents results of numerical experiments conducted on the high explosive PBX 9502 to investigate how recently observed grain-scale damage mechanisms of ratchet growth affect uniaxial compression measurements. Simulations are multiscale in the sense of directly resolving grains, pores, cracks, and grain-interfaces based upon scanning electron microscope (SEM) images of damaged and undamaged samples. The combined finite-discrete element method (FDEM) is utilized to resolve both grain-scale microfracture and elastoplastic deformation of solid grains. Pristine (undamaged) and damaged microstructures are compared in simulation of unconfined compression tests of the same material from the literature. Here, the simulation results show the observed microscale mechanisms of damage, specifically microfracture predominantly around and sometimes through grains and crack-associated pore growth, can well-explain the effective degradation of strength and stiffness observed in the laboratory measurements.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Optimum power analysis of a self-reactive wave energy point absorber with mechanically-driven power take-offs

This paper studies a self-reactive ocean wave energy converter (WEC) where energy is absorbed through the relative motion between a floating buoy on the ocean surface and a submerged body. Two types of direct-drive power take off (PTO) systems are examined. One adopts a ball screw system in inverse driving mode with mechanical motion rectifier (MMR) that converts the oscillating bi-directional input heave motion into unidirectional rotation to drive an electromagnetic generator, and another one uses a ball screw system to directly drive the generator (called non-MMR). Dynamic models for both types of PTOs are established and integrated with the overall WEC system model in both the time and frequency domains. Numerical simulation is used to investigate the dynamic performance of both PTOs. The influences of the PTO inerter, which is from the rotational inertia of the flywheel and generator, and the WEC drag damping coefficient are discussed and explored. The analytical closed-form solution for the optimum condition of the non-MMR system is derived based on the criterion of maximum power absorption, and numerical simulation yields the optimum condition for the MMR system. Overall, the results under regular and irregular waves show that the PTO inerter in the MMR system can improve the power absorption for small wave periods and maintain the same performance for large wave periods. The PTO inerter in the non-MMR system shifts the peak frequency. The WEC drag damping negatively influences the power absorption in both PTO systems, especially at large wave periods.

16 TIDAL AND WAVE POWER↗

DOE OSTI2023

Staggered scheme for the compressible fluctuating hydrodynamics of multispecies fluid mixtures

Here, we present a numerical formulation for the solution of nonisothermal, compressible Navier-Stokes equations with thermal fluctuations to describe mesoscale transport phenomena in multispecies fluid mixtures. The novelty of our numerical method is the use of staggered grid momenta along with a finite volume discretization of the thermodynamic variables to solve the resulting stochastic partial differential equations. The key advantages of the numerical scheme are that it significantly simplifies the discretization of diffusive and stochastic momentum fluxes into a more compact form, and it provides an unambiguous prescription of boundary conditions involving pressure. The staggered grid scheme more accurately reproduces the equilibrium static structure factor of hydrodynamic fluctuations in gas mixtures compared to a collocated scheme described previously by Balakrishnan et al. [Phys. Rev. E 89, 013017 (2014)1539-375510.1103/PhysRevE.89.013017]. The numerical method is tested for ideal noble gases mixtures under various nonequilibrium conditions, such as applied thermal and concentration gradients, to assess the role of cross-diffusion effects, such as Soret and Dufour, on the long-ranged correlations of hydrodynamic fluctuations, which are also more accurately reproduced compared to the collocated scheme. We numerically study giant nonequilibrium fluctuations driven by concentration gradients and fluctuation-driven Rayleigh-Taylor instability in gas mixtures. Wherever applicable, excellent agreement is observed with theory and measurements from the direct simulation Monte Carlo method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Electromagnetic Signatures from the Tidal Tail of a Black Hole—Neutron Star Merger

Black hole - neutron star (BH-NS) mergers are a major target for ground-based gravitational wave observatories. A merger can also produce an electromagnetic counterpart (a kilonova) if it ejects neutron-rich matter that assembles into heavy elements through r-process nucleosynthesis. We study the kilonova signatures of the unbound dynamical ejecta of a BH-NS merger. We take as our initial state the results from a numerical relativity simulation and then use a general relativistic hydrodynamics code to study the evolution of the ejecta with parameterized r-process heating models. The unbound dynamical ejecta is initially a flattened, directed tidal tail largely confined to a plane. Heating from the r-process inflates the ejecta into a more spherical shape and smooths its small-scale structure, though the ejecta retains its bulk directed motion. We calculate the electromagnetic signatures using a 3D radiative transfer code and a parameterized opacity model for lanthanide-rich matter. The light curve varies with viewing angle because of two effects: asphericity results in brighter emission for orientations with larger projected areas, while Doppler boosting results in brighter emission for viewing angles more aligned with the direction of bulk motion. For typical r-process heating rates, the peak bolometric luminosity varies by a factor of ~3 with orientation while the peak in the optical bands varies by ~3 magnitudes. The spectrum is blueshifted at viewing angles along the bulk motion, which increases the V-band peak magnitude to ~-14 despite the lanthanide-rich composition.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1991

Numerical simulation of icing, deicing, and shedding

An algorithm has been developed to numerically model the concurrent phenomena of two-dimensional transient heat transfer, ice accretion, ice shedding and ice trajectory which arise from the use of electrothermal pad. The Alternating Direction Implicit method is used to simultaneously solve the heat transfer and accretion equations occurring in the multilayered body covered with ice. In order to model the phase change between ice and water, a technique was used which assumes a phase for each node. This allows the equations to be linearized such that a direct solution is possible. This technique requires an iterative procedure to find the correct phase at each node. The computer program developed to find this solution has been integrated with the NASA-Lewis flow/trajectory code LEWICE.

Wright, W. B.↗

NASA NTRS1991

Charge collected by diffusion from an ion track under mixed boundary conditions

Charge-carrier diffusion from an ion track in a silicon substrate at least a few hundred microns thick is analyzed. The substrate upper surface is treated as reflective except for a small section, intended to represent a reverse-biased junction, which is treated as a sink. Total charge collected by the sink is calculated by assuming transport to be governed by an ambipolar diffusion equation with temporally constant and spatially uniform carrier lifetime and diffusion coefficient. Present results apply to a normally incident track but could easily be generalized to arbitrary track direction. The collected charge is found to depend on track length and on the electrostatic capacitance, rather than the area, of the sink. Theoretical predictions are compared to the results of a numerical simulation called the Poisson and Continuity Equation Solver (PISCES) for three cases and are found to agree within a factor of two in the worst case.

Edmonds, Larry D.↗

NASA NTRS2000

High Resolution Experiments of Compressible Turbulence Interacting With Shock Waves

Most of the second year of our research program focused on exploring the potentially favorable the effects of expansion waves on homogeneous and isotropic turbulence, which is formed downstream of a grid. Expansion waves are associated with compressible flows and may reduce the drag over airfoils by suppressing turbulence. In the very few previous investigations of interactions of turbulence with expansion waves the effects due to stabilizing streamline curvature substantially masked the effects of turbulence suppression due to flow expansion though the waves. In the present flow configuration planar expansion waves interact with grid generated turbulence in our high-resolution shock tube research facility. This approach will assess directly the effects of the interaction on turbulence. The first objective of our study was to identify the nature of expansion waves present in our shock tube facility. Our time-dependent numerical simulations of the flow in our facility indicated the existence of two regions of traveling expansion waves. The system of expansion waves utilized in this investigation is generated by the exiting shock wave and the induced flow behind it at the end of the driver. Several new measuring techniques are being developed which are capable of providing velocity-gradient-related quantities in compressible flows for the first time.

Agui, Juan↗

NASA NTRS2004

HiPEP Ion Optics System Evaluation Using Gridlets

Experimental measurements are presented for sub-scale ion optics systems comprised of 7 and 19 aperture pairs with geometrical features that are similar to the HiPEP ion optics system. Effects of hole diameter and grid-to-grid spacing are presented as functions of applied voltage and beamlet current. Recommendations are made for the beamlet current range where the ion optics system can be safely operated without experiencing direct impingement of high energy ions on the accelerator grid surface. Measurements are also presented of the accelerator grid voltage where beam plasma electrons backstream through the ion optics system. Results of numerical simulations obtained with the ffx code are compared to both the impingement limit and backstreaming measurements. An emphasis is placed on identifying differences between measurements and simulation predictions to highlight areas where more research is needed. Relatively large effects are observed in simulations when the discharge chamber plasma properties and ion optics geometry are varied. Parameters investigated using simulations include the applied voltages, grid spacing, hole-to-hole spacing, doubles-to-singles ratio, plasma potential, and electron temperature; and estimates are provided for the sensitivity of impingement limits on these parameters.

Willliams, John D.↗

NASA NTRS2020

Cascade Model for Planetesimal Formation by Turbulent Clustering

We use a newly developed cascade model of turbulent concentration of particles in protoplanetary nebulae to calculate several properties of interest to the formation of primitive planetesimals and to the meteorite record. The model follows, and corrects, calculations of the primary planetesimal Initial Mass Function (IMF) by Cuzzi et al. (2010), in which an incorrect cascade model was used. Here we use the model of Hartlep et al. (2017), which has been validated against several published numerical simulations of particle concentration in turbulence. We find that, for a range of nebula and particle properties, planetesimals may be “born big," formed as sandpiles with diameters in the 10 - 100km range, directly from freely floating particles. The IMFs have a modal nature, with a well-defined peak rather than a powerlaw size dependence. Predictions for the inner and outer nebula behave similarly in these regards, and observations of inner and outer nebula primitive bodies support such modal IMFs. Also, we present predictions of local particle concentrations on several lengthscales in which particles “commonly" find themselves, which have significance for meteoritical observations of the redox state and isotopic fractionation in regions of chondrule formation. An important difference between these results, and those of Cuzzi et al. (2010), is that particle growth-by-sticking must proceed to at least the 1-few cm radius range for the IMF and meteoritical properties to be most plausibly satisfied. That is, as far as the inner nebula goes, the predominant “particles" must be aggregates of chondrules (or chondrule-size precursors) rather than individual chondrules themselves.

Planetesimals↗

DOE OSTI2023

Large deformation and brittle failure calculated using the dual-domain material point method

The dual domain material point (DDMP) method is explored as a candidate to be implemented in a general purpose code to perform simulations of materials with complex geometry that undergo large history-dependent deformation and failure. To test its candidacy, we study its mesh convergence, its sensitivity to mesh orientation, and its ability to handle softening and failure of a material. Simulations of large deformation and simulations of mechanical failure are performed using both DDMP and the material point method (MPM). When cell-crossing of material points is not an issue and when there are a sufficient number of material points in each computation cell, the numerical error decreases with the square of the cell size as expected for both MPM and DDMP. DDMP has reduced error compared with MPM when there are many instances of material points crossing cell boundaries due to the continuous nature of the modified gradient of the shape functions. Simulations of a specimen under tension are also performed where the background mesh is aligned and misaligned with the tension direction. MPM displays a significant mesh-dependent stress field, DDMP shows negligible mesh dependency. Despite a mesh orientation-dependent stress field from MPM, the critical tension and failure mode from both MPM and DDMP calculations have negligible mesh dependency when using a non-local failure model. If only the failure mode is important (i.e., local stresses are unimportant), MPM with a non-local failure model is a suitable method for modeling failure with small deformations. However, if local stresses are also important or if there are large deformations with many cell-crossings before failure, DDMP should be the method that is used. A needed improvement for DDMP is identified from our numerical simulations.

36 MATERIALS SCIENCE↗

DOE OSTI2020

The TRANSCEND University Consortium: Theme 2 - Site Decommissioning, Deconstruction and Remediation - 20329

This paper will discuss the aims, objectives and progress to-date under Theme 2 of the TRANSCEND consortium project: Site decommissioning, deconstruction and remediation. Decommissioning nuclear sites involves waste retrieval, decontamination, deconstruction and, where necessary, containment and/or remediation of the remaining structure and surrounding land. Critical to management of these processes is limiting radiation exposure for the workforce, restricting the spread of radionuclides in groundwater, surface water and airborne particulates, and minimising the volume of contaminated waste for disposal. The aim of Theme 2 research is to develop new technologies for monitoring, remediation and containment that serve to minimise the volume of radioactively contaminated waste for disposal, for application prior to, during and after retrieval, deconstruction and decontamination operations. Prior research, conducted under the previous DISTINCTIVE project, demonstrated that colloidal silica grout can penetrate low permeability materials (including cement) for hydraulic barrier formation, and improved sorption capacity. The silica grout can be injected at surface using extremely small (potentially gravity-driven) fluid pressures, without the need for borehole drilling. Current work is investigating the erodibility of silica-grouted soils for inhibition of airborne and water-borne particulates as well as enhancing the grout's sorption capacity by addition of other materials to provide a chemical, as well as hydraulic, barrier to subsurface migration. EK remediation uses low voltage DC current to control migration of contaminants in porous media as well as to remove or degrade them. Researchers in the consortium have already demonstrated that low-energy ex-situ EK techniques can be used to provide remediation and volume minimisation for AWE legacy wastes in the UK. Under the current project researchers are building EK test cells containing simulated site materials at laboratory and intermediate(m)-scales to: remove, focus or degrade contaminants (remediation or waste minimisation); and direct subsurface water, chemical and colloid flow (fencing/containment or forced migration). This lab-based research is being informed by numerical models of EK processes that can subsequently be used to design full-scale on-site applications by nuclear site holders. Research under theme 2 of TRANSCEND will ultimately combine novel non-invasive detection technologies with EK techniques and colloidal silica grout barriers, to optimise the containment of radionuclide contamination in soils. Thus, allowing us to detect in-situ contamination, mobilise it to a selected location and grout it in-situ, prior to the initiation of decommissioning and deconstruction operations. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2025

Numerical demonstration of Biermann-generated magnetic fields in a laser-driven cylindrical implosion platform

We present two-dimensional radiation-magnetohydrodynamics FLASH simulations studying Biermann battery generated magnetic fields in a convergent geometry based on modifying an existing laser direct drive cylindrical implosion platform at Los Alamos National Laboratory. During the laser driven implosion, instabilities can arise as the target converges, and in a plasma, such instabilities can generate magnetic fields via the Biermann battery. After significant redesign (lowering the inner foam fill density, deceasing the ablator mass, and modifying the laser pulse), we demonstrate that magnetic fields reach strengths of many hundred kilogauss and persist over a nanosecond time scale. Such results pave the way for an experimental campaign for further code verification and validation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2003

The Formation of Asteroid Satellites in Catastrophic Impacts: Results from Numerical Simulations

We have performed new simulations of the formation of asteroid satellites by collisions, using a combination of hydrodynamical and gravitational dynamical codes. This initial work shows that both small satellites and ejected, co-orbiting pairs are produced most favorably by moderate-energy collisions at more direct, rather than oblique, impact angles. Simulations so far seem to be able to produce systems qualitatively similar to known binaries. Asteroid satellites provide vital clues that can help us understand the physics of hypervelocity impacts, the dominant geologic process affecting large main belt asteroids. Moreover, models of satellite formation may provide constraints on the internal structures of asteroids beyond those possible from observations of satellite orbital properties alone. It is probable that most observed main-belt asteroid satellites are by-products of cratering and/or catastrophic disruption events. Several possible formation mechanisms related to collisions have been identified: (i) mutual capture following catastrophic disruption, (ii) rotational fission due to glancing impact and spin-up, and (iii) re-accretion in orbit of ejecta from large, non-catastrophic impacts. Here we present results from a systematic investigation directed toward mapping out the parameter space of the first and third of these three collisional mechanisms.

Durda, D. D.↗