• Computers and Concrete
• /
• v.7 no.2
• /
• pp.159-167
• /
• 2010

This paper demonstrates numerical techniques for complex large-scale modeling with microplane constitutive theories for reinforced high strength concrete, which for these applications, is defined to be around the 7000 psi (48 MPa) strength as frequently found in protective structural design. Applications involve highly impulsive loads, such as an explosive detonation or impact-penetration event. These capabilities were implemented into the authors' finite element code, ParaAble and the PRONTO 3D code from Sandia National Laboratories. All materials are explicitly modeled with eight-noded hexahedral elements. The concrete is modeled with a microplane constitutive theory, the reinforcing steel is modeled with the Johnson-Cook model, and the high explosive material is modeled with a JWL equation of state and a programmed burn model. Damage evolution, which can be used for erosion of elements and/or for post-analysis examination of damage, is extracted from the microplane predictions and computed by a modified Holmquist-Johnson-Cook approach that relates damage to levels of inelastic strain increment and pressure. Computation is performed with MPI on parallel processors. Several practical analyses demonstrate that large-scale analyses of this type can be reasonably run on large parallel computing systems.

• The Bulletin of The Korean Astronomical Society
• /
• v.46 no.2
• /
• pp.47.3-48
• /
• 2021

We studied the large scale dynamo process in a system forced by helical magnetic field. The dynamo process is basically nonlinear, but can be linearized with 𝛼&𝛽 coefficients and large scale magnetic field $\bar{B}$. This is very useful to the investigation of solar (stellar) dynamo. A coupled semi-analytic equations based on statistical mechanics are used to investigate the exact evolution of 𝛼&𝛽. This equation set needs only magnetic helicity ${\bar{H}}_M({\equiv}{\langle}{\bar{A}}{\cdot}{\bar{B}}{\rangle},\;{\bar{B}}={\nabla}{\times}{\bar{A}})$ and magnetic energy ${\bar{E}}_M({\equiv}{\langle}{\bar{B}}^2{\rangle}/2)$. They are fundamental physics quantities that can be obtained from the dynamo simulation or observation without any artificial modification or assumption. 𝛼 effect is thought to be related to magnetic field amplification. However, in reality the averaged 𝛼 effect decreases very quickly without a significant contribution to ${\bar{B}}$ field amplification. Conversely, 𝛽 effect contributing to the magnetic diffusion maintains a negative value, which plays a key role in the amplification with Laplacian ∇²(= - k²) for the large scale regime. In addition, negative magnetic diffusion accounts for the attenuation of plasma kinetic energy E_V(= 〈 U² 〉/2) (U: plasma velocity) when the system is saturated. The negative magnetic diffusion is from the interaction of advective term - U • ∇ B from magnetic induction equation and the helical velocity field. In more detail, when 'U' is divided into the poloidal component U_pol and toroidal one U_tor in the absence of reflection symmetry, they interact with - B • ∇ U and - U • ∇ B from ∇ × 〈 U × B 〉 leading to 𝛼 effect and (negative) 𝛽 effect, respectively. We discussed this process using the theoretical method and intuitive field structure model supported by the simulation result.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.32 no.2
• /
• pp.106-121
• /
• 2020

In order to make harbor outskirt facilities robust using the reliability-based design, probabilistic models of wave heights at varying stage of shoaling process optimized for Korean sea waves are prerequisite. In this rationale, we numerically simulate the nonlinear shoaling process of random waves over the beach with a sandbar at its foreshore. In doing so, comprehensive numerical models made of spatially filtered Navier-Stokes Eq., LES [Large Eddy Simulation], dynamic Smagorinsky turbulence closure were used. Considering the characteristics of swells observed at the east coast of Korean Peninsula, random waves were simulated using JONSWAP wave spectrum of various peak enhancement coefficients and random phase method. The coefficients of probabilistic models proposed in this study are estimated from the results of frequency analysis of wave crests and its associated trough detected by Wave by Wave Analysis of the time series of numerically simulated free surface displacements based on the threshold crossing method. Numerical results show that Modified Glukhovskiy wave height distribution, the most referred probabilistic models at finite water depth in the literature, over-predicts the occurring probability of relatively large and small wave heights, and under predicts the occurrence rate of waves of moderate heights. On the other hand, probabilistic models developed in this study show vary encouraging agreements. In addition, the discrepancy of the Modified Glukhovskiy distribution from the measured one are most visible over the surf zone, and as a result, the Modified Glukhovskiy distribution should be applied with caution for the reliability-based design of harbor outskirt facilities deployed near the surf-zone.