Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “first wall”

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 235 records · Page 13

Heat transfer in rotating serpentine passages with trips skewed to the flow

Experiments were conducted to determine the effects of buoyancy and Coriolis forces on heat transfer in turbine blade internal coolant passages. The experiments were conducted with a large scale, multi-pass, heat transfer model with both radially inward and outward flow. Trip strips, skewed at 45 deg to the flow direction, were machined on the leading and trailing surfaces of the radial coolant passages. An analysis of the governing flow equations showed that four parameters influence the heat transfer in rotating passages: coolant-to-wall temperature, rotation number, Reynolds number, and radius-to-passage hydraulic diameter ratio. The first three of these four parameters were varied over ranges which are typical of advanced gas turbine engine operating conditions. Results were correlated and compared to previous results from similar stationary and rotating models with smooth walls and with trip strips normal to the flow direction. The heat transfer coefficients on surfaces, where the heat transfer decreased with rotation and buoyancy, decreased to as low as 40 percent of the value without rotation. However, the maximum values of the heat transfer coefficients with high rotation were only slightly above the highest levels previously obtained with the smooth wall models. It was concluded that (1) both Coriolis and buoyancy effects must be considered in turbine blade cooling designs with trip strips, (2) the effects of rotation are markedly different depending upon the flow direction, and (3) the heat transfer with skewed trip strips is less sensitive to buoyancy than the heat transfer in models with either smooth or normal trips. Therefore, skewed trip strips rather than normal trip strips are recommended and geometry-specific tests are required for accurate design information.

Johnson, B. V.↗

Bubble Formation on a Wall in Cross-Flowing Liquid and Surrounding Fluid Motion,With and Without Heating

Application of gas-liquid two-phase flow systems for space-based thermal management and for the HEDS program demands a precise control of bubble size distribution in liquid. The necessity of bulk liquid motion for controlling bubble size and frequency in the space environment has been suggested by recent studies on pool, forced convection boiling and bubble formation in flowing liquid. The present work, consisting of two parts, explores bubble generation at wall in a cross-flowing liquid, i.e., in a forced convection boiling configuration. A schematic is shown. The first part looks into the bubble formation process under isothermal conditions in a reduced gravity environment, by injecting gas through a hole in the wall of a flowing liquid channel. In the latter part with channel wall heating, flow and temperature fields near a single bubble are studied under normal (1-g) and micro-gravity (mu-g) conditions. The objective of the isothermal experiments is to experimentally investigate the effects of liquid cross-flow velocity, gas flow rate, and orifice diameter on bubble formation. Data were taken mainly under reduced gravity conditions but some data were taken in normal gravity for comparison. The reduced gravity experiment was conducted aboard the NASA DC-9 Reduced Gravity Aircraft. The results show that the process of bubble formation and detachment depends on gravity, the orifice diameter (D(sub N)), the gas flow rate (Q(sub g)), and the liquid cross-flow velocity (U(sub L)). The reduced gravity data are shown. The data are analyzed based on a force balance, and two different detachment mechanisms are identified. When the gas momentum is large, the bubble detaches from the injection orifice as the gas momentum overcomes the attaching effects of liquid drag and inertia. The surface tension force is much reduced because a large part of the bubble pinning edge at the orifice is lost as the bubble axis is tilted by the liquid flow. When the gas momentum is small, the force balance in the liquid flow direction is important, and the bubble detaches when the bubble axis inclination exceeds a certain angle. With wall heating, liquid motion around an air bubble in cross-flowing 2cs silicone oil is experimentally investigated in 1-g. A spectral element based steady 2D numerical model is also developed. The traces of particles from experimental flow visualization and the corresponding computed streamlines are shown. At the upstream side of the bubble facing the cross-flow, thermocapillary and forced convection create liquid motion away from the wall, up along the surface. At the downstream side, a competing interaction between the two creates a recirculation cell, causing the bulk liquid to stagnate on the surface and separate thereafter. The important dimensionless parameters are - Surface tension and local cross-flow Reynolds numbers R(sub sigma)(U(sub ref)a/v) and Re(sub loc)(U(sub L)a/v), respectively based on reference thermocapillary U(sub ref)(sigma(sub T)Delta.T/mu, Delta T=T(sub wall)-T(sub liquid)) and local cross-flow velocity U(sub L), Prandtl number P(sub r) and Grashoff number Gr(rho.g.beta.DeltaTa(sup 3)/mu.v). Variation of the stagnation point with R(sub sigma) and Re(sub loc) is shown. Also shown are good agreement between experimental and numerical results in 1-g. The computational model is extended to mu-g condition to investigate temperature and velocity on the bubble surface, stagnation and reattachment points of the recirculation cell and wall heat transfer. It is observed that wall heating significantly alters the flow field around the bubble and thus the forces acting on the bubble, which control its detachment. Thus a combination of heating and liquid cross-flow can be utilized to precisely control bubble formation in a mu-g environment.

Bhunia, Avijit↗

On the Goertler instability in hypersonic flows: Sutherland law fluids and real gas effects

The Goertler vortex instability mechanism in a hypersonic boundary layer on a curved wall is investigated. The precise roles of the effects of boundary layer growth, wall cooling, and gas dissociation is clarified in the determination of stability properties. It is first assumed that the fluid is an ideal gas with viscosity given by Sutherland's law. It is shown that when the free stream Mach number M is large, the boundary layer divides into two sublayers: a wall layer of O(M sup 3/2) thickness over which the basic state temperature is O(M squared) and a temperature adjustment layer of O(1) thickness over which the basic state temperature decreases monotonically to its free stream value. Goertler vortices which have wavelengths comparable with the boundary layer thickness are referred to as wall modes. It is shown that their downstream evolution is governed by a set of parabolic partial differential equations and that they have the usual features of Goertler vortices in incompressible boundary layers. As the local wavenumber increases, the neutral Goertler number decreases and the center of vortex activity moves towards the temperature adjustment layer. Goertler vortices with wavenumbers of order one or larger must necessarily be trapped in the temperature adjustment layer and it is this mode which is most dangerous. For this mode, it was found that the leading order term in the Goertler number expansion is independent of the wavenumber and is due to the curvature of the basic state. This term is also the asymptotic limit of the neutral Goertler numbers of the wall mode. To determine the higher order corrections terms in the Goertler number expansion, two wall curvature cases are distinguished. Real gas effects were investigated by assuming that the fluid is an ideal dissociating gas. It was found that both gas dissociation and wall cooling are destabilizing for the mode trapped in the temperature adjustment layer, but for the wall mode trapped near the wall the effect of gas dissociation can be either destabilizing or stabilizing.

Fu, Yibin B.↗

Analytical Modeling of Pressure Wall Hole Size and Maximum Tip-to-Tip Crack Length for Perforating Normal and Oblique Orbital Debris Impacts

This report presents the results of a study whose objective was to develop first-principles-based models of hole size and maximum tip-to-tip crack length for a spacecraft module pressure wall that has been perforated in an orbital debris particle impact. The hole size and crack length models are developed by sequentially characterizing the phenomena comprising the orbital debris impact event, including the initial impact, the creation and motion of a debris cloud within the dual-wall system, the impact of the debris cloud on the pressure wall, the deformation of the pressure wall due to debris cloud impact loading prior to crack formation, pressure wall crack initiation, propagation, and arrest, and finally pressure wall deformation following crack initiation and growth. The model development has been accomplished through the application of elementary shock physics and thermodynamic theory, as well as the principles of mass, momentum, and energy conservation. The predictions of the model developed herein are compared against the predictions of empirically-based equations for hole diameters and maximum tip-to-tip crack length for three International Space Station wall configurations. The ISS wall systems considered are the baseline U.S. Lab Cylinder, the enhanced U.S. Lab Cylinder, and the U.S. Lab Endcone. The empirical predictor equations were derived from experimentally obtained hole diameters and crack length data. The original model predictions did not compare favorably with the experimental data, especially for cases in which pressure wall petalling did not occur. Several modifications were made to the original model to bring its predictions closer in line with the experimental results. Following the adjustment of several empirical constants, the predictions of the modified analytical model were in much closer agreement with the experimental results.

Schonberg, William P.↗

Computational Fluid Dynamics Simulation of Dual Bell Nozzle Film Cooling

Marshall Space Flight Center (MSFC) and Armstrong Flight Research Center (AFRC) are working together to advance the technology readiness level (TRL) of the dual bell nozzle concept. Dual bell nozzles are a form of altitude compensating nozzle that consists of two connecting bell contours. At low altitude the nozzle flows fully in the first, relatively lower area ratio, nozzle. The nozzle flow separates from the wall at the inflection point which joins the two bell contours. This relatively low expansion results in higher nozzle efficiency during the low altitude portion of the launch. As ambient pressure decreases with increasing altitude, the nozzle flow will expand to fill the relatively large area ratio second nozzle. The larger area ratio of the second bell enables higher Isp during the high altitude and vacuum portions of the launch. Despite a long history of theoretical consideration and promise towards improving rocket performance, dual bell nozzles have yet to be developed for practical use and have seen only limited testing. One barrier to use of dual bell nozzles is the lack of control over the nozzle flow transition from the first bell to the second bell during operation. A method that this team is pursuing to enhance the controllability of the nozzle flow transition is manipulation of the film coolant that is injected near the inflection between the two bell contours. Computational fluid dynamics (CFD) analysis is being run to assess the degree of control over nozzle flow transition generated via manipulation of the film injection. A cold flow dual bell nozzle, without film coolant, was tested over a range of simulated altitudes in 2004 in MSFC's nozzle test facility. Both NASA centers have performed a series of simulations of that dual bell to validate their computational models. Those CFD results are compared to the experimental results within this paper. MSFC then proceeded to add film injection to the CFD grid of the dual bell nozzle. A series of nozzle pressure ratios and film coolant flow rates are investigated to determine the effect of the film injection on the nozzle flow transition behavior. The results of this CFD study of a dual bell with film injection are presented in this paper.

Braman, Kalen↗

Flow control in a diffusing S-Duct

Accurate measurements have been made of secondary flow in a 1.51 area ratio diffusing 30 deg - 30 deg S-Duct with circulair cross section. Turbulent flow was entering the duct at Mach number of 0.6, the boundary layer thickness at the duct entrance was ten percent of the duct inlet diameter. Through measurements made, local flow velocity vector as well as static and total pressures mapping of the flow at several stations were obtained. Strong secondary flow was measured in the first bend which continued into the second bend with new vorticity produced in there in the opposite direction. Surface oil flow visualization and wall pressures indicated a region of separated flow starting at theta approximately equal to 22 deg on the inside of the first bend up to theta approximately equal to 44 deg on the outside of the second bend. The flow separated in 'cyclone' form and never reattached in the duct. As a result of the secondary flow and the flow separation, significant total pressure distortion was observed at the exit of the duct. Using flow control devices the separation was eliminated while the exit distortion was improved.

Vakili, A. D.↗

Boiling inception in trichlorotrifluoroethane during forced convection at high pressures

The inception of bubbles during forced convection was studied experimentally by using trichlorotrifluoroethane (R-113 or Freon-113). The experiments were performed in a rectangular channel, 12.7 x 9.5 mm in cross section. Heating was from a 3.2 mm wide strip embedded in the longer side of the channel. The pressures studied ranged from 3.6 to 20.7 bar, mass velocities from 700 to 600 kg/sq m/sec, and inlet subcoolings from 26 to 97 C. Photographs of the flow were used to determine when bubbles first appeared on the heated surface. These data were compared with wall temperature measurements and inception theories. A reasonable method for calculating the complete boiling curve was found to agree with these results.

Dougall, R. S.↗

A new method of optimal design for a two-dimensional diffuser by using dynamic programming

A new method for predicting the optimal velocity distribution on the wall of a two dimensional diffuser is presented. The method uses dynamic programming to solve the optimal control problem with inequality constraints of state variables. The physical model of optimization is designed to prevent the separation of the boundary layer while approaching the maximum pressure ratio in a diffuser of a specified length. The computational results are in fair agreement with the experimental ones. Optimal velocity distribution on a diffuser wall is said to occur when the flow decelerates quickly at first and then smoothly, while the flow is near separation, but always protected from it. The optimal velocity distribution can be used to design the contour of the diffuser.

Gu, Chuangang↗

Supported transition metal catalysts for para- to ortho-hydrogen conversion

The main goal of this study was to develop and improve on existing catalysts for the conversion of ortho- to para-hydrogen. Starting with a commercially available Air Products nickel silicate, which had a beta value of 20, we were trying to synthesize catalysts that would be an improvement to AP. This was accomplished by preparing silicates with various metals as well as different preparation methods. We also prepared supported ruthenium catalysts by various techniques using several metal precursors to improve present technology. What was also found was that the activation conditions prior to catalytic testing was highly important for both the silicates and the supported ruthenium catalysts. While not the initial focus of the research, we made some interesting observations into the adsorption of H2 on ruthenium. This helped us to get a better understanding of how ortho- to para-H2 conversion takes place, and what features in a catalyst are important to optimize activity. Reactor design was the final area in which some interesting conclusions were drawn. As discussed earlier, the reactor catalyst bed must be constructed using straight 1/8 feet OD stainless steel tubing. It was determined that the use of 1/4 feet OD tubing caused two problems. First, the radius from the center of the bed to the wall was too great for thermal equilibrium. Since the reaction of ortho- to para-H2 is exothermic, the catalyst bed center was warmer than the edges. Second, the catalyst bed was too shallow using a 1/4 feet tube. This caused reactant blow-by which was thought to decrease the measured activity when the flow rate was increased. The 1/8 feet tube corrected both of these concerns.

Brooks, Christopher J.↗

Computational Hemodynamics Involving Artificial Devices

This paper reports the progress being made towards developing complete blood flow simulation capability in human, especially, in the presence of artificial devices such as valves and ventricular assist devices. Devices modeling poses unique challenges different from computing the blood flow in natural hearts and arteries. There are many elements needed such as flow solvers, geometry modeling including flexible walls, moving boundary procedures and physiological characterization of blood. As a first step, computational technology developed for aerospace applications was extended in the recent past to the analysis and development of mechanical devices. The blood flow in these devices is practically incompressible and Newtonian, and thus various incompressible Navier-Stokes solution procedures can be selected depending on the choice of formulations, variables and numerical schemes. Two primitive variable formulations used are discussed as well as the overset grid approach to handle complex moving geometry. This procedure has been applied to several artificial devices. Among these, recent progress made in developing DeBakey axial flow blood pump will be presented from computational point of view. Computational and clinical issues will be discussed in detail as well as additional work needed.

Kwak, Dochan↗

Blood Pump Development Using Rocket Engine Flow Simulation Technology

This paper reports the progress made towards developing complete blood flow simulation capability in humans, especially in the presence of artificial devices such as valves and ventricular assist devices. Devices modeling poses unique challenges different from computing the blood flow in natural hearts and arteries. There are many elements needed to quantify the flow in these devices such as flow solvers, geometry modeling including flexible walls, moving boundary procedures and physiological characterization of blood. As a first step, computational technology developed for aerospace applications was extended to the analysis and development of a ventricular assist device (VAD), i.e., a blood pump. The blood flow in a VAD is practically incompressible and Newtonian, and thus an incompressible Navier-Stokes solution procedure can be applied. A primitive variable formulation is used in conjunction with the overset grid approach to handle complex moving geometry. The primary purpose of developing the incompressible flow analysis capability was to quantify the flow in advanced turbopump for space propulsion system. The same procedure has been extended to the development of NASA-DeBakey VAD that is based on an axial blood pump. Due to massive computing requirements, high-end computing is necessary for simulating three-dimensional flow in these pumps. Computational, experimental, and clinical results are presented.

Kwak, Dochan↗

Benchmark Wall Heat Flux Data for a GO2/GH2 Single Element Combustor

Wall heat flux measurements in a 1.5 in. diameter circular cross-section rocket chamber for a uni-element shear coaxial injector element operating on gaseous oxygen (GOz)/gaseous hydrogen (GH,) propellants are presented. The wall heat flux measurements were made using arrays of Gardon type heat flux gauges and coaxial thermocouple instrumentation. Wall heat flux measurements were made for two cases. For the first case, GOZ/GHz oxidizer-rich (O/F=l65) and fuel-rich preburners (O/F=1.09) integrated with the main chamber were utilized to provide vitiated hot fuel and oxidizer to the study shear coaxial injector element. For the second case, the preburners were removed and ambient temperature gaseous oxygen/gaseous hydrogen propellants were supplied to the study injector. Experiments were conducted at four chamber pressures of 750, 600, 450 and 300psia for each case. The overall mixture ratio for the preburner case was 6.6, whereas for the ambient propellant case, the mixture ratio was 6.0. Total propellant flow was nominally 0.27-0.29 Ibm/s for the 750 psia case with flowrates scaled down linearly for lower chamber pressures. The axial heat flux profile results for both the preburner and ambient propellant cases show peak heat flux levels a t axial locations between 2.0 and 3.0 in. from the injector face. The maximum heat flux level was about two times greater for the preburner case. This is attributed to the higher injector fuel-to-oxidizer momentum flux ratio that promotes mixing and higher initial propellant temperature for the preburner case which results in a shorter reaction zone. The axial heat flux profiles were also scaled with respect to the chamber pressure to the power 0.8. The results at the four chamber pressures for both cases collapsed to a single profile indicating that at least to first approximation, the basic fluid dynamic structures in the flow field are pressure independent as long as the chamber/njector/nozzle geometry and injection velocities remain the same.

Marshall, William M.↗

Validation of High-Resolution CFD Method for Slosh Damping Extraction of Baffled Cryogenic Propellant Tanks

Propellant slosh is a potential source of disturbance critical to the stability of space vehicles. The slosh dynamics are typically represented by a mechanical model of a spring-mass-damper. This mechanical model is then included in the equation of motion of the entire vehicle for Guidance, Navigation and Control analysis. A Volume-Of-Fluid (VOF) based Computational Fluid Dynamics (CFD) program developed at MSFC was applied to extract slosh damping in the baffled tank from the first principle. First the experimental data using water with sub-scale smooth wall tank were used as the baseline validation. It is demonstrated that CFD can indeed accurately predict low damping values from the smooth wall at different fill levels. The damping due to a ring baffles at different depths from the free surface was then simulated, and fairly good agreement with experimental measurement was observed. Comparison with an empirical correlation of Miles equation is also made.

Yang, H. Q.↗

Nonlinear Buckling Analysis of Tow-Steered Composite Cylinders with Cutouts

The buckling and postbuckling behavior of two composite tow-steered shells with cutouts of different sizes is assessed using nonlinear finite element (FE) analysis and compared to experimental measurements. The cylindrical shells are manufactured using an automated fiber placement system, where the shells’ fiber orientation angles vary continuously around the shell circumference from ±10 degrees on the crown and keel to ±45 degrees on the sides. One shell features thickness variations due to tow overlaps that result from application of all 24 tows during each pass of the fiber placement system. The second shell uses the system’s tow drop/add capability to achieve a more uniform wall thickness without overlaps. Unreinforced cutouts of two different sizes—the first smaller cutout representing a passenger door on a commercial aircraft and the second larger cutout a cargo door—were machined into each of the two cylinders resulting in a total of four test cases. These cylinders were tested in axial compression and buckled elastically in previous work and are now analyzed using nonlinear FE models to compare bifurcation buckling loads as well as the load-displacement response in the prebuckling and postbuckling regimes. For all four shells analyzed, the prebuckling stiffness, buckling load, and deformation mode sequence throughout the loading-unloading cycle is accurately reproduced by the models. In particular, the shells first buckle locally around the cutouts in a stable (super-critical) manner with only a slight decrease in axial stiffness, which occurs due to the favorable load redistribution facilitated by tow steering. The shells then buckle globally in an unstable (sub-critical) manner with diamond-shaped buckles forming to the left and right of the cutouts. The buckling load of all shells with cutouts is at least 82% of the buckling load of the pristine shells without cutouts. Overall, the ability to sustain local buckling phenomena, and the relatively small reductions in global buckling load compared to pristine shells without cutouts, demonstrates the great potential of using tow steering to mitigate the adverse effects of cutouts in axially-compressed shell structures.

composites↗

Nonlinear Buckling Analysis of Tow-Steered Composite Cylinders with Cutouts

The buckling and postbuckling behavior of two composite tow-steered shells with cutouts of different sizes is assessed using nonlinear finite element (FE) analysis and compared to experimental measurements. The cylindrical shells are manufactured using an automated fiber placement system, where the shells’ fiber orientation angles vary continuously around the shell circumference from 10 degrees on the crown and keel to 45 degrees on the sides. One shell features thickness variations due to tow overlaps that result from application of all 24 tows during each pass of the fiber placement system. The second shell uses the system’s tow drop/add capability to achieve a more uniform wall thickness without overlaps. Unreinforced cutouts of two different sizes—the first smaller cutout representing a passenger door on a commercial aircraft and the second larger cutout a cargo door—were machined into each of the two cylinders resulting in a total of four test cases. These cylinders were tested in axial compression and buckled elastically in previous work and are now analyzed using nonlinear FE models to compare bifurcation buckling loads as well as the load-displacement response in the prebuckling and postbuckling regimes. For all four shells analyzed, the prebuckling stiffness, buckling load, and deformation mode sequence throughout the loading-unloading cycle is accurately reproduced by the models. In particular, the shells first buckle locally around the cutouts in a stable (super-critical) manner with only a slight decrease in axial stiffness, which occurs due to the favorable load redistribution facilitated by tow steering. The shells then buckle globally in an unstable (sub-critical) manner with diamond-shaped buckles forming to the left and right of the cutouts. The buckling load of all shells with cutouts is at least 82% of the buckling load of the pristine shells without cutouts. Overall, the ability to sustain local buckling phenomena, and the relatively small reductions in global buckling load compared to pristine shells without cutouts, demonstrates the great potential of using tow steering to mitigate the adverse effects of cutouts in axially-compressed shell structures.

Buckling↗

Development of Advanced Manufacturing Approaches for Single-Piece Launch Vehicle Structures

Advanced, near-net shape manufacturing methods have the potential to enable production of structures with fewer welds and reduced machining requirements. Two such methods that are suitable for the manufacturing of single-piece, stiffened barrel-shaped structures are presented. The first solution deployed existing manufacturing technology to produce a thick-walled barrel with integrally machined stiffeners. In this study, an 8-ft. diameter barrel was produced and subjected to a buckling test. Data on the manufacturing and testing of this barrel are provided and compared against the traditional multi-piece weld construction approach. The single-piece barrel resulted in a 28% greater load carrying capacity than the welded barrel. A second solution utilized a novel flow-forming technique to produce the barrel and stiffeners in one process without the need for welding or machining. This is an emerging manufacturing method known as the Integrally Stiffened Cylinder (ISC) process. This innovative process is being evaluated for launch vehicle and commercial aircraft manufacturing. The one-piece, stiffened ISC barrels, which have been successfully fabricated at 10 ft. in diameter, offer a direct replacement for conventional multi-piece, welded or riveted structures. A cost-benefit analysis for launch vehicle cryogenic propellant tanks estimated that the ISC process offers up to a 50% reduction in manufacturing costs and a 10% reduction in mass.

integrally stiffened cylinder (ISC)↗

Sonic wind tunnel of the Institute of Fluid Mechanics of Lille

A 65 hp wind tunnel with a 40 mm by 240 mm airstream is described. This wind tunnel can achieve speeds in the neighborhood of the speed of sound, both subsonic and supersonic. It is useful in studying the transonic bump technique. The test section is 600 mm long. The side walls are made of transparent glass, and both the upper and lower walls are deformable, each through the use of nine jacks with elastic sleeves. So as to avoid condensation, the airstream's temperature is stabilized by an air exchanger at the temperature of the outside air. The first results for supersonic operation, the distribution of Mach numbers within the airstream between the parallel walls, the value of the use factor, and the diffuser's efficiency are all given.

Gontier, G.↗

The effects of compliant walls on transition and turbulence

Two theoretical approaches to the phenomenon of drag reduction by a compliant wall are proposed. One is to study the limit cycle oscillations of the Navier-Stokes equations for a boundary layer flow over a compliant wall or a Poiseuille flow in a channel with a compliant wall. The other is a complementary study of a simpler, sequential model, which partially decouples the fluid and wall motions. In this approach, a 'Total Model' is constructed from suitable component models. While the first study is more rigorously based, it will be initially limited, by computational constraints, to relatively low Reynolds numbers. The second, more approximate approach, will allow consideration of a fully developed turbulent flow. The Total Model will predict the compliant wall motion accurately, but will not per se be sufficient to predict the skin friction drag. However, once the wall motion is known, the skin friction drag and the pressure drag can be computed separately, using one of the existing models.

Dowell, E. H.↗