• Proceedings of the Korean Society of Computational and Applied Mathematics Conference
• /
• 2003.09a
• /
• pp.1-1
• /
• 2003

In this talk, a reduced-order modeling methodology based on centroidal Voronoi tessellations (CVT's)is introduced. CVT's are special Voronoi tessellations for which the generators of the Voronoi diagram are also the centers of mass (means) of the corresponding Voronoi cells. The discrete data sets, CVT's are closely related to the h-means clustering techniques. Even with the use of good mesh generators, discretization schemes, and solution algorithms, the computational simulation of complex, turbulent, or chaotic systems still remains a formidable endeavor. For example, typical finite element codes may require many thousands of degrees of freedom for the accurate simulation of fluid flows. The situation is even worse for optimization problems for which multiple solutions of the complex state system are usually required or in feedback control problems for which real-time solutions of the complex state system are needed. There hava been many studies devoted to the development, testing, and use of reduced-order models for complex systems such as unsteady fluid flows. The types of reduced-ordered models that we study are those attempt to determine accurate approximate solutions of a complex system using very few degrees of freedom. To do so, such models have to use basis functions that are in some way intimately connected to the problem being approximated. Once a very low-dimensional reduced basis has been determined, one can employ it to solve the complex system by applying, e.g., a Galerkin method. In general, reduced bases are globally supported so that the discrete systems are dense; however, if the reduced basis is of very low dimension, one does not care about the lack of sparsity in the discrete system. A discussion of reduced-ordering modeling for complex systems such as fluid flows is given to provide a context for the application of reduced-order bases. Then, detailed descriptions of CVT-based reduced-order bases and how they can be constructed of complex systems are given. Subsequently, some concrete incompressible flow examples are used to illustrate the construction and use of CVT-based reduced-order bases. The CVT-based reduced-order modeling methodology is shown to be effective for these examples and is also shown to be inexpensive to apply compared to other reduced-order methods.

• Journal of computational fluids engineering
• /
• v.21 no.2
• /
• pp.112-119
• /
• 2016

NEXTSat-1 is the next-generation small-size artificial satellite system planed by the Satellite Technology Research Center(SatTReC) in Korea Advanced Institute of Science and Technology(KAIST). For the control of attitude and transition of the orbit, the system has adopted a RHM(Resisto-jet Head Module), which has a very simple geometry with a reasonable efficiency. An axisymmetric model is devised with two coil-resistance heaters using xenon(Xe) gas, and the minimum required specific impulse is 60 seconds under the thrust more than 30 milli-Newton. To design the module, seven basic parameters should be decided: the nozzle shape, the power distribution of heater, the pressure drop of filter, the diameter of nozzle throat, the slant length and the angle of nozzle, and the size of reservoir, etc. After quasi one-dimensional analysis, a theoretical value of specific impulse is calculated, and the optima of parameters are found out from the baseline with a series of multi-physical numerical simulations based on the compressible Navier-Stokes equations for gas and the heat conduction energy equation for solid. A commercial code, COMSOL Multiphysics is used for the computation with a FEM (finite element method) based numerical scheme. The final values of design parameters indicate 5.8% better performance than those of baseline design after the verification with all the tuned parameters. The present method should be effective to reduce the time cost of trial and error in the development of RHM, the thruster of NEXTSat-1.

• International Journal of Aeronautical and Space Sciences
• /
• v.18 no.2
• /
• pp.186-196
• /
• 2017

In order to resolve the trimming difficulty in rotor CFD calculations, a high-efficiency and improved "delta trim method" is established to compute the blade control settings that are necessary to identify the blade motion. In this method, a simplified model which combines the blade element theory and different inflow models is employed to calculate the control settings according to the target aerodynamic forces, then it is coupled into a CFD solver with unsteady Navier-Stokes equations by the delta methodology, which makes the control settings and aerodynamics calculated and updated in the meantime at every trim cycle. Different from the previous work, the current research combines the inflow model based on prescribed wake theory. Using the method established, the control settings and aerodynamic characteristics of Helishape 7A, AH-1G and Caradonna-Tung rotors are calculated. The influence of different inflow models on trimming calculations is analyzed and the computational efficiency of the current "delta trim method" is compared with that of the "CFD-based trim method". Furthermore, for the sake of improving the calculation efficiency, a novel acceleration factor method is introduced to accelerate the trimming process. From the numerical cases, it is demonstrated that the current "delta trim method" has higher computational efficiency than "CFD-based trim method" in both hover and forward flight, and up to 70% of the amount of calculation can be saved by current "delta trim method" which turns out to be satisfactory for engineering applications. In addition, the proposed acceleration factor shows a good ability to accelerate the trim procedure, and the prescribed wake inflow model is always of better stability than other simple inflow models whether the acceleration factor is utilized in trimming calculations.

• Wind and Structures
• /
• v.15 no.6
• /
• pp.481-494
• /
• 2012

The greenhouse type metal structures are increasingly used in modern construction of livestock farms because they are less laborious to construct and they provide a more favorable microclimate for the growth of animals compared to conventional livestock structures. A key stress factor for metal structures is the wind. The external pressure coefficient ($c_{pe}$) is used for the calculation of the wind effect on the structures. A high pressure coefficient value leads to an increase of the construction weight and subsequently to an increase in the construction cost. The EC1 in conjunction with EN 13031-1:2001, which is specialized for greenhouses, gives values for this coefficient. This value must satisfy two requirements: the safety of the structure and a reduced construction cost. In this paper, the Navier - Stokes and continuity equations are solved numerically with the finite element method (Galerkin Method) in order to simulate the two dimensional, incompressible, viscous air flow over the vaulted roofs of single span and twin-span with eaves livestock greenhouses' structures, with a height of 4.5 meters and with length of span of 9.6 and 14 m. The simulation was carried out in a wind tunnel. The numerical results of pressure coefficients, as well as, the distribution of them are presented and compared with data from Eurocodes for wind actions (EC1, EN 13031-1:2001). The results of the numerical experiment were close to the values given by the Eurocodes mainly on the leeward area of the roof while on the windward area a further segmentation is suggested.

• Ocean Systems Engineering
• /
• v.6 no.1
• /
• pp.61-97
• /
• 2016

This is the second of two papers on the 3D numerical modeling of nearshore hydro- and morphodynamics. In Part I, the focus was on surf and swash zone hydrodynamics in the cross-shore and longshore directions. Here, we consider nearshore processes with an emphasis on the effects of oceanic forcing and beach characteristics on sediment transport in the cross- and longshore directions, as well as on foreshore bathymetry changes. The Delft3D and XBeach models were used with four turbulence closures (viz., ${\kappa}-{\varepsilon}$, ${\kappa}-L$, ATM and H-LES) to solve the 3D Navier-Stokes equations for incompressible flow as well as the beach morphology. The sediment transport module simulates both bed load and suspended load transport of non-cohesive sediments. Twenty sets of numerical experiments combining nine control parameters under a range of bed characteristics and incident wave and tidal conditions were simulated. For each case, the general morphological response in shore-normal and shore-parallel directions was presented. Numerical results showed that the ${\kappa}-{\varepsilon}$ and H-LES closure models yield similar results that are in better agreement with existing morphodynamic observations than the results of the other turbulence models. The simulations showed that wave forcing drives a sediment circulation pattern that results in bar and berm formation. However, together with wave forcing, tides modulate the predicted nearshore sediment dynamics. The combination of tides and wave action has a notable effect on longshore suspended sediment transport fluxes, relative to wave action alone. The model's ability to predict sediment transport under propagation of obliquely incident wave conditions underscores its potential for understanding the evolution of beach morphology at field scale. For example, the results of the model confirmed that the wave characteristics have a considerable effect on the cumulative erosion/deposition, cross-shore distribution of longshore sediment transport and transport rate across and along the beach face. In addition, for the same type of oceanic forcing, the beach morphology exhibits different erosive characteristics depending on grain size (e.g., foreshore profile evolution is erosive or accretive on fine or coarse sand beaches, respectively). Decreasing wave height increases the proportion of onshore to offshore fluxes, almost reaching a neutral net balance. The sediment movement increases with wave height, which is the dominant factor controlling the beach face shape.