• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.26 no.4
• /
• pp.223-231
• /
• 2013

This paper presents a time-domain Gauss-Newton full-waveform inversion method for the material profile reconstruction in heterogeneous semi-infinite solid media. To implement the inverse problem in a finite computational domain, perfectly-matchedlayers( PMLs) are introduced as wave-absorbing boundaries within which the domain's wave velocity profile is to be reconstructed. The inverse problem is formulated in a partial-differential-equations(PDE)-constrained optimization framework, where a least-squares misfit between measured and calculated surface responses is minimized under the constraint of PML-endowed wave equations. A Gauss-Newton-Krylov optimization algorithm is utilized to iteratively update the unknown wave velocity profile with the aid of a specialized regularization scheme. Through a series of one-dimensional examples, the solution of the Gauss-Newton inversion was close enough to the target profile, and showed superior convergence behavior with reduced wall-clock time of implementation compared to a conventional inversion using Fletcher-Reeves optimization algorithm.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.26 no.4
• /
• pp.213-221
• /
• 2013

This paper introduces a mesh continuation scheme for a one-dimensional inverse medium problem to reconstruct the spatial distribution of elastic wave velocities in heterogeneous semi-infinite solid domains. To formulate the inverse problem, perfectly-matched-layers(PMLs) are introduced as wave-absorbing boundaries that surround the finite computational domain truncated from the originally semi-infinite extent. To tackle the inverse problem in the PML-truncated domain, a partial-differential-equations(PDE)-constrained optimization approach is utilized, where a least-squares misfit between calculated and measured surface responses is minimized under the constraint of PML-endowed wave equations. The optimization problem iteratively solves for the unknown wave velocities with their updates calculated by Fletcher-Reeves conjugate gradient algorithms. The optimization is performed using a mesh continuation scheme through which the wave velocity profile is reconstructed in successively denser mesh conditions. Numerical results showed the robust performance of the mesh continuation scheme in reconstructing target wave velocity profile in a layered heterogeneous solid domain.

• Journal of applied mathematics & informatics
• /
• v.8 no.3
• /
• pp.811.1-818
• /
• 2001

• Wind and Structures
• /
• v.17 no.1
• /
• pp.87-105
• /
• 2013

Modeling an equilibrium atmospheric boundary layer (ABL) in an empty computational domain has routinely been performed with the k-${\varepsilon}$ turbulence model. However, the research objects of structural wind engineering are bluff bodies, and the SST k-${\omega}$ turbulence model is more widely used in the numerical simulation of flow around bluff bodies than the k-${\varepsilon}$ turbulence model. Therefore, to simulate an equilibrium ABL based on the SST k-${\omega}$ turbulence model, the inlet profiles of the mean wind speed U, turbulence kinetic energy k, and specific dissipation rate ${\omega}$ are proposed, and the source terms for the U, k and ${\omega}$ are derived by satisfying their corresponding transport equations. Based on the proposed inlet profiles, numerical comparative studies with and without considering the source terms are carried out in an empty computational domain, and an actual numerical simulation with a trapezoidal hill is further conducted. It shows that when the source terms are considered, the profiles of U, k and ${\omega}$ are all maintained well along the empty computational domain and the accuracy of the actual numerical simulation is greatly improved. The present study could provide a new methodology for modeling the equilibrium ABL problem and for further CFD simulations with practical value.

• Wind and Structures
• /
• v.15 no.4
• /
• pp.313-329
• /
• 2012

There have existed for a number of years good practice guidelines for the use of Computational Fluid Dynamics (CFD) in the field of wind engineering. As part of those guidelines, details are given for the size of flow domain that should be used around a building of height, H. For low-rise buildings, the domain sizes produced by following the guidelines are reasonable and produce results that are largely free from blockage effects. However, when high-rise or tall buildings are considered, the domain size based solely on the building height produces very large domains. A large domain, in most cases, leads to a large cell count, with many of the cells in the grid being used up in regions far from the building/wake region. This paper challenges this domain size guidance by looking at the effects of changing the domain size around a tall building. The RNG ${\kappa}-{\varepsilon}$ turbulence model is used in a series of steady-state solutions where the only parameter varied is the domain size, with the mesh resolution in the building/wake region left unchanged. Comparisons between the velocity fields in the near-field of the building and pressure coefficients on the building are used to inform the assessment. The findings of the work for this case suggest that a domain of approximately 10% the volume of that suggested by the existing guidelines could be used with a loss in accuracy of less than 10%.

• Journal of applied mathematics & informatics
• /
• v.16 no.1_2
• /
• pp.383-405
• /
• 2004

We study the incomprssible Navier Stokes equations for the flow inside contraction geometry. The governing equations are expressed in the vorticity-stream function formulations. A rectangular computational domain is arised by elliptic grid generation technique. The numerical solution is based on a technique of automatic numerical generation of acurvilinear coordinate system by transforming the governing equation into computational plane. The transformed equations are approximated using central differences and solved simultaneously by successive over relaxation iteration. The time dependent of the vorticity equation solved by using explicit marching procedure. We will apply the technique on several irregular-shapes.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.32 no.5
• /
• pp.323-331
• /
• 2019

In this study, a computational model approach for the modeling of hydrodynamic pressures acting on radial gates during strong earthquakes is proposed. The use of the dynamic layering method with the Arbitrary Lagrangian Eulerian (ALE) algorithm and the SIMPLE method for simulating free reservoir surface flow in addition to moving boundary interfaces between the fluid domain and a structure due to earthquake excitation are suggested. The verification and validation of the proposed approach are realized by comparisons performed using the renowned formulation derived by the experimental results for vertical and inclined dam surfaces subjected to earthquake excitation. A parameter study for the truncated lengths of the two-dimensional fluid domain demonstrates that twice the water level leads to efficient and converged computational results. Finally, numerical simulations for large radial gates with different curvatures subjected to two strong earthquakes are successfully performed using the suggested computational model.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.24 no.4
• /
• pp.365-373
• /
• 2011

Damage of reinforced concrete panel subjected to blast load using parallel and domain decomposition is analyzed. The numerical results are sensitive to the mesh size because blast waves are generated during the extremely short term. In order to investigate the effect of mesh size on the blast wave, the analysis results from various wave mesh size using AUTODYN, the explicit finite element analysis program, were compared with existing experimental results. The smaller mesh size was, the higher accuracy was. However, in this case, the analysis was inefficient. Therefore, in order to increase numerical efficiency, the parallel analysis using decomposed method based on Euler and Lagrangian description was performed. Finally, the decomposed method using both the structure domain based on Lagrange description and the blast wave domain based on Euler description was more efficient than the decomposed method using only the Lagrange mesh on structure domain.

• Journal of the Korean Institute of Telematics and Electronics
• /
• v.22 no.2
• /
• pp.18-24
• /
• 1985

Above certain filter order, the frequency domain -adaptive filter is superior to the time domain-adaptive filter in computational complexity. In this paper a new type algorithm, $\mu$-FLMs algorithm, is proposed for the frequency domain- adaptive filter and the characteristics of the proposed algorithm is compared with that of the time domain- adaptive filter algorithm($\mu$-FLMS algorithm). The simulation results showed that under the same convergence rate , the frequency domain-adaptive filter is efficient in compu tational burden.

• International Journal of CAD/CAM
• /
• v.11 no.1
• /
• pp.1-10
• /
• 2011

Product lifecycle management (PLM) is a new business strategy for enterprise's product R&D. A PLM system holds and maintaining the integrity of the product data produced throughout its entire lifecycle. There is, therefore, a need to build a safe and effective product data model to support PLM system. The paper proposes a domain-based product data model for PLM. The domain modeling method is introduced, including the domain concept and its defining standard along the product evolution process. The product data model in every domain is explained, and the mapping rules among these models are discussed. Mapped successively among these models, product data can be successfully realized the dynamic evolution and the historical traceability in PLM system.