Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Coupled Simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 937 records · Page 52

MFiX: Fractional-Step Method Implementation

A comprehensive, multiphase computational fluid dynamics (CFD) simulation solves several coupled transport equations including continuity, momentum, species, and energy. Chemical reactions further couple these equations through heats of reaction and rates of formation of products and rates of destruction of reactants. A fractional-step method separates changes attributed to chemical reactions from transport phenomena like convection and diffusion. When the governing equations are split into the transport and reacting components, efficient and independent methodologies can be exploited to solve the different systems. Specifically, discretization of field variable transport equations results in large, sparse matrices which are loosely coupled. These systems are solved in succession using iterative techniques that take advantage of the matrix structure. In contrast, chemical reactions tightly couple field variables locally within the domain (e.g., within a single computational cell) resulting in low dimensional but dense, nonlinear systems that are better solved using direct integration techniques.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Plasma-Assisted Pre-Chamber Ignition System for Highly Dilute Stoichiometric Heavy-Duty Natural Gas Engines (Final Technical Report)

This project explored advanced ignition technologies to significantly enhance efficiency and reduce operating costs for heavy-duty natural gas engines operating at stoichiometric conditions, while meeting ultra-low NOx emission standards. The main goal was to develop and validate a plasma-assisted pre-chamber ignition system that could deliver at least a 2% increase in brake thermal efficiency (BTE) and a 4% decrease in total cost of ownership (TCO) compared to a typical multi-cylinder engine with three-way catalyst aftertreatment, ensuring compatibility with the expected 2027 EPA/CARB regulations. In the first half of the project, the research team concentrated on developing and testing plasma-assisted pre-chamber ignition using nanosecond pulsed discharges. Extensive experiments were conducted in an optically accessible rapid-compression and expansion machine, a constant-volume chamber, and an optical single-cylinder engine. Experiments were coupled with CFD simulations. The work produced unique insights into pre-chamber flame formation, jet ignition, dilution effects, and flame quenching at pressures, temperatures, and dilution levels relevant to engines. Although plasma-assisted ignition showed promise in controlled lab settings, the research also identified fundamental and practical challenges when applying this technology to real engine conditions. Midway through the project, a crucial pivot was made, guided by three key findings. First, the power electronics required for nanosecond plasma discharges were found to be too costly for commercial use, undermining the project’s cost-of-ownership goals. Second, nanosecond plasma ignition was highly sensitive to turbulent flow in the pre-chamber, resulting in lower ignition reliability than traditional spark under engine-like conditions. Third, achieving a truly diffuse low-temperature plasma at high pressures near top dead center was not possible, reducing the anticipated chemical enhancement benefits. These results collectively suggested that continuing with plasma-assisted ignition was unlikely to meet both efficiency and cost objectives. In response, the project shifted focus to a more realistic approach: enhancing traditional spark-based pre-chamber ignition with significantly less spark energy. Using insights gained earlier in the project, the team redesigned the pre-chamber to maintain high dilution tolerance and quick combustion, even with lower ignition energy. Testing confirmed that with optimized pre-chamber design and combustion timing, a lower-energy spark could reliably ignite highly diluted stoichiometric mixtures, reduce burn time, and boost thermal efficiency. Final engine testing and techno-economic analysis verified that this revised approach successfully achieved the project goals. The optimized pre-chamber ignition system provided over a 2% increase in calculated brake thermal efficiency compared to the baseline engine. Notably, the lower ignition energy and simplified hardware reduced component stress, extended maintenance intervals, and lowered the total cost of ownership. When used with stoichiometric operation and traditional three-way aftertreatment, the system remained compatible with near-zero NOx emissions targets without increasing cost or complexity in the emissions control system. In summary, although the project deviated from its initial plasma-assisted ignition idea, the work produced a more practical and commercially viable solution. The results show that precisely optimized, low-energy pre-chamber spark ignition can significantly improve efficiency and reduce overall ownership costs for heavy-duty natural gas engines. This directly aligns with DOE goals for cleaner, more efficient, and cost-effective transportation technologies.

03 NATURAL GAS↗

Optimizing Price-Informed Operation of a Battery Storage System in an Office Building

New prescriptive efficiency requirements amended to US model energy codes historically have been evaluated using an average, blended electricity rate, which obscures demand charges. Post 2019, new measures can be evaluated using a representative time-of-use (TOU) tariff yet more sophisticated analysis methods are needed to consider a variety of TOU tariffs and assess price-informed control. To address these needs, this study couples building prototype simulation model with varying-in-sophistication battery storage operating strategies for different electricity TOU tariffs. The analysis compares the impact of operating a battery storage system following simpler rule-of-thumb methods versus a semi-optimized priced-informed heuristic approach. The investigation demonstrates that the heuristic approach results in greater electricity cost savings and is practical to implement.

Lei, Xuechen↗

Efficient Ensemble-Based Stochastic Gradient Methods for Optimization Under Geological Uncertainty

Ensemble-based stochastic gradient methods, such as the ensemble optimization (EnOpt) method, the simplex gradient (SG) method, and the stochastic simplex approximate gradient (StoSAG) method, approximate the gradient of an objective function using an ensemble of perturbed control vectors. These methods are increasingly used in solving reservoir optimization problems because they are not only easy to parallelize and couple with any simulator but also computationally more efficient than the conventional finite-difference method for gradient calculations. In this work, we show that EnOpt may fail to achieve sufficient improvement of the objective function when the differences between the objective function values of perturbed control variables and their ensemble mean are large. On the basis of the comparison of EnOpt and SG, we propose a hybrid gradient of EnOpt and SG to save on the computational cost of SG. We also suggest practical ways to reduce the computational cost of EnOpt and StoSAG by approximating the objective function values of unperturbed control variables using the values of perturbed ones. We first demonstrate the performance of our improved ensemble schemes using a benchmark problem. Results show that the proposed gradients saved about 30–50% of the computational cost of the same optimization by using EnOpt, SG, and StoSAG. As a real application, we consider pressure management in carbon storage reservoirs, for which brine extraction wells need to be optimally placed to reduce reservoir pressure buildup while maximizing the net present value. Results show that our improved schemes reduce the computational cost significantly.

58 GEOSCIENCES↗

Cabin Thermal Management Analysis for SuperTruck II Next-Generation Hybrid Electric Truck Design

In this article, we present a multistage, coupled thermal management simulation approach, informed by physical testing where available, to aid design decisions for PACCAR's SuperTruck II hybrid truck cabin concept. Focus areas include cabin insulation, battery sizing, and sleeper curtain position, as well as heating, ventilating, and air-conditioning (HVAC) component and accessory configurations, to maintain or improve thermal comfort while saving energy. The authors analyzed weather data and determined the national vehicle miles traveled weighted temperature and solar conditions for long-haul trucks. Example weather day profiles were selected to approximate the 5th and 95th percentile weighted conditions. A daylong drive cycle was developed to impose appropriate external wind conditions during rest and driving periods. Using the National Renewable Energy Laboratory's vehicle HVAC modeling and simulation tool VTCab, HVAC load design trade-off studies for the new truck geometry concept were completed. Parameters analyzed included effects of paint color, insulation, glass transmissivity, and curtain location. Simulation results helped with early design material selections for efficient cabin climate control. A detailed three-dimensional computer-aided engineering (CAE), computational fluid dynamics (CFD), radiation, and human physiology co-simulation, referred to in this article as CAE Thermal-CFD, was used to evaluate thermal comfort and energy impacts of diffuser configurations and air supply settings in driving and hoteling modes. Analysis revealed that it is more difficult to heat the cabin in hoteling mode during the winter than to cool the space in the summer. This seasonal load profile drives the requirement of additional energy storage for heating comfort. To determine the battery capacity requirement, multiday HVAC operation drive cycle simulations were then completed, showing that a 15-kWh battery would be required for HVAC operation during hoteling. Results helped reduce cabin thermal loads, determine component sizing requirements, and improve occupant comfort to save fuel and contribute to the economic viability of the hybrid system.

33 ADVANCED PROPULSION SYSTEMS↗

Optimizing radar scan strategies for tracking isolated deep convection using observing system simulation experiments

Abstract. Optimizing radar observation strategies is one of the most important considerations in pre-field campaign periods. This is especially true for isolated convective clouds that typically evolve faster than the observations captured by operational radar networks. This study investigates uncertainties in radar observations of the evolution of the microphysical and dynamical properties of isolated deep convective clouds developing in clean and polluted environments. It aims to optimize the radar observation strategy for deep convection through the use of high-spatiotemporal cloud-resolving model simulations, which resolve the evolution of individual convective cells every 1 min, coupled with a radar simulator and a cell tracking algorithm. The radar simulation settings are based on the Tracking Aerosol Convection Interactions ExpeRiment (TRACER) and Experiment of Sea Breeze Convection, Aerosols, Precipitation and Environment (ESCAPE) field campaigns held in the Houston, TX, area but are generalizable to other field campaigns focusing on isolated deep convection. Our analysis produces the following four outcomes. First, a 5–7 m s−1 median difference in maximum updrafts of tracked cells is shown between the clean and polluted simulations in the early stages of the cloud lifetimes. This demonstrates the importance of obtaining accurate estimates of vertical velocity from observations if aerosol impacts are to be properly resolved. Second, tracking of individual cells and using vertical cross section scanning every minute capture the evolution of precipitation particle number concentration and size represented by polarimetric observables better than the operational radar observations that update the volume scan every 5 min. This approach also improves multi-Doppler radar updraft retrievals above 5 km above ground level for regions with updraft velocities greater than 10 m s−1. Third, we propose an optimized strategy composed of cell tracking by quick (1–2 min) vertical cross section scans from more than one radar in addition to the operational volume scans. We also propose the use of a single-RHI (range height indicator) updraft retrieval technique for cells close to the radars, for which multi-Doppler radar retrievals are still challenging. Finally, increasing the number of deep convective cells sampled by such observations better represents the median maximum updraft evolution with sample sizes of more than 10 deep cells, which decreases the error associated with sampling the true population to less than 3 m s−1.

54 ENVIRONMENTAL SCIENCES↗

Microwave resonance thermomagnetic analysis - A new method for characterizing fine-grained ferromagnetic constituents in lunar materials

Microwave resonance thermomagnetic analysis (MRTA) is the name given to a newly evolved technique for inferring the natures of fine-grained ferromagnetic constituents in lunar materials. Based on standard ferromagnetic resonance (FMR) procedures, the method makes use of the microwave skin effect for diagnosing the presence of metallic iron. Modelling experiments carried out on well-characterized iron and magnetitelike precipitates produced independently in simulated lunar glasses, coupled with published data for magnetite, provide a potential basis for detecting and discriminating between iron metal and ferric iron spinel, even when both are present in an unknown sample. Application of the technique to the lunar samples indicates the possible existence of magnetitelike phases in amounts up to about 0.3 wt% in soils from seven samples regions of the moon. These findings do not require any special geologic processes for their explanation, although some evidence supports the suggestion that fumarolic activity may have occurred in the lunar highlands.

Griscom, D. L.↗

Effect of combined roll rate and sideslip angle on aircraft flight stability

The combined destabilizing effects of roll rate and sideslip angle are studied for a high angle-of-attack flight condition using a linear simulation of the coupled motions of a Space Shuttle-like configuration. Variations in the eigenvalues for a flight condition nominally identified by an angle of attack of 33.29 deg and a Mach number of 4.9 are examined with illustrations of the separate effects of body-axis roll rate and yaw rate, stability-axis roll rate, dynamic pressure, and sideslip angle. It is shown that the stability of the free motion of the vehicle is sensitive to roll rate and that this sensitivity is enhanced by nonzero sideslip angle. The results suggest that fully coupled linearized equations can be of value both for the study of Space Shuttle stability and for a better understanding of post-stall gyrations, incipient spin, and departure prevention for high-performance aircraft.

Stengel, R. F.↗

Recovery of spinning satellites

The behavior of a space tug and a spinning satellite in a coupled configuration was simulated and analyzed. A docking concept was developed to investigate the requirements pertaining to the design of a docking interface. Sensing techniques and control requirements for the chase vehicle were studied to assess the feasibility of an automatic docking. The effects of nutation dampers and liquid propellant slosh motion upon the docking transient were investigated.

Coppey, J. M.↗

On the supposed anticorrelation of solar polar and equatorial rotation rates

The anticorrelation between two of the three parameters used to calculate the angular velocity of the sun's differential rotation is here said to be due to numerical coupling. A computer simulation technique shows that the relationship between the two parameters is caused by the effect of noise on the least-squares analysis used to obtain the three parameters used to determine the angular velocity in terms of the heliographic latitude. The computer simulation technique is described. The supposed anticorrelation had been used to infer that variations of the sun's polar and equatorial rotation rates are anticorrelated.

Duvall, T. L., Jr.↗

Effects of magnification and visual accommodation on aimpoint estimation in simulated landings with real and virtual image displays

Twenty professional pilots observed a computer-generated airport scene during simulated autopilot-coupled night landing approaches and at two points (20 sec and 10 sec before touchdown) judged whether the airplane would undershoot or overshoot the aimpoint. Visual accommodation was continuously measured using an automatic infrared optometer. Experimental variables included approach slope angle, display magnification, visual focus demand (using ophthalmic lenses), and presentation of the display as either a real (direct view) or a virtual (collimated) image. Aimpoint judgments shifted predictably with actual approach slope and display magnification. Both pilot judgments and measured accommodation interacted with focus demand with real-image displays but not with virtual-image displays. With either type of display, measured accommodation lagged far behind focus demand and was reliably less responsive to the virtual images. Pilot judgments shifted dramatically from an overwhelming perceived-overshoot bias 20 sec before touchdown to a reliable undershoot bias 10 sec later.

Randle, R. J.↗

Jet production in super-Eddington accretion disks

A two-dimensional, radiation-coupled, Newtonian hydrodynamic simulation is reported for a super-Eddington, mass accretion rate, M = 4 M(E) disk accretion flow onto a 3-solar mass pseudoblack hole. Near the disk midplane, convection cells effectively block the accretion flow, even though viscous heating maximizes there. Accretion predominantly occurs in a supersonic inflow which follows streamlines of approximately constant angular momentum. The optically thick inflow traps radiation so that 80 percent of the luminosity is absorbed by the black hole; the emergent power is sub-Eddington. An axial jet self consistently forms just outside a conical photosphere which bounds the accretion zone; radiation pressure accelerates the jet to about 10 to the 10th cm/s. The jet's mass efflux is only 0.4 percent of the total mass accretion rate.

Eggum, G. E.↗

Analysis of coaxial spray combustion flames and related numerical issues

An approach to the simulation of strongly coupled multiphase flows in combustion hardware is sketched and its unique requirements highlighted. An example of a successful application to a coaxial injector flame is presented. Furthermore, several numerical issues that tend to interact with the physics of the problem are discussed with special regard to their potential impact on the choices of numerical parameters by the analyst. These include the issues of stability, numerical diffusivity, stiffness, and boundary conditions. The theme of this paper focuses on the intriguing relationships among the grid, the solution algorithm, and the actual physical mechanisms themselves.

Liang, P. Y.↗

Surface brightness profiles of 10 comets

CCD photometric observations of the comae of 10 comets, obtained at the 4-m and 2.1-m telescopes at KPNO during 1985-1986 using filters centered at 700.5, 650.0, or 546.0 nm, are reported. The data are presented in extensive tables and graphs and characterized in detail. The radial surface brightness profiles are shown to be steeper than predicted by an idealized spherically symmetric steady-state comet model, the steepness increasing with the projected distance from the nucleus. These profiles are attributed, on the basis of Monte Carlo simulations, to imperfect coupling between the sublimated gas and the optically dominant grains of the coma.

Jewitt, D. C.↗

Adaptive remeshing for transient problems with moving bodies

The combination of adaptive remeshing techniques, flow solvers for transient problems with moving grids, and integrators for rigid body motion is presented. The resulting scheme allows the economic simulation of fully coupled fluid-rigid body interaction problems of arbitrary geometric complexity. Several results are given to demonstrate the capabilities developed.

Lohner, Rainald↗

An adaptive finite element solver for transient problems with moving bodies

The combination of adaptive remeshing techniques, flow solvers for transient problems with moving grids, and integrators for rigid body motion is presented. The resulting scheme allows the economic simulation of fully coupled fluid-rigid body interaction problems of arbitrary geometric complexity. Several results are given to demonstrate the capabilities developed.

Lohner, Rainald↗