For the aerodynamic design of aircraft, CFD has made a remarkable progress to provide lift, drag and surface pressure distribution collaborating with wing tunnel experiment. However, not much information is known about the flow phenomena in wakes of airplanes and wings. In the article, we concentrate on wake flow analysis as a post work of computations and experiments done for APC-I workshops. Through the analysis, we have found that wing tip vortices and boundary layers of a wing affect velocity distribution on a wake, while only the tip vortices do pressure. The precise prediction of wake physics deeply depends on grid resolution quality in subspaces such as cross-sectional planes in a wake away from an airplane surface, for both of experiment and computation.
The paper reports numerical simulations of aerodynamic characteristics of a commercial aircraft model known as the NASA-CRM. Specific test cases and detailed experimental results have been provided by JAXA for the CFD workshop “Aerodynamic Prediction Challenge (APC).” We employed the cell-vertex finite-volume method on unstructured grids and carried out the RANS simulations with the k-ω SST turbulence model to present results in the APC-I. After the verification of the model in a benchmark problem, both alpha-sweep simulations from low- to high-angles of attack and the grid-convergence study at a cruise conditions are carried out. The simulations well reproduced the experimental results at moderate angles of attack by both types of computational grids used in the study but the results showed the dependency on the grid type and grid resolutions at the highest angle of attack of 5.72 degrees due to different predictions of flow separation near the root chord of the main wing.
Thanks to recent lunar exploration missions, high-resolution lunar surface observation data was obtained. In future lunar exploration, it is requested to land at a specific point in which scientific interest is higher than other areas. SLIM project is demonstrating pinpoint landing technology, which is achieved, by combination of "autonomous image-based high-precision navigation technology" and "autonomous guidance technology aims to generate fuel optimum landing trajectory". This paper presents powered descending trajectory design in terms of trajectory optimization. As usually considered in general space mission development, an optimal solution in terms of minimum fuel consumption is the basis of investigation. This study addresses trajectory optimization considering specific objective functions derived from practical constraints regarding mission design, e.g., altitude, length of downrange, and visibility from ground stations. In this paper, nominal trajectory design considering minimum fuel consumption is presented, followed by parametric studies to identify sensitivity to change of initial conditions with which powered descending starts. Finally, trajectory optimization results with various type of objective function are presented.
The computational grid dependency is an important problem for CFD. We have computed aerodynamics on NASA-Common Research Model (CRM) with FaSTAR and various grids to investigate the grid dependency. We employed four grids: two Cartesian-based unstructured grids, a tetrahedral unstructured grid, and a hexahedral structured grid. The computational conditions are based on the test cases of Aerodynamics Prediction Challenge (APC). First, the grid convergence at a fixed angle of attack and the trend of an angle-of-attack sweep are compared between the four grids. The lift coefficients computed with the two similar Cartesian-based grids are different, and this is caused by the grid difference around the leading edge. However, the overall trend of angle-of-attack sweep is almost same between the four grids. Next, we computed aerodynamics on NASA-CRM with a support device to investigate the support interference. It is found that the support interference on the drag and pitching moment is large and should be considered.
Ceramic/metal brazing was investigated to produce light-weight and highly-efficient ceramic thrusters. Silicon nitride ceramic and metal bars were brazed using an Ag-based brazing material. Four-point bend tests were conducted at room and high temperatures to evaluate the strength of the brazed joints. Computational fluid dynamics (CFD) and finite element method (FEM) analyses were also performed to investigate the effect of the construction and shape of the joints on the stress distribution around them. It was demonstrated that brazing was a great candidate as the joining technique, and a 20 N ceramics/metal brazed thruster was successfully produced.
The conventional spectral volume (SV) method for three-dimensional tetrahedral unstructured mesh is extended to use hexahedral mesh. In the test calculations of linear scalar advection problem and diffusion problem, the formal spatial order of accuracy is achieved even for skewed computational meshes. The Spalart Allmaras turbulence model implemented in the present code is verified by calculating the grid-converged skin friction of turbulent boundary layer on a flat plate. Then verified code is applied to compute the flowfield around the NASA-CRM. In this study, we examine the agreement of the computed aerodynamic coefficients with the corresponding experimental data for different angles of attack. We also examine the change of aerodynamic coefficients with varying Reynolds number, although the experimental data is not available. It is shown that the present SV code can predict aerodynamic coefficients around the cruise angle of attack conditions fairly well. When higher Reynolds number is assumed, the computed lift increases while the viscous drag is reduced, as is expected. On the other hand, the pressure drag is increased with increasing Reynolds number due to the shock wave on the wing which moves toward downstream.
In this paper, the authors propose a novel landing method named “Two-step Landing Method” for small lunar lander which is needed to be designed considering constrains from envelope area of rocket and the weight of the lander. The proposed method enforces intentional body tumbling at the contact of main leg. We analyzed its dynamics by three-dimensional simulations which consider lander’s attitude and lateral velocity and landing site’s slope angle. Numerical simulation models have been designed on Mechanical Dynamics Software “ADAMS”, and lander models refer to “SLIM” which is a small lander proposed by ISAS/JAXA. It is found that the proposed landing method can land on steep slope by tilting body attitude toward inclination direction of landing site. Especially in the case of landing with lateral residual velocity, the proposed method has higher landing stability than conventional landing method.
This paper focuses on the Evolutional Triangle Similarity Matching (ETSM) method for estimating spacecraft location in Smart Lander for Investigating Moon (SLIM) mission and improves it by adding the functions of elimination of line symmetric triangles between crater map and camera shot image, comparison of rotation relationship of triangles and triangle formation method using Delaunay triangulation and introducing point group matching as a coordinate calculation function. To evaluate the robustness of the improved method, we conduct simulation experiments using the crater map and camera shot images in six situations. This experiments have revealed the following implications: (1) this method improved accuracy of location estimation within 5.1 pixels by the functions of elimination of line symmetric triangles between crater map and camera shot image, (2) this method slight got worse accuracy at low or high altitude of spacecraft, however, this method successfully reduced incorrect spacecraft location estimation by comparison of rotation relationship of triangles, (3) this method improved accuracy of location estimation by triangle formation method using Delaunay triangulation, but possibility of incorrect spacecraft location estimation is slight increased, and (4) integration method of these three mechanism can estimate spacecraft location within 5 pixels without being affected altitude difference and rotation of camera shot image.
Aerodynamics coefficients of NASA common research model (NASA-CRM) are computationally investigated by revisiting some test cases appeared in aerodynamics prediction challenge (APC) using high-order discontinuous Galerkin (DG) methods. While the employed high-order DG methods reasonably well predict the overall aerodynamics for the NASA-CRM, some discrepancies appear between the CFD and experimental data, especially in a lift and a pitching moment at low angle of attack. Although those discrepancies still exist, installing a sting enhances the reproduction of experimental data of pitching moment. Additionally, fluid-structure coupled simulations used to consider aero-elastic deformations give better agreements of the pitching moment with experimental data.
In response to the First and Second Aerodynamics Prediction Challenges, held in Tokyo, July 2015 and in Kanazawa, July 2016, respectively, computational fluid dynamics simulations were performed for the NASA Common Research Model using the Tohoku University Aerodynamic Simulation (TAS) Code. Our results were summarized in this manuscript, with an emphasis on key computational techniques and mesh generation methods. Unstructured hybrid meshes were generated using the Mixed-Element Grid Generator in 3 Dimensions (MEGG3D), and were deformed based on wing deformation data obtained during wind tunnel testing. The effects of support system interference, of mesh density and of laminar to turbulent boundary layer transition are shown to discuss the validity of computational results. Aerodynamic coefficients were well predicted at low angles of attack when the support system interference effect was considered, while an accurate prediction of pitching moment at high angles of attack was challenging because mesh density affected the shock location on the wing and the size of side-of-body separation.
In order to validate CFD results using wind tunnel test (WTT) data, it is necessary to consider what types of corrections are applied to the WTT data, such as wall interference correction, near-field support interference correction and so on. In this paper, a CFD tool called “Cflow” developed by Kawasaki Heavy Industries, Ltd. is validated using results of wind tunnel test conducted by JAXA with NASA-CRM (Common Research Model). At first, effects of near-field support interference and wing deformation on aerodynamic performance are discussed. Then it is confirmed that Cflow results are well consistent with experimental results, taking these two effects into account.
SLIM (Smart Lander for Investigating Moon) is the Lunar Landing Demonstrator which is under development at ISAS/JAXA. SLIM demonstrates not only so-called Pin-Point Landing Technique to the lunar surface, but also demonstrates the design to make the explorer small and lightweight. Realizing the compact explorer is one of the key points to achieve the frequent lunar and planetary explorations. This paper summarizes the preliminary system design of SLIM, especially the way to reduce the size.
SLIM project which aim for pin-point landing on the moon surface. For achieving this plan, it is necessary to estimate the flight position of the space probe. The estimation is performed by matching the detected craters with database. This paper introduces a crater detection method using Principal Component Analysis (PCA) and its evaluation. This method is capable of real-time processing under low computational resources such as Field-Programmable Gate Array (FPGA). In this research, we report improvement of robustness at detection and high accuracy of crater size measurement.
Energy absorbing system for landing gears is an important on the SLIM project. Open cell porous aluminum manufactured through 3D selective laser melting (SLM) process has been applied on the energy absorbing system. Compressive tests for cylindrical and hemispherical shaped porous aluminum with different porosities revealed the high potential as an energy absorbing component. Heat treatment after SLM processing is effective to increase the energy absorbing potential of the porous aluminum.
The landing radar employs a pulse-type radar using 4.3 GHz C-band microwave radiation. It has a wide beam for measuring the altitude in vertical direction, as well as four narrow tilted beams for measuring the velocity in horizontal direction. In this paper, development of the Bread Board Model (BBM), a field experiment, and the design of SLIM loading Flight Model (FM) are introduced. Furthermore, the radar simulator required for FM development of a radar is explained.
Next generation moon landing mission will require autonomous pinpoint landing capability because of requirements for landing on specific terrains in a limited area. This capability requires precise absolute self-localization of the lander during braking descent phase. The purpose of this paper is to propose an algorithm to estimate the lander position and to evaluate its mountability to a space-grade FPGA. In this method, the position estimation is performed by matching crater point patterns with database point patterns by finding topological correspondences using crater-based linear features. In addition, we confirmed the resource amount and the calculation time when this algorithm is implemented on the FPGA using high-level synthesis.
Inviscid and Reynolds-averaged Navier-Stokes (RANS) simulations of transonic flows around the NASA Common Research Model are conducted using the Cartesian flow solver UTCart. The immersed boundary method is used to represent the smooth geometry surfaces on the Cartesian grids. The wall function is combined with the immersed boundary method to reproduce the turbulent boundary layer on the geometry surface in the RANS simulations. In the inviscid calculations, the qualitative flow feature including the position on the shock-wave on the wing shows agreement with the reference result a body-fitted grid. In the RANS calculations, the trend of pitching moment and drag shows fair agreement with the reference result, while prediction of the flow separation at high angle of attack is still difficult. Compared with the reference result, the differences in the total drag coefficient at a moderate angle of attack on the medium grid (33 million cells) and the fine grid (99 million cells) are 31 drag counts (10%) and 20 drag counts (6.5%), respectively. Furthermore, each of the calculated aerodynamic coefficients shows a consistent trend of grid convergence toward the reference result.
Summary of First Aerodynamics Prediction Challenge (APC-I) is presented. The APC-I is a domestic CFD prediction workshop that was held on July 3, 2015. The test cases include aerodynamic prediction of NASA-CRM with and without aeroelastic effects, and its wake flow prediction. We compare the CFD results with JAXA’s wind tunnel measurements. There are 15 participants from government, academia, industry, and commercial. The CFD results submitted from the participants are compared and discussed.
Summary of Second Aerodynamics Prediction Challenge (APC-II) is presented. The APC-II is a domestic CFD prediction workshop that was held on July 6, 2016. The test cases include aerodynamic prediction of NASA-CRM with and without support effects, and buffet prediction. We compare the CFD results with JAXA’s wind tunnel measurements. There are 9 participants from national research agency, academia, industry, and commercial software vendor. The CFD results submitted from the participants are compared and discussed in the presentation.