In this talk, I will introduce briefly a moving mesh method based on harmonic mapping. As a rare character, the unique existence of the harmonic mapping is the basic motivation for us to develop this method. The method is implemented in finite element and an iterative procedure is adopted to avoid mesh tangling caused by numerical factors. Our method can move the mesh interior the domain and mesh on the boundary in coupling for both 2D and 3D problems. The moving mesh module can be a black box added on the whole solver to the PDE under consideration that it is very convenient for coding - no modifications to the solver of the PDE are required. The inter-mesh mesh updation is implemented by a linear convection equation instead of generally adopted interpolation methods, thus the method can be easily to applied to problem as incompressible Navier-Stokes equation where the divergence free interpolation can be a big problem, and problem as conservation laws where the conservative interpolation is not trival to be implemented. Numerical results including viscos Burgers equation, reaction-diffusion equation, incompressible Navier-Stokes equation and its coupling with level set method will be shown.
Indiana University-Purdue University in Indianapolis
Time:
Thursday, May 12, 2005 - 2:00pm
Location:
MSTB 254
Probability of a gap in the spectrum of a matrix from the
Gaussian Unitary Ensemble is given by a Fredholm determinant. Its asymptotics when the gap becomes large is an interesting problem related to Painleve equations, random permutations, etc. These asymptotics for the gap in the bulk of the spectrum were conjectured by Dyson. The proof was given over the years by Widom, Deift, Its , Zhou, and the speaker. In particular, the proof of the multiplicative constant in the asymptotics
was the last difficulty recently resolved. I will explain the method of determining this constant and the rest of the asymptotics (applicable also to other important Fredholm, and also Hankel, and Toeplitz determinants where the corresponding constant is not yet determined). The method uses the Riemann-Hilbert approach. This part of the talk will be based on the works of Deift, Its, Zhou, and the speaker.
Partial differential equations were discretized using the mortar finite element method, where the mortar space contains piecewise quadratic and cubic functions. We first proved the wellposedness of the saddle point system. To solve this saddle point system, the existing domain decomposition (DD) algorithm of the linear system solver requires the communication of both the interface solution and the interface residual of every local problem. We have developed a new algorithm where only the interface solution is communicated, to accommodate globally nonconforming meshes. The resulting communication complexity is reduced. The scalability and parallel efficiency of this new algorithm were tested with a highly adaptive mesh. As the number of processors increases, linear and logarithmic speed-up of the solving time were observed with this new linear system solver, for the convection-diffusion equation and the Poisson equation respectively.
Three-dimensional dynamic simulation in computational biology provides an emerging field for the application of efficient distributed linear system solvers. We developed a 3D continuum model to investigate the role of structural and functional cellular components in regulating synchronized calcium signaling (SCS), characterized by high gradient near the t-tubule membrane and low gradient in the cytoplasm along the transverse direction, which enables ventricular myocyte to respond rapidly and forcefully to electrical andchemical stimuli. The distributed linear system solver improved the simulation speed by ~10 folds. Simulation results suggest that both t-tubule structure and the spatially heterogeneous distribution of calcium-handling-proteins are important for SCS. The model also predicts that two aspects of heterogeneous distribution are required: the concentration of calcium-handling-proteins in the t-tubule membrane to be ~6 times of that in the surface membrane; and the concentration of L-type calcium channels, in the cytoplasmic end of the t-tubule, to be ~2.3 times of that in the surface membrane end. These results have provided a foundation for further studies on the effects of three- dimensional t-tubule geometry and ion channel distribution on calcium dynamics.
We will discuss two approaches for modeling cancer growth. One approach is to employ a deterministic space-independent model that includes separate components to represent specific and non-specific immune function. The other approach employs a spatially dependent hybrid cellular-automata (HCA) model whose rules are driven by cellular metabolic function. For the determinstic model, numerical simulations of mixed chemo-immuno therapy and vaccine therapy using both mouse and human parameters are presented. We illustrate situations for which neither chemotherapy nor immunotherapy alone are sufficient to control tumor growth, but in combination the therapies are able to eliminate the entire tumor burden. The HCA model is in its initial stages of development, and preliminary results will be presented.