4.3NAMar 14, 2011
Variational Multiscale Proper Orthogonal Decomposition: Convection-Dominated Convection-Diffusion EquationsTraian iliescu, Zhu Wang
We introduce a variational multiscale closure modeling strategy for the numerical stabilization of proper orthogonal decomposition reduced-order models of convection-dominated equations. As a first step, the new model is analyzed and tested for convection-dominated convection-diffusion equations. The numerical analysis of the finite element discretization of the model is presented. Numerical tests show the increased numerical accuracy over the standard reduced-order model and illustrate the theoretical convergence rates.
4.3NAOct 17, 2017
An ensemble algorithm for numerical solutions to deterministic and random parabolic PDEsYan Luo, Zhu Wang
In this paper, we develop an ensemble-based time-stepping algorithm to efficiently find numerical solutions to a group of linear, second-order parabolic partial differential equations (PDEs). Particularly, the PDE models in the group could be subject to different diffusion coefficients, initial conditions, boundary conditions, and body forces. The proposed algorithm leads to a single discrete system for the group with multiple right-hand-side vectors by introducing an ensemble average of the diffusion coefficient functions and using a new semi-implicit time integration method. The system could be solved more efficiently than multiple linear systems with a single right-hand-side vector. We first apply the algorithm to deterministic parabolic PDEs and derive a rigorous error estimate that shows the scheme is first-order accurate in time and is optimally accurate in space. We then extend it to find stochastic solutions of parabolic PDEs with random coefficients and put forth an ensemble-based Monte Carlo method. The effectiveness of the new approach is demonstrated through theoretical analysis. Several numerical experiments are presented to illustrate our theoretical results.
3.3HCMay 31, 2020
A Virtual Obstacle Course within Diverse Sensory EnvironmentsZhu Wang, Anat Lubetzky, Charles Hendee et al.
We developed a novel assessment platform with untethered virtual reality, 3-dimensional sounds, and pressure sensing floor mat to help assess the walking balance and negotiation of obstacles given diverse sensory load and/or cognitive load. The platform provides an immersive 3D city-like scene with anticipated/unanticipated virtual obstacles. Participants negotiate the obstacles with perturbations of: auditory load by spatial audio, cognitive load by a memory task, and visual flow by generated by avatars movements at various amounts and speeds. A VR headset displays the scenes while providing real-time position and orientation of the participant's head. A pressure-sensing walkway senses foot pressure and visualizes it in a heatmap. The system helps to assess walking balance via pressure dynamics per foot, success rate of crossing obstacles, available response time as well as head kinematics in response to obstacles and multitasking. Based on the assessment, specific balance training and fall prevention program can be prescribed.
1.2CODec 23, 2019
MM for Penalized EstimationZhu Wang
Penalized estimation can conduct variable selection and parameter estimation simultaneously. The general framework is to minimize a loss function subject to a penalty designed to generate sparse variable selection. The majorization-minimization (MM) algorithm is a computational scheme for stability and simplicity, and the MM algorithm has been widely applied in penalized estimation. Much of the previous work have focused on convex loss functions such as generalized linear models. When data are contaminated with outliers, robust loss functions can generate more reliable estimates. Recent literature has witnessed a growing impact of nonconvex loss-based methods, which can generate robust estimation for data contaminated with outliers. This article investigates MM algorithm for penalized estimation, provide innovative optimality conditions and establish convergence theory with both convex and nonconvex loss functions. With respect to applications, we focus on several nonconvex loss functions, which were formerly studied in machine learning for regression and classification problems. Performance of the proposed algorithms are evaluated on simulated and real data including healthcare costs and cancer clinical status. Efficient implementations of the algorithms are available in the R package mpath in CRAN.
5.6HCFeb 8, 2019
Virtual Environments for Rehabilitation of Postural Control DysfunctionZhu Wang, Anat Lubetzky, Marta Gospodarek et al.
We developed a novel virtual reality [VR] platform with 3-dimensional sounds to help improve sensory integration and visuomotor processing for postural control and fall prevention in individuals with balance problems related to sensory deficits, such as vestibular dysfunction (disease of the inner ear). The system has scenes that simulate scenario-based environments. We can adjust the intensity of the visual and audio stimuli in the virtual scenes by controlling the user interface (UI) settings. A VR headset (HTC Vive or Oculus Rift) delivers stereo display while providing real-time position and orientation of the participants' head. The 3D game-like scenes make participants feel immersed and gradually exposes them to situations that may induce dizziness, anxiety or imbalance in their daily-living.
1.2NAMay 25, 2017
An efficient algorithm for simulating ensembles of parameterized flow problemsMax Gunzburger, Nan Jiang, Zhu Wang
Many applications of computational fluid dynamics require multiple simulations of a flow under different input conditions. In this paper, a numerical algorithm is developed to efficiently determine a set of such simulations in which the individually independent members of the set are subject to different viscosity coefficients, initial conditions, and/or body forces. The proposed scheme applied to the flow ensemble leads to need to solve a single linear system with multiple right-hand sides, and thus is computationally more efficient than solving for all the simulations separately. We show that the scheme is nonlinearly and long-term stable under certain conditions on the time-step size and a parameter deviation ratio. Rigorous numerical error estimate shows the scheme is of first-order accuracy in time and optimally accurate in space. Several numerical experiments are presented to illustrate the theoretical results.