Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Compression Ratio”

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 199 records · Page 11

Autoignition study of iso -cetane/tetralin blends at low temperature

iso-Cetane and tetralin are the two fuel components commonly considered in literature diesel surrogates, and are hydrocarbon classes representative of iso-alkanes and naphthoaromatics, respectively. Since both surrogate components are involved in the key reactions controlling the autoignition reactivity of diesel surrogates, autoignition studies of iso-cetane, tetralin, and their blending behavior at low temperatures can be helpful in developing/refining chemical kinetic models of diesel surrogates, in order to better predict diesel ignition response under low temperature combustion conditions. In this investigation, autoignition experiments of iso-cetane, tetralin, and their binary mixtures at varying blending ratios have been conducted in a rapid compression machine (RCM) in the temperature/pressure range of 630–930 K and 10–20 bar. For the neat fuel components, the newly-acquired RCM data have been compared with the literature experimental results at some overlapping conditions. In addition, a chemical kinetic model describing the oxidation of iso-cetane and tetralin compiled/updated from the diesel surrogate model developed by Lawrence Livermore National Laboratory is used to simulate the experimental datasets. It is shown that this model exhibits good agreement with the experimental results of tetralin and also predicts the total ignition delay times of iso-cetane at lower temperatures reasonably well. However, the current model fails to fully capture the increase of total ignition delay time with increasing amount of tetralin in the binary blends observed in the present RCM experiments. Therefore, chemical kinetic analyses of the current model are conducted to help identify possible reasons leading to the discrepancies between simulations and experiments. Furthermore, the chemical kinetic interactions between iso-cetane and tetralin are explored by including additional cross-reactions among iso-cetane, tetralin, and their fuel fragments in the current model to assess their effects on model predictions. As a result, based on the chemical kinetic analysis results, future direction for model refinements is discussed.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Sputtered ferroelectric aluminum scandium boron nitride (Al 1−x−y B x Sc y N)/ n -GaN heterostructures

This work demonstrates ferroelectric switching in magnetron sputtered Al 1−x−y B x Sc y N/ n -GaN heterostructures. Using high power impulse magnetron sputtering, a silicon doped n -GaN bottom electrode with an electron concentration of 6.0 × 10 19 cm −3 is grown on c-plane sapphire. Al 1−x−y B x Sc y N films are prepared on the GaN surface with Al:B:Sc ratios that produce tensile, lattice matched, or compressive epitaxial strains. X-ray diffraction shows that lattice matched and compressively strained Al 1−x−y B x Sc y N compositions are pseudomorphic, while partial relaxation is observed for tensilely strained Al 1−x−y B x Sc y N/ n -GaN heterostructures. Electrically, the Al 1−x−y B x Sc y N/ n -GaN stacks show robust hysteresis; the P–E loops are fully saturated with both lattice matched and compressively strained Al 1−x−y B x Sc y N compositions exhibiting remanent polarization values of 135 μC/cm 2 . For comparison, Al 1−x−y B x Sc y N films are also prepared on metal organic chemical vapor deposition and single crystal GaN substrates to extend strain and morphology trends to more common substrate types. This report validates that sputter deposition is a feasible technique for fabricating strain-tunable ferroelectric III–N heterostructures with high crystalline fidelity and smooth surface morphologies.

42 ENGINEERING↗

Spin-density-wave order controlled by uniaxial stress in CeAuSb 2

We report the tetragonal heavy-fermion compound CeAuSb 2 (space group P4/nmm ) exhibits incommensurate spin-density wave (SDW) order below T N ≈ 6.5 K with the propagation vector q A = (δ A , δ A , 1/2 ). The application of uniaxial stress along the [010] direction induces a sudden change in the resistivity ratio ρ a /ρ b at a compressive strain of ε ≈ - 0.5 %. Here we use neutron scattering to show that the uniaxial stress induces a first-order transition to a SDW state with a different propagation vector (0, δ B ,1/2) with δ B = 0.25.The magnetic structure of the new (B) phase consists of Ce layers with ordered moments alternating with layers with zero moment stacked along the c axis. The ordered layers have an up-up-down-down configuration along the b axis. This is an unusual situation in which the loss of spatial inversion in a metallic system is driven by the magnetic order. We argue that the change in SDW wave vector leads to Fermi-surface reconstruction and a concomitant change in the transport properties.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Mechanical Properties of Functionally Graded Concrete Lining for Deep Underground Structures

With the mining depth of coal resources increasing, the thickness of traditional lining for deep mines becomes large. The bearing capacity of the outer lining cannot be fully utilized, so a new type of functionally graded lining structure (FGL) with radial Young’s modulus varying in gradient is proposed. In this study, through theoretical analysis and numerical simulation, the mechanical properties of the functionally graded lining were studied, including the characteristics of the elastic working state and the ultimate bearing state. The influences of the structural parameters of the functionally graded lining, including the inner radius, the thickness, Young’s modulus, the compressive strength of the concrete, Poisson’s ratio, and the number of layers on the mechanical properties were analyzed. The calculation formula of the ultimate bearing capacity of the lining and the calculation formula of the maximum tangential strain at the time of lining failure were put forward. The accuracy of the formula was verified by comparing with the numerical simulation results. The research results provide the basis for the construction of the design theory of a functionally graded lining structure and have great engineering significance for the construction of the kilometer-deep mines in the future.

Li, De-Chun↗

Experimental Study of the Fundamental Properties of Warm Dense Mixtures

The aim of the proposed research was to provide a multi-scale study of the properties of warm dense hydrocarbons, by studying the thermodynamic properties through the equation of state and the microscopic properties by x-ray scattering. Understanding the fundamental properties of warm dense mixtures is an intellectual challenge due to the complexity of the system. Unlike liquids or gases, where constituent particle interactions occur through collisions or bonding between valence electrons, in these strongly coupled systems the atoms will be partially ionized and compressed so tightly together that interactions between the inner core electrons can play a role in the systems chemistry. Advances in computational capabilities and development of new theoretical models have been used to predict the properties of mixtures but there is currently no experimental data of the fundamental interaction between the particles in mixtures to test these predictions against. The proposed research was to experimentally investigate the interaction of the different species in hydrocarbons by measuring the compressibility of substances with different carbon and hydrogen ratios and the complexity of the microscopic interactions through elastic and inelastic x-ray scattering. To study the bulk properties of warm dense hydrocarbons we established the Warm Dense Matter Research Laboratory (WDMRL) in the Institute for Shock Physics at Washington State University. The goal was for experiments in the WDMRL to determine shock loading conditions of interest in the hydrocarbon mixtures. We planned on using Hugoniot EOS measurements and a range of carbon and hydrogen concentrations to determine conditions when the EOS of the mixture varied significantly from that of the classical mixing model. Even though shock transit measurements through aluminum foils suggested pressures upto 400GPa, the experiments in the WDMRL were unsuccessful in getting usable shockwave compression data above 100GPa in polystyrene which was below the pressure of interest for hydrocarbon mixtures (>200GPa). To complete the project, we tested a technique using x-ray phase contrast imaging to map the location of tracer layers in a test sample of polycarbonate to record the material motion in dynamically compressed samples. A technique that will be useful for future warm dense matter experiments. These experiments used <200nm gold layers in polycarbonate samples to measure the material velocity and shock speed using x-ray phase contrast imaging. These results were compared to continuum surface measurements performed at ISP and show that the tracer layer technique can measure hydrodynamic properties accurately in dynamically compressed materials.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Computationally efficient evaluation of optimum homogeneous charge compression ignition operating range with accelerated multizone engine cycle simulation

The very intensive calculations necessary to define a performance map requiring evaluation of over a hundred individual operating points can be efficiently conducted with accelerated multizone for engine cycle simulation, leading to a definition of regions of acceptable and optimum homogeneous charge compression ignition operation. Accelerated multizone for engine cycle simulation has the virtue of enabling accurate evaluation of many operating conditions based on thermal stratification data from a single fluid mechanics run at motored conditions. This is possible because thermal stratification is more sensitive to engine geometry than to operating conditions. In this article, accuracy of accelerated multizone for engine cycle simulation is demonstrated by comparison with experimental data for iso-octane homogeneous charge compression ignition operation over a broad range of lean equivalence ratios (0.14–0.28). The validated accelerated multizone for engine cycle simulation model is then applied to generating a performance map for an engine controlled by appropriately adjusting equivalence ratio and internal exhaust gas recirculation. Regions of acceptable and optimum combustion are identified. It is finally demonstrated that while indicated mean effective pressure remains low for optimum homogeneous charge compression ignition operation (1–4 bar), this is sufficient for a large fraction of typical driving in light-duty vehicles. Much driving including idle can therefore be done in homogeneous charge compression ignition mode at high efficiency and low (essentially zero) NOx and particulate matter emissions.

30 DIRECT ENERGY CONVERSION↗

Asynchronicity in opposed-piston RCMs: Does it matter?

Rapid Compression Machines (RCMs) are widely utilized to study combustion phenomena at engine-relevant conditions, and significant efforts are typically made to create a quiescent environment, particularly for investigations of autoignition chemistry. Opposed-piston configurations can be advantageous due to shorter compression times and reduced surface area to volume ratios. Each side must be actuated simultaneously, but this can be challenging in practice. These devices, like most RCMs, utilize hydraulics for actuation, speed control and arrestation of the piston at the end of the stroke; there is no mechanical control or linkage of the two piston trajectories. To quantify the magnitudes and effects of piston asynchronous behavior, this work employs both detailed experimental measurements and, for the first time, high-fidelity, Direct Numerical Simulation (DNS). The boundary conditions are carefully considered applying insight from high-resolution linear variable differential transformer (LVDT) measurements of the piston trajectory and a zero-dimensional kinematics model of the piston-shaft assembly. Sufficient resolution in the piston crevice region is used. The complicated fluid dynamical behavior that can evolve during piston compression and the ensuing delay processes due to offset timings from t offset = 0-10 ms is elucidated. It is found that near t offset = 6 ms and beyond, the boundary layer on the face of the first-seating piston can be sufficiently perturbed, due initially to reemergence of gas from the crevice of the firstseating piston, so that the adiabatic core can become degraded at long ignition delay times. Substantial mixing of colder gas into the interior of the reaction chamber can alter the measurements, similar to effects previously observed for improper piston crevice configuration. In conclusion, experimental techniques to mitigate asynchronous behavior are discussed and demonstrated.

33 ADVANCED PROPULSION SYSTEMS↗

Discrepancy in the Critical State Void Ratio of Poorly Graded Sand due to Shear Strain Localization

The critical state (CS) concept is a theoretical framework that models the constitutive behavior of soils, including sand and other granular materials. It supports the notion of a unique postfailure state, where the soil ultimately experiences continuous shearing with no change in the plastic volumetric strain. However, the published literature has frequently noted the nonconvergence of sand specimens with different initial densities to a unique CS in the compression plane due to many factors such as specimen fabric, particle morphology, breakage, and grain size distribution. This paper examines the CS for poorly graded (uniform) glass beads and 3 different types of silica sands using 50 conventional triaxial compression (CTC) experiments, 12 oedometer tests, and in situ synchrotron microcomputed tomography (SMT) scans for 10 CTC experiments. The results of the 50 CTC experiments revealed a diffused CS zone in the compression plane, which was further examined using the in situ SMT scans. A thorough three-dimensional image analysis of the SMT scans accurately quantified the evolution of the local void ratio (elocal) versus axial compression within zones of intensive shearing toward the center of the specimen. The evolution of the void ratio was also measured using the entire volume of the specimen (eglobal). At the CS, the elocal/eglobal ratio was assessed to be ~1.25 when a single shear band developed within the scanned specimens and ~1.1-1.15 for specimens that failed via external bulging that was internally manifested by the development of multiple shear bands. This finding suggests that the CS zone in the compression plane can be attributed to the common wrong consideration of eglobal evolution in lieu of elocal within the developing shear bands. Furthermore, the lack of shear band development in uniaxial compression has made the results of the oedometer test reliable in quantifying the CS parameters in the compression plane.

58 GEOSCIENCES↗

Assessing continuum plasticity postulates with grain stress and local strain measurements in triaxially compressed sand

Critical state and continuum plasticity theories have been used in research and engineering practice in soil and rock mechanics for decades. These theories rely on postulated relationships between material stresses and strains. Some classical postulates include coaxiality between stress and strain rates, stress–dilatancy relationships, and kinematic assumptions in shear bands. Although numerical and experimental data have quantified the strains and grain kinematics in such experiments, little data quantifying grain stresses are available. Here, we report the first-known grain stress and local strain measurements in triaxial compression tests on synthetic quartz sands using synchrotron X-ray tomography and 3D X-ray diffraction. We use these data to examine the micromechanics of shear banding, with a focus on coaxiality, stress-dilatancy, and kinematics within bands. Our results indicate the following: 1) elevated deviatoric stress, strain, and stress ratios in shear bands throughout experiments; 2) coaxial principal compressive stresses and strains throughout samples; 3) significant contraction along shear bands; 4) vanishing volumetric strain but nonvanishing stress fluctuations throughout samples at all stages of deformation. Our results provide some of the first-known in situ stress and strain measurements able to aid in critically evaluating postulates employed in continuum plasticity and strain localization theories for sands.

42 ENGINEERING↗

In situ compression artifact removal in scientific data using deep transfer learning and experience replay

The massive amount of data produced during simulation on high-performance computers has grown exponentially over the past decade, exacerbating the need for streaming compression and decompression methods for efficient storage and transfer of this data---key to realizing the full potential of large-scale computational science. Lossy compression approaches such as JPEG when applied to scientific simulation data realized as a stream of images can achieve good compression rates but at the cost of introducing compression artifacts and loss of information. This paper develops a unified framework for in situ compression artifact removal in which the fully convolutional neural network architectures are combined with scalable training, transfer learning, and experience replay to achieve superior accuracy and efficiency while significantly decreasing the storage footprint as compared with the traditional optimization-based approaches. We demonstrate the proposed approach and compare it with compressed sensing postprocessing and other baseline deep learning models using climate simulations and nuclear reactor simulations, both of which are driven by hyperbolic partial differential equations. Our approach when applied to remove the compression artifacts on the JPEG-compressed nuclear reactor simulation data (using a transfer-trained model that was pretrained on the climate simulation data and updated incrementally as the nuclear reactor simulation progressed), achieved a significant improvement---mean peak signal-to-noise ratio of 42.438 as compared with 27.725 obtained with the compressed sensing approach.

97 MATHEMATICS AND COMPUTING↗

Ultrafast X-ray Diffraction Study of a Shock-Compressed Iron Meteorite above 100 GPa

Natural kamacite samples (Fe92.5Ni7.5) from a fragment of the Gibeon meteorite were studied as a proxy material for terrestrial cores to examine phase transition kinetics under shock compression for a range of different pressures up to 140 GPa. In situ time-resolved X-ray diffraction (XRD) data were collected of a body-centered cubic (bcc) kamacite section that transforms to the high-pressure hexagonal close-packed (hcp) phase with sub-nanosecond temporal resolution. The coarse-grained crystal of kamacite rapidly transformed to highly oriented crystallites of the hcp phase at maximum compression. The hcp phase persisted for as long as 9.5 ns following shock release. Comparing the c/a ratio with previous static and dynamic work on Fe and Fe-rich Fe-Ni alloys, it was found that some shots exhibit a larger than ideal c/a ratio, up to nearly 1.65. This work represents the first time-resolved laser shock compression structural study of a natural iron meteorite, relevant for understanding the dynamic material properties of metallic planetary bodies during impact events and Earth’s core elasticity.

36 MATERIALS SCIENCE↗

Effect of fiber sizing and glass fiber laminate hybridization on vibration damping and mechanical properties of banana fiber reinforced polypropylene composites

Modern automotive applications demand lightweight, multifunctional materials to reach mileage goals and natural fiber reinforced composites (NFRCs) are one of the classes of materials proposed as a solution. NFRCs exhibit good vibration damping properties and have low density, but are often limited by processing challenges, poor-fiber matrix compatibility and variable performance. Herein, we investigate non-woven wet-lay of comingled banana fiber (BF), recycled glass fiber (rGF), and polypropylene (PP) fibers to in situ sizing and preparation of composite feedstocks for compression molding. BF and rGF hybrids were prepared by stacking rGF layers during compression molding to produce composites with various fiber ratios. The effect of fiber content, in-situ sizing and ratio of BF to rGF on tensile, flexural and vibration damping performance are investigated. Key results are the significant increase in tensile strength by in situ sizing (40 % sized at 60 wt% BF) and in flexural modulus (+58 % sized at 60 wt% BF) and flexural strength (+41 % sized 60 wt% BF) compared to the unsized equivalent. For BF-rGF hybrid composites with40 wt% total fiber content, flexural strength and modulus were improved by 51 % and 231 % respectively for a 1:1 ratio BF:rGF compared to BF reinforced system. Lastly, identifying the cross-over point where damping and stiffness are optimized for a hybrid composite. These findings demonstrate that these composites can be used as alternative to synthetic fiber or mineral filled composites in automotive applications, particularly where weight reduction, vibration damping and stiffness are desired.

Banana fiber↗

Strong gravity extruding peaks in speed of sound profiles of massive neutron stars

The speed of sound squared (SSS) s 2 in massive neutron stars (NSs) characterizes not only the stiffness of supradense neutron-rich matter within but also equivalently properties of the curved geometry due to the strong-field gravity and matter-geometry coupling. A peaked density or radius profile of s 2 has been predicted for massive NSs using various NS equation of state (EOS) models. However, the nature, cause, location and size of the peak in s 2 profiles are still very EOS model dependent. In this work, we investigate systematically s 2 profiles in massive NSs in a new approach that is independent of the nuclear EOS model and without any presumption about the NS structure and/or composition. In terms of the small quantities (reduced radius, the energy density and pressure scaled by their central values), we perform double-element perturbative expansions in solving perturbatively the scaled Tolman-Oppenheimer-Volkoff (TOV) equations and analyzing s 2 profiles from the Newtonian limit to the general relativistic (GR) case. The GR term in the TOV equations plays a twofold role: it compresses NS matter and modifies the pressure/energy density ratio from small values in Newtonian stars showing no s 2 peak to large ones for massive NSs possessing a peak in their s 2 profiles, and eventually takes away the peak in extremely compact/massive NSs approaching the causality limit. In particular, the peaked behavior in s 2 is expected to emerge near the center of massive NSs like PSR J⁢0740+6620, while a sharp phase transition is unlikely to occur there. Furthermore, these features revealed from our analyses are universal as they are intrinsic properties of the GR stellar structure equations independent of the still very uncertain EOS of supradense neutron-rich matter in NSs.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Lift Equivalence and Cancellation for Airfoil Surge–Pitch–Plunge Oscillations

A NACA 0018 airfoil in freestream velocity is oscillated in longitudinal, transverse, and angle-of-attack directions with respect to the freestream velocity, known as surge, plunge, and pitch. The lift-based equivalence method introduces phase shifts between these three motions to construct in-phase sinusoidal components for maximum lift, waveform construction. Lift cancellation is also determined with the exact negative pitch and plunge motion amplitudes found from the equivalence method to achieve out-of-phase wave destruction. Lift cancellation occurs when a combination of these motions is sought to obtain a constant lift magnitude throughout the oscillation cycle. To achieve both equivalence and cancellation of lift, a prescribed pure pitch amplitude through the Theodorsen theory equates the corresponding equivalent plunge amplitude and pitch–plunge phase shift. These Theodorsen, linear superposition findings of pitch–plunge are leveraged toward the Greenberg theory to determine a closed-form, surge–pitch–plunge solution through the addition of a surge–plunge phase shift and optimal surge amplitude for lift cancellation. The lift cancellation surge–pitch–plunge amplitudes define the equivalence amplitude investigated here and theoretically limit the experiment to combinations of the first lift harmonic of the Greenberg theory. The analytical results are then compared with experimental lift force measurements and dye visualization. Finally, the normalized lift differences due to unsteady wake and boundary-layer behavior are examined to explore the extents of the Greenberg theory for these cases of lift-based equivalence and cancellation.

42 ENGINEERING↗

Simple Statistical Models for Predicting Overpressure Due to CO2 and Low-Salinity Waste-Fluid Injection into Deep Saline Formations

Deep saline aquifers have been used for waste-fluid disposal for decades and are the proposed targets for large-scale CO2 storage to mitigate CO2 concentration in the atmosphere. Due to relatively limited experience with CO2 injection in deep saline formations and given that the injection targets for CO2 sometimes are the same as waste-fluid disposal formations, it could be beneficial to model and compare both practices and learn from the waste-fluid disposal industry. In this paper, we model CO2 injection in the Patterson Field, which has been proposed as a site for storage of 50 Mt of industrial CO2 over 25 years. We propose general models that quickly screen the reservoir properties and calculate pressure changes near and far from the injection wellbore, accounting for variable reservoir properties. The reservoir properties we investigated were rock compressibility, injection rate, vertical-to-horizontal permeability ratio, average reservoir permeability and porosity, reservoir temperature and pressure, and the injectant total dissolved solids (TDS) in cases of waste-fluid injection. We used experimental design to select and perform simulation runs, performed a sensitivity analysis to identify the important variables on pressure build-up, and then fit a regression model to the simulation runs to obtain simple proxy models for changes in average reservoir pressure and bottomhole pressure. The CO2 injection created more pressure compared to saline waste-fluids, when similar mass was injected. However, we found a more significant pressure buildup at the caprock-reservoir interface and lower pressure buildup at the bottom of the reservoir when injecting CO2 compared with waste-fluid injection.

Ansari, Esmail↗

Methods for the separation of CO2 from a gas stream

Disclosed herein are 2-stage membrane separation methods for capturing CO2 from a feed gas. The methods can employ two selectively permeable membranes, which may be the same or different. The selectively permeable membrane can have a carbon dioxide permeance of from 500 to 3000 GPU at 57° C. and 1 atm feed pressure and a carbon dioxide:nitrogen selectivity of from 10 to 1000 at 57° C. and 1 atm feed pressure. High pressure ratios across the membranes can be achieved by compressing the feed gas to a high pressure, by using vacuum pumps to create a lowered pressure on the permeate side of the membrane, by using a sweep stream, or a combination thereof. When a sweep stream is used, the sweep stream may include a portion of the retentate gas stream obtained from the retentate side of one or more of the membranes used.

Ho, W. S. Winston↗

Saturation of Vortex Rings Ejected from Shock-Accelerated Interfaces

Structures evoking vortex rings can be discerned in shock-accelerated flows ranging from astrophysics to inertial confinement fusion. By constructing an analogy between vortex rings produced in conventional propulsion systems and rings generated by a shock impinging upon a high-aspect-ratio protrusion along a material interface, we extend classical, constant-density vortex-ring theory to compressible multifluid flows. Here, we further demonstrate saturation of such vortex rings as the protrusion aspect ratio is increased, thus explaining morphological differences observed in practice.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Autoignition of Premixed Liquefied Petroleum Gas in a Rapid Compression Machine: Experimental Results and Chemical Kinetic Mechanism Reduction

Liquefied petroleum gas (LPG) has many properties that make it an attractive alternative fuel such as lower cost than conventional fuels and an established distribution infrastructure. The development of high efficiency, spark ignited LPG engines is currently limited by engine knock and misfire. The knock and misfire limits are further complicated by the wide range of chemical reactivity in LPG, particularly in international markets. In this study, a rapid compression machine (RCM) was used to characterize the effects of variation in LPG fuel reactivity, equivalence ratio, and exhaust gas recirculation (EGR) on autoignition of LPG/oxidizer/inert/EGR blends. Experiments were conducted with 100% propane (C3H8) and blends of propane with propene, ethane, isobutane, and n-butane. EGR was simulated with mixtures of Ar, CO2, CO, and NO at substitution percentages from 0 to 30 mass percent. Equivalence ratio was varied from 0.75 to 1.5. Ignition delay period under homogeneous autoignition conditions was measured at compressed pressures and temperatures of 23 to 25 bar and 701 to 921 K, respectively. Zero-dimensional simulations of the RCM experiments were performed using Chemkin with several detailed chemical kinetic mechanisms to determine their suitability at predicting ignition delay periods. Multiple reduced chemical kinetic mechanisms were created from the NUIGMech1.1 mechanism to determine the optimal balance between accuracy and computational efficiency for future three-dimensional, time-dependent spark-ignited engine computations.

42 ENGINEERING↗