Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “drag estimation”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Particle kinetic simulation of high altitude hypervelocity flight

Rarefied flows about hypersonic vehicles entering the upper atmosphere or through nozzles expanding into a near vacuum may only be simulated accurately with a direct simulation Monte Carlo (DSMC) method. Under this grant, researchers enhanced the models employed in the DSMC method and performed simulations in support of existing NASA projects or missions. DSMC models were developed and validated for simulating rotational, vibrational, and chemical relaxation in high-temperature flows, including effects of quantized anharmonic oscillators and temperature-dependent relaxation rates. State-of-the-art advancements were made in simulating coupled vibration-dissociation recombination for post-shock flows. Models were also developed to compute vehicle surface temperatures directly in the code rather than requiring isothermal estimates. These codes were instrumental in simulating aerobraking of NASA's Magellan spacecraft during orbital maneuvers to assess heat transfer and aerodynamic properties of the delicate satellite. NASA also depended upon simulations of entry of the Galileo probe into the atmosphere of Jupiter to provide drag and flow field information essential for accurate interpretation of an onboard experiment. Finally, the codes have been used extensively to simulate expanding nozzle flows in low-power thrusters in support of propulsion activities at NASA-Lewis. Detailed comparisons between continuum calculations and DSMC results helped to quantify the limitations of continuum CFD codes in rarefied applications.

Boyd, Iain↗
NASA NTRS1994

Particle kinetic simulation of high altitude hypervelocity flight

Rarefied flows about hypersonic vehicles entering the upper atmosphere or through nozzles expanding into a near vacuum may only be simulated accurately with a direct simulation Monte Carlo (DSMC) method. Under this grant, researchers enhanced the models employed in the DSMC method and performed simulations in support of existing NASA projects or missions. DSMC models were developed and validated for simulating rotational, vibrational, and chemical relaxation in high-temperature flows, including effects of quantized anharmonic oscillators and temperature-dependent relaxation rates. State-of-the-art advancements were made in simulating coupled vibration - dissociation - recombination for post-shock flows. Models were also developed to compute vehicle surface temperatures directly in the code rather than requiring isothermal estimates. These codes were instrumental in simulating aerobraking of NASA's Magellan spacecraft during orbital maneuvers to assess heat transfer and aerodynamic properties of the delicate satellite. NASA also depended upon simulations of entry of the Galileo probe into the atmosphere of Jupiter to provide drag and flow field information essential for accurate interpretation of an onboard experiment. Finally, the codes have been used extensively to simulate expanding nozzle flows in low-power thrusters in support of propulsion activities at NASA-Lewis. Detailed comparisons between continuum calculations and DSMC results helped to quantify the limitations of continuum CFD codes in rarefied applications.

Boyd, Iain↗
NASA NTRS2005

ESTIMATING THE BENEFIT OF TRMM TROPICAL CYCLONE DATA IN SAVING LIVES

The Tropical Rainfall Measuring Mission (TRMM) is a joint NASA/JAXA research mission launched in late 1997 to improve our knowledge of tropical rainfall processes and climatology (Kummerow et ai., 2000; Adler et ai., 2003). In addition to being a highly successful research mission, its data are available in real time and operational weather agencies in the U.S. and internationally are using TRMM data and images to monitor and forecast hazardous weather (tropical cyclones, floods, etc.). For example, in 2004 TRMM data were used 669 times for determining tropical cyclone location fixes (National Research Council, 2004). TRMM flies at a relatively low altitude, 400 km, and requires orbit adjustment maneuvers to maintain altitude against the small drag of the atmosphere. There is enough fuel used for these maneuvers remaining on TRMM for the satellite to continue flying until 2011-12. However, most of the remaining fuel may be used to perform a controlled re-entry of the satellite into the Pacific Ocean. The fuel threshold for this operation will be reached in the summer of 2005, although the maneuver would actually occur in late 2006 or 2007. The full science mission would end in 2005 under the controlled re-entry option. This re-entry option is related to the estimated probability of injury (1/5,000) that might occur during an uncontrolled re-entry of the satellite. If the estimated probability of injury exceeds 1/10,000 a satellite is a candidate for a possible controlled re-entry. In the TRMM case the NASA Safety Office examined the related issues and concluded that, although TRMM exceeded the formal threshold, the use of TRMM data in the monitoring and forecasting of hazardous weather gave a public safety benefit that compensated for TRMM slightly exceeding the orbital debris threshold (Martin, 2002). This conclusion was based in part on results of an independent panel during a workshop on benefits of TRMM data in concluded that the benefit of TRMM data in saving lives through its use in operational forecasting could not be quantified. The objective of this paper is to describe a possible technique to estimate the number of lives saved per year and apply it to the TRMM case and the use of its data in monitoring and forecasting tropical cyclones.

Adler, Robert F.↗
NASA NTRS1989

Lifetime predictions for the Solar Maximum Mission (SMM) and San Marco spacecraft

Lifetime prediction techniques developed by the Goddard Space Flight Center (GSFC) Flight Dynamics Division (FDD) are described. These techniques were developed to predict the Solar Maximum Mission (SMM) spacecraft orbit, which is decaying due to atmospheric drag, with reentry predicted to occur before the end of 1989. Lifetime predictions were also performed for the Long Duration Exposure Facility (LDEF), which was deployed on the 1984 SMM repair mission and is scheduled for retrieval on another Space Transportation System (STS) mission later this year. Concepts used in the lifetime predictions were tested on the San Marco spacecraft, which reentered the Earth's atmosphere on December 6, 1988. Ephemerides predicting the orbit evolution of the San Marco spacecraft until reentry were generated over the final 90 days of the mission when the altitude was less than 380 kilometers. The errors in the predicted ephemerides are due to errors in the prediction of atmospheric density variations over the lifetime of the satellite. To model the time dependence of the atmospheric densities, predictions of the solar flux at the 10.7-centimeter wavelength were used in conjunction with Harris-Priester (HP) atmospheric density tables. Orbital state vectors, together with the spacecraft mass and area, are used as input to the Goddard Trajectory Determination System (GTDS). Propagations proceed in monthly segments, with the nominal atmospheric drag model scaled for each month according to the predicted monthly average value of F10.7. Calibration propagations are performed over a period of known orbital decay to obtain the effective ballistic coefficient. Progagations using plus or minus 2 sigma solar flux predictions are also generated to estimate the despersion in expected reentry dates. Definitive orbits are compared with these predictions as time expases. As updated vectors are received, these are also propagated to reentryto continually update the lifetime predictions.

Smith, E. A.↗
NASA NTRS2003

Aeroheating Thermal Model Correlation for Mars Global Surveyor (MGS) Solar Array

The Mars Global Surveyor (MGS) Spacecraft made use of aerobraking to gradually reduce its orbit period from a highly elliptical insertion orbit to its final science orbit. Aerobraking produces a high heat load on the solar arrays, which have a large surface area exposed to the airflow and relatively low mass. To accurately model the complex behavior during aerobraking, the thermal analysis needed to be tightly coupled to the spatially varying, time dependent aerodynamic heating. Also, the thermal model itself needed to accurately capture the behavior of the solar array and its response to changing heat load conditions. The correlation of the thermal model to flight data allowed a validation of the modeling process, as well as information on what processes dominate the thermal behavior. Correlation in this case primarily involved detailing the thermal sensor nodes, using as-built mass to modify material property estimates, refining solar cell assembly properties, and adding detail to radiation and heat flux boundary conditions. This paper describes the methods used to develop finite element thermal models of the MGS solar array and the correlation of the thermal model to flight data from the spacecraft drag passes. Correlation was made to data from four flight thermal sensors over three of the early drag passes. Good correlation of the model was achieved, with a maximum difference between the predicted model maximum and the observed flight maximum temperature of less than 5%. Lessons learned in the correlation of this model assisted in validating a similar model and method used for the Mars Odyssey solar array aeroheating analysis, which were used during onorbit operations.

Amundsen, Ruth M.↗
NASA NTRSDate not supplied

Machine Learning Approach for Aircraft Performance Model Parameter Estimation for Trajectory Prediction Applications

Inaccurate prediction of aircraft trajectory by ground-based decision support tools (DST) is a major concern in air traffic management (ATM). Aircraft trajectory prediction tools rely on a simplified point-mass aircraft performance model (APM) to make their predictions. Even though the performance coefficients and weight of an aircraft are a vital part of the APM’s predictions and accuracy, these coefficients are proprietary in nature and therefore, unavailable to DSTs. Current ATM research focuses on improving the estimate of some APM parameters by freezing all other coefficients. This simplified approach introduces unwanted sources of bias and negatively impacts the accuracy of the performance model. In this paper, we apply machine learning (ML) techniques for the simultaneous prediction of three key APM parameters (two drag coefficients and the initial aircraft weight). To accomplish this, we employ an ordinary differential equation (ODE) fitting approach to generate optimized APM parameter labels customized to each individual flight record. Subsequently, we train ML models to capture the relationship between the historical data and the optimized APM parameters. Two different ML model solutions are applied and APM coefficients are predicted for unseen flights. The results indicate that the ML models are able to capture the relationship between APM parameters and flight-related features with good accuracy.

trajectory prediction, machine learning, aircraft ↗
NASA NTRS1980

A low-order global gravity field of Venus and dynamical implications

An estimate of sixth-degree and sixth-order harmonic coefficients of the global gravity field of Venus is obtained by processing the long periodic variations of the mean orbital elements of the Pioneer Venus orbiter. Approximately 220 days of data are included in this reduction, which provides almost complete longitudinal coverage. Oblateness is estimated to be -5.97 + or - 3.2 x 10 to the -6. It is noted that the amplitudes of other coefficients are similar to the predicted coefficients using Kaula's rule under equal stress assumption. Atmospheric density values as a function of altitude are obtained to help model drag perturbation. A radial acceleration map at 100 km above the Venus surface is generated, and correlation between gravity anomalies and major topographic features is observed. Spectral analysis of the harmonic model suggests that the interior density anomalies are like those on earth. The orientation angles of the principal axes of the moments of inertia are computed, and deviation of the maximum moment of inertia axis from the spin axis is observed to be small (less than 5 deg). It is found that the minimum moment of inertia axis passes through the Aphrodite and Beta regions.

Ananda, M. P.↗
NASA NTRS1986

Forward-swept wing configuration designed for high maneuverability by use of a transonic computational method

A transonic computational analysis method and a transonic design procedure have been used to design the wing and the canard of a forward-swept-wing fighter configuration for good transonic maneuver performance. A model of this configuration was tested in the Langley 16-Foot Transonic Tunnel. Oil-flow photographs were obtained to examine the wind flow patterns at Mach numbers from 0.60 to 0.90. The transonic theory gave a reasonably good estimate of the wing pressure distributions at transonic maneuver conditions. Comparison of the forward-swept-wing configuration with an equivalent aft-swept-wing-configuration showed that, at a Mach number of 0.90 and a lift coefficient of 0.9, the two configurations have the same trimmed drag. The forward-swept wing configuration was also found to have trimmed drag levels at transonic maneuver conditions which are comparable to those of the HiMAT (highly maneuverable aircraft technology) configuration and the X-29 forward-swept-wing research configuration. The configuration of this study was also tested with a forebody strake.

Mann, M. J.↗
NASA NTRS1999

A Study of Nucleate Boiling with Forced Convection in Microgravity

The ultimate objective of basic studies of flow boiling in microgravity is to improve the understanding of the processes involved, as manifested by the ability to predict its behavior. This is not yet the case for boiling heat transfer even in earth gravity, despite the considerable research activity over the past 30 years. The elements that constitute the nucleate boiling process - nucleation, growth, motion, and collapse of the vapor bubbles (if the bulk liquid is subcooled) - are common to both pool and flow boiling. It is well known that the imposition of bulk liquid motion affects the vapor bubble behavior relative to pool boiling, but does not appear to significantly influence the heat transfer. Indeed, it has been recommended in the past that empirical correlations or experimental data of pool boiling be used for design purposes with forced convection nucleate boiling. It is anticipated that such will most certainly not be possible for boiling in microgravity, based on observations made with pool boiling in microgravity. In earth gravity buoyancy will act to remove the vapor bubbles from the vicinity of the heater surface regardless of how much the imposed bulk velocity is reduced, depending, of course, on the geometry of the system. Vapor bubbles have been observed to dramatically increase in size in pool boiling in microgravity, and the heat flux at which dryout took place was reduced considerably below what is generally termed the critical heat flux (CHF) in earth gravity, depending on the bulk liquid subcooling. However, at heat flux levels below dryout, the nucleate pool boiling process was enhanced considerably over that in earth gravity, in spite of the large vapor bubbles formed in microgravity and perhaps as a consequence. These large vapor bubbles tended to remain in the vicinity of the heater surface, and the enhanced heat transfer appeared to be associated with the presence of what variously has been referred to as a liquid microlayer between the bubble and the heater surface. The enhancement of the boiling process with low velocities in earth gravity for those orientations producing the formation of a liquid macrolayer described above, accompanied by "sliding" vapor bubbles, has been demonstrated. The enhancement was presented as a function of orientation, subcooling, and heated length, while a criterion for the heat transfer for mixed natural/forced convection nucleate boiling was given previously. A major unknown in the prediction and application of flow boiling heat transfer in microgravity is the upper limit of the heat flux for the onset of dryout (or critical heat flux - CHF), for given conditions of fluid-heater surfaces, including geometry, system pressure and bulk liquid subcooling. It is clearly understood that the behavior in microgravity will be no different than on earth with sufficiently high flow velocities, and would require no space experimentation. However, the boundary at which this takes place is still an unknown. Previous results of CHF measurements were presented for low velocity flow boiling at various orientations in earth gravity as a function of flow velocity and bulk liquid subcooling, along with preliminary measurements of bubble residence times on a flat heater surface. This showed promise as a parameter to be used in modeling the CHF, both in earth gravity and in microgravity. The objective of the work here is to draw attention to and show results of current modeling efforts for the CHF, with low velocities in earth gravity at different orientations and subcoolings. Many geometrical possibilities for a heater surface exist in flowing boiling, with boiling on the inner and outer surfaces of tubes perhaps being the most common. If the vapor bubble residence time on and departure size from the heater surface bear a relationship to the CHF, as results to be given indicate, it is important that visualization of and access to vapor bubble growth be conveniently available for research purposes. In addition, it is desirable to reduce the number of variables as much as possible in a fundamental study. These considerations dictated the use of a flat heater surface, which is rectangular in shape, 1.91 cm by 3.81 cm (0.75 x 1.5 inches), consisting either of a 400 Angstrom thick semi-transparent gold film sputtered on a quartz substrate which serves simultaneously as a heater and a resistance thermometer, or a copper substrate of the same size. The heater substrate is a disc which can be rotated so that the heated length in the flow direction can be changed from 1.91 to 3.81 cm (0.75 to 1.5 inches). The fluid is R-113, and the velocities can be varied between 0.5 cm/s and 60 cm/s. For a sufficiently low velocity the CHF can be modeled reasonably well at various orientations by the correlation for pool boiling corrected for the influence of bulk liquid subcooling, multiplied by the square root of q, the angle relative to horizontal. This arises from equating buoyancy and drag forces in the inverted positions where the vapor bubbles are held against the heater surface as they slide. A distortion of the measurements relative to pool boiling occurs as the flow velocity increases. In modeling this effect at different levels of subcooling it appeared appropriate to estimate the volumetric rate of vapor generation, using measurements of bubble frequency (or residence time), void fraction and average bubble boundary layer thickness. These were determined with the use of a platinum hot wire probe 0.025 mm in diameter by 1.3 mm long, applying a constant current to distinguish between contact with liquid or vapor. Two-dimensional spatial variations are obtained with a special mechanism to resolve displacements in increments of 0.025 mm. From such measurements it was determined that the fraction of the surface heat transfer resulting in evaporation varies inversely with the subcooling correction factor for the CHF. The measured inverse bubble residence time is normalized relative to that predicted for an infinite horizontal flat plate at the CHF, and is correlated well with the CHF normalized relative to that for pool boiling, for various orientation angles and subcooling levels. This correspondence is then combined with a normalizing factor for the energy flux leaving the heater surface at the CHF and the computed bubble radius at departure, determined from the balance between the outward velocity of the interface due to evaporation and the buoyance induced velocity of the center of mass of the bubble. The product of the CHF and the corresponding residence time was determined to be a constant for all orientations at a given bulk flow velocity and liquid subcooling, and must be determined empirically for each velocity and subcooling at present. It then becomes possible to predict the CHF for the different orientations, velocities, and subcoolings. These are compared with normalized measurements of the CHF for velocities ranging from 4 cm/s to 55 cm/s, subcoolings from 2.8 to 22.2 K, over orientations angles of 360 degrees.

Merte, Herman, Jr.↗
NASA NTRS1975

Sensitivity analysis of short-arc station coordinate determinations from range data

The accurate determination of the geocentric coordinates of a tracking station is essential for most geodetic and geophysical satellite applications. Since most of these satellites are close to the earth, the geopotential model is a dominant source of error which significantly influences station coordinate determinations. Other sources, such as GM error and drag, also influence the accuracy of the station coordinate determination. One technique for reducing the effect of these errors is to use short-arcs consisting of a few passes of the satellite over the tracking station. This paper analyzes the sensitivity of short-arc station coordinate estimates to various errors in the physical model, to the number of observations, and to the station-satellite geometry using simulated as well as real data.

Schutz, B. E.↗
NASA NTRS1976

Sensitivity analysis of short-arc station coordinate determinations from range data

The accurate determination of the geocentric coordinates of a tracking station is essential for most geodetic and geophysical satellite applications. Since most of these satellites are close to the earth, the geopotential model is a dominant source of error which significantly influences station coordinate determinations. Other sources, such as GM error and drag, also influence the accuracy of the station coordinate determination. One technique for reducing the effect of these errors is to use short-arcs consisting of a few passes of the satellite over the tracking station. This paper analyzes the sensitivity of short-arc station coordinate estimates to various errors in the physical model, to the number of observations, and to the station-satellite geometry using simulated as well as real data.

Schutz, B. E.↗
NASA NTRS1987

Effect of a round airfoil nose on leading-edge suction

Kulfan (1979) assumed that the angle of attack for initial vortex separation on a slender wing with rounded leading edges could be obtained by equating the leading-edge suction (LES) and nose drag coefficients. In the present study, this assumption is examined and is shown to predict reasonably well the initial angle of attack at which laminar separation occurs near the airfoil nose. However, the assumption is shown to be slightly less accurate for thick or cambered airfoils. Attainable LES estimated by Kulfan's method seemed to agree well with that obtained from an airfoil aerodynamics code and experimental data on a NACA 64A009 airfoil at M = 0.4 and Re = 0.86 x 10 to the 6th.

Lan, C. Edward↗
NASA NTRS1999

Piloted Simulation Study of the Effects of High-Lift Aerodynamics on the Takeoff Noise of a Representative High-Speed Civil Transport

As part of an effort between NASA and private industry to reduce airport-community noise for high-speed civil transport (HSCT) concepts, a piloted simulation study was initiated for the purpose of predicting the noise reduction benefits that could result from improved low-speed high-lift aerodynamic performance for a typical HSCT configuration during takeoff and initial climb. Flight profile and engine information from the piloted simulation were coupled with the NASA Langley Aircraft Noise Prediction Program (ANOPP) to estimate jet engine noise and to propagate the resulting source noise to ground observer stations. A baseline aircraft configuration, which also incorporated different levels of projected improvements in low-speed high-lift aerodynamic performance, was simulated to investigate effects of increased lift and lift-to-drag ratio on takeoff noise levels. Simulated takeoff flights were performed with the pilots following a specified procedure in which either a single thrust cutback was performed at selected altitudes ranging from 400 to 2000 ft, or a multiple-cutback procedure was performed where thrust was reduced by a two-step process. Results show that improved low-speed high-lift aerodynamic performance provides at least a 4 to 6 dB reduction in effective perceived noise level at the FAA downrange flyover measurement station for either cutback procedure. However, improved low-speed high-lift aerodynamic performance reduced maximum sideline noise levels only when using the multiple-cutback procedures.

Glaab, Louis J.↗
NASA NTRS1973

The slab thickness of the mid-latitude ionosphere.

The thickness of the peak of the ionosphere depends primarily on the temperature T sub n of the neutral gas, and corresponds approximately to an alpha-Chapman layer at a temperature of 0.87T sub n. The overall slab thickness, as given by Faraday rotation measurements, is then tau = 0.22T sub n + 7 km. Expansion of the topside ionosphere, and changes in the E- and F1-regions increase tau by about 20 km during the day in summer. Near solar minimum, tau is increased by a lowering of the O(+)/H(+) transition height; if the neutral temperature T sub n is estimated, this height can be obtained from observed values of tau. Hourly values of slab thickness were determined over a period of 6 yr at 34 and 42 S. Near solar maximum the nighttime values were about 260 km in all seasons. The corresponding neutral temperatures agree with satellite drag values; they show a semiannual variation of 14% and a seasonal change of 5%. Daytime values of tau were about 230 km in winter and 320 km in summer, implying a seasonal change of 30% in T sub n.

Titheridge, J. E.↗
NASA NTRS1984

Long-lifetime Martian orbit selection using a time-dependent model of the Martian atmosphere

A mathematical model of the time-dependent Martian atmosphere has been developed in order to accurately calculate the effects of aerodynamic drag on a low altitude satellite. The time-dependent properties of the model include solar activity effects, dust storm effects, seasonal and diurnal variations, and annual motion effects. Position effects are accounted for through Martian latitude and longitude. Expected values of mass density, temperature, scale height, and the estimated standard deviation of the mass density are provided. An example of the use of the model in selecting an orbit for the Mars Geochemical/Climatology Orbiter is given.

Culp, R. D.↗
NASA NTRS1995

A case for bone canaliculi as the anatomical site of strain generated potentials

We address the question of determining the anatomical site that is the source of the experimentally observed strain generated potentials (SGPs) in bone tissue. There are two candidates for the anatomical site that is the SGP source, the collagen-hydroxyapatite porosity and the larger size lacunar-canalicular porosity. In the past it has been argued, on the basis of experimental data and a reasonable model, that the site of the SGPs in bone is the collagen-hydroxyapatite porosity. The theoretically predicted pore radius necessary for the SGPs to reside in this porosity is 16 nm, which is somewhat larger than the pore radii estimated from gas adsorption data where the preponderance of the pores were estimated to be in the range 5-12.5 nm. However, this pore size is significantly larger than the 2 nm size of the small tracer, microperoxidase, which appears to be excluded from the mineralized matrix. In this work a similar model, but one in which the effects of fluid dynamic drag of the cell surface matrix in the bone canaliculi are included, is used to show that it is possible for the generation of SGPs to be associated with the larger size lacunar-canalicular porosity when the hydraulic drag and electrokinetic contribution of the bone fluid passage through the cell coat (glycocalyx) is considered. The consistency of the SGP data with this model is demonstrated. A general boundary condition is introduced to allow for current leakage at the bone surface. The results suggest that the current leakage is small for the in vitro studies in which the strain generated potentials have been measured.

Non-NASA Center↗
NASA NTRS2004

Adjoint Algorithm for CAD-Based Shape Optimization Using a Cartesian Method

Adjoint solutions of the governing flow equations are becoming increasingly important for the development of efficient analysis and optimization algorithms. A well-known use of the adjoint method is gradient-based shape optimization. Given an objective function that defines some measure of performance, such as the lift and drag functionals, its gradient is computed at a cost that is essentially independent of the number of design variables (geometric parameters that control the shape). More recently, emerging adjoint applications focus on the analysis problem, where the adjoint solution is used to drive mesh adaptation, as well as to provide estimates of functional error bounds and corrections. The attractive feature of this approach is that the mesh-adaptation procedure targets a specific functional, thereby localizing the mesh refinement and reducing computational cost. Our focus is on the development of adjoint-based optimization techniques for a Cartesian method with embedded boundaries.12 In contrast t o implementations on structured and unstructured grids, Cartesian methods decouple the surface discretization from the volume mesh. This feature makes Cartesian methods well suited for the automated analysis of complex geometry problems, and consequently a promising approach to aerodynamic optimization. Melvin et developed an adjoint formulation for the TRANAIR code, which is based on the full-potential equation with viscous corrections. More recently, Dadone and Grossman presented an adjoint formulation for the Euler equations. In both approaches, a boundary condition is introduced to approximate the effects of the evolving surface shape that results in accurate gradient computation. Central to automated shape optimization algorithms is the issue of geometry modeling and control. The need to optimize complex, "real-life" geometry provides a strong incentive for the use of parametric-CAD systems within the optimization procedure. In previous work, we presented an effective optimization framework that incorporates a direct-CAD interface. In this work, we enhance the capabilities of this framework with efficient gradient computations using the discrete adjoint method. We present details of the adjoint numerical implementation, which reuses the domain decomposition, multigrid, and time-marching schemes of the flow solver. Furthermore, we explain and demonstrate the use of CAD in conjunction with the Cartesian adjoint approach. The final paper will contain a number of complex geometry, industrially relevant examples with many design variables to demonstrate the effectiveness of the adjoint method on Cartesian meshes.

Nemec, Marian↗