Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Constrained control”

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

Optimizing LSST observing strategy for weak lensing systematics

ABSTRACT The Legacy Survey of Space and Time (LSST) survey will provide unprecedented statistical power for measurements of dark energy. Consequently, controlling systematic uncertainties is becoming more important than ever. The LSST observing strategy will affect the statistical uncertainty and systematics control for many science cases; here, we focus on weak lensing (WL) systematics. The fact that the LSST observing strategy involves hundreds of visits to the same sky area provides new opportunities for systematics mitigation. We explore these opportunities by testing how different dithering strategies (pointing offsets and rotational angle of the camera in different exposures) affect additive WL shear systematics on a baseline operational simulation, using the ρ-statistics formalism. Some dithering strategies improve systematics control at the end of the survey by a factor of up to ∼3–4 better than others. We find that a random translational dithering strategy, applied with random rotational dithering at every filter change, is the most effective of those strategies tested in this work at averaging down systematics. Adopting this dithering algorithm, we explore the effect of varying the area of the survey footprint, exposure time, number of exposures in a visit, and exposure to the Galactic plane. We find that any change that increases the average number of exposures (in filters relevant to WL) reduces the additive shear systematics. Some ways to achieve this increase may not be favorable for the WL statistical constraining power or for other probes, and we explore the relative trade-offs between these options given constraints on the overall survey parameters.

79 ASTRONOMY AND ASTROPHYSICS↗

Spatial and Temporal Organization of Chromosome Duplication and Segregation in the Cyanobacterium Synechococcus elongatus PCC 7942

The spatial and temporal control of chromosome duplication and segregation is crucial for proper cell division. While this process is well studied in eukaryotic and some prokaryotic organisms, relatively little is known about it in prokaryotic polyploids such as Synechococcus elongatus PCC 7942, which is known to possess one to eight copies of its single chromosome. Using a fluorescent repressor-operator system, S. elongatus chromosomes and chromosome replication forks were tagged and visualized. We found that chromosomal duplication is asynchronous and that the total number of chromosomes is correlated with cell length. Thus, replication is independent of cell cycle and coupled to cell growth. Replication events occur in a spatially random fashion. However, once assembled, replisomes move in a constrained manner. On the other hand, we found that segregation displays a striking spatial organization in some cells. Chromosomes transiently align along the major axis of the cell and timing of alignment was correlated to cell division. This mechanism likely contributes to the non-random segregation of chromosome copies to daughter cells.

59 BASIC BIOLOGICAL SCIENCES↗

The expanded Lagrangian system for constrained optimization problems

Smooth penalty functions can be combined with numerical continuation/bifurcation techniques to produce a class of robust and fast algorithms for constrained optimization problems. The key to the development of these algorithms is the Expanded Lagrangian System which is derived and analyzed in this work. This parameterized system of nonlinear equations contains the penalty path as a solution, provides a smooth homotopy into the first-order necessary conditions, and yields a global optimization technique. Furthermore, the inevitable ill-conditioning present in a sequential optimization algorithm is removed for three penalty methods: the quadratic penalty function for equality constraints, and the logarithmic barrier function (an interior method) and the quadratic loss function (an interior method) for inequality constraints. Although these techniques apply to optimization in general and to linear and nonlinear programming, calculus of variations, optimal control and parameter identification in particular, the development is primarily within the context of nonlinear programming.

Poore, A. B.↗

Aeroelastic Model Structure Computation for Envelope Expansion

Structure detection is a procedure for selecting a subset of candidate terms, from a full model description, that best describes the observed output. This is a necessary procedure to compute an efficient system description which may afford greater insight into the functionality of the system or a simpler controller design. Structure computation as a tool for black-box modeling may be of critical importance in the development of robust, parsimonious models for the flight-test community. Moreover, this approach may lead to efficient strategies for rapid envelope expansion that may save significant development time and costs. In this study, a least absolute shrinkage and selection operator (LASSO) technique is investigated for computing efficient model descriptions of non-linear aeroelastic systems. The LASSO minimises the residual sum of squares with the addition of an l(Sub 1) penalty term on the parameter vector of the traditional l(sub 2) minimisation problem. Its use for structure detection is a natural extension of this constrained minimisation approach to pseudo-linear regression problems which produces some model parameters that are exactly zero and, therefore, yields a parsimonious system description. Applicability of this technique for model structure computation for the F/A-18 (McDonnell Douglas, now The Boeing Company, Chicago, Illinois) Active Aeroelastic Wing project using flight test data is shown for several flight conditions (Mach numbers) by identifying a parsimonious system description with a high percent fit for cross-validated data.

Kukreja, Sunil L.↗

The Primordial Inflation Explorer (PIXIE): A Nulling Polarimeter for Cosmic Microwave Background Observations

The Primordial Inflation Explorer (PIXIE) is a concept for an Explorer-class mission to measure the gravity-wave signature of primordial inflation through its distinctive imprint on the linear polarization of the cosmic microwave background. The instrument consists of a polarizing Michelson interferometer configured as a nulling polarimeter to measure the difference spectrum between orthogonal linear polarizations from two co-aligned beams. Either input can view the sky or a temperature-controlled absolute reference blackbody calibrator. Rhe proposed instrument can map the absolute intensity and linear polarization (Stokes I, Q, and U parameters) over the full sky in 400 spectral channels spanning 2.5 decades in frequency from 30 GHz to 6 THz (1 cm to 50 micron wavelength). Multi-moded optics provide background-limited sensitivity using only 4 detectors, while the highly symmetric design and multiple signal modulations provide robust rejection of potential systematic errors. The principal science goal is the detection and characterization of linear polarization from an inflationary epoch in the early universe, with tensor-to-scalar ratio r < 10..3 at 5 standard deviations. The rich PIXIE data set can also constrain physical processes ranging from Big Bang cosmology to the nature of the first stars to physical conditions within the interstellar medium of the Galaxy.

Kogut, Alan J.↗

Kaon BSM B -parameters using improved staggered fermions from N f = 2 + 1 unquenched QCD

In this paper, we present results for the matrix elements of the additional ΔS = 2 operators that appear in models of physics beyond the Standard Model (BSM), expressed in terms of four BSM B -parameters. Combined with experimental results for ΔM K and ε K , these constrain the parameters of BSM models. We use improved staggered fermions, with valence hypercubic blocking transfromation (HYP)-smeared quarks and N f = 2 + 1 flavors of “asqtad” sea quarks. The configurations have been generated by the MILC Collaboration. The matching between lattice and continuum four-fermion operators and bilinears is done perturbatively at one-loop order. We use three lattice spacings for the continuum extrapolation: a ≈ 0.09 , 0.06 and 0.045 fm. Valence light-quark masses range down to ≈ m$phys\atop{s}$ /13 while the light sea-quark masses range down to ≈ m$phys\atop{s}$ / 20 . Compared to our previous published work, we have added four additional lattice ensembles, leading to better controlled extrapolations in the lattice spacing and sea-quark masses. We report final results for two renormalization scales, μ = 2 and 3 GeV, and compare them to those obtained by other collaborations. Agreement is found for two of the four BSM B-parameters (B 2 and B$SUSY\atop{3}$ ). The other two (B 4 and B 5 ) differ significantly from those obtained using regularization independent momentum subtraction (RI-MOM) renormalization as an intermediate scheme, but are in agreement with recent preliminary results obtained by the RBC-UKQCD Collaboration using regularization independent symmetric momentum subtraction (RI-SMOM) intermediate schemes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Baryon fraction from the BAO amplitude: a consistent approach to parameterizing perturbation growth

Galaxy clustering constrains the baryon fraction Omega_b/Omega_m through the amplitude of baryon acoustic oscillations and the suppression of perturbations entering the horizon before recombination. This produces a different pre-recombination distribution of baryons and dark matter. After recombination, the gravitational potential responds to both components in proportion to their mass, allowing robust measurement of the baryon fraction. This is independent of new-physics scenarios altering the recombination background (e.g. Early Dark Energy). The accuracy of such measurements does, however, depend on how baryons and CDM are modeled in the power spectrum. Previous template-based splitting relied on approximate transfer functions that neglected part of information. We present a new method that embeds an extra parameter controlling the balance between baryons and dark matter in the growth terms of the perturbation equations in the CAMB Boltzmann solver. This approach captures the baryonic suppression of CDM prior to recombination, avoids inconsistencies, and yields a clean parametrization of the baryon fraction in the linear power spectrum, separating out the simple physics of growth due to the combined matter potential. We implement this framework in an analysis pipeline using Effective Field Theory of Large-Scale Structure with HOD-informed priors and validate it against noiseless LCDM and EDE cosmologies with DESI-like errors. The new scheme achieves comparable precision to previous splitting while reducing systematic biases, providing a more robust way to baryon-fraction measurements. In combination with BBN constraints on the baryon density and Alcock-Paczynski estimates of the matter density, these results strengthen the use of baryon fraction measurements to derive a Hubble constant from energy densities, with future DESI and Euclid data expected to deliver competitive constraints.

Crespi, Andrea [U. Waterloo (main); Waterloo U., I↗

An acoustic charge transport imager for high definition television applications

The primary goal of this research is to develop a solid-state television (HDTV) imager chip operating at a frame rate of about 170 frames/sec at 2 Megapixels/frame. This imager will offer an order of magnitude improvements in speed over CCD designs and will allow for monolithic imagers operating from the IR to UV. The technical approach of the project focuses on the development of the three basic components of the imager and their subsequent integration. The camera chip can be divided into three distinct functions: (1) image capture via an array of avalanche photodiodes (APD's); (2) charge collection, storage, and overflow control via a charge transfer transistor device (CTD); and (3) charge readout via an array of acoustic charge transport (ACT) channels. The use of APD's allows for front end gain at low noise and low operating voltages while the ACT readout enables concomitant high speed and high charge transfer efficiency. Currently work is progressing towards the optimization of each of these component devices. In addition to the development of each of the three distinct components, work towards their integration and manufacturability is also progressing. The component designs are considered not only to meet individual specifications but to provide overall system level performance suitable for HDTV operation upon integration. The ultimate manufacturability and reliability of the chip constrains the design as well. The progress made during this period is described in detail.

Hunt, W. D.↗

The Primordial Inflation Explorer (PIXIE) Mission

The Primordial Inflation Explorer (PIXIE) is an Explorer-class mission to map the absolute intensity and linear polarization of the cosmic microwave background and diffuse astrophysical foregrounds over the full sky from frequencies 30 GHz to 6 THz (I cm to 50 I-tm wavelength). PIXIE uses a polarizing Michelson interferometer with 2.7 K optics to measure the difference spectrum between two orthogonal linear polarizations from two co-aligned beams. Either input can view either the sky or a temperature-controlled absolute reference blackbody calibrator. The multimoded optics and high etendu provide sensitivity comparable to kilo-pixel focal plane arrays, but with greatly expanded frequency coverage while using only 4 detectors total. PIXIE builds on the highly successful COBEIFIRAS design by adding large-area polarization-sensitive detectors whose fully symmetric optics are maintained in thermal equilibrium with the CMB. The highly symmetric nulled design provides redundant rejection of major sources of systematic uncertainty. The principal science goal is the detection and characterization of linear polarization from an inflationary epoch in the early universe, with tensor-to-scalar ratio r much less than 10(exp -3). PIXIE will also return a rich data set constraining physical processes ranging from Big Bang cosmology, reionization, and large-scale structure to the local interstellar medium. Keywords: cosmic microwave background, polarization, FTS, bolometer

Kogut, Alan J.↗

Isolating the Source Region of Infrasound Travel Time Variability Using Acoustic Sensors on High-Altitude Balloons

High-altitude balloons carrying infrasound sensor payloads can be leveraged toward monitoring efforts to provide some advantages over other sensing modalities. On 10 July 2020, three sets of controlled surface explosions generated infrasound waves detected by a high-altitude floating sensor. One of the signal arrivals, detected when the balloon was in the acoustic shadow zone, could not be predicted via propagation modeling using a model atmosphere. Considering that the balloon’s horizontal motion showed direct evidence of gravity waves, we examined their role in infrasound propagation. Implementation of gravity wave perturbations to the wind field explained the signal detection and aided in correctly predicting infrasound travel times. Our results show that the impact of gravity waves is negligible below 20 km altitude; however, their effect is important above that height. The results presented here demonstrate the utility of balloon-borne acoustic sensing toward constraining the source region of variability, as well as the relevance of complexities surrounding infrasound wave propagation at short ranges for elevated sensing platforms.

47 OTHER INSTRUMENTATION↗

An acoustic charge transport imager for high definition television applications

The primary goal of this research is to develop a solid-state high definition television (HDTV) imager chip operating at a frame rate of about 170 frames/sec at 2 Megapixels per frame. This imager offers an order of magnitude improvement in speed over CCD designs and will allow for monolithic imagers operating from the IR to the UV. The technical approach of the project focuses on the development of the three basic components of the imager and their integration. The imager chip can be divided into three distinct components: (1) image capture via an array of avalanche photodiodes (APD's), (2) charge collection, storage and overflow control via a charge transfer transistor device (CTD), and (3) charge readout via an array of acoustic charge transport (ACT) channels. The use of APD's allows for front end gain at low noise and low operating voltages while the ACT readout enables concomitant high speed and high charge transfer efficiency. Currently work is progressing towards the development of manufacturable designs for each of these component devices. In addition to the development of each of the three distinct components, work towards their integration is also progressing. The component designs are considered not only to meet individual specifications but to provide overall system level performance suitable for HDTV operation upon integration. The ultimate manufacturability and reliability of the chip constrains the design as well. The progress made during this period is described in detail in Sections 2-4.

Hunt, W. D.↗

The Role of Lightning in Controlling Interannual Variability of Tropical Tropospheric Ozone and OH and its Implications for Climate

Nitrogen oxides (NO(x) = NO + NO2) produced by lightning make a major contribution to the production of the dominant tropospheric oxidants (OH and ozone). These oxidants control the lifetime of many trace gases including long-lived greenhouse gases, and control the source-receptor relationship of inter-hemispheric pollutant transport. Lightning is affected by meteorological variability, and therefore represents a potentially important tropospheric chemistry-climate feedback. Understanding how interannual variability (IAV) in lightning affects IAV in ozone and OH in the recent past is important if we are to predict how oxidant levels may change in a future warmer climate. However, lightning parameterizations for chemical transport models (CTMs) show low skill in reproducing even climatological distributions of flash rates from the Lightning Imaging Sensor (LIS) and the Optical Transient Detector (OTD) satellite instruments. We present an optimized regional scaling algorithm for CTMs that enables sufficient sampling of spatiotemporally sparse satellite lightning data from LIS to constrain the spatial, seasonal, and interannual variability of tropical lightning. We construct a monthly time series of lightning flash rates for 1998-2010 and 35degS-35degN, and find a correlation of IAV in total tropical lightning with El Nino. We use the IAV-constraint to drive a 9-year hindcast (1998-2006) of the GEOS-Chem 3D chemical transport model, and find the increased IAV in LNO(x) drives increased IAV in ozone and OH, improving the model fs ability to simulate both. Although lightning contributes more than any other emission source to IAV in ozone, we find ozone more sensitive to meteorology, particularly convective transport. However, we find IAV in OH to be highly sensitive to lightning NO(x), and the constraint improves the ability of the model to capture the temporal behavior of OH anomalies inferred from observations of methyl chloroform and other gases. The sensitivity of OH is explained using photochemical reaction rates which show a "magnification" effect of the initial lightning NO perturbation on OH primary production, HO(x) recycling, and OH loss frequencies. This influence on OH may represent a negative feedback, if lightning increases in a warming world..

Murray, Lee T.↗

SpF: Enabling Petascale Performance for Pseudospectral Dynamo Models

Pseudospectral (PS) methods possess a number of characteristics (e.g., efficiency, accuracy, natural boundary conditions) that are extremely desirable for dynamo models. Unfortunately, dynamo models based upon PS methods face a number of daunting challenges, which include exposing additional parallelism, leveraging hardware accelerators, exploiting hybrid parallelism, and improving the scalability of global memory transposes. Although these issues are a concern for most models, solutions for PS methods tend to require far more pervasive changes to underlying data and control structures. Further, improvements in performance in one model are difficult to transfer to other models, resulting in significant duplication of effort across the research community.We have developed an extensible software framework for pseudospectral methods called SpF that is intended to enable extreme scalability and optimal performance. High-level abstractions provided by SpF unburden applications of the responsibility of managing domain decomposition and load balance while reducing the changes in code required to adapt to new computing architectures. The key design concept in SpF is that each phase of the numerical calculation is partitioned into disjoint numerical kernels that can be performed entirely in-processor. The granularity of domain-decomposition provided by SpF is only constrained by the data-locality requirements of these kernels. SpF builds on top of optimized vendor libraries for common numerical operations such as transforms, matrix solvers, etc., but can also be configured to use open source alternatives for portability. SpF includes several alternative schemes for global data redistribution and is expected to serve as an ideal testbed for further research into optimal approaches for different network architectures.In this presentation, we will describe the basic architecture of SpF as well as preliminary performance data and experience with adapting legacy dynamo codes. We will conclude with a discussion of planned extensions to SpF that will provide pseudospectral applications with additional flexibility with regard to time integration, linear solvers, and discretization in the radial direction.

Pseudospectral (PS)↗

Evaluating turbulence models at high energy densities

The RESHOCK campaign at Lawrence Livermore National Laboratory has been working on increasing our understanding of the evolution of turbulent, unstable plasma interfaces which are applicable in various HED applications including ICF implosions. A common approach to model such interfaces involves the use of Renolds-Averaged-Navier-Stokes (RANS) models whose parameters have been constrained by theory and experiments. These parameter sets are not unique and the applicability of one set of tuned model parameters to systems with differing shock strengths or material densities has not been well tested. Our new NIF experiments, along with related studies of turbulent shear flows, provide data and model validation for plasma interfaces in the HED regime. Our experiments measure mixing-layer width at an unstable interface in the plasma regime and are specifically designed to challenge mix models of the Reynolds-Averaged-Navier-Stokes (RANS) type where turbulence is assumed to be fully developed, i.e. the flow exhibits a broad spectrum of length scales without memory of the initial condition.4 We utilize precise control of the initial interface conditions (densities on both sides of the interface, an initial rippled perturbation pattern consisting of equally weighted wavelengths between 10 and 20 microns, and initial material compositions) along with a repeatable drive history

36 MATERIALS SCIENCE↗

Large Space Optics: From Hubble to JWST and Beyond

If necessity truly is the mother of invention, then advances in lightweight space mirror technology have been driven by launch vehicle mass and volume constraints. In the late 1970 s, at the start of Hubble development, the state of the art in ground based telescopes was 3 to 4 meter monolithic primary mirrors with masses of 6000 to 10,000 kg - clearly too massive for the planned space shuttle 25,000 kg capability to LEO. Necessity led Hubble to a different solution. Launch vehicle mass constraints (and cost) resulted in the development of a 2.4 meter lightweight eggcrate mirror. At 810 kg (180 kg/m2), this mirror was approximately 7.4% of HST s total 11,110 kg mass. And, the total observatory structure at 4.3 m x 13.2 m fit snuggly inside the space shuttle 4.6 m x 18.3 m payload bay. In the early 1990 s, at the start of JWST development, the state of the art in ground based telescopes was 8 meter class monolithic primary mirrors (16,000 to 23,000 kg) and 10 meter segmented mirrors (14,400 kg). Unfortunately, launch vehicles were still constrained to 4.5 meter payloads and 25,000 kg to LEO or 6,600 kg to L2. Furthermore, science now demanded a space telescope with 6 to 8 meter aperture operating at L2. Mirror technology was identified as a critical capability necessary to enable the next generation of large aperture space telescopes. Specific telescope architectures were explored via three independent design concept studies conducted during the summer of 1996 (1). These studies identified two significant architectural constraints: segmentation and areal density. Because the launch vehicle fairing payload dynamic envelop diameter is approximately 4.5 meters, the only way to launch an 8 meter class mirror is to segment it, fold it and deploy it on orbit - resulting in actuation and control requirements. And, because of launch vehicle mass limits, the primary mirror allocation was only 1000 kg - resulting in a maximum areal density of 20 kg/m2. At the inception of JWST in 1996, such a capability did not exist. A highly successful technology development program was initiated resulting in matured and demonstrated mirror technology for JWST (2, 3). Today, the JWST 6.5 meter primary mirror has an areal density of 25 kg/m2 for a total mass of 625 kg or 9.6% of the total JWST observatory mass of 6,500 kg. Looking into the future, science requires increasing larger collecting apertures. Ground based telescopes are already moving towards 30+ meter mirrors. The only way to meet this challenge for space telescopes is via even lower areal density mirrors or on-orbit assembly or larger launch vehicles (4). The planned NASA Ares V with its 10 meter fairing and 55,000 kg payload to L2 eliminates this constraint (5).

Stahl, H. Philip↗

Subtask 2.18 - Advancing CO 2 Capture Technology: Partnership for CO 2 Capture (PCO 2 C) Phase III

Industries and utilities continue to investigate ways to decrease their carbon footprint. Carbon capture and storage (CCS) can enable existing power generation facilities to meet the current national CO 2 reduction goals. The Partnership for CO2 Capture Phase III focused on several important research areas in an effort to find ways to decrease the cost of capture across both precombustion and postcombustion platforms. Two flue gas pretreatment technologies for postcombustion capture, an SO 2 reduction scrubbing technology from Cansolv Technologies Inc. and the Tri-Mer filtration technology that combines particulate, NOx, and SO 2 control, were evaluated on the Energy & Environmental Research Center’s (EERC’s) pilot-scale test system. Pretreating the flue gas should enable more efficient, and therefore less expensive, CO 2 capture. Both technologies were found to be effective in pretreating flue gas prior to CO 2 capture. Two new postcombustion capture solvents were tested, one from the Korea Carbon Capture and Sequestration R&D Center (KCRC) and one from CO 2 Solutions Incorporated. Both of these solvents showed the ability to capture CO 2 while requiring less regeneration energy, which would reduce the cost of capture. Hydrogen separation membranes from Commonwealth Scientific and Industrial Research Organisation were evaluated through precombustion testing. They are composed of vanadium alloy, which is less expensive than the palladium alloys that are typically used. Their performance was comparable to that of other membranes that have been tested at the EERC. Aspen Plus® software was used to model the KCRC and CO 2 Solutions solvents and found that they would result in significantly improved overall plant performance. The modeling effort also showed that the parasitic steam load at partial capture of 45% is less than half that of 90% overall capture, indicating savings that could be accrued if 90% capture is not required. Modeling of three regional power plants using the Carnegie Mellon Integrated Environmental Control Model showed that, among other things, the use of a bypass during partial capture may minimize the size of the capture tower(s) and result in a slight reduction in the revenue required to operate the capture facility. The results reinforced that a one-size-fits-all approach cannot be taken to adding capture to a power plant. Laboratory testing indicated that Fourier transform infrared spectroscopy could be used to continuously sample stack emissions at CO 2 capture facilities to detect and quantify any residual amine or its degradation products, particularly nitrosamines. The information gathered during Phase III is important for utility stakeholders as they determine how to reduce their CO 2 emissions in a carbon-constrained world. This subtask was funded through the EERC–U.S. Department of Energy (DOE) Joint Program on Research and Development for Fossil Energy-Related Resources Cooperative Agreement No. DE-FC26-08NT43291. Nonfederal funding was provided by the North Dakota Industrial Commission, PPL Montana, Nebraska Public Power District, Tri-Mer Corporation, Montana–Dakota Utilities Co., Basin Electric Power Cooperative, KCRC/Korean Institute of Energy Research, Cansolv Technologies, and CO 2 Solutions, Inc.

01 COAL, LIGNITE, AND PEAT↗

Inferred rheological structure and mantle conditions from postseismic deformation following the 2010 M w 7.2 El Mayor-Cucapah Earthquake

The 2010M w 7.2 El Mayor-Cucapah earthquake provides a unique target of postseismic study as deformation extends across several distinct geological provinces, including the cold Mesozoic arc crust of the Peninsular Ranges and newly formed, hot, extending lithosphere within the Salton Trough. We use five years of global positioning system measurements to invert for afterslip and constrain a 3-D finite-element model that simulates viscoelastic relaxation. We find that afterslip cannot readily explain far-field displacements (more than 50 km from the epicentre). These displacements are best explained by viscoelastic relaxation of a horizontally and vertically heterogeneous lower crust and upper mantle. Lower viscosities beneath the Salton Trough compared to the Peninsular Ranges and other surrounding regions are consistent with inferred differences in the respective geotherms. Our inferred viscosity structure suggests that the depth of the Lithosphere/Asthenosphere Boundary (LAB) is ~65 km below the Peninsular Ranges and ~32 km beneath the Salton Trough. These depths are shallower than the corresponding seismic LAB. This suggests that the onset of partial melting in peridotite may control the depth to the base of the mechanical lithosphere. In contrast, the seismic LAB may correspond to an increase in the partial melt percentage associated with the change from a conductive to an adiabatic geotherm.

58 GEOSCIENCES↗

Development of a two-zone zero-dimensional mixing model for actively fueled prechamber engines

In this numerical work, a zero-dimensional, two-zone mixing model was developed and validated to capture the mixture preparation process of an actively fueled prechamber engine. The model was developed for ethanol-gasoline fuel blends but is applicable to various fuels of interest relevant to prechamber applications. The model was exercised to demonstrate the effects of boundary conditions, fuel composition, and residual content on the mixture stoichiometry predictions. The 0D model showed reasonable agreement in prechamber equivalence ratio predictions with higher fidelity computational fluid dynamics modeling, but the accuracy of model was constrained to the governing assumptions of the injection process. The modeling assumptions surrounding the non-dimensional injection of mass into the prechamber were scrutinized to identify key areas of model improvement. The results indicate heat transfer effects play a crucial role in the evaporation of the fuel within the confined auxiliary chamber, especially when considering fuels with high latent heat of vaporization. Additionally, the presence of fuel films that collect on the prechamber surfaces complicated the modeled evaporative cooling loses as the modeling suggests the films evaporate primarily by means of conduction from the boundary. Furthermore, the developed model serves as a simplified framework to generate calibration maps of prechamber fueling or run real time as an open loop control strategy such that desired prechamber mixture proportions are effectively targeted.

Zero-dimensional↗