Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic equivalent”

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 613 records · Page 34

Experimental Study of an Axisymmetric Dual Throat Fluidic Thrust Vectoring Nozzle for Supersonic Aircraft Application

An axisymmetric version of the Dual Throat Nozzle concept with a variable expansion ratio has been studied to determine the impacts on thrust vectoring and nozzle performance. The nozzle design, applicable to a supersonic aircraft, was guided using the unsteady Reynolds-averaged Navier-Stokes computational fluid dynamics code, PAB3D. The axisymmetric Dual Throat Nozzle concept was tested statically in the Jet Exit Test Facility at the NASA Langley Research Center. The nozzle geometric design variables included circumferential span of injection, cavity length, cavity convergence angle, and nozzle expansion ratio for conditions corresponding to take-off and landing, mid climb and cruise. Internal nozzle performance and thrust vectoring performance was determined for nozzle pressure ratios up to 10 with secondary injection rates up to 10 percent of the primary flow rate. The 60 degree span of injection generally performed better than the 90 degree span of injection using an equivalent injection area and number of holes, in agreement with computational results. For injection rates less than 7 percent, thrust vector angle for the 60 degree span of injection was 1.5 to 2 degrees higher than the 90 degree span of injection. Decreasing cavity length improved thrust ratio and discharge coefficient, but decreased thrust vector angle and thrust vectoring efficiency. Increasing cavity convergence angle from 20 to 30 degrees increased thrust vector angle by 1 degree over the range of injection rates tested, but adversely affected system thrust ratio and discharge coefficient. The dual throat nozzle concept generated the best thrust vectoring performance with an expansion ratio of 1.0 (a cavity in between two equal minimum areas). The variable expansion ratio geometry did not provide the expected improvements in discharge coefficient and system thrust ratio throughout the flight envelope of typical a supersonic aircraft. At mid-climb and cruise conditions, the variable geometry design compromised thrust vector angle achieved, but some thrust vector control would be available, potentially for aircraft trim. The fixed area, expansion ratio of 1.0, Dual Throat Nozzle provided the best overall compromise for thrust vectoring and nozzle internal performance over the range of NPR tested compared to the variable geometry Dual Throat Nozzle.

Flamm, Jeffrey D.↗

A 3D High-Order Unstructured Finite-Volume Algorithm for Solving Maxwell's Equations

A three-dimensional finite-volume algorithm based on arbitrary basis functions for time-dependent problems on general unstructured grids is developed. The method is applied to the time-domain Maxwell equations. Discrete unknowns are volume integrals or cell averages of the electric and magnetic field variables. Spatial terms are converted to surface integrals using the Gauss curl theorem. Polynomial basis functions are introduced in constructing local representations of the fields and evaluating the volume and surface integrals. Electric and magnetic fields are approximated by linear combinations of these basis functions. Unlike other unstructured formulations used in Computational Fluid Dynamics, the new formulation actually does not reconstruct the field variables at each time step. Instead, the spatial terms are calculated in terms of unknowns by precomputing weights at the beginning of the computation as functions of cell geometry and basis functions to retain efficiency. Since no assumption is made for cell geometry, this new formulation is suitable for arbitrarily defined grids, either smooth or unsmooth. However, to facilitate the volume and surface integrations, arbitrary polyhedral cells with polygonal faces are used in constructing grids. Both centered and upwind schemes are formulated. It is shown that conventional schemes (second order in Cartesian grids) are equivalent to the new schemes using first degree polynomials as the basis functions and the midpoint quadrature for the integrations. In the new formulation, higher orders of accuracy are achieved by using higher degree polynomial basis functions. Furthermore, all the surface and volume integrations are carried out exactly. Several model electromagnetic scattering problems are calculated and compared with analytical solutions. Examples are given for cases based on 0th to 3rd degree polynomial basis functions. In all calculations, a centered scheme is applied in the interior, while an upwind matching scheme is employed at material interfaces and the Engquist-Majda non-reflecting boundary condition is implemented at the numerical outer boundaries. The staggered leapfrog scheme and the Runge-Kutta methods are utilized for the time integration. Excellent agreements are found between the numerical and analytical solutions.

Liu, Yen↗

Transients in Pacific/North American Plate Boundary Deformation: Synthesis and Modeling of GPS and Borehole Strain Observations

This is the Final Technical Report on research conducted between 1 June 1997 and 14 September 2001 entitled "Transients in Pacific/North American plate boundary deformation: Synthesis and modeling of GPS and borehole strain observations." As the project title implies, our effort involved a geodetic study of strain transients, i.e., temporal variations in deformation rates, that occur within plate boundary zones and their relationship to earthquakes and plate motions. Important transients occur during and following large earthquakes, and there are also strain transients not apparently associated with earthquakes. A particularly intriguing class of transients, for which there is a modest but growing list of examples, are preseismic anomalies. Such earthquake precursors, if further documented and understood, would have obvious importance for earthquake hazard mitigation. Because the timescales for these diverse transients range over at least 6 orders of magnitude (minutes to years), no single geodetic technique is optimum. We therefore undertook a systematic synthesis of Global Positioning Satellite (GPS) and borehole strainmeter data in three areas in California where there are adequate numbers of both types of instruments (or their equivalent): the San Francisco Bay region (within the Bay Area Regional Deformation network), southern California (within the Southern California Integrated GPS Network), and Parkfield (where a two-color laser system provides a proxy for continuous GPS measurements). An integral component of our study was the elucidation of the physical mechanisms by which such transients occur and propagate. We therefore initiated the development of multiple forward models, using two independent approaches. In the first, we explored the response to specified earthquake slip in viscoelastic models that incorporated failure criteria and the geometry of major faults in California. In the second approach, we examined the dynamical response of a complex rheological medium to the application of a far-field stress imposed by plate motions. The forward models were used both to gain insight into the range of strain transients to be expected under different assumed mechanical conditions and to develop representations for strain fields that allow GPS, borehole, and other strain data to be combined in a self-consistent, yet well-determined, manner. The models also provided a basis for hypothesis testing, by which data from a strain transient well characterized by GPS and borehole observations were utilized to distinguish among competing candidates for the causative physical mechanism and the governing physical characteristics. During the three years of this project, continued to a fourth year through a no-cost extension of the grant, we published 14 papers and presented or co-authored 37 papers at national scientific meetings.

Solomon, Sean C.↗

Verification and Validation of a Three-Dimensional Composite Impact Model with Tabulated Input

A material model which incorporates several key capabilities which have been identified by the aerospace community as lacking in the composite impact models currently available in the commercial transient dynamic finite element code LS-DYNA® has been developed. The material model utilizes experimentally based tabulated input to define the evolution of plasticity and damage as opposed to specifying discrete input parameters (such as modulus and strength. The plasticity portion of the orthotropic, three-dimensional, macroscopic composite constitutive model is based on an extension of the Tsai-Wu composite failure model into a generalized yield function with a non-associative flow rule. The capability to account for the rate and temperature dependent deformation response of composites has also been incorporated into the material model. For the damage model, a strain equivalent formulation is utilized to allow for the uncoupling of the deformation and damage analyses. In the damage model, a diagonal damage tensor is defined to account for the directionally dependent variation of damage. However, the terms in the damage matrix are semi-coupled such that the damage in a particular coordinate direction is a function of the stresses and plastic strains in all of the coordinate directions. For the failure model, a tabulated approach is utilized in which a stress or strain based invariant is defined as a function of the location of the current stress state in stress space to define the initiation of failure, which allows an arbitrarily shaped failure surface to be defined. A systematic series of validation and verification studies, at a variety of length scales ranging from single element simulations to simulations of a flat panel impact test, have been performed to fully exercise and evaluate the capabilities of the developed model.

Impact↗

Dynamic Mode Decomposition of Unsteady Pressure-Sensitive Paint Measurements for the NASA Unitary Plan Wind Tunnel Tests

This paper describes the Dynamic Mode Decomposition (DMD) of the pressures on the scale model of the Space Launch System (SLS) Block 1 cargo vehicle with the Unsteady Pressure-Sensitive Paint (uPSP) measurements, which were collected in the Ascent Transient Aerodynamics Tests with the Unitary Plan Wind Tunnel 11-by-11-foot Transonic Wind Tunnel in September 2019 at NASA Ames Research Center. The work described in this paper is a part of NASA’s development of a new state-of-the-art uPSP capability in production wind tunnels. The conventional DMD algorithm is based on the Singular Value Decomposition (SVD) of the data matrix. For the matrix of the uPSP measurements of the SLS ATAT, the number of rows is equal to the number of nodes in the grid of the scale model, and the number of columns is equal to the number of frames in the videos taken with 4 Phantom high-speed cameras. In this paper, it is verified that, for the time series with zero mean value, the DMD is equivalent to the decomposition with the Discrete Fourier Transform (DFT). Considering the uPSP is mainly used in the assessment of the unsteady, aerodynamic phenomena, the DMD of the uPSP measurements can be implemented in two steps: (1) subtract the mean value from the uPSP measurement on each of the grid nodes; (2) apply the Fast Fourier Transform (FFT) on the resulting zero-mean time series. The DMD of the uPSP measurements with FFT has two advantages: (1) the computational complexity of FFT is O(N*logN), where N is the length of the time series; (2) compared to the SVD-based DMD algorithm, the DMD with FFT can be easily implemented in parallel processing. A sample matrix of uPSP measurements, at the size of 341 grid nodes and 128 frames, is generated. Figures 1 and 2 show the eigenvalues and the ratios of the eigenvectors, respectively, of the sample matrix, without and with the mean value removed on each of the grid nodes, computed with the SVD-based DMD and the FFT. The figures demonstrate the equivalence of the SVD-based DMD and the decomposition with DFT/FFT for the time series with zero mean value. The results of DMD of the uPSP measurements of the SLS ATAT in September 2019 are presented in the paper. The DMD modes at different frequencies are shown, the aerodynamic phenomena (e.g. shockwave and vortex shedding) are demonstrated and the correlation of the DMD modes with the test configuration parameter (e.g., the Mach Number) is discussed. Figure 3 shows a software tool to visualize the DMD modes. The code to implement the algorithm described in this paper was written in C, with libraries of FFTW for FFT and MPI/OpenMP for parallel processing, and executed on the NASA Pleiades supercomputer. Funding for this research was provided by the NASA Aerosciences Evaluation and Test Capabilities Project.

Pressure-Sensitive Paint↗

A Comparison of Multivariable Control Design Techniques for a Turbofan Engine Control

This paper compares two previously published design procedures for two different multivariable control design techniques for application to a linear engine model of a jet engine. The two multivariable control design techniques compared were the Linear Quadratic Gaussian with Loop Transfer Recovery (LQG/LTR) and the H-Infinity synthesis. The two control design techniques were used with specific previously published design procedures to synthesize controls which would provide equivalent closed loop frequency response for the primary control loops while assuring adequate loop decoupling. The resulting controllers were then reduced in order to minimize the programming and data storage requirements for a typical implementation. The reduced order linear controllers designed by each method were combined with the linear model of an advanced turbofan engine and the system performance was evaluated for the continuous linear system. Included in the performance analysis are the resulting frequency and transient responses as well as actuator usage and rate capability for each design method. The controls were also analyzed for robustness with respect to structured uncertainties in the unmodeled system dynamics. The two controls were then compared for performance capability and hardware implementation issues.

Garg, Sanjay↗

GMS-based"Future Time" Rainfall Data Assimilation for Mesoscale Weather Prediction over Korean Peninsula and Future Prospects with GPM Satellite Measurements

This study examines the use of satellite-derived nowcasted (short-term forecasted) rainfall over 3-hour time periods to gain an equivalent time increment in initializing a nonhydrostatic mesoscale model used for predicting convective rainfall events over the Korean peninsula. Infrared (IR) window measurements from the Japanese Geostationary Meteorological Satellite (GMS) are used to specify latent heating for a spinup period of the model - but in future time -- thus initializing in advance of actual time in the framework of a prediction scenario. The main scientific objective of the study is to investigate the strengths and weaknesses of this approach insofar as data assimilation, in which the nowcasted assimilation data are derived independently of the prognostic model itself. Although there have been various recent improvements in formulating the dynamics, thermodynamics, and microphysics of mesoscale models, as well as computer advances which allow the use of high resolution cloud-resolving grids and explicit latent heating over regional domains, spinup remains at the forefront of unresolved mesoscale modeling problems. In general, non-realistic spinup limits the skill in predicting the spatial-temporal distribution of convection and precipitation, primarily in the early hours of a. forecast, stemming from standard prognostic variables not representing the initial diabatic heating field produced by the ambient convection and cloud fields. The long-term goal of this research is to improve short-range (12-hour) quantitative precipitation forecasting (QPF) over the Korean peninsula through the use of innovative data assimilation methods based on geosynchronous satellite measurements. As a step in th~s direction, a non-standard data assimilation experiment in conjunction with GMS-retrieved nowcasted rainfall information introduced to the mesoscale model is conducted. The 3-hourly precipitation forecast information is assimilated through nudging the associated diabatic heating during the early stages of a forecast period. This procedure is expected to enhance details in the moisture field during model integration, and thus improve spinup performance, assuming the errors in the future time latent heating data ate less than intrinsic model background errors.

Smith, Eric A.↗

Examination of time-variable input effects in a nonlinear analogue magnetosphere model

The plasma physical analog model (an extension of the damped, harmonic-oscillator dripping faucet model) is employed to consider explicitly the effect of time-varying the inputs. This work is equivalent to considering the effects of northward and southward turnings of the interplanetary magnetic field for various periods of time. It is found that relatively extended episodes (not less than 2 hours) of turned-on input with shorter (about 1 hour) periods of turned-off input lead to model behavior much like the continuously driven case. Going to short input intervals with longer periods of zero input leads to highly irregular and dramatically fluctuating episodes of magnetotail unloading. These results give an insight into the diversity of apparent magnetospheric responses during relatively isolated substorm conditions. This work shows the absolutely critical interdependence (in a nonlinear dynamical system) of input phasing and internal magnetospheric response cycles.

Baker, D. N.↗

CFD Analysis for Assessing the Effect of Wind on the Thermal Control of the Mars Science Laboratory Curiosity Rover

The challenging range of landing sites for which the Mars Science Laboratory Rover was designed, requires a rover thermal management system that is capable of keeping temperatures controlled across a wide variety of environmental conditions. On the Martian surface where temperatures can be as cold as -123 C and as warm as 38 C, the rover relies upon a Mechanically Pumped Fluid Loop (MPFL) Rover Heat Rejection System (RHRS) and external radiators to maintain the temperature of sensitive electronics and science instruments within a -40 C to 50 C range. The RHRS harnesses some of the waste heat generated from the rover power source, known as the Multi Mission Radioisotope Thermoelectric Generator (MMRTG), for use as survival heat for the rover during cold conditions. The MMRTG produces 110 W of electrical power while generating waste heat equivalent to approximately 2000 W. Heat exchanger plates (hot plates) positioned close to the MMRTG pick up this survival heat from it by radiative heat transfer. Winds on Mars can be as fast as 15 m/s for extended periods. They can lead to significant heat loss from the MMRTG and the hot plates due to convective heat pick up from these surfaces. Estimation of this convective heat loss cannot be accurately and adequately achieved by simple textbook based calculations because of the very complicated flow fields around these surfaces, which are a function of wind direction and speed. Accurate calculations necessitated the employment of sophisticated Computational Fluid Dynamics (CFD) computer codes. This paper describes the methodology and results of these CFD calculations. Additionally, these results are compared to simple textbook based calculations that served as benchmarks and sanity checks for them. And finally, the overall RHRS system performance predictions will be shared to show how these results affected the overall rover thermal performance.

wind↗

Ignition and combustion characteristics of metallized propellants

Shakedown and calibration of the experimental apparatus designed to study secondary atomization and ignition characteristics of aluminum slurry propellants was begun and nears completion. The burner was tested over an array of equivalence ratios and oxygen mass fractions to determine the permissible range of operating conditions. A burner flow control system was designed and fabricated which permits independent control of flame stoichiometry and oxidant O2 concentrations. The Sun slurry formulation was selected as the base for initial testing. Laser beam waists at the focal volume of the single particle sizing system were determined to be 396 microns for the Ar-ion and 81 microns for the He-Ne laser. Size calibrations using pinholes were completed and compared to Mie theory for light scattering from particles. A vibrating office droplet generator was used to test the single particle sizing system dynamically by producing a water droplet laden gas flow. The data acquisition and analysis programming was completed and testing is in progress. From initial results, the particle sizing, atomization and ignition detection systems appear to be performing as expected.

Turns, S. R.↗

Assimilation of Altimeter Data into a Quasigeostrophic Model of the Gulf Stream System: Dynamical Considerations - Part 1

The dynamical consequences of constraining a numerical model with sea surface height data have been investigated. The model used for this study is a quasigeostrophic model of the Gulf Stream region. The data that have been assimilated are maps of sea surface height obtained as the superposition of sea surface height variability deduced from the Geosat altimeter measurements and a mean field constructed from historical hydrographic data. The method used for assimilating the data is the nudging technique. Nudging has been implemented in such a way as to achieve a high degree of convergence of the surface model fields toward the observations. The assimilation of the surface data is thus equivalent to the prescription of a surface pressure boundary condition. The authors analyzed the mechanisms of the model adjustment and the characteristics of the resultant equilibrium state when the surface data are assimilated. Since the surface data are the superposition of a mean component and an eddy component, in order to understand the relative role of these two components in determining the characteristics of the final equilibrium state, two different experiments have been considered: in the first experiment only the climatological mean field is assimilated, while in the second experiment the total surface streamfunction field (mean plus eddies) has been used. It is shown that the model behavior in the presence of the surface data constraint can be conveniently described in terms of baroclinic Fofonoff modes. The prescribed mean component of the surface data acts as a 'surface topography' in this problem. Its presence determines a distortion of the geostrophic contours in the subsurface layers, thus constraining the mean circulation in those layers. The intensity of the mean flow is determined by the inflow/outflow conditions at the open boundaries, as well as by eddy forcing and dissipation.

Capotondi, Antonietta↗

Dynamics and Structure of Weakly-Strained Flames In Normal- and Micro-Gravity

Strained laminar flames have been systematically studied, as the understanding of their structure and dynamic behavior is of relevance to turbulent combustion. Most of these studies have been conducted in opposed-jet, stagnation-type flow configurations. Flame studies in stagnation flows also allow for the determination of fundamental flame properties such as laminar flame speeds and extinction states that can be also conveniently modeled in detail. Studies at high strain rates are important in quantifying and understanding the response of vigorously-burning flames under conditions in which the transport time scales become comparable to the chemical time scales. Studies of weakly-strained flames can be of particular interest for all stoichiometries. For example, the laminar flame speeds for any equivalence ratio, phi, can be accurately determined by using the counterflow technique only if measurements are obtained at very low strain rates. Furthermore, near-limit flames can be only stabilized by weak strain rates. Previous studies have shown that weakly-burning flames are particularly sensitive to chain mechanisms, thermal radiation, and unsteadiness. The stabilization and study of weakly-strained flames is complicated by the presence of buoyancy that can render the flames unstable to the point of extinction. Such instabilities are caused either by the induced natural convection or the fact that higher density fluid can find itself on top of a lower density fluid. Thus, the use of microgravity (mu-g) becomes essential in order to provide meaningful insight into this important combustion regime. In view of the foregoing considerations, the main objectives of the program are to: (1) Experimentally determine the laminar flame speed at near-zero strain rates; (2) Experimentally determine the extinction limits of near-limit flames; (3) Experimentally determine the response of near-limit flames to unsteadiness and heat loss; (4) Introduce Digital Particle Image Velocimetry (DPIV) in microgravity; (5) Conduct detailed numerical simulations of the experiments; and (6) Provide physical insight into the underlying physico-chemical mechanisms.

Zhang, Hai↗

Synthetic Jets in Cross-flow: Jets From Orifices of Different Geometry - Part 2

The flow fields of synthetic jets in a cross-flow from orifices of different geometry are investigated. The geometries include a straight, a tapered, a pitched and a cluster of nine orifices, all having the same cross-sectional area through which the perturbation is discharged into the cross-flow. The strength of the jet from the tapered orifice in comparison to that from the straight one is found to be only slightly enhanced. The flow field from the cluster of orifices, when viewed a few equivalent diameters downstream, is similar to that from the single orifice. However, the penetration is somewhat lower in the former case due to the increased mixing of the distributed jets with the cross-flow. The penetration for the pitched configuration is the lowest, as expected. The jet trajectories for the straight and pitched orifices are well represented by correlation equations available for steady jets-in-cross-flow. Distributions of streamwise velocity, vorticity as well as turbulence intensity are documented for various cases. In addition, distributions of phase-averaged velocity and vorticity for the cylindrical and the clustered orifices are presented providing an insight into the flow dynamics.

Milanovic, Ivana M.↗

Inductive Pulsed Plasma Thruster Model with Time-Evolution of Energy and State Properties

A model for pulsed inductive plasma acceleration is presented that consists of a set of circuit equations coupled to both a one-dimensional equation of motion and an equation governing the partitioning of energy. The latter two equations are obtained for the plasma current sheet by treating it as a single element of finite volume and integrating the governing equations over that volume. The integrated terms are replaced where necessary by physically-equivalent quantities that are calculated through the solution of other parts of the governing equation set. The model improves upon previous one-dimensional performance models by permitting the time-evolution of the energy and state properties of the plasma, the latter allowing for the tailoring of the model to different gases that may be chosen as propellants. The time evolution of the various energy modes in the system and the associated plasma properties, calculated for argon propellant, are presented to demonstrate the efficacy of the model. The model produces a result where efficiency is maximized at a given value of the electrodynamic scaling term known as the dynamic impedance parameter. Qualitatively and quantitatively, the model compares favorably with performance measured for two separate inductive pulsed plasma thrusters, with disagreements attributable to simplifying assumptions employed in the generation of the model solution.

Polzin, Kurt A.↗

An adaptive control scheme for a flexible manipulator

The problem of controlling a single link flexible manipulator is considered. A self-tuning adaptive control scheme is proposed which consists of a least squares on-line parameter identification of an equivalent linear model followed by a tuning of the gains of a pole placement controller using the parameter estimates. Since the initial parameter values for this model are assumed unknown, the use of arbitrarily chosen initial parameter estimates in the adaptive controller would result in undesirable transient effects. Hence, the initial stage control is carried out with a PID controller. Once the identified parameters have converged, control is transferred to the adaptive controller. Naturally, the relevant issues in this scheme are tests for parameter convergence and minimization of overshoots during control switch-over. To demonstrate the effectiveness of the proposed scheme, simulation results are presented with an analytical nonlinear dynamic model of a single link flexible manipulator.

Yang, T. C.↗

Numerical Analysis of a Rotating Detonation Engine in the Relative Reference Frame

A two-dimensional, computational fluid dynamic (CFD) simulation of a semi-idealized rotating detonation engine (RDE) is described. The simulation operates in the detonation frame of reference and utilizes a relatively coarse grid such that only the essential primary flow field structure is captured. This construction yields rapidly converging, steady solutions. Results from the simulation are compared to those from a more complex and refined code, and found to be in reasonable agreement. The performance impacts of several RDE design parameters are then examined. Finally, for a particular RDE configuration, it is found that direct performance comparison can be made with a straight-tube pulse detonation engine (PDE). Results show that they are essentially equivalent.

pulse detonation engines↗

Numerical Analysis of a Rotating Detonation Engine in the Relative Reference Frame

A two-dimensional, computational fluid dynamic (CFD) simulation of a semi-idealized rotating detonation engine (RDE) is described. The simulation operates in the detonation frame of reference and utilizes a relatively coarse grid such that only the essential primary flow field structure is captured. This construction yields rapidly converging, steady solutions. Results from the simulation are compared to those from a more complex and refined code, and found to be in reasonable agreement. The performance impacts of several RDE design parameters are then examined. Finally, for a particular RDE configuration, it is found that direct performance comparison can be made with a straight-tube pulse detonation engine (PDE). Results show that they are essentially equivalent.

propulsion↗

Multiobjective Multidisciplinary Optimization of Low-Boom Supersonic Transports Using Multifidelity Models

A multidisciplinary optimization (MDO) method has been developed to design a computational fluid dynamics (CFD) based low-boom configuration that can be obtained from a Pareto solution of a low-fidelity multiobjective MDO problem with mission constraints. This paper refines the developed MDO method by using multifidelity models for CFD-based multiobjective MDO and a better method for the system-level trade between the target low boom ground noise level and the overland range. The refined MDO method can generate a low-boom configuration that satisfies the mission requirements, has the lowest takeoff gross weight and the longest range for the low-boom mission, trims the low-boom cruise flight with fuel redistributions, and has a reversed equivalent area distribution closely matching a low-boom target with ground noise level below 70 PLdB. The validity of the refined MDO method is demonstrated by a design study of a low-boom supersonic transport that carries 40 passengers, flies a low-boom mission with cruise Mach of 1.7 and range of 3500 nm, and cruises overwater at Mach 1.8 with range of 3882 nm. Moreover, the refined MDO method eliminates the difference between the assumed cruise weight for CFD-based low-boom inverse design optimization and the estimated cruise weight of the optimal inverse design solution with respect to the mission requirements.

MDO↗