Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical continuation”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
NASA NTRS1961

Calculation of Flow Fields from Bow-Wave Profiles for the Downstream Region of Blunt-Nosed Circular Cylinders in Axial Hypersonic Flight

A method by which known bow-wave profiles may be analyzed to give the flow fields around blunt-nosed cylinders in axial hypersonic flow is presented. In the method, the assumption is made that the pressure distribution curve in a transverse plane is similar to that given by blast- wave theory. Numerical analysis based on the one-dimensional energy and continuity equations then leads to distributions of all the flow variables in the cross section, for either a perfect gas or a real gas. The entire flow field need not be solved. Attention can be confined to any desired station. The critical question is the validity of the above assumption. It is tested for the case of a hemisphere cylinder in flight at 20,000 ft/sec. The flow is analyzed for three stations along the cylindrical afterbody, and found to compare very closely with the results of an exact (inviscid) solution. The assumed form of the pressure distribution occurs at stations as close as 1.2 diameters to the body nose. However, it is suggested that the assumption may not apply this far forward in general, particularly when bodies of nonsmooth contour are considered.

Seiff, Alvin↗
NASA NTRS1999

In Situ Monitoring of Crystal Growth Using MEPHISTO, Mission STS 87-Program USMP-4: Experimental Analysis

This report summarizes the results of the In situ Monitoring of Crystal Growth Using MEPHISTO (Material por l'Etude des Phenomenes Interessant de la Solidification sur Terre et en Orbite) experiment on USMP-4. The report includes microstructural and compositional data obtained during the first year of the post flight analysis, as well as numerical simulation of the flight experiment. Additional analyses are being continued and will be reported in the near future. The experiments utilized MEPHISTO hardware to study the solidification and melting behavior of bismuth alloyed with 1 at% tin. The experiments involved repeated melting and solidification of three samples, each approximately 90 cm long and 6mm in diameter. Half of each sample also included a 2 mm. diameter growth capillary, to assist in the formation of single grain inside. One sample provided the Seebeck voltage generated during melting and freezing processes. Another one provided temperature data and Peltier pulsed demarcation of the interface shape for post flight analysis. The third sample provided resistance and velocity measurements, as well as additional thermal data. The third sample was also quenched at the end of the mission to preserve the interface composition for post flight determination. A total of more than 45cm of directionally solidified alloy were directionally solidified at the end of the flight for post mission structural and compositional characterization. Metallurgical analysis of the samples has shown that the interfacial kinetics play a key role in controlling the morphological stability of faceted alloys. Substantial differences were observed in the Seebeck signal between the ground-based experiments and the space-based experiments. The temperature gradient in the liquid for the ground-based experiments was also significantly lower than the temperature gradient in the liquid for the space-based experiments. Both of these observations indicate significant influence of liquid convection for the ground-based experiments.

Abbaschian, Reza↗
NASA NTRS1999

Bubble Dynamics on a Heated Surface

In this work, we study steady and oscillatory thermocapillary and natural convective flows generated by a bubble on a heated solid surface. The interaction between gas and vapor bubbles with the surrounding fluid is of interest for both space and ground-based processing. A combined numerical-experimental approach is adopted here. The temperature field is visualized using Mach-Zehnder and/or Wollaston Prism Interferometry and the flow field is observed by a laser sheet flow visualization technique. A finite element numerical model is developed which solves the transient two-dimensional continuity, momentum, and energy equations and includes the effects of temperature-dependent surface tension and bubble surface deformation. Below the critical Marangoni number, the steady state low-g and 1-g temperature and velocity fields predicted by the finite element model are in excellent agreement with both the visualization experiments in our laboratory and recently published experimental results in the literature. Above the critical Marangoni number, the model predicts an oscillatory flow which is also closely confirmed by experiments. It is shown that the dynamics of the oscillatory flow are directly controlled by the thermal and hydrodynamic interactions brought about by combined natural and thermocapillary convection. Therefore, as numerical simulations show, there are considerable differences between the 1-g and low-g temperature and flow fields at both low and high Marangoni numbers. This has serious implications for both materials processing and fluid management in space.

Kassemi, M.↗
NASA NTRS2008

Ocean Carbon and Biogeochemistry Scoping Workshop on Terrestrial and Coastal Carbon Fluxes in the Gulf of Mexico, St. Petersburg, FL

Despite their relatively small surface area, ocean margins may have a significant impact on global biogeochemical cycles and, potentially, the global air-sea fluxes of carbon dioxide. Margins are characterized by intense geochemical and biological processing of carbon and other elements and exchange large amounts of matter and energy with the open ocean. The area-specific rates of productivity, biogeochemical cycling, and organic/inorganic matter sequestration are high in coastal margins, with as much as half of the global integrated new production occurring over the continental shelves and slopes (Walsh, 1991; Doney and Hood, 2002; Jahnke, in press). However, the current lack of knowledge and understanding of biogeochemical processes occurring at the ocean margins has left them largely ignored in most of the previous global assessments of the oceanic carbon cycle (Doney and Hood, 2002). A major source of North American and global uncertainty is the Gulf of Mexico, a large semi-enclosed subtropical basin bordered by the United States, Mexico, and Cuba. Like many of the marginal oceans worldwide, the Gulf of Mexico remains largely unsampled and poorly characterized in terms of its air-sea exchange of carbon dioxide and other carbon fluxes. The goal of the workshop was to bring together researchers from multiple disciplines studying terrestrial, aquatic, and marine ecosystems to discuss the state of knowledge in carbon fluxes in the Gulf of Mexico, data gaps, and overarching questions in the Gulf of Mexico system. The discussions at the workshop were intended to stimulate integrated studies of marine and terrestrial biogeochemical cycles and associated ecosystems that will help to establish the role of the Gulf of Mexico in the carbon cycle and how it might evolve in the face of environmental change.

Air sea exchanges↗
NASA NTRS2011

Impacts of Aerosol Direct Effects on the South Asian Climate: Assessment of Radiative Feedback Processes Using Model Simulations and Satellite/Surface Measurements

Current assessment of aerosol radiative effect is hindered by our incomplete knowledge of aerosol optical properties, especially absorption, and our current inability to quantify physical and microphysical processes. In this research, we investigate direct aerosol radiative effect over heavy aerosol loading areas (e.g., Indo-Gangetic Plains, South/East Asia) and its feedbacks on the South Asian climate during the pre-monsoon season (March-June) using the Purdue Regional Climate Model (PRCM) with prescribed aerosol data derived by the NASA Goddard Earth Observing System Model (GEOS-5). Our modeling domain covers South and East Asia (60-140E and 0-50N) with spatial resolutions of 45 km in horizontal and 28 layers in vertical. The model is integrated from 15 February to 30 June 2008 continuously without nudging (i.e., only forced by initial/boundary conditions). Two numerical experiments are conducted with and without the aerosol-radiation effects. Both simulations are successful in reproducing the synoptic patterns on seasonal-to-interannual time scales and capturing a pre-monsoon feature of the northward rainfall propagation over Indian region in early June which shown in Tropical Rainfall Measuring Mission (TRMM) observation. Preliminary result suggests aerosol-radiation interactions mainly alter surface-atmosphere energetics and further result in an adjustment of the vertical temperature distribution in lower atmosphere (below 700 hPa). The modifications of temperature and associated rainfall and circulation feedbacks on the regional climate will be discussed in the presentation.

Wang, Sheng-Hsiang↗
DOE OSTI2023

Tools and Methods for Optimization of Nuclear Plant Outages

Refueling outages are one of the most challenging phases in a nuclear power plants (NPP) operating cycle. Refueling outages are extremely costly for a NPP due to the large amount of required resources and because of lost revenue due to plant being off the grid. Outage durations have steadily decreased across the industry over that last few decades primarily due to improved planning and coordination, but there are still many plants that struggle to meet the performance metrics accomplished by other utilities. Schedule resilience is one of the issues. NPP outages require scheduling thousands of activities in a duration of around 30 days on average. Outage staff begin working on the schedule more than a year ahead of the outage start and make every effort to build a robust schedule. Despite the robust and detailed planning, once the outage starts, numerous emergent issues typically appear along with schedule delays requiring continuous replanning and adjusting. When schedule disruption occurs during an outage, plant staff make urgent efforts to recover but often not able to maintain the planned outage duration. These outage delays can cost a utility several million dollars per day. Tools that could help outage schedulers create a more resilient schedule and allow them to optimally reschedule emergent work could significantly reduce outage delays. One key aspect of creating a resilient schedule is to have accurate estimates for activity duration. Another important outage scheduling capability is the ability to schedule emergent work with minimal disruption. The Optimization of Outage Activities project under the Risk-Informed Systems Analysis Pathway (RISA) sponsored by Department of Energy (DOE) Light Water Reactor Sustainability (LWRS) program focuses on developing tools and methods to support nuclear power plants with optimization of outage schedules. The goal of the outage optimization is the completion of all planned and emergent outage activities as fast as possible while maintaining highest level of safety. This report describes the initial development of tools to support outage management that leverage computational and machine learning methods developed in other RISA and LWRS projects.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Methods and Tools To Assess Robustness of Nuclear Plant Outages

Refueling outages are one of the most challenging phases in a nuclear power plant (NPP) operating cycle. Refueling outages are extremely costly for an NPP due to the large amount of required resources and because of lost revenue due to plant being off the grid. Outage durations have steadily decreased across the industry over that last few decades primarily due to improved planning and coordination, but there are still many plants that struggle to meet the performance metrics accomplished by other utilities. Schedule resilience is one of the issues. NPP outages require scheduling thousands of activities within 30 days on average. Despite detailed planning, once the outage starts, numerous emergent issues typically appear along with schedule delays requiring continuous replanning and adjustment. When schedule disruption occurs during an outage, plant staff make urgent efforts to recover but are often not able to maintain the planned outage duration. These outage delays can cost a utility several million dollars per day. Tools that could help outage schedulers create a more resilient schedule and allow them to optimally reschedule emergent work could significantly reduce outage delays. One key aspect of creating a resilient schedule is to have accurate estimates for activity duration. Another important outage scheduling capability is the ability to schedule emergent work with minimal disruption. This paper focuses on developing tools and methods to support NPPs with outage schedule optimization and it describes the initial development of tools to support outage management that leverage computational and machine learning methods.

97 - MATHEMATICS AND COMPUTING↗
NASA NTRS1995

The method of space-time conservation element and solution element-applications to one-dimensional and two-dimensional time-marching flow problems

A nontraditional numerical method for solving conservation laws is being developed. The new method is designed from a physicist's perspective, i.e., its development is based more on physics than numerics. Even though it uses only the simplest approximation techniques, a 2D time-marching Euler solver developed recently using the new method is capable of generating nearly perfect solutions for a 2D shock reflection problem used by Helen Yee and others. Moreover, a recent application of this solver to computational aeroacoustics (CAA) problems reveals that: (1) accuracy of its results is comparable to that of a 6th order compact difference scheme even though nominally the current solver is only of 2nd-order accuracy; (2) generally, the non-reflecting boundary condition can be implemented in a simple way without involving characteristic variables; and (3) most importantly, the current solver is capable of handling both continuous and discontinuous flows very well and thus provides a unique numerical tool for solving those flow problems where the interactions between sound waves and shocks are important, such as the noise field around a supersonic over- or under-expansion jet.

Chang, Sin-Chung↗
NASA NTRS1990

Multi-blocking strategies for the INS3D incompressible Navier-Stokes code

With the continuing development of bigger and faster supercomputers, computational fluid dynamics (CFD) has become a useful tool for real-world engineering design and analysis. However, the number of grid points necessary to resolve realistic flow fields numerically can easily exceed the memory capacity of available computers. In addition, geometric shapes of flow fields, such as those in the Space Shuttle Main Engine (SSME) power head, may be impossible to fill with continuous grids upon which to obtain numerical solutions to the equations of fluid motion. The solution to this dilemma is simply to decompose the computational domain into subblocks of manageable size. Computer codes that are single-block by construction can be modified to handle multiple blocks, but ad-hoc changes in the FORTRAN have to be made for each geometry treated. For engineering design and analysis, what is needed is generalization so that the blocking arrangement can be specified by the user. INS3D is a computer program for the solution of steady, incompressible flow problems. It is used frequently to solve engineering problems in the CFD Branch at Marshall Space Flight Center. INS3D uses an implicit solution algorithm and the concept of artificial compressibility to provide the necessary coupling between the pressure field and the velocity field. The development of generalized multi-block capability in INS3D is described.

Gatlin, Boyd↗
NASA NTRS1949

Response of a Rotating Propeller to Aerodynamic Excitation

The flexural vibration of a rotating propeller blade with clamped shank is analyzed with the object of presenting, in matrix form, equations for the elastic bending moments in forced vibration resulting from aerodynamic forces applied at a fixed multiple of rotational speed. Matrix equations are also derived which define the critical speeds end mode shapes for any excitation order and the relation between critical speed and blade angle. Reference is given to standard works on the numerical solution of matrix equations of the forms derived. The use of a segmented blade as an approximation to a continuous blade provides a simple means for obtaining the matrix solution from the integral equation of equilibrium, so that, in the numerical application of the method presented, the several matrix arrays of the basic physical characteristics of the propeller blade are of simple form, end their simplicity is preserved until, with the solution in sight, numerical manipulations well-known in matrix algebra yield the desired critical speeds and mode shapes frame which the vibration at any operating condition may be synthesized. A close correspondence between the familiar Stodola method and the matrix method is pointed out, indicating that any features of novelty are characteristic not of the analytical procedure but only of the abbreviation, condensation, and efficient organization of the numerical procedure made possible by the use of classical matrix theory.

Arnoldi, Walter E.↗
NASA NTRS1978

A composition and thickness model for lunar impact crater and basin deposits

A model was developed for determining the composition and thicknesses of continuous deposits of large lunar impact craters and basins. Results of a photogeologic study and topographic analysis of continuous deposits of the lunar crater Delisle, which support the model, show that numerous secondary craters and crater chains with concentric dunes on their uprange rims occur well within the mapped deposits. At any given radius from Delisle, the secondary craters are of equal freshness and not vastly different in size, but with increasing radius they become better defined. The upper surface of the continuous deposits reaches the level of the pre-existing mare surface well within the mapped boundary of the deposits.

Morrison, R. H.↗
NASA NTRS1998

Studying Turbulence Using Numerical Simulation Databases: Proceedings of the Summer Program - No. 7

The Seventh Summer Program of the Center for Turbulence Research took place in the four-week period, July 5 to July 31, 1998. This was the largest CTR Summer Program to date, involving thirty-six participants from the U. S. and nine other countries. Thirty-one Stanford and NASA-Ames staff members facilitated and contributed to most of the Summer projects. A new feature, and perhaps a preview of the future programs, was that many of the projects were executed on non-NASA computers. These included supercomputers located in Europe as well as those operated by the Departments of Defense and Energy in the United States. In addition, several simulation programs developed by the visiting participants at their home institutions were used. Another new feature was the prevalence of lap-top personal computers which were used by several participants to carry out some of the work that in the past were performed on desk-top workstations. We expect these trends to continue as computing power is enhanced and as more researchers (many of whom CTR alumni) use numerical simulations to study turbulent flows. CTR's main role continues to be in providing a forum for the study of turbulence for engineering analysis and in facilitating intellectual exchange among the leading researchers in the field. Once again the combustion group was the largest. Turbulent combustion has enjoyed remarkable progress in using simulations to address increasingly complex and practically more relevant questions. The combustion group's studies included such challenging topics as fuel evaporation, soot chemistry, and thermonuclear reactions. The latter study was one of three projects related to the Department of Energy's ASCI Program (www.llnl.gov/asci); the other two (rocket propulsion and fire safety) were carried out in the turbulence modeling group. The flow control and acoustics group demonstrated a successful application of the so-called evolution algorithms which actually led to a previously unknown forcing strategy for jets yielding increased spreading rate. A very efficient algorithm for flow in complex geometries with moving boundaries based on the immersed boundary forcing technique was tested with very encouraging results. Also a new strategy for the destruction of aircraft trailing vortices was introduced and tested. The Reynolds Averaged Modeling (RANS) group demonstrated that the elliptic relaxation concept for RANS calculations is also applicable to transonic flows with shocks; however, prediction of laminar/turbulent transition remains an important pacing item. A large fraction of the LES effort was devoted to the development and testing of a new algorithmic procedure (as opposed to phenomenological model) for subgrid scale modeling based on regularized de-filtering of the flow variables. This appears to be a very promising approach, and a significant effort is currently underway to assess its robustness in high Reynolds number flows and in conjunction with numerical methods for complex flows. As part of the Summer Program two review tutorials were given on Turbulent structures in hydrocarbon pool fires (Sheldon Tieszen), and Turbulent combustion modeling: from RANS to LES via DNS (Luc Vervisch); and two seminars entitled Assessment of turbulence models for engineering applications (Paul Durbin) and Subgrid-scale modeling for non-premixed, turbulent reacting flows (James Riley) were presented. A number of colleagues from universities, government agencies, and industry attended the final presentations of the participants on July 31 and participated in the discussions. There are twenty-six papers in this volume grouped in five areas. Each group is preceded with an overview by its coordinator.

Source record↗
DOE OSTI2023

Demonstrating Computational Equivalence Between Continuous and Discrete Adjoint Methods by Calculating Time-Dependent Adjoint Solutions with Neutron Diffusion Models

The continuous adjoint method and the discrete adjoint method are two alternative approaches used to calculate adjoint solutions for adjoint systems. The continuous adjoint method derives adjoint equations analytically from continuous forward equations and then solves the adjoint equations either analytically or numerically in a discretized form whereas the discrete adjoint method calculates the adjoint solutions directly from the discretized forward equations. With regard to the methodology development and calculation procedure, distinct differences are well recognized between the two methods. For certain reasons, both methods are exclusively preferred and commonly used by different computational communities, but limited studies clarify the connections between the two adjoint methods from either of the communities. Herein, this paper demonstrates the computational equivalence between the continuous and discrete adjoint methods by investigating time-dependent adjoint solutions to the two-group neutron diffusion model in nuclear reactor analysis problems using both methods. Adjoint solutions can be used to estimate system parameters for reactor safety analysis. Appropriate final state conditions for the adjoint systems are specified in both of the methods, and the conditions are clarified with proper physical explanations. With the help of an event-based case study on neutron diffusion models, the accuracy of the time-dependent adjoint fluxes obtained from both methods is verified, and the pros and cons of both adjoint methods are examined. More importantly, the computational equivalence of both methods is demonstrated when they are applied to multigroup neutron diffusion systems. The advantage of calculating time-dependent adjoint fluxes by directly solving time-dependent adjoint systems rather than taking steady-state approximations as in common practice is also demonstrated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Multi-scale Coupled Processes Modeling of Fractures as Porous, Interfacial and Granular Systems from Rock Images with the Numerical Manifold Method

Abstract The greatest challenges of rigorously modeling coupled hydro-mechanical processes in fractured rocks at different scales are associated with computational geometry. In addition, selections of continuous or discontinuous models, physical laws, and coupling priorities at different scales based on different geometric features determine the applicability of a numerical model for a certain type of problem. In this study, we present our multi-scale modeling capabilities that have been developed based on the numerical manifold method for analyzing coupled hydro-mechanical processes in fractured rocks. Based on their geometric features, the fractures are modeled as continua—finite-thickness porous zones, and discontinua—discontinuous interfaces and microscale asperities and granular systems. Different governing equations, physical laws, coupling priorities, and approaches for addressing fracture intersections and shearing are then applied to describe these. We applied these models to simulate coupled processes in fractured rocks using realistic geometry obtained from rock images at different scales. We first calculated shearing of a single fracture with different models and demonstrated the impacts of asperities on shearing. We then applied the continuous and discontinuous models to simulate a network of rough fractures, demonstrating that contact dynamics contribute significantly to the geometric, multi-physical evolution of systems where rough fractures are not mineral filled. For a discrete fracture network, our coupled processes modeling demonstrates that shearing of the discrete fractures can have a major impact on stress and pore pressure distribution. Lastly, we applied the discontinuous granular model to simulate evolution of a complex granular system with a deformation band, demonstrating that the deformation band can dominate contact dynamics, the structural and the stress evolution of the granular system.

58 GEOSCIENCES↗
DOE OSTI2022

Magnon heat transport in a two-dimensional Mott insulator

Whether or not anomalies in the thermal conductivity in insulating cuprates can be attributed to antiferromagnetic order and magnons in a 2D Mott insulator remains an intriguing open question. To shed light on this issue, we investigate the thermal conductivity $\kappa$ and its relationship with the specific heat $c_v$ in the half-filled 2D single-band Hubbard model, using the numerically exact determinant quantum Monte Carlo algorithm and maximum entropy analytic continuation. At low temperatures where the charge degrees of freedom are gapped-out and $c_v$ exhibits a clear magnon peak, we observe that thermal conductivity $\kappa$ also tends to form a peak at similar temperatures. Reducing temperature further produces a sharp upturn in $\kappa$, associated with an increasing mean-free path. Here we identify this as the high-temperature side of the anomalous peak in insulating cuprates, where the mean-free path eventually is cut-off by other scattering effects, including phonons, disorder, and physical size. Different scattering effects in our model are identified and analyzed in the thermal diffusivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1989

Parallel and vector computation for stochastic optimal control applications

A general method for parallel and vector numerical solutions of stochastic dynamic programming problems is described for optimal control of general nonlinear, continuous time, multibody dynamical systems, perturbed by Poisson as well as Gaussian random white noise. Possible applications include lumped flight dynamics models for uncertain environments, such as large scale and background random atmospheric fluctuations. The numerical formulation is highly suitable for a vector multiprocessor or vectorizing supercomputer, and results exhibit high processor efficiency and numerical stability. Advanced computing techniques, data structures, and hardware help alleviate Bellman's curse of dimensionality in dynamic programming computations.

Hanson, F. B.↗
NASA NTRS1994

Prediction of Undsteady Flows in Turbomachinery Using the Linearized Euler Equations on Deforming Grids

A linearized Euler solver for calculating unsteady flows in turbomachinery blade rows due to both incident gusts and blade motion is presented. The model accounts for blade loading, blade geometry, shock motion, and wake motion. Assuming that the unsteadiness in the flow is small relative to the nonlinear mean solution, the unsteady Euler equations can be linearized about the mean flow. This yields a set of linear variable coefficient equations that describe the small amplitude harmonic motion of the fluid. These linear equations are then discretized on a computational grid and solved using standard numerical techniques. For transonic flows, however, one must use a linear discretization which is a conservative linearization of the non-linear discretized Euler equations to ensure that shock impulse loads are accurately captured. Other important features of this analysis include a continuously deforming grid which eliminates extrapolation errors and hence, increases accuracy, and a new numerically exact, nonreflecting far-field boundary condition treatment based on an eigenanalysis of the discretized equations. Computational results are presented which demonstrate the computational accuracy and efficiency of the method and demonstrate the effectiveness of the deforming grid, far-field nonreflecting boundary conditions, and shock capturing techniques. A comparison of the present unsteady flow predictions to other numerical, semi-analytical, and experimental methods shows excellent agreement. In addition, the linearized Euler method presented requires one or two orders-of-magnitude less computational time than traditional time marching techniques making the present method a viable design tool for aeroelastic analyses.

Clark, William S.↗
NASA NTRS1993

Numerical prediction of transition of the F-16 wing at supersonic speeds

This work is part of the high speed research program currently underway at NASA. This project has the goal of gaining understanding of the technical requirements for supersonic-hypersonic flight. Specifically, this research is part of a continuing project to study the laminar flow over swept wings at high speeds and involves the numerical prediction of the flow about the F-16XL wing. The research uses the CNS/ARC3D codes and the resulting crossflow velocity components in order to estimate transition locations on the wing. Effects of angle of attack on the extent of laminar flows was found to be minimal. This result can be attributed to the fact that a laminar flow airfoil was used in this study, which has a continuous favorable pressure gradient over approximately the first 20 percent of the chord for angles of attacks up to 10 degrees. It should also be noted that even after 20 percent chord the pressure gradient either slowly continued to increase, but never decreased before 90 percent chord, except for the higher swept cases when separation occurs. Angles of attack greater than 10 degrees were not considered since this study assumes natural laminar flow for normal supersonic cruise flight conditions.

Cummings, Russell M.↗