• Interaction and multiscale mechanics
• /
• v.5 no.1
• /
• pp.27-40
• /
• 2012

The two-dimensional incompressible flow of a linear viscoelastic fluid we considered in this research has rapidly oscillating initial conditions which contain both the large scale and small scale information. In order to grasp this double-scale phenomenon of the complex flow, a multiscale analysis method is developed based on the mathematical homogenization theory. For the incompressible flow of a linear viscoelastic Maxwell fluid, a well-posed multiscale system, including averaged equations and cell problems, is derived by employing the appropriate multiple scale asymptotic expansions to approximate the velocity, pressure and stress fields. And then, this multiscale system is solved numerically using the pseudospectral algorithm based on a time-splitting semi-implicit influence matrix method. The comparisons between the multiscale solutions and the direct numerical simulations demonstrate that the multiscale model not only captures large scale features accurately, but also reflects kinetic interactions between the large and small scale of the incompressible flow of a linear viscoelastic fluid.

• Interaction and multiscale mechanics
• /
• v.2 no.4
• /
• pp.353-374
• /
• 2009

A new seamless multiscale simulation was developed for coupling the continuum model with its molecular dynamics. Kriging-based Finite Element Method (K-FEM) is employed to model the continuum base of the entire domain, while the molecular dynamics (MD) is confined in a localized domain of interest. In the coupling zone, where the MD domain overlaps the continuum model, the overall Hamiltonian is postulated by contributions from the continuum and the molecular overlays, based on a quartic spline scaling parameter. The displacement compatibility in this coupling zone is then enforced by the Lagrange multiplier technique. A multiple-time-step velocity Verlet algorithm is adopted for its time integration. The validation of the present method is reported through numerical tests of one dimensional atomic lattice. The results reveal that at the continuum/MD interface, the commonly reported spurious waves in the literature are effectively eliminated in this study. In addition, the smoothness of the transition from MD to the continuum can be significantly improved by either increasing the size of the coupling zone or expanding the nodal domain of influence associated with K-FEM.

• Journal of Mechanical Science and Technology
• /
• v.20 no.1
• /
• pp.110-124
• /
• 2006

The multi scale wavelet-Galerkin method implemented in an adaptive manner has an advantage of obtaining accurate solutions with a substantially reduced number of interpolation points. The method is becoming popular, but its numerical efficiency still needs improvement. The objectives of this investigation are to present a new numerical scheme to improve the performance of the multi scale adaptive wavelet-Galerkin method and to give detailed implementation procedure. Specifically, the subdomain technique suitable for multiscale methods is developed and implemented. When the standard wavelet-Galerkin method is implemented without domain subdivision, the interaction between very long scale wavelets and very short scale wavelets leads to a poorly-sparse system matrix, which considerably worsens numerical efficiency for large-sized problems. The performance of the developed strategy is checked in terms of numerical costs such as the CPU time and memory size. Since the detailed implementation procedure including preprocessing and stiffness matrix construction is given, researchers having experiences in standard finite element implementation may be able to extend the multi scale method further or utilize some features of the multiscale method in their own applications.

• Nuclear Engineering and Technology
• /
• v.36 no.3
• /
• pp.229-236
• /
• 2004

Radiation hardening is a multiscale phenomenon involving various processes over a wide range of time and length. We present a multiscale model for estimating the amount of radiation hardening in pressure vessel steel in the environment of a light water reactor. The model comprises two main parts: molecular dynamics (MD) simulation and a point defect cluster (PDC) model. The MD simulation was used to investigate the primary damage caused by displacement cascades. The PDC model mathematically formulates interactions between point defects and their clusters, which explains the evolution of microstructures. We then used a dislocation barrier model to calculate the hardening due to the PDCs. The key input for this multiscale model is a neutron spectrum at the inner surface of reactor pressure vessel steel of the Younggwang Nuclear Power Plant No.5. A combined calculation from the MD simulation and the PDC model provides a convenient tool for estimating the amount of radiation hardening.

• Korean Chemical Engineering Research
• /
• v.45 no.6
• /
• pp.582-590
• /
• 2007

In view of the analysis of the phenomena and the prediction of the performance, mathematical modelling and simulation of a high temperature pilot reactor for water gas shift reaction (WGSR) has been carried out. Multiscale simulation incorporated computational fluid dynamics (CFD) technique, which has the capability to deal with the reactor shape, fluid and energy transport with extensive degree of accuracy, and process modeling technique, which, in turn is responsible for reaction kinetics and mass transport. This research employed multiscale simulation and the results were compared with those from process simulation. From multiscale simulation, the maximum conversion of was predicted approximately 0.85 and the maximum temperature at the reactor was calculated 720 K, resulting from the heat of reaction. Dynamic simulation was also performed for the time transient profile of temperature, conversion, etc. Considering the results, it is concluded that multiscale simulation is a safe and accurate technique to predict reactor behaviors, and consequently will be available for the design of commercial size chemical reactors as well as other commercial unit operations.

• The Korean Journal of Applied Statistics
• /
• v.24 no.6
• /
• pp.995-1006
• /
• 2011

A time series can be decomposed into simple components with a multiscale method. Empirical mode decomposition(EMD) is a recently invented multiscale method in Huang et al. (1998). It is natural to apply a classical prediction method such a vector autoregressive(AR) model to the obtained simple components instead of the original time series; in addition, a prediction procedure combining a classical prediction model to EMD and Hilbert spectrum is proposed in Kim et al. (2008). In this paper, we suggest to adopt principal component analysis(PCA) to the prediction procedure that enables the efficient selection of input variables among obtained components by EMD. We discuss the utility of adopting PCA in the prediction procedure based on EMD and Hilbert spectrum and analyze the daily worm account data by the proposed PCA adopted prediction method.

• Computers and Concrete
• /
• v.33 no.3
• /
• pp.325-340
• /
• 2024

The growing demands for sustainable and high-performance construction materials necessitate a deep understanding of their physicochemical properties by that of these heterogeneities. This paper presents a comprehensive review of the state-of-the-art hierarchical multiscale modeling approach aimed at predicting the intricate physicochemical characteristics of construction materials. Emphasizing the heterogeneity inherent in these materials, the review briefly introduces single-scale analyses, including the ab initio method, molecular dynamics, and micromechanics, through a scale-bridging technique. Herein, the limitations of these models are also overviewed by that of effectively scale-bridging methods of length or time scales. The hierarchical multiscale model demonstrates these physicochemical properties considering chemical reactions, material defects from nano to macro scale, microscopic properties, and their influence on macroscopic events. Thereby, hierarchical multiscale modeling can facilitate the efficient design and development of next-generation construction.

• Composites Research
• /
• v.36 no.5
• /
• pp.329-334
• /
• 2023

The classical multiscale finite element (FE² ) method involves iterative calculations of micro-boundary value problems for representative volume elements at every integration point in macro scale, making it a computationally time and data storage space. To overcome this, we developed the data-driven multiscale analysis method based on the mean-field homogenization (MFH). Data-driven computational mechanics (DDCM) analysis is a model-free approach that directly utilizes strain-stress datasets. For performing multiscale analysis, we efficiently construct a strain-stress database for the microstructure of composite materials using mean-field homogenization and conduct data-driven computational mechanics simulations based on this database. In this paper, we apply the developed multiscale analysis framework to an example, confirming the results of data-driven computational mechanics simulations considering the microstructure of a hyperelastic composite material. Therefore, the application of data-driven computational mechanics approach in multiscale analysis can be applied to various materials and structures, opening up new possibilities for multiscale analysis research and applications.

• Nuclear Engineering and Technology
• /
• v.41 no.1
• /
• pp.11-20
• /
• 2009

Radiation effects on materials are inherently multiscale phenomena in view of the fact that various processes spanning a broad range of time and length scales are involved. A multiscale modeling approach to embrittlement of pressure vessel steels is presented here. The approach includes an investigation of the mechanisms of defect accumulation, microstructure evolution and the corresponding effects on mechanical properties. An understanding of these phenomena is required to predict the behavior of structural materials under irradiation. We used molecular dynamics (MD) simulations at an atomic scale to study the evolution of high-energy displacement cascade reactions. The MD simulations yield quantitative information on primary damage. Using a database of displacement cascades generated by the MD simulations, we can estimate the accumulation of defects over diffusional length and time scales by applying kinetic Monte Carlo simulations. The evolution of the local microstructure under irradiation is responsible for changes in the physical and mechanical properties of materials. Mechanical property changes in irradiated materials are modeled by dislocation dynamics simulations, which simulate a collective motion of dislocations that interact with the defects. In this paper, we present a multi scale modeling methodology that describes reactor pressure vessel embrittlement in a light water reactor environment.

• Interaction and multiscale mechanics
• /
• v.1 no.1
• /
• pp.1-31
• /
• 2008

The present paper presents multiscale modelling via coupling of the discrete and finite element methods. Theoretical formulation of the discrete element method using spherical or cylindrical particles has been briefly reviewed. Basic equations of the finite element method using the explicit time integration have been given. The micr-macro transition for the discrete element method has been discussed. Theoretical formulations for macroscopic stress and strain tensors have been given. Determination of macroscopic constitutive properties using dimensionless micro-macro relationships has been proposed. The formulation of the multiscale DEM/FEM model employing the DEM and FEM in different subdomains of the same body has been presented. The coupling allows the use of partially overlapping DEM and FEM subdomains. The overlap zone in the two coupling algorithms is introduced in order to provide a smooth transition from one discretization method to the other. Coupling between the DEM and FEM subdomains is provided by additional kinematic constraints imposed by means of either the Lagrange multipliers or penalty function method. The coupled DEM/FEM formulation has been implemented in the authors' own numerical program. Good performance of the numerical algorithms has been demonstrated in a number of examples.