Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic model reduction”

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,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Cross-Linker Selection Controls Glass Transition Elevation or Reduction in Dynamic Covalently Bonded Polymer Networks

Introducing cross-links is a powerful approach to improve polymeric material performance relevant to controlling viscoelasticity, thermal and creep resistance, degradability, and efficient membrane separations. The chemically specific glass transition temperature T g is of fundamental importance in determining the time scales of key dynamical processes and physical state of the material in such applications. Here, we study experimentally how the introduction of relatively large cross-linking molecules in slowly exchanging dynamic bond-forming polymers (vitrimers) impacts vitrification for diverse polymer chemistries and a wide range of cross-link fractions. We find T g can increase, decrease, or even remain essentially unchanged, in qualitative contrast to the generic elevation of T g in traditional permanent polymer networks. We formulate an effective terpolymer network model to understand this rich behavior, which emerges as a consequence of a competition between pure cross-linking and generalized plasticization effects. The latter is associated with the tunable cross-linker size and intrinsic dynamic mobility that can offset slowing down due to traditional permanent cross-linking constraints. Here, a new strategy for functional polymer network design is suggested based on adjusting the relative importance of the two competing physical effects, which potentially can significantly enhance energy savings in applications while retaining other intrinsic properties germane to advanced materials performance.

Copolymers↗
DOE OSTI2021

Separating Physics and Dynamics Grids for Improved Computational Efficiency in Spectral Element Earth System Models

Previous studies have shown that atmospheric models with a spectral element grid can benefit from putting physics calculations on a relatively coarse finite volume grid. Here we demonstrate an alternative high-order, element-based mapping approach used to implement a quasi-equal-area, finite volume physics grid in E3SM. Unlike similar methods, the new method in E3SM requires topology data purely local to each spectral element, which trivially allows for regional mesh refinement. Simulations with physics grids defined by 2 × 2, 3 × 3, and 4 × 4 divisions of each element are shown to verify that the alternative physics grid does not qualitatively alter the model solution. The model performance is substantially affected by the reduction of physics columns when using the 2 × 2 grid, which can increase the throughput of physics calculations by roughly 60%–120% depending on whether the computational resources are configured to maximize throughput or efficiency. A pair of regionally refined cases are also shown to highlight the refinement capability.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2021

Impact of Regional Northern Hemisphere Mid-latitude Anthropogenic Sulfur Dioxide Emissions on Local and Remote Tropospheric Oxidants

The unintended consequences of reductions in regional anthropogenic sulfur dioxide (SO2) emissions implemented to protect human health are poorly understood. SO2 decreases began in the 1970s in the US and Europe and are expected to continue into the future, while recent emissions decreases in China are also projected to continue. In addition to the well-documented climate effects (warming) from reducing aerosols, tropospheric oxidation is impacted via aerosol modification of photolysis rates and radical sinks. Impacts on the hydroxyl radical and other trace constituents directly affect climate and air quality metrics such as surface ozone levels. We use the Geophysical Fluid Dynamics Laboratory Atmospheric Model version 3 nudged towards National Centers for Environmental Prediction (NCEP) reanalysis wind velocities to estimate the impact of SO2 emissions from the US, Europe, and China by differencing a control simulation with an otherwise identical simulation in which 2015 anthropogenic SO2 emissions are set to zero over one of the regions. Springtime sulfate aerosol changes occur both locally to the emission region and also throughout the Northern Hemispheric troposphere, including remote oceanic regions and the Arctic. Hydroperoxy (HO2) radicals are directly removed via heterogeneous chemistry on aerosol surfaces, including sulfate, in the model, and we find that sulfate aerosol produced by SO2 emissions from the three individual northern mid-latitude regions strongly reduces both HO2 and hydroxyl (OH) by up to 10 % year-round throughout most of the troposphere north of 30∘ N latitude. Regional SO2 emissions significantly increase nitrogen oxides (NOx) by about 5 %–8 % throughout most of the free troposphere in the Northern Hemisphere by increasing the NOx lifetime as the heterogeneous sink of HO2 on sulfate aerosol declines. Despite the NOx increases, tropospheric ozone decreases at northern mid-latitudes by 1 %–4 % zonally averaged and by up to 5 ppbv in summertime surface air over China, where the decreases in HO2 and OH suppress O3 production. Since 2015 anthropogenic SO2 emissions in China exceed those in the US or Europe, the oxidative response is greatest for the China perturbation simulation. Chemical effects of aerosols on oxidation (reactive uptake) dominate over radiative effects (photolysis rates), the latter of which are only statistically significant locally for the large perturbation over China. We find that the SO2 emissions decrease in China, which has yet to be fully realized, will have the largest impact on oxidants and related species in the Northern Hemisphere free troposphere compared to future decreases in Europe or the US. Our results bolster previous calls for a multipollutant strategy for air pollution mitigation to avoid the unintended consequence of aerosol removal leading to surface ozone increases that offset or mask surface ozone gains achieved by regulation of other pollutants, especially in countries where current usage of high-sulfur emitting fuels may be phased out in the future.

anthropogenic sulfur dioxide emissions↗
NASA NTRSDate not supplied

Isolator Shock Train Response to Dual-Mode Scramjet Throttling-Like Forcing

Typical dual-mode scramjet designs employ an isolator component situated between the inlet and combustor. The isolator fluidically buffers the inlet from the combustor via the formation of a shock train, an interconnected series of shock-boundary layer interactions. The isolator is required during ramjet-engine operation and at the lower edges of the scramjet operational envelope. In an operational engine, transient behavior imposes fluidic dynamics due to attitude adjustments, combustor throttling, aberrant resonances within the combustor, etc. To accommodate these transients, the shock train adjusts on both temporal and spatial scales. Sufficiently large downstream pressure increases in the combustor play a critical role during unstart, whereby strong compression waves propagate upstream in the isolator to ultimately force the shock train out the inlet, if it cannot stabilize to the pressure increase. This critically reduces the engine mass flow rate, potentially causing a catastrophic reduction in thrust and loss of the vehicle. A nearly instantaneous unstart can, in principle, be predicted by a ratio of the normal shock pressure jump in an inviscid isolator. However, lower-magnitude perturbations, below this bound, also induce unstarts. Understanding the transient behavior inherent in unsteady isolator operation is critical to reducing safety margin in isolator length and refining engine control algorithms to improve engine performance. The Isolator Dynamics Research Lab (IDRL) is a cold-flow, direct-connect isolator facility ideal for assessing the fundamental fluid dynamics of shock trains. The internal flow of a scramjet isolator is modeled using a 2D, symmetric Mach 2.5 nozzle followed by a 695 mm long, 25.4x50.8 mm rectangular isolator test section. A blockage cone on a high-speed ball screw actuator downstream of the isolator is pulsed to simulate the throttling behavior in an operational combustor. Wall flush high-speed static pressure transducers detect the streamwise shock train within the isolator. Time delays between the onset of the throttling and the response motion of the shock train are analyzed to assess the controlling scales and parameters. From this analysis, the response time of the isolator shock train system is found to depend primarily on the isolator back pressure time-rate-of-change but quickly saturates to a maximum propagation speed along the isolator. These results inform dual-mode scramjet engine designers of the critical parameters and limits of an unstart resistant engine.

scramjet↗
NASA NTRS2008

Spiking Neurons for Analysis of Patterns

Artificial neural networks comprising spiking neurons of a novel type have been conceived as improved pattern-analysis and pattern-recognition computational systems. These neurons are represented by a mathematical model denoted the state-variable model (SVM), which among other things, exploits a computational parallelism inherent in spiking-neuron geometry. Networks of SVM neurons offer advantages of speed and computational efficiency, relative to traditional artificial neural networks. The SVM also overcomes some of the limitations of prior spiking-neuron models. There are numerous potential pattern-recognition, tracking, and data-reduction (data preprocessing) applications for these SVM neural networks on Earth and in exploration of remote planets. Spiking neurons imitate biological neurons more closely than do the neurons of traditional artificial neural networks. A spiking neuron includes a central cell body (soma) surrounded by a tree-like interconnection network (dendrites). Spiking neurons are so named because they generate trains of output pulses (spikes) in response to inputs received from sensors or from other neurons. They gain their speed advantage over traditional neural networks by using the timing of individual spikes for computation, whereas traditional artificial neurons use averages of activity levels over time. Moreover, spiking neurons use the delays inherent in dendritic processing in order to efficiently encode the information content of incoming signals. Because traditional artificial neurons fail to capture this encoding, they have less processing capability, and so it is necessary to use more gates when implementing traditional artificial neurons in electronic circuitry. Such higher-order functions as dynamic tasking are effected by use of pools (collections) of spiking neurons interconnected by spike-transmitting fibers. The SVM includes adaptive thresholds and submodels of transport of ions (in imitation of such transport in biological neurons). These features enable the neurons to adapt their responses to high-rate inputs from sensors, and to adapt their firing thresholds to mitigate noise or effects of potential sensor failure. The mathematical derivation of the SVM starts from a prior model, known in the art as the point soma model, which captures all of the salient properties of neuronal response while keeping the computational cost low. The point-soma latency time is modified to be an exponentially decaying function of the strength of the applied potential. Choosing computational efficiency over biological fidelity, the dendrites surrounding a neuron are represented by simplified compartmental submodels and there are no dendritic spines. Updates to the dendritic potential, calcium-ion concentrations and conductances, and potassium-ion conductances are done by use of equations similar to those of the point soma. Diffusion processes in dendrites are modeled by averaging among nearest-neighbor compartments. Inputs to each of the dendritic compartments come from sensors. Alternatively or in addition, when an affected neuron is part of a pool, inputs can come from other spiking neurons. At present, SVM neural networks are implemented by computational simulation, using algorithms that encode the SVM and its submodels. However, it should be possible to implement these neural networks in hardware: The differential equations for the dendritic and cellular processes in the SVM model of spiking neurons map to equivalent circuits that can be implemented directly in analog very-large-scale integrated (VLSI) circuits.

Huntsberger, Terrance↗
DOE OSTI2024

Controlling Ion Uptake in Carboxylated Mixed Conductors

Organic mixed ionic‐electronic conductors (OMIECs) have garnered significant attention due to their capacity to transport both ions and electrons, making them ideal for applications in energy storage, neuromorphics, and bioelectronics. However, charge compensation mechanisms during the polymer redox process remain poorly understood, and are often oversimplified as single‐ion injection with little attention to counterion effects. To advance understanding and design strategies toward next‐generation OMIEC systems, a series of p‐channel carboxylated mixed conductors is investigated. Varying side‐chain functionality, distinctive swelling character is uncovered during electrochemical doping/dedoping with model chao‐/kosmotropic electrolytes. Carboxylic acid functionalized polymers demonstrate strong deswelling and mass reduction during doping, indicating cation expulsion, while ethoxycarbonyl counterparts exhibit prominent mass increase, pointing to an anion‐driven doping mechanism. By employing operando grazing incidence X‐ray fluorescence (GIXRF), it is revealed that the carboxyl functionalized polymer engages in robust cation interaction, whereas ester functionalization shifts the mechanism towards no cation involvement. It is demonstrated that cations are pivotal in mitigating swelling by counterbalancing anions, enabling efficient anion uptake without compromising performance. These findings underscore the transformative influence of functionality‐driven factors and side‐chain chemistry in governing ion dynamics and conduction, providing new frameworks for designing OMIECs with enhanced performance and reduced swelling.

carboxylated polythiophenes↗
DOE OSTI2023

Exploring the Meteorological Impacts of Surface and Rooftop Heat Mitigation Strategies Over a Tropical City

Different heat mitigation technologies have been developed to improve the thermal environment in cities. However, the regional impacts of such technologies, especially in the context of a tropical city, remain unclear. The deployment of heat mitigation technologies at city-scale can change the radiation balance, advective flow, and energy balance between urban areas and the overlying atmosphere. We used the mesoscale Weather Research and Forecasting model coupled with a physically based single-layer urban canopy model to assess the impacts of five different heat mitigation technologies on surface energy balance, standard surface meteorological fields, and planetary boundary layer (PBL) dynamics for premonsoon typical hot summer days over a tropical coastal city in the month of April in 2018, 2019, and 2020. Results indicate that the regional impacts of cool materials (CMs), super-cool broadband radiative coolers, green roofs (GRs), vegetation fraction change, and a combination of CMs and GRs (i.e., “Cool city (CC)”) on the lower atmosphere are different at diurnal scale. Results showed that super-cool materials have the maximum potential of ambient temperature reduction of 1.6°C during peak hour (14:00 LT) compared to other technologies in the study. During the daytime hours, the PBL height was considerably lower than the reference scenario with no implementation of strategies by 700 m for super-cool materials and 500 m for both CMs and CC cases; however, the green roofing system underwent nominal changes over the urban area. During the nighttime hours, the PBL height increased by CMs and the CC strategies compared to the reference scenario, but minimal changes were evident for super-cool materials. The changes of temperature on the vertical profile of the heat mitigation implemented city reveal a stable PBL over the urban domain and a reduction of the vertical mixing associated with a pollution dome. This would lead to crossover phenomena above the PBL due to the decrease in vertical wind speed. Therefore, assessing the coupled regional impact of urban heat mitigation over the lower atmosphere at city-scale is urgent for sustainable urban planning.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2001

Summary Report of the Orbital X-34 Wing Static Aeroelastic Study

This report documents the results of a computational study conducted on the Orbital Sciences X-34 vehicle to compute its inviscid aerodynamic characteristics taking into account the wing structural flexibility. This was a joint exercise between LaRC and SDRC of California. SDRC modeled the structural details of the wing, and provided the structural deformation for a given pressure distribution on its surfaces. This study was done for a Mach number of 1.35 and an angle of attack of 9 deg.; the freestream dynamic pressure was assumed to be 607 lb/sq ft. Only the wing and the body were simulated in the CFD computations. Two wing configurations were examined. The first had the elevons in the undeflected position and the second had the elevons deflected 20 deg. up. The results indicated that with elevon undeflected, the wing twists by about 1.5 deg. resulting in a reduction in the angle of attack at the wing tip to by 1.5 deg. The maximum vertical deflection of the wing is about 3.71 inches at the wing tip. For the wing with the undeflected elevons, the effect of this wing deformation is to reduce the normal force coefficient (C(sub N)) by 0.012 and introduce a noise up pitching moment coefficient (C(sub m)) of 0.042.

Prabhn, Ramadas K.↗
DOE OSTI2022

Optimal Sizing and Operation of a Hybrid Clean Energy Center

In this paper, we consider the deployment of battery energy storage systems (BESSs) paired with photovoltaics (PV) to form a hybrid clean energy center (HCEC), where the BESS is equipped for two main purposes: i) reducing the volatility of PV output and ii) shifting energy to meet an export limit. We develop a linear dynamical model to characterize the output of an HCEC and express the volatility and export power as functions of the battery's dispatch operations. We also propose an optimal sizing method via second-order cone programming (SOCP) to maximize the net benefits of the HCEC. Numerical studies are performed to illustrate and validate the proposed method.

Battery energy storage systems, photovoltaics, sec↗
DOE OSTI2024

Modeling oxygen reduction activity loss mechanisms in atomically dispersed Fe–N–C electrocatalysts

Materials degradation is a major factor that limits the wider adoption of renewable and clean energy technologies. This is particularly true for the Pt group metal-free (PGM-free) atomically dispersed metal-nitrogen-carbon (M-N-C) catalysts. Here, while many experimental studies have investigated and reported the phenomenological aspects of M-N-C degradation, only a few modeling studies have considered degradation mechanisms at the atomic level. Understanding the mechanisms responsible for activity loss occurring in atomically dispersed M-N-C’s is crucial towards rationally designing active, durable, and less expensive Earth-abundant catalysts. Towards this end, we have surveyed recent literature concerning the modeling of corrosion mechanisms that impact M-N-C catalysts (Fe–N–C, in particular) and offer our own perspectives on the future direction of this field.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Scalable training of graph convolutional neural networks for fast and accurate predictions of HOMO-LUMO gap in molecules

Abstract Graph Convolutional Neural Network (GCNN) is a popular class of deep learning (DL) models in material science to predict material properties from the graph representation of molecular structures. Training an accurate and comprehensive GCNN surrogate for molecular design requires large-scale graph datasets and is usually a time-consuming process. Recent advances in GPUs and distributed computing open a path to reduce the computational cost for GCNN training effectively. However, efficient utilization of high performance computing (HPC) resources for training requires simultaneously optimizing large-scale data management and scalable stochastic batched optimization techniques. In this work, we focus on building GCNN models on HPC systems to predict material properties of millions of molecules. We use HydraGNN, our in-house library for large-scale GCNN training, leveraging distributed data parallelism in PyTorch. We use ADIOS, a high-performance data management framework for efficient storage and reading of large molecular graph data. We perform parallel training on two open-source large-scale graph datasets to build a GCNN predictor for an important quantum property known as the HOMO-LUMO gap. We measure the scalability, accuracy, and convergence of our approach on two DOE supercomputers: the Summit supercomputer at the Oak Ridge Leadership Computing Facility (OLCF) and the Perlmutter system at the National Energy Research Scientific Computing Center (NERSC). We present our experimental results with HydraGNN showing (i) reduction of data loading time up to 4.2 times compared with a conventional method and (ii) linear scaling performance for training up to 1024 GPUs on both Summit and Perlmutter.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Scale-up of microwave-assisted, continuous flow, liquid phase reactors: Application to 5-Hydroxymethylfurfural production

Microwave (MW) technology can be powerful for the electrification and intensification of chemical manufacturing. However, very few examples of scaled MW processes exist. In this work, we build a scaled, MW-assisted, continuous flow reactor for the processing of liquid phase chemistries with specific application to 5-hydroxymethylfurfural (HMF) production. Here, we demonstrate a corresponding computational fluid dynamics model to simulate the reactor’s temperature profile and performance, both of which are in good agreement with experiments. We construct a surrogate model to relate operating parameters to performance for guiding experiments, without demanding simulations, via active learning. Heat recirculation is demonstrated in this scaled reactor, further extending its scale and reducing the energy demand. An HMF yield of ~55 % and a productivity of 0.1 kg/hr, 8x higher than any other reactor, at a flowrate >20x than prior work, are demonstrated while maintaining energy efficiency of >98 %. A basic economic analysis estimates a $1.85/kg cost of HMF from fructose and >60 % reduction in CO 2 emissions than a conventional system without considering the impact of green electricity. We demonstrate that MW technology is well-applied to systems larger than the laboratory scale, and strategies to further reduction in cost and CO 2 are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Atomistic Origins of Reversible Noncatalytic Gas–Solid Interfacial Reactions

Noncatalytic gas-solid reactions are a large group of heterogeneous reactions that are usually assumed to occur irreversibly because of the strong driving force to favor the forward direction toward the product formation. Using the example of Ni oxidation into NiO with CO 2 , herein we demonstrate the existence of the reverse element that results in the NiO reduction from the countering effect of the gaseous product of CO. Using in-situ electron microscopy observations and atomistic modeling, we show that the oxidation process occurs via preferential CO 2 adsorption along step edges that results in step-flow growth of NiO layers, and the presence of Ni atoms on the flat NiO surface promotes the nucleation of NiO layers. Simultaneously, the NiO reduction happens via preferential step-edge adsorption of CO that leads to the receding motion of atomic steps, and the presence of Ni vacancies in the NiO surface facilitates the CO adsorption induced surface pitting. Based on the in-situ TEM observations, temperature and CO 2 pressure effect maps are constructed to illustrate the spatiotemporal dynamics of the competing NiO redox reactions. We report, these results demonstrate the rich gas-solid surface reaction dynamics induced by the co-existing forward and reverse reaction elements and have practical applicability in manipulating gas-solid reactions via controlling the gas environment or atomic structure of the solid surface to steer the reaction toward the desired direction.

25 ENERGY STORAGE↗
DOE OSTI2021

Structural insights into photoactivation of plant Cryptochrome-2

Cryptochromes (CRYs) are evolutionarily conserved photoreceptors that mediate various light-induced responses in bacteria, plants, and animals. Plant cryptochromes govern a variety of critical growth and developmental processes including seed germination, flowering time and entrainment of the circadian clock. CRY’s photocycle involves reduction of their flavin adenine dinucleotide (FAD)-bound chromophore, which is completely oxidized in the dark and semi to fully reduced in the light signaling-active state. Despite the progress in characterizing cryptochromes, important aspects of their photochemistry, regulation, and light-induced structural changes remain to be addressed. In this study, we determine the crystal structure of the photosensory domain of Arabidopsis CRY2 in a tetrameric active state. Systematic structure-based analyses of photo-activated and inactive plant CRYs elucidate distinct structural elements and critical residues that dynamically partake in photo-induced oligomerization. Our study offers an updated model of CRYs photoactivation mechanism as well as the mode of its regulation by interacting proteins.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Transformer Neural Networks with Spatiotemporal Attention for Predictive Control and Optimization of Industrial Processes

In the context of real-time optimization and model predictive control of industrial systems, machine learning, and neural networks represent cutting-edge tools that hold promise for enhancing dynamic modeling. This work presents a novel transformer neural network architecture for real-time optimization and model predictive control. This network design includes a modified attention mechanism inspired by positional embedding attention from vision transformers and task-specific modifications to the input-output structure of the transformer’s decoder stack. Experiments were conducted using data from a 450 MW coal-fired power plant to evaluate this approach's effectiveness. The transformer neural network was compared with conventional recurrent models, including GRU and LSTM. The transformer exhibited a 6% increase in the R-squared (R2) value of predictions and an 83% reduction in mean squared error (MSE). Computation time was also reduced by 84% compared to conventional recurrent models.

Gallup, Ethan R.↗
DOE OSTI2020

Reduced Models for Chemical Kinetics derived from Parallel Ensemble Simulations of Stirred Reactors

Detailed chemical kinetic models are essential to improve the predictive capability of computational fluid dynamics (CFD) simulations and to impact the development of next generation combustion technologies. However, such detailed kinetic models are beyond the reach of CFD codes due to the computational cost of transporting additional chemical degrees of freedom. Reduced order models for combustion kinetics exploit the presence of lower dimensional manifolds in the state space and require far fewer number of scalars to be transported as part of CFD simulations. In this paper, we describe the software infrastructure developed to allow large ensembles of stirred reactor configurations with detailed chemical kinetics. The canonical reactor simulations are incorporated in a workflow with sparse grid and data reduction methodologies aimed at development of reduced order models. We present results from analyzing the results of the stirred reactor calculations for syngas mixture using principal component analysis (PCA) and neural network methodologies. The results show that nonlinear autoencoder model can surpass PCA in representing the original dataset with far fewer degrees of freedom.

Zhang, Pei↗
DOE OSTI2021

Material flows and GHG emissions from housing stock evolution in US counties, 2020–60

The evolution of housing stocks determines demand for construction materials and energy, and associated emissions of greenhouse gasses (GHGs). The contribution of construction to building life-cycle emissions is growing as buildings become more energy efficient and the energy supply decarbonizes. A housing stock model is developed for counties in the United States using dynamic vacancy rates which endogenously influence stock out- and inflows. Stocks of three house types and 10 construction cohorts are projected for all contiguous US counties for the period 2020–60. Inflows and outflows of construction materials are then estimated along with GHG emissions associated with material production and construction activities in scenarios defined by stock turnover rates, population share by house type, and floor area characteristics of new houses. The results provide new insights into the drivers of construction-related emissions at local and national levels, and identify opportunities for their reduction. Demolition material flows grow in relation to construction material flows over the analysis period. Increasing the stock turnover rate increases future floor area per person, material requirements, and emissions from construction. Scenarios with reduced floor area and more multifamily homes in new construction have lower floor area growth, material requirements, and emissions from construction.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS1988

The use of multigrid techniques in the solution of the Elrod algorithm for a dynamically loaded journal bearing

A numerical solution to a theoretical model of vapor cavitation in a dynamically loaded journal bearing is developed, utilizing a multigrid iterative technique. The code is compared with a presently existing direct solution in terms of computational time and accuracy. The model is based on the Elrod algorithm, a control volume approach to the Reynolds equation which mimics the Jakobssen-Floberg and Olsson cavitation theory. Besides accounting for a moving cavitation boundary and conservation of mass at the boundary, it also conserves mass within the cavitated region via liquid striations. The mixed nature of the equations (elliptic in the full film zone and nonelliptic in the cavitated zone) coupled with the dynamic aspects of the problem create interesting difficulties for the present solution approach. Emphasis is placed on the methods found to eliminate solution instabilities. Excellent results are obtained for both accuracy and reduction of computational time.

Woods, Claudia M.↗