Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical diffusion”

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,405 records · Page 78

A Model of Charge-Transfer Excitons: Diffusion, Spin Dynamics, and Magnetic Field Effects

In this Letter, we investigate how the microscopic dynamics of charge-transfer (CT) excitons are influenced by the presence of an external magnetic field in disordered molecular semiconductors. This influence is driven by the dynamic interplay between the spin and spatial degrees of freedom of the electron–hole pair. To account for this interplay, we have developed a numerical framework that combines a traditional model of quantum spin dynamics with a stochastic coarse-grained model of charge transport. This combination gives a general and efficient methodology for simulating the effects of magnetic field on CT state dynamics, therefore providing a basis for revealing the microscopic origin of experimentally observed magnetic field effects. We reveal that simulations carried out on our model are capable of reproducing experimental results as well as generating theoretical predictions related to the efficiency of organic electronic materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Design of Graded Cathode Catalyst Layers with Various Ionomers for Fuel Cell Application

Proton exchange membrane fuel cells (PEMFCs) powered by green hydrogen (H2) have become a promising alternative to conventional hydrocarbon-fueled power generators. Despite technological advancements, further improvements in efficiency, durability, and low-cost production are required for the widespread adoption of PEMFCs. Though numerous approaches to improve PEMFC electrodes have been reported, most strategies utilize a single material set (e.g., one combination of catalyst and ionomer) to improve performance. Alternatively, anisotropic (graded) electrode structures with locally tunable properties may yield superior electrode performance due to improved ionic and gas phase transport. In this work, graded cathode catalyst layers (CCLs) incorporating different ionomers (Nafion D2020, Aquivion D79-25BS, and HOPI) were designed and prepared. Performance screening shows that some of these graded electrode structures have comparable performance to optimized single-ionomer electrode structures (D2020) suggesting some synergistic benefit. Additionally, we show that electrodes with lower equivalent weight (EW) D79 ionomer near the membrane and D2020 ionomer near the gas diffusion media outperformed electrodes with the inverted configuration. Finally, EIS analysis shows some graded ionomer structures (e.g. D79/D2020) have better than expected H+ conductivity, generally leading to better electrode performance. However, further optimization of ionomer content and catalyst ink formulations is needed to improve overall PEMFC performance.

DIRECT ENERGY CONVERSION↗

Neutron noise calculations in a hexagonal geometry and comparison with analytical solutions

This paper presents the development of a neutronic and kinetic solver for hexagonal geometries. The tool is developed based on the diffusion theory with multi-energy groups and multi-groups of delayed neutron precursors allowing the solutions of forward and adjoint problems of static and dynamic states, and is applicable to both thermal and fast systems with hexagonal geometries. In the dynamic problems, the small stationary fluctuations of macroscopic cross sections are considered as noise sources, and then the induced first order noise is calculated fully in the frequency domain. Numerical algorithms for solving the static and noise equations are implemented with a spatial discretization based on finite differences and a power iterative solution. A coarse mesh finite difference method has been adopted for speeding up the convergence. Since no other numerical tool could calculate frequency-dependent noise in hexagonal geometry, validation calculations have been performed and benchmarked to analytical solutions based on a 2-D homogeneous system with two-energy groups and one-group of delayed neutron precursor, in which point-like perturbations of thermal absorption cross section at central and non-central positions are considered as noise sources. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Validation Data Acquisition in HTTF during PCC Events

A validation experimental campaign was conducted in an Integral Effect Test (IET) facility of High Temperature Gas Reactors (HTGR), the High-Temperature Test Facility (HTTF) at Oregon State University (OSU). The HTTF simulates Depressurized and Pressurized Conduction Cooldowns (DCC and PCC). This campaign required the development of a new laser spectroscopic diagnostic to measure velocity in the challenging conditions of the HTTF: low speed (~1 m/s) gas flows at HTGR prototypical temperature and 1/10th pressure. This was a collaborative effort between co-PIs at The George Washington University (GW), Oregon State University (OSU), and NASA Langley Research Center. The main accomplishments of this project include the record for dynamic range for velocimetry, highest accuracy obtained with this technique, successful deployment in an IET leading to new validation matrix for CFD. These are detailed below and in manuscript appended to this executive summary. For this project, we introduced a new class of laser spectroscopic diagnostics to Thermal-Hydraulics to measure velocity of gases; furthermore, the diagnostic was demonstrated in-situ in an IET during DCC events. In such scenarios, particles used in mainstream techniques, like Particle Image Velocimetry (PIV) are not appropriate as they settle down too rapidly and also contaminate the experimental facility. Molecular tracers stay mixed within the working gas and can seed the flow in a technique called Molecular Tagging Velocimetry (MTV). In MTV a molecular tracer is photo-dissociated by a first (write) laser into radicals or molecules. The pattern created by the write laser is then interrogated with planar laser-induced fluorescence (PLIF), the read pulse(s), which are recorded with a camera. The pattern is probed and matched at two times (interval or probe time, dt), resulting in a time-of-flight velocimetry technique. This project demonstrated a new application range for MTV in gases. MTV has been extensively used for high-speed gas dynamics: molecules do not suffer from inertia associated with particles in regions of high accelerations (such as shock waves). For high-speed flows, the probe time is on the order of 0.1-10 µs. In time-of-flight techniques, the velocity is reconstructed from the recorded displacement, dx, between time interval dt: U = dx/dt. To first order, the resolution of the pattern matching is typically a (fixed) fraction of a pixel; therefore the resolution of the velocity measurement increases with dx and, for a given flow velocity, dt. For the low-speed flows of interest here, this requires long probe time, and therefore finding tracers that are stable for a long time. After an extensive literature survey we identified 6 tracers that would be suitable. We also worked with Spectra Physics to design a custom laser system that is flexible enough to probe these various tracers. The final instrument is made of 4 lasers that are integrated into a compact and movable cart to facilitate the deployment of the instrument off-site. One unique aspect of the laser system is the tunable dye laser (needed to do PLIF of the radicals). By being pumped by two Nd:YAG lasers, this laser is able to provide two pulses with any probe time (dt) for wavelengths between 200 and 900 nm with a linewidth of 3 pm. Accommodating the dual-pulse capability required customizing the pump lasers and the dye laser, making this system unique. Additionally, a custom UV-intensified CCD camera capable of frame-straddling and exposure time as short as 10 ns was acquired. More than two-third of the instrument cost (~$200k) was covered by a cost sharing from the PI’s startup funds. Of the 6 potential tracers identified, two were tested extensively (H2O and N2O) at high-temperature and elevated pressure at GW prior to conducting the tests at OSU. For both tracers, an Ar-F laser (193 nm, 10 ns, 10 mJ/pulse) photo-dissociates trace amounts (0.1 – 5% molar fraction) of the gases into OH and NO. OH is probed at 281.905 nm and NO at 226.186 nm (4-30mJ, 10 ns pulse). NO is very stable in inert environment and we could reconstruct velocities precisely with probe delay times as long at 40 ms, compared to 2 ms for OH, resulting in increased precision. Therefore, for the tests in the HTTF, N2O was selected as the final tracer and its optimal concentration was determined to be 0.5% molar. In the GW laboratory, we obtained a velocity measurement precision of 4 mm/s at conditions mimicking the HTTF tests. We also demonstrated the benefit of molecular tracers over particles and the flexibility of our instrument by continuously measuring velocity from 1000 to 0.1 m/s in a blow down test. This represents a dynamic range on velocity of 80 dB, which is unprecedented for velocimetry! Once the diagnostic was optimized, two two-week campaigns were conducted in the HTTF. This experimental campaign required numerous practical solutions to successfully and rapidly deploy this very sensitive table-top diagnostic in-situ. This was accomplished thanks to very extensive planning from all parties involved in the project. Due to constraints on availability of the HTTF, only adiabatic DCC tests were conducted. Specifically, the velocity was measured at the exit of the hot leg into the Reactor Cavity Simulation Tank (RCST). We were able to record the gas velocity associated with the air ingress for long times (up to one hour). Depending on the initial conditions (density of the gases in the Reactor Pressure Vessel and the RCST) the lock exchange lasted 30-60 s past the valve openings. Large vortices were identified within the shear layer at the interface of the cross flow. These vortices are an indication of shear-driven instability between the counter streams of different density. After the lock-exchange, a small residual flow was still present and its strength decayed exponentially. It finds its source in a buoyancy difference between the gases in both vessels due to the mixing and diffusion of the gases in the cross-over duct and the RPV lower plenum. We did not expect to capture this part of the transient, which would have been less significant (if not negligible) in the presence of a heated core; the precision of our measurement in the HTTF was 6 mm/s, which exceeded our expectations for such a complex test! Using commercial software (Fluent) we simulated the experiment with Computational Fluid Dynamic (CFD) using Reynolds Averaged Navier-Stokes (RANS) turbulence models. We were not able to accurately reproduce the shear layer instability, mixing, and diffusion captured in our experimental data, even with large 3D meshes (>10 M points). We attributed this to the inadequacy of the RANS models for the flow of interest: our experimental Reynolds number was too low and was outside the applicability range of the RANS models. Instead, we should have employed CFD with less modeling, such as Large Eddy (LES) or even Direct Numerical (DNS) simulations. Unfortunately, we did not have the computational resources to do so. Finally, it should be noted that this experimental technique holds great promises for measuring velocity in thermalhydraulics. MTV has an application domain (highlighted in the narrative) that is complimentary to more established (and simpler) techniques, such as PIV. It has attracted the attention of several thermalhydraulics researchers and sponsors to measure velocity in harsh environments. The spectroscopic nature of this technique also enables to expand it to measure temperature and concentration simultaneously to velocity. In lieu of a narrative, 4 detailed journal publications are appended to this executive summary to describe the various development and accomplishments during this project. They include 2 published manuscripts (Experiments in Fluids and Measurement Science and Technology), and two invited contributions following the NURETH-17 conference (Nuclear Science and Technology). Our current Nuclear Engineering and Design draft is not included at this time.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A NUMERICAL SIMULATION OF COSMIC-RAY MODULATION NEAR THE HELIOPAUSE

Based on a hybrid galactic cosmic-ray transport model, which incorporated MHD global heliospheric data into Parker’s cosmic-ray transport equation, we studied the behavior of the transport of galactic cosmic rays and the corresponding gradients in their flux near the heliopause (HP). We found that, (1) by increasing the ratio of the parallel diffusion coefficient to the perpendicular diffusion coefficient in the interstellar magnetic field of the outer heliosheath, the simulated radial flux near the HP increases as well. As the ratio multiplying factor reached 10{sup 10}, the radial flux experienced a sudden jump near the HP, similar to what Voyager 1 observed in 2012. (2) The effect of changing the diffusion coefficients’ ratio on the radial flux variation depends on the energy of the cosmic rays, the lower the energy, the more pronounced the effect is. (3) The magnitude of the diffusion coefficients also affect the radial flux near the HP, the modulation beyond the HP varies by adjusting the magnitude multiplying factor.

79 ASTRONOMY AND ASTROPHYSICS↗

Laboratory Study of Collisionless Magnetic Reconnection

A concise review is given on the past two decades’ results from laboratory experiments on collisionless magnetic reconnection in direct relation with space measurements, especially by the Magnetospheric Multiscale (MMS) mission. Highlights include spatial structures of electromagnetic fields in ion and electron diffusion regions as a function of upstream symmetry and guide field strength, energy conversion and partitioning from magnetic field to ions and electrons including particle acceleration, electrostatic and electromagnetic kinetic plasma waves with various wavelengths, and plasmoid-mediated multiscale reconnection. Combined with the progress in theoretical, numerical, and observational studies, the physics foundation of fast reconnection in collisionless plasmas has been largely established, at least within the parameter ranges and spatial scales that were studied. Immediate and long-term future opportunities based on multiscale experiments and space missions supported by exascale computation are discussed, including dissipation by kinetic plasma waves, particle heating and acceleration, and multiscale physics across fluid and kinetic scales.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multiscale dendritic needle network model of alloy solidification with fluid flow

We present a mathematical formulation of a multiscale model for solidification with convective flow in the liquid phase. The model is an extension of the dendritic needle network approach for crystal growth in a binary alloy. We propose a simple numerical implementation based on finite differences and step-wise approximations of parabolic dendritic branches of arbitrary orientation. Results of the two-dimensional model are verified against reference benchmark solutions for steady, unsteady, and buoyant flow, as well as steady-state dendritic growth in the diffusive regime. Simulations of equiaxed growth under forced flow yield dendrite tip velocities within 10% of quantitative phase-field results from the literature. Finally, we perform illustrative simulations of polycrystalline solidification using physical parameters for an aluminum-10wt% copper alloy. Resulting microstructures show notable differences when taking into account natural buoyancy in comparison to a purely diffusive transport regime. The resulting model opens new avenues for computationally and quantitatively investigating the influence of fluid flow and gravity-induced buoyancy upon the selection of dendritic microstructures. Further ongoing developments include an equivalent formulation for directional solidification conditions and the implementation of the model in three dimensions, which is critical for quantitative comparison to experimental measurements.

36 MATERIALS SCIENCE↗

Axisymmetric global gravitational equilibrium for magnetized, rotating hot plasma

In this work, we present analytic solutions for three-dimensional magnetized axisymmetric equilibria confining rotating hot plasma in a gravitational field. Our up–down symmetric solution to the full Grad–Shafranov equation can exhibit equatorial plane localization of the plasma density and current, resulting in disk equilibria for the plasma density. For very weak magnetic fields and high plasma pressure, we find strongly rotating thin plasma disk gravitational equilibria that satisfy strict Keplerian motion provided the gravitational energy is much larger than the plasma pressure, which must be large compared to the magnetic energy of the poloidal magnetic field. When the rotational energy exceeds the gravitational energy and it is larger than the plasma pressure, diffuse disk equilibrium solutions continue to exist provided the poloidal magnetic energy remains small. For stronger magnetic fields and lower plasma pressure and rotation, we can also find gravitational equilibria with strong localization to the equatorial plane. However, a toroidal magnetic field is almost always necessary to numerically verify these equilibria are valid solutions in the presence of gravity for the cases considered in Catto & Krasheninnikov (J. Plasma Phys., vol. 81, 2015, 105810301). In all cases both analytic and numerical results are presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A coarse-grained model for the simulations of biomolecular interactions in cellular environments

The interactions of bio-molecules constitute the key steps of cellular functions. However, in vivo binding properties differ significantly from their in vitro measurements due to the heterogeneity of cellular environments. Here we introduce a coarse-grained model based on rigid-body representation to study how factors such as cellular crowding and membrane confinement affect molecular binding. The macroscopic parameters such as the equilibrium constant and the kinetic rate constant are calibrated by adjusting the microscopic coefficients used in the numerical simulations. By changing these model parameters that are experimentally approachable, we are able to study the kinetic and thermodynamic properties of molecular binding, as well as the effects caused by specific cellular environments. We investigate the volumetric effects of crowded intracellular space on bio-molecular diffusion and diffusion-limited reactions. Furthermore, the binding constants of membrane proteins are currently difficult to measure. We provide quantitative estimations about how the binding of membrane proteins deviates from soluble proteins under different degrees of membrane confinements. The simulation results provide biological insights to the functions of membrane receptors on cell surfaces. Overall, our studies establish a connection between the details of molecular interactions and the heterogeneity of cellular environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Evolution of the marker distribution in gyrokinetic $δf$ particle-in-cell simulations

The evolution of the particle weight in a δf particle-in-cell simulation depends on the marker distribution that can evolve in a turbulent field due to turbulent diffusion. When Monte Carlo methods are used to implement the test particle collision operator, or when the particle motion is not strictly Hamiltonian in a collisionless simulation, the marker distribution will evolve along the particle trajectory and, in general, cannot be known exactly. A two-dimensional numerical marker distribution is proposed as an approximation. It is shown to be advantageous over other common methods for evaluating the marker distribution in long-time turbulence simulations. A generalized two-weight δf-method is proposed to mitigate the marker evolution problem.

, Monte Carlo methods↗

Multilevel Methods for Elliptic Problems with Highly Varying Coefficients on Nonaligned Coarse Grids

We generalize the analysis of classical multigrid and two-level overlapping Schwarz methods for 2nd order elliptic boundary value problems to problems with large discontinuities in the coefficients that are not resolved by the coarse grids or the subdomain partition. The theoretical results provide a recipe for designing hierarchies of standard piecewise linear coarse spaces such that the multigrid convergence rate and the condition number of the Schwarz preconditioned system do not depend on the coefficient variation or on any mesh parameters. One assumption we have to make is that the coarse grids are sufficiently fine in the vicinity of cross points or where regions with large diffusion coefficients are separated by a narrow region where the coefficient is small. We do not need to align them with possible discontinuities in the coefficients. The proofs make use of novel stable splittings based on weighted quasi-interpolants and weighted Poincaré-type inequalities. Finally, numerical experiments are included that illustrate the sharpness of the theoretical bounds and the necessity of the technical assumptions.

97 MATHEMATICS AND COMPUTING↗

Double Diffusion Maps and their Latent Harmonics for scientific computations in latent space

In this work, we introduce a data-driven approach to building reduced dynamical models through manifold learning; the reduced latent space is discovered using Diffusion Maps (a manifold learning technique) on time series data. A second round of Diffusion Maps on those latent coordinates allows the approximation of the reduced dynamical models. This second round enables mapping the latent space coordinates back to the full ambient space (what is called lifting); it also enables the approximation of full state functions of interest in terms of the reduced coordinates. In our work, we develop and test three different reduced numerical simulation methodologies, either through pre-tabulation in the latent space and integration on the fly or by going back and forth between the ambient space and the latent space. The data-driven latent space simulation results, based on the three different approaches, are validated through (a) the latent space observation of the full simulation through the Nyström Extension formula, or through (b) lifting the reduced trajectory back to the full ambient space, via Latent Harmonics. Latent space modeling often involves additional regularization to favor certain properties of the space over others, and the mapping back to the ambient space is then constructed mostly independently from these properties; here, we use the same data-driven approach to construct the latent space and then map back to the ambient space.

97 MATHEMATICS AND COMPUTING↗

Faster ablative Kelvin–Helmholtz instability growth in a magnetic field

Shear flows along a plasma interface will quickly grow unstable due to the Kelvin–Helmholtz instability. If there is a concurrent temperature gradient across the interface, higher modes are stabilized by the thermal diffusion. These ablative effects must be considered in, for example, jet features in inertial confinement fusion hot-spots, or plasma plumes in young supernovae. We show that magnetization of the plasma can greatly affect the instability, even if magnetic pressure is small. This is because electrons are localized by their gyromotion, reducing the heat flux and material ablation. We use a two-dimensional numerical extended-magnetohydrodynamics approach to assess this effect for dense fusion conditions. In comparison with the unmagnetized case, self-generated Biermann fields make only a minor difference to growth rates. However, simulations with a large 50 kT external field found that the growth rate of the least stable mode increased by 40%. In conclusion, this has implications for mix processes in Z-pinches and magnetized inertial confinement fusion concepts.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Magnetically self-regulated formation of early protoplanetary disks

The formation of protoplanetary disks during the collapse of molecular dense cores is significantly influenced by angular momentum transport, notably by the magnetic torque. In turn, the evolution of the magnetic field is determined by dynamical processes and non-ideal MHD effects such as ambipolar diffusion. Considering simple relations between various timescales characteristic of the magnetized collapse, we derive an expression for the early disk radius, r≃18 au (η{sub AD}/0.1s){sup 2/9}(B{sub z}/0.1G){sup −4/9}(M/0.1M{sub ⊙}){sup 1/3}, where M is the total disk plus protostar mass, η{sub AD} is the ambipolar diffusion coefficient, and B {sub z} is the magnetic field in the inner part of the core. This is significantly smaller than the disks that would form if angular momentum was conserved. The analytical predictions are confronted against a large sample of 3D, non-ideal MHD collapse calculations covering variations of a factor 100 in core mass, a factor 10 in the level of turbulence, a factor 5 in rotation, and magnetic mass-to-flux over critical mass-to-flux ratios 2 and 5. The disk radius estimates are found to agree with the numerical simulations within less than a factor 2. A striking prediction of our analysis is the weak dependence of circumstellar disk radii upon the various relevant quantities, suggesting weak variations among class-0 disk sizes. In some cases, we note the onset of large spiral arms beyond this radius.

79 ASTRONOMY AND ASTROPHYSICS↗

The influence of vacancy generation at the initial stage of ion implantation

The paper presents a coupled isothermal model at the initial stage of a solid surface treatment with particle beams. Mechanical stresses arising due to the interaction of particles with the surface affect the redistribution of the implanted impurity. Vacancies in the metal surface and their generation under stress are also taken into account. The kinetic law is formulated on the basis of thermodynamics of irreversible processes. The authors used numerical investigation methods. As a result, they have obtained the distributions of impurity concentration and deformations for various time moments. The authors also compare the concentration and deformation profiles with and without vacancies and study the influences of some model parameters. The effect of vacancy generation on the diffusion has been established to lead to an increase in the depth of penetration, as well as in the concentration of impurities.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Two-dimensional hybrid Monte Carlo–fluid modelling of dc glow discharges: Comparison with fluid models, reliability, and accuracy

Two-dimensional hybrid Monte Carlo–fluid numerical code is developed and applied to model the dc glow discharge. The model is based on the separation of electrons into two parts: the low energetic (slow) and high energetic (fast) electron groups. Ions and slow electrons are described within the fluid model using the drift-diffusion approximation for particle fluxes. Fast electrons, represented by suitable number of super particles emitted from the cathode, are responsible for ionization processes in the discharge volume, which are simulated by the Monte Carlo collision method. Electrostatic field is obtained from the solution of Poisson equation. The test calculations were carried out for an argon plasma. Main properties of the glow discharge are considered. Current-voltage curves, electric field reversal phenomenon, and the vortex current formation are developed and discussed. The results are compared to those obtained from the simple and extended fluid models. Contrary to reports in the literature, the analysis does not reveal significant advantages of existing hybrid methods over the extended fluid model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

X-RAY POINT-SOURCE POPULATIONS CONSTITUTING THE GALACTIC RIDGE X-RAY EMISSION

Apparently diffuse X-ray emission has been known to exist along the central quarter of the Galactic Plane since the beginning of X-ray astronomy; this is referred to as the Galactic Ridge X-ray emission (GRXE). Recent deep X-ray observations have shown that numerous X-ray point sources account for a large fraction of the GRXE in the hard band (2-8 keV). However, the nature of these sources is poorly understood. Using the deepest X-ray observations made in the Chandra bulge field, we present the result of a coherent photometric and spectroscopic analysis of individual X-ray point sources for the purpose of constraining their nature and deriving their fractional contributions to the hard-band continuum and Fe K line emission of the GRXE. Based on the X-ray color-color diagram, we divided the point sources into three groups: A (hard), B (soft and broad spectrum), and C (soft and peaked spectrum). The group A sources are further decomposed spectrally into thermal and non-thermal sources with different fractions in different flux ranges. From their X-ray properties, we speculate that the group A non-thermal sources are mostly active galactic nuclei and the thermal sources are mostly white dwarf (WD) binaries such as magnetic and non-magnetic cataclysmic variables (CVs), pre-CVs, and symbiotic stars, whereas the group B and C sources are X-ray active stars in flares and quiescence, respectively. In the log N-log S curve of the 2-8 keV band, the group A non-thermal sources are dominant above Almost-Equal-To 10{sup -14} erg cm{sup -2} s{sup -1}, which is gradually taken over by Galactic sources in the fainter flux ranges. The Fe K{alpha} emission is mostly from the group A thermal (WD binaries) and the group B (X-ray active stars) sources.

79 ASTRONOMY AND ASTROPHYSICS↗

A Novel Multi-scale Simulation Strategy for Turbulent Reacting Flows

In this project, a new methodology was proposed to bridge the gap between Direct Numerical Simulation (DNS) and Large Eddy Simulation (LES). This novel methodology, titled Lattice-Based Multiscale Simulation (LBMS), creates a lattice structure of One-Dimensional Turbulence (ODT) models. This model has been shown to capture turbulent combustion with high fidelity by fully resolving interactions between turbulence and diffusion. By creating a lattice of ODT models, which are then coupled, LBMS overcomes the shortcomings of ODT, which are its inability to capture large scale three dimensional flow structures. However, by spacing these lattices significantly apart, LBMS can avoid the curse of dimensionality that creates untenable computational costs associated with DNS. This project has shown that LBMS is capable of reproducing statistics of isotropic turbulent flows while coarsening the spacing between lines significantly. It also investigates and resolves issues that arise when coupling ODT lines, such as flux reconstruction perpendicular to a given ODT line, preservation of conserved quantities when eddies cross a course cell volume and boundary condition application. Robust parallelization is also investigated.

97 MATHEMATICS AND COMPUTING↗