• 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.

• Structural Engineering and Mechanics
• /
• v.15 no.6
• /
• pp.669-683
• /
• 2003

The stochastic optimal nonlinear control of coupled adjacent building structures is studied based on the stochastic dynamical programming principle and the stochastic averaging method. The coupled structures with control devices under random seismic excitation are first condensed to form a reduced-order structural model for the control analysis. The stochastic averaging method is applied to the reduced model to yield stochastic differential equations for structural modal energies as controlled diffusion processes. Then a dynamical programming equation for the energy processes is established based on the stochastic dynamical programming principle, and solved to determine the optimal nonlinear control law. The seismic response mitigation of the coupled structures is achieved through the structural energy control and the dimension of the optimal control problem is reduced. The seismic excitation spectrum is taken into account according to the stochastic dynamical programming principle. Finally, the nonlinear controlled structural response is predicted by using the stochastic averaging method and compared with the uncontrolled structural response to evaluate the control efficacy. Numerical results are given to demonstrate the response mitigation capabilities of the proposed stochastic optimal control method for coupled adjacent building structures.

• Computers and Concrete
• /
• v.8 no.2
• /
• pp.125-142
• /
• 2011

In the present paper, a numerical simulation method based on mesoscopic composite structure of concrete, the truss network model, is developed to evaluate the diffusivity of concrete in order to account for the microstructure of concrete, the binding effect of chloride ions and the chloride concentration dependence. In the model, concrete is described as a three-phase composite, consisting of mortar, coarse aggregates and the interfacial transition zones (ITZs) between them. The advantage of the current model is that it can easily represent the movement of mass (e.g. water or chloride ions) through ITZs or the potential cracks within concrete. An analytical method to estimate the chloride diffusivity of mortar and ITZ, which are both treated as homogenious materials in the model, is introduced in terms of water-to-cement ratio (w/c) and sand volume fraction. Using the newly developed approaches, the effect of cracking of concrete on chloride diffusion is reflected by means of the similar process as that in the test. The results of calculation give close match with experimental observations. Furthermore, with consideration of the binding capacity of chloride ions to cement paste and the concentration dependence for diffusivity, the one-dimensional nonlinear diffusion equation is established, as well as its finite difference form in terms of the truss network model. A series of numerical analysises performed on the model find that the chloride diffusion is substantially influenced by the binding capacity and concentration dependence, which is same as that revealed in some experimental investigations. This indicates the necessity to take into account the binding capacity and chloride concentration dependence in the durability analysis and service life prediction of concrete structures.

• Structural Engineering and Mechanics
• /
• v.81 no.4
• /
• pp.523-527
• /
• 2022

The objective of this study is to simulate the flow in pipes with various boundary conditions. Free-pressure fluid model, is used in the pipe based on Navier-Stokes equation. The models are solved by using the numerical method. A problem called "stability of pipes" is used in order to compare frequency and critical fluid velocity. When the initial conditions of problem satisfied the instability conditions, the free-pressure model could accurately predict discontinuities in the solution field. Employing nonlinear strains-displacements, stress-strain energy method the governing equations were derived using Hamilton's principal. Differential quadrature method (DQM) is used for obtaining the frequency and critical fluid velocity. The results of this paper are analyzed by hyperbolic numerical method. Results show that the level of numerical diffusion in the solution field and the range of well-posedness are two important criteria for selecting the two-fluid models. The solutions for predicting the flow variables is approximately equal to the two-pressure model 2. Therefore, the predicted pressure changes profile in the two-pressure model is more consistent with actual physics. Therefore, in numerical modeling of gas-liquid two-phase flows in the vertical pipe, the present model can be applied.

• Wind and Structures
• /
• v.36 no.2
• /
• pp.133-147
• /
• 2023

The random vibration of saddle membrane structures under wind load is studied theoretically and experimentally. First, the nonlinear random vibration differential equations of saddle membrane structures under wind loads are established based on von Karman's large deflection theory, thin shell theory and potential flow theory. The probabilistic density function (PDF) and its corresponding statistical parameters of the displacement response of membrane structure are obtained by using the diffusion process theory and the Fokker Planck Kolmogorov equation method (FPK) to solve the equation. Furthermore, a wind tunnel test is carried out to obtain the displacement time history data of the test model under wind load, and the statistical characteristics of the displacement time history of the prototype model are obtained by similarity theory and probability statistics method. Finally, the rationality of the theoretical model is verified by comparing the experimental model with the theoretical model. The results show that the theoretical model agrees with the experimental model, and the random vibration response can be effectively reduced by increasing the initial pretension force and the rise-span ratio within a certain range. The research methods can provide a theoretical reference for the random vibration of the membrane structure, and also be the foundation of structural reliability of membrane structure based on wind-induced response.

• Fire Science and Engineering
• /
• v.10 no.2
• /
• pp.28-39
• /
• 1996

A recently developed electro-thermal simulation methodology is used to analyze the behavior of a PWM(Pulse-Width-Modulated) voltage source inverter which uses IGBT(Insulated Gate Bipolar Transistor) as the switching devices. In the electro-thermal network simulation methdology, the simulator solves for the temperature distribution within the power semiconductor devices(IGBT electro-thermal model), control logic circuitry, the IGBT gate drivers, the thermal network component models for the power silicon chips, package, and heat sinks as well as the current and voltage within the electrical network. The thermal network describes the flow of heat form the chip surface through the package and heat sink and thus determines the evolution of the chip surface temperature used by the power semiconductor device models. The thermal component model for the device silicon chip, packages, and heat sink are developed by discretizing the nonlinear heat diffusion equation and are represented in component from so that the thermal component models for various package and heat sink can be readily connected to on another to form the thermal network.

• 한국신재생에너지학회:학술대회논문집
• /
• 2009.06a
• /
• pp.373-373
• /
• 2009

Relative humidity, membrane conductivity and water activity are critical parameters of polymer electrolyte membrane fuel cells (PEMFC) for high performance and reliability. These parameters are closely related with temperature. Moreover, the ideal values of these parameters are not always identical along the channels. Therefore, the cooling channel design and its operating condition should be well optimized along the all location of the channels. In the present study, we have performed a numerical investigation on the effects of cooling channels on performance of a PEMFC. Three-dimensional Navier-Stokes equations are solved with the energy equation including heat generated by the electrochemical reactions in the fuel cell. The present numerical model includes the gas diffusion layers (GDL) and serpentine channels for both anode and cathode gas flows, as well as cooling channels. To accurately predict the water transport across the membrane, the distribution of water content in the membrane is calculated by solving a nonlinear differential equation with a nonlinear coefficient, i.e., the water diffusivity which is a function of water content as well as temperature. Main emphasis is placed on the heat transfer between the solid bipolar plate and coolant flow. The present results show that local current density is affected by cooling channels due to the change of the oxygen concentration and the membrane conductivity as well as the water content. It is also found that the relative humidity is influenced by the generated water and the gas temperature and thus it affects the distribution of fuel concentration and the conductivity of the membrane, ultimately fuel cell performance. Unit-cell experiments are also carried out to validate the numerical models. The performance curves between the models and experiments show reasonable results.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.36 no.11
• /
• pp.1297-1303
• /
• 2012

In this paper, the standard procedures for measuring the moisture sorption properties of thin polymeric materials such as polyethylene terephthalate (PET) by using the thermo-gravimetric method to characterize the moisture diffusion in the polymer are presented, and the sorption properties are quantified. The moisture diffusivity and solubility are characterized to investigate the effect of temperature and humidity on the moisture sorption properties according to the Arrhenius equation. The validation of the obtained sorption properties using thermogravimetry is discussed with the measured permeability based on Fickian diffusion. The nonlinear behavior of the concentration dependent moisture diffusion is investigated experimentally, and the nonlinearity is characterized numerically for the case of having an interface with an inorganic material such as a metal. The Fickian/Non-Fickian model based on the obtained moisture sorption properties is compared and discussed.

• Water for future
• /
• v.27 no.2
• /
• pp.73-83
• /
• 1994

The unproper development of the nearshore zone can enhance the diffusion of pollutant in the nearshore zone resulting in unbalanced sediment budget of beach which causes alteration of beach topography. Therefore, it is required to predict the effects of the envirnmental change quantitatively. In this paper, the depth-averaged and time-averaged energy balance equation is selected to acount for the wave transformation such as refraction, shoaling effect, the surf zone energy disipation, wave breaking index and bore, due to wave breaking in the shore region.(Numerical solutions are obtained by a finite difference method, ADI and Upwind. For the calculation of the wave-induced current, the unsteady nonlinear depth-averaged and time-averaged governing equation is derived based on the continuity and momentum equation for imcompressible fluid.) Numerical solutions are obtained by finite difference method considering influences of factors such as lateral mixing coefficient, bed shear stress, wave direction angle, wave steepness, wave period and bottom slope. The model is applied to the computation of wave transformation, wave-induced current and variation of mean water leel on a uniformly sloping beach.

• Journal of the Korean Society of Safety
• /
• v.11 no.4
• /
• pp.143-150
• /
• 1996

This is an explicit-Implicit, finite element analysis for linear as well as nonlinear hygrothermal stress problems. Additional features, such as moisture diffusion equation, crack element and virtual crack extension(VCE ) method for evaluating J-integral are implemented in this program. The Linear Elastic Fracture Mechanics(LEFM) Theory is employed to estimate the crack driving force under the transient condition for and existing crack. Pores in materials are assumed to be saturated with moisture in the liquid form at the room temperature, which may vaporize as the temperature increases. The vaporization effects on the crack driving force are also studied. The Ideal gas equation is employed to estimate the thermodynamic pressure due to vaporization at each time step after solving basic nodal values. A set of field equations governing the time dependent response of porous media are derived from balance laws based on the mixture theory Darcy's law Is assumed for the fluid flow through the porous media. Perzyna's viscoplastic model incorporating the Von-Mises yield criterion are implemented. The Green-Naghdi stress rate is used for the invariant of stress tensor under superposed rigid body motion. Isotropic elements are used for the spatial discretization and an iterative scheme based on the full newton-Raphson method is used for solving the nonlinear governing equations.